Linear reversible second-order cellular automata and their first-order matrix equivalents

International Nuclear Information System (INIS)

Macfarlane, A J

2004-01-01

Linear or one-dimensional reversible second-order cellular automata, exemplified by three cases named as RCA1-3, are introduced. Displays of their evolution in discrete time steps, t=0, 1, 2, ..., from their simplest initial states and on the basis of updating rules in modulo 2 arithmetic, are presented. In these, shaded and unshaded squares denote cells whose cell variables are equal to one and zero respectively. This paper is devoted to finding general formulas for, and explicit numerical evaluations of, the weights N(t) of the states or configurations of RCA1-3, i.e. the total number of shaded cells in tth line of their displays. This is achieved by means of the replacement of RCA1-3 by the equivalent linear first-order matrix automata MCA1-3, for which the cell variables are 2x2 matrices, instead of just numbers (element of Z 2 ) as for RCA1-3. MCA1-3 are tractable because it has been possible to generalize to them the heavy duty methods already well-developed for ordinary first-order cellular automata like those of Wolfram's Rules 90 and 150. While the automata MCA1-3 are thought to be of genuine interest in their own right, with untapped further mathematical potential, their treatment has been applied here to expediting derivation of a large body of general and explicit results for N(t) for RCA1-3. Amongst explicit results obtained are formulas also for each of RCA1-3 for the total weight of the configurations of the first 2 M times, M=0, 1, 2, ..

Symmetry Classification of First Integrals for Scalar Linearizable Second-Order ODEs

Directory of Open Access Journals (Sweden)

K. S. Mahomed

2012-01-01

Full Text Available Symmetries of the fundamental first integrals for scalar second-order ordinary differential equations (ODEs which are linear or linearizable by point transformations have already been obtained. Firstly we show how one can determine the relationship between the symmetries and the first integrals of linear or linearizable scalar ODEs of order two. Secondly, a complete classification of point symmetries of first integrals of such linear ODEs is studied. As a consequence, we provide a counting theorem for the point symmetries of first integrals of scalar linearizable second-order ODEs. We show that there exists the 0-, 1-, 2-, or 3-point symmetry cases. It is shown that the maximal algebra case is unique.

On the classical theory of ordinary linear differential equations of the second order and the Schroedinger equation for power law potentials

International Nuclear Information System (INIS)

Lima, M.L.; Mignaco, J.A.

1983-01-01

The power law potentials in the Schroedinger equation solved recently are shown to come from the classical treatment of the singularities of a linear, second order differential equation. This allows to enlarge the class of solvable power law potentials. (Author) [pt

Schwarz maps of algebraic linear ordinary differential equations

Science.gov (United States)

Sanabria Malagón, Camilo

2017-12-01

A linear ordinary differential equation is called algebraic if all its solution are algebraic over its field of definition. In this paper we solve the problem of finding closed form solution to algebraic linear ordinary differential equations in terms of standard equations. Furthermore, we obtain a method to compute all algebraic linear ordinary differential equations with rational coefficients by studying their associated Schwarz map through the Picard-Vessiot Theory.

Convolution of second order linear recursive sequences II.

Directory of Open Access Journals (Sweden)

Szakács Tamás

2017-12-01

Full Text Available We continue the investigation of convolutions of second order linear recursive sequences (see the first part in [1]. In this paper, we focus on the case when the characteristic polynomials of the sequences have common root.

Runge-Kutta Methods for Linear Ordinary Differential Equations

Science.gov (United States)

Zingg, David W.; Chisholm, Todd T.

1997-01-01

Three new Runge-Kutta methods are presented for numerical integration of systems of linear inhomogeneous ordinary differential equations (ODES) with constant coefficients. Such ODEs arise in the numerical solution of the partial differential equations governing linear wave phenomena. The restriction to linear ODEs with constant coefficients reduces the number of conditions which the coefficients of the Runge-Kutta method must satisfy. This freedom is used to develop methods which are more efficient than conventional Runge-Kutta methods. A fourth-order method is presented which uses only two memory locations per dependent variable, while the classical fourth-order Runge-Kutta method uses three. This method is an excellent choice for simulations of linear wave phenomena if memory is a primary concern. In addition, fifth- and sixth-order methods are presented which require five and six stages, respectively, one fewer than their conventional counterparts, and are therefore more efficient. These methods are an excellent option for use with high-order spatial discretizations.

Algebraic properties of first integrals for systems of second-order ...

African Journals Online (AJOL)

Symmetries of the rst integrals for scalar linear or linearizable second- order ordinary differential equations (ODEs) have already been derived and shown to exhibit interesting properties. One of these is that the symmetry algebra sl(3; R ) is generated by the three triplets of symmetries of the functionally independent first ...

Nonlocal symmetries of a class of scalar and coupled nonlinear ordinary differential equations of any order

International Nuclear Information System (INIS)

Pradeep, R Gladwin; Chandrasekar, V K; Senthilvelan, M; Lakshmanan, M

2011-01-01

In this paper, we devise a systematic procedure to obtain nonlocal symmetries of a class of scalar nonlinear ordinary differential equations (ODEs) of arbitrary order related to linear ODEs through nonlocal relations. The procedure makes use of the Lie point symmetries of the linear ODEs and the nonlocal connection to deduce the nonlocal symmetries of the corresponding nonlinear ODEs. Using these nonlocal symmetries, we obtain reduction transformations and reduced equations to specific examples. We find that the reduced equations can be explicitly integrated to deduce the general solutions for these cases. We also extend this procedure to coupled higher order nonlinear ODEs with specific reference to second-order nonlinear ODEs. (paper)

ADM For Solving Linear Second-Order Fredholm Integro-Differential Equations

Science.gov (United States)

Karim, Mohd F.; Mohamad, Mahathir; Saifullah Rusiman, Mohd; Che-Him, Norziha; Roslan, Rozaini; Khalid, Kamil

2018-04-01

In this paper, we apply Adomian Decomposition Method (ADM) as numerically analyse linear second-order Fredholm Integro-differential Equations. The approximate solutions of the problems are calculated by Maple package. Some numerical examples have been considered to illustrate the ADM for solving this equation. The results are compared with the existing exact solution. Thus, the Adomian decomposition method can be the best alternative method for solving linear second-order Fredholm Integro-Differential equation. It converges to the exact solution quickly and in the same time reduces computational work for solving the equation. The result obtained by ADM shows the ability and efficiency for solving these equations.

Second-order kinetic model for the sorption of cadmium onto tree fern: a comparison of linear and non-linear methods.

Science.gov (United States)

Ho, Yuh-Shan

2006-01-01

A comparison was made of the linear least-squares method and a trial-and-error non-linear method of the widely used pseudo-second-order kinetic model for the sorption of cadmium onto ground-up tree fern. Four pseudo-second-order kinetic linear equations are discussed. Kinetic parameters obtained from the four kinetic linear equations using the linear method differed but they were the same when using the non-linear method. A type 1 pseudo-second-order linear kinetic model has the highest coefficient of determination. Results show that the non-linear method may be a better way to obtain the desired parameters.

Focal decompositions for linear differential equations of the second order

Directory of Open Access Journals (Sweden)

L. Birbrair

2003-01-01

two-points problems to itself such that the image of the focal decomposition associated to the first equation is a focal decomposition associated to the second one. In this paper, we present a complete classification for linear second-order equations with respect to this equivalence relation.

Lie group classification of first-order delay ordinary differential equations

Science.gov (United States)

Dorodnitsyn, Vladimir A.; Kozlov, Roman; Meleshko, Sergey V.; Winternitz, Pavel

2018-05-01

A group classification of first-order delay ordinary differential equations (DODEs) accompanied by an equation for the delay parameter (delay relation) is presented. A subset of such systems (delay ordinary differential systems or DODSs), which consists of linear DODEs and solution-independent delay relations, have infinite-dimensional symmetry algebras—as do nonlinear ones that are linearizable by an invertible transformation of variables. Genuinely nonlinear DODSs have symmetry algebras of dimension n, . It is shown how exact analytical solutions of invariant DODSs can be obtained using symmetry reduction.

New results for exponential synchronization of linearly coupled ordinary differential systems

International Nuclear Information System (INIS)

Tong Ping; Chen Shi-Hua

2017-01-01

This paper investigates the exponential synchronization of linearly coupled ordinary differential systems. The intrinsic nonlinear dynamics may not satisfy the QUAD condition or weak-QUAD condition. First, it gives a new method to analyze the exponential synchronization of the systems. Second, two theorems and their corollaries are proposed for the local or global exponential synchronization of the coupled systems. Finally, an application to the linearly coupled Hopfield neural networks and several simulations are provided for verifying the effectiveness of the theoretical results. (paper)

On nonnegative solutions of second order linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander; Vodstrčil, Petr

2004-01-01

Roč. 32, č. 1 (2004), s. 59-88 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z1019905 Keywords : second order linear functional differential equations * nonnegative solution * two-point boundary value problem Subject RIV: BA - General Mathematics

Higher derivative discontinuous solutions to linear ordinary differential equations: a new route to complexity?

International Nuclear Information System (INIS)

Datta, Dhurjati Prasad; Bose, Manoj Kumar

2004-01-01

We present a new one parameter family of second derivative discontinuous solutions to the simplest scale invariant linear ordinary differential equation. We also point out how the construction could be extended to generate families of higher derivative discontinuous solutions as well. The discontinuity can occur only for a subset of even order derivatives, viz., 2nd, 4th, 8th, 16th,.... The solutions are shown to break the discrete parity (reflection) symmetry of the underlying equation. These results are expected to gain significance in the contemporary search of a new dynamical principle for understanding complex phenomena in nature

An introduction to linear ordinary differential equations using the impulsive response method and factorization

CERN Document Server

Camporesi, Roberto

2016-01-01

This book presents a method for solving linear ordinary differential equations based on the factorization of the differential operator. The approach for the case of constant coefficients is elementary, and only requires a basic knowledge of calculus and linear algebra. In particular, the book avoids the use of distribution theory, as well as the other more advanced approaches: Laplace transform, linear systems, the general theory of linear equations with variable coefficients and variation of parameters. The case of variable coefficients is addressed using Mammana’s result for the factorization of a real linear ordinary differential operator into a product of first-order (complex) factors, as well as a recent generalization of this result to the case of complex-valued coefficients.

On the Liouvillian solution of second-order linear differential equations and algebraic invariant curves

International Nuclear Information System (INIS)

Man, Yiu-Kwong

2010-01-01

In this communication, we present a method for computing the Liouvillian solution of second-order linear differential equations via algebraic invariant curves. The main idea is to integrate Kovacic's results on second-order linear differential equations with the Prelle-Singer method for computing first integrals of differential equations. Some examples on using this approach are provided. (fast track communication)

Numerical solution of second-order stochastic differential equations with Gaussian random parameters

Directory of Open Access Journals (Sweden)

Rahman Farnoosh

2014-07-01

Full Text Available In this paper, we present the numerical solution of ordinary differential equations (or SDEs, from each orderespecially second-order with time-varying and Gaussian random coefficients. We indicate a complete analysisfor second-order equations in specially case of scalar linear second-order equations (damped harmonicoscillators with additive or multiplicative noises. Making stochastic differential equations system from thisequation, it could be approximated or solved numerically by different numerical methods. In the case oflinear stochastic differential equations system by Computing fundamental matrix of this system, it could becalculated based on the exact solution of this system. Finally, this stochastic equation is solved by numericallymethod like E.M. and Milstein. Also its Asymptotic stability and statistical concepts like expectationand variance of solutions are discussed.

Order-constrained linear optimization.

Science.gov (United States)

Tidwell, Joe W; Dougherty, Michael R; Chrabaszcz, Jeffrey S; Thomas, Rick P

2017-11-01

Despite the fact that data and theories in the social, behavioural, and health sciences are often represented on an ordinal scale, there has been relatively little emphasis on modelling ordinal properties. The most common analytic framework used in psychological science is the general linear model, whose variants include ANOVA, MANOVA, and ordinary linear regression. While these methods are designed to provide the best fit to the metric properties of the data, they are not designed to maximally model ordinal properties. In this paper, we develop an order-constrained linear least-squares (OCLO) optimization algorithm that maximizes the linear least-squares fit to the data conditional on maximizing the ordinal fit based on Kendall's τ. The algorithm builds on the maximum rank correlation estimator (Han, 1987, Journal of Econometrics, 35, 303) and the general monotone model (Dougherty & Thomas, 2012, Psychological Review, 119, 321). Analyses of simulated data indicate that when modelling data that adhere to the assumptions of ordinary least squares, OCLO shows minimal bias, little increase in variance, and almost no loss in out-of-sample predictive accuracy. In contrast, under conditions in which data include a small number of extreme scores (fat-tailed distributions), OCLO shows less bias and variance, and substantially better out-of-sample predictive accuracy, even when the outliers are removed. We show that the advantages of OCLO over ordinary least squares in predicting new observations hold across a variety of scenarios in which researchers must decide to retain or eliminate extreme scores when fitting data. © 2017 The British Psychological Society.

Factorization of a class of almost linear second-order differential equations

International Nuclear Information System (INIS)

Estevez, P G; Kuru, S; Negro, J; Nieto, L M

2007-01-01

A general type of almost linear second-order differential equations, which are directly related to several interesting physical problems, is characterized. The solutions of these equations are obtained using the factorization technique, and their non-autonomous invariants are also found by means of scale transformations

Some oscillation criteria for the second-order linear delay differential equation

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2011-01-01

Roč. 136, č. 2 (2011), s. 195-204 ISSN 0862-7959 Institutional research plan: CEZ:AV0Z10190503 Keywords : second-order linear differential equation with a delay * oscillatory solution Subject RIV: BA - General Mathematics http://www.dml.cz/handle/10338.dmlcz/141582

Myshkis type oscillation criteria for second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2015-01-01

Roč. 178, č. 1 (2015), s. 143-161 ISSN 0026-9255 Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillation criteria Subject RIV: BA - General Mathematics Impact factor: 0.664, year: 2015 http://link.springer.com/article/10.1007%2Fs00605-014-0719-y

Exact Solutions for Certain Nonlinear Autonomous Ordinary Differential Equations of the Second Order and Families of Two-Dimensional Autonomous Systems

Directory of Open Access Journals (Sweden)

M. P. Markakis

2010-01-01

Full Text Available Certain nonlinear autonomous ordinary differential equations of the second order are reduced to Abel equations of the first kind ((Ab-1 equations. Based on the results of a previous work, concerning a closed-form solution of a general (Ab-1 equation, and introducing an arbitrary function, exact one-parameter families of solutions are derived for the original autonomous equations, for the most of which only first integrals (in closed or parametric form have been obtained so far. Two-dimensional autonomous systems of differential equations of the first order, equivalent to the considered herein autonomous forms, are constructed and solved by means of the developed analysis.

Second order oscillations of a Vlasov-Poisson plasma in the Fourier transformed space

International Nuclear Information System (INIS)

Sedlacek, Z.; Nocera, L.

1991-05-01

The Vlasov-Poisson system of equations in the Fourier-transformed velocity space is studied. At first some results of the linear theory are reformulated: in the new representation the Van Kampen eigenmodes and their adjoint are found to be ordinary functions with convenient piece-wise continuity properties. A transparent derivation is given of the free-streaming temporal echo in terms of the kinematics of wave packets in the Fourier-transformed velocity space. This analysis is further extended to include Coulomb interactions which allows to establish a connection between the echo theory, the second order oscillations of Best and the phenomenon of linear sidebands. The calculation of the time evolution of the global second order electric field is performed in detail in the case of a Maxwellian equilibrium distribution function. It is concluded that the phenomenon of linear sidebands may be properly explained in terms of the intrinsic features of the equilibrium distribution function. (author) 5 figs., 32 refs

On oscillations of solutions to second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2013-01-01

Roč. 20, č. 1 (2013), s. 65-94 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillatory solution Subject RIV: BA - General Mathematics Impact factor: 0.340, year: 2013 http://www.degruyter.com/view/j/gmj.2013.20.issue-1/gmj-2013-0001/gmj-2013-0001.xml?format=INT

On oscillations of solutions to second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2013-01-01

Roč. 20, č. 1 (2013), s. 65-94 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillatory solution Subject RIV: BA - General Mathematics Impact factor: 0.340, year: 2013 http://www.degruyter.com/view/j/gmj.2013.20.issue-1/gmj-2013-0001/gmj-2013-0001. xml ?format=INT

Ordinary differential equations principles and applications

CERN Document Server

Nandakumaran, A K; George, Raju K

2017-01-01

Written in a clear, logical and concise manner, this comprehensive resource allows students to quickly understand the key principles, techniques and applications of ordinary differential equations. Important topics including first and second order linear equations, initial value problems and qualitative theory are presented in separate chapters. The concepts of two point boundary value problems, physical models and first order partial differential equations are discussed in detail. The text uses tools of calculus and real analysis to get solutions in explicit form. While discussing first order linear systems, linear algebra techniques are used. The real-life applications are interspersed throughout the book to invoke reader's interest. The methods and tricks to solve numerous mathematical problems with sufficient derivations and explanation are provided. The proofs of theorems are explained for the benefit of the readers.

Bounded solutions of self-adjoint second order linear difference equations with periodic coeffients

Directory of Open Access Journals (Sweden)

Encinas A.M.

2018-02-01

Full Text Available In this work we obtain easy characterizations for the boundedness of the solutions of the discrete, self–adjoint, second order and linear unidimensional equations with periodic coefficients, including the analysis of the so-called discrete Mathieu equations as particular cases.

Hyers-Ulam stability for second-order linear differential equations with boundary conditions

Directory of Open Access Journals (Sweden)

Pasc Gavruta

2011-06-01

Full Text Available We prove the Hyers-Ulam stability of linear differential equations of second-order with boundary conditions or with initial conditions. That is, if y is an approximate solution of the differential equation $y''+ eta (x y = 0$ with $y(a = y(b =0$, then there exists an exact solution of the differential equation, near y.

Maximum principles for boundary-degenerate second-order linear elliptic differential operators

OpenAIRE

Feehan, Paul M. N.

2012-01-01

We prove weak and strong maximum principles, including a Hopf lemma, for smooth subsolutions to equations defined by linear, second-order, partial differential operators whose principal symbols vanish along a portion of the domain boundary. The boundary regularity property of the smooth subsolutions along this boundary vanishing locus ensures that these maximum principles hold irrespective of the sign of the Fichera function. Boundary conditions need only be prescribed on the complement in th...

A Fresh Look at Linear Ordinary Differential Equations with Constant Coefficients. Revisiting the Impulsive Response Method Using Factorization

Science.gov (United States)

Camporesi, Roberto

2016-01-01

We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations of any order based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as…

Solving Ordinary Differential Equations

Science.gov (United States)

Krogh, F. T.

1987-01-01

Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Growth and Zeros of Meromorphic Solutions to Second-Order Linear Differential Equations

Directory of Open Access Journals (Sweden)

Maamar Andasmas

2016-04-01

Full Text Available The main purpose of this article is to investigate the growth of meromorphic solutions to homogeneous and non-homogeneous second order linear differential equations f00+Af0+Bf = F, where A(z, B (z and F (z are meromorphic functions with finite order having only finitely many poles. We show that, if there exist a positive constants σ > 0, α > 0 such that |A(z| ≥ eα|z|σ as |z| → +∞, z ∈ H, where dens{|z| : z ∈ H} > 0 and ρ = max{ρ(B, ρ(F} < σ, then every transcendental meromorphic solution f has an infinite order. Further, we give some estimates of their hyper-order, exponent and hyper-exponent of convergence of distinct zeros.

Internal crisis in a second-order non-linear non-autonomous electronic oscillator

International Nuclear Information System (INIS)

Stavrinides, S.G.; Deliolanis, N.C.; Miliou, A.N.; Laopoulos, Th.; Anagnostopoulos, A.N.

2008-01-01

The internal crisis of a second-order non-linear non-autonomous chaotic electronic circuit is studied. The phase portraits consist of two interacting sub-attractors, a chaotic and a periodic one. Maximal Lyapunov exponents were calculated, for both the periodic and the chaotic waveforms, in order to confirm their nature. Transitions between the chaotic and the periodic sub-attractors become more frequent by increasing the circuit driving frequency. The frequency distribution of the corresponding laminar lengths and their average values indicate that an internal crisis takes place in this circuit, manifested in the intermittent behaviour of the corresponding orbits

Hyers-Ulam stability of linear second-order differential equations in complex Banach spaces

Directory of Open Access Journals (Sweden)

Yongjin Li

2013-08-01

Full Text Available We prove the Hyers-Ulam stability of linear second-order differential equations in complex Banach spaces. That is, if y is an approximate solution of the differential equation $y''+ alpha y'(t +eta y = 0$ or $y''+ alpha y'(t +eta y = f(t$, then there exists an exact solution of the differential equation near to y.

Solution of second order linear fuzzy difference equation by Lagrange's multiplier method

Directory of Open Access Journals (Sweden)

Sankar Prasad Mondal

2016-06-01

Full Text Available In this paper we execute the solution procedure for second order linear fuzzy difference equation by Lagrange's multiplier method. In crisp sense the difference equation are easy to solve, but when we take in fuzzy sense it forms a system of difference equation which is not so easy to solve. By the help of Lagrange's multiplier we can solved it easily. The results are illustrated by two different numerical examples and followed by two applications.

OSCILLATION OF A SECOND-ORDER HALF-LINEAR NEUTRAL DAMPED DIFFERENTIAL EQUATION WITH TIME-DELAY

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

In this paper,the oscillation for a class of second-order half-linear neutral damped differential equation with time-delay is studied.By means of Yang-inequality,the generalized Riccati transformation and a certain function,some new sufficient conditions for the oscillation are given for all solutions to the equation.

Linear Matrix Inequalities for Analysis and Control of Linear Vector Second-Order Systems

DEFF Research Database (Denmark)

Adegas, Fabiano Daher; Stoustrup, Jakob

2015-01-01

the Lyapunov matrix and the system matrices by introducing matrix multipliers, which potentially reduce conservativeness in hard control problems. Multipliers facilitate the usage of parameter-dependent Lyapunov functions as certificates of stability of uncertain and time-varying vector second-order systems......SUMMARY Many dynamical systems are modeled as vector second-order differential equations. This paper presents analysis and synthesis conditions in terms of LMI with explicit dependence in the coefficient matrices of vector second-order systems. These conditions benefit from the separation between....... The conditions introduced in this work have the potential to increase the practice of analyzing and controlling systems directly in vector second-order form. Copyright © 2014 John Wiley & Sons, Ltd....

The mixed BVP for second order nonlinear ordinary differential equation at resonance

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan

2017-01-01

Roč. 290, 2-3 (2017), s. 393-400 ISSN 0025-584X Institutional support: RVO:67985840 Keywords : mixed problem at resonance * nonlinear ordinary differencial equation Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.742, year: 2016

An accurate scheme by block method for third order ordinary ...

African Journals Online (AJOL)

problems of ordinary differential equations is presented in this paper. The approach of collocation approximation is adopted in the derivation of the scheme and then the scheme is applied as simultaneous integrator to special third order initial value problem of ordinary differential equations. This implementation strategy is ...

A fresh look at linear ordinary differential equations with constant coefficients. Revisiting the impulsive response method using factorization

Science.gov (United States)

Camporesi, Roberto

2016-01-01

We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations of any order based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as well as of the other more advanced approaches: Laplace transform, linear systems, the general theory of linear equations with variable coefficients and variation of parameters. The approach presented here can be used in a first course on differential equations for science and engineering majors.

Non-monotone positive solutions of second-order linear differential equations: existence, nonexistence and criteria

Directory of Open Access Journals (Sweden)

Mervan Pašić

2016-10-01

Full Text Available We study non-monotone positive solutions of the second-order linear differential equations: $(p(tx'' + q(t x = e(t$, with positive $p(t$ and $q(t$. For the first time, some criteria as well as the existence and nonexistence of non-monotone positive solutions are proved in the framework of some properties of solutions $\\theta (t$ of the corresponding integrable linear equation: $(p(t\\theta''=e(t$. The main results are illustrated by many examples dealing with equations which allow exact non-monotone positive solutions not necessarily periodic. Finally, we pose some open questions.

Comparison of second and third orders Runge-Kutta methods for ...

African Journals Online (AJOL)

This work is concerned with the analysis of second and third orders Runge- Kutta formulae capable of solving initial value problems in Ordinary Differential Equations of the form: y1 = f(x, y), y(x0) = y0, a £ x £ b. The intention is to find out which of these two orders can improve the performance of results when implemented ...

Passive impedance-based second-order sliding mode control for non-linear teleoperators

Directory of Open Access Journals (Sweden)

Luis G García-Valdovinos

2017-02-01

Full Text Available Bilateral teleoperation systems have attracted significant attention in the last decade mainly because of technological advancements in both the communication channel and computers performance. In addition, non-linear multi-degree-of-freedom bilateral teleoperators along with state observers have become an open research area. In this article, a model-free exact differentiator is used to estimate the full state along with a chattering-free second-order sliding mode controller to guarantee a robust impedance tracking under both constant and an unknown time delay of non-linear multi-degree-of-freedom robots. The robustness of the proposed controller is improved by introducing a change of coordinates in terms of a new nominal reference similar to that used in adaptive control theory. Experimental results that validate the predicted behaviour are presented and discussed using a Phantom Premium 1.0 as the master robot and a Catalyst-5 virtual model as the slave robot. The dynamics of the Catalyst-5 system is solved online.

Binocular Combination of Second-Order Stimuli

Science.gov (United States)

Zhou, Jiawei; Liu, Rong; Zhou, Yifeng; Hess, Robert F.

2014-01-01

Phase information is a fundamental aspect of visual stimuli. However, the nature of the binocular combination of stimuli defined by modulations in contrast, so-called second-order stimuli, is presently not clear. To address this issue, we measured binocular combination for first- (luminance modulated) and second-order (contrast modulated) stimuli using a binocular phase combination paradigm in seven normal adults. We found that the binocular perceived phase of second-order gratings depends on the interocular signal ratio as has been previously shown for their first order counterparts; the interocular signal ratios when the two eyes were balanced was close to 1 in both first- and second-order phase combinations. However, second-order combination is more linear than previously found for first-order combination. Furthermore, binocular combination of second-order stimuli was similar regardless of whether the carriers in the two eyes were correlated, anti-correlated, or uncorrelated. This suggests that, in normal adults, the binocular phase combination of second-order stimuli occurs after the monocular extracting of the second-order modulations. The sensory balance associated with this second-order combination can be obtained from binocular phase combination measurements. PMID:24404180

First and second order vortex dynamics

International Nuclear Information System (INIS)

Kim, Yoonbai; Lee, Kimyeong

2002-01-01

The low energy dynamics of vortices in self-dual Abelian Higgs theory in (2+1)-dimensional spacetime is of second order in vortex velocity and characterized by the moduli space metric. When the Chern-Simons term with a small coefficient is added to the theory, we show that a term linear in vortex velocity appears and can be consistently added to the second order expression. We provide an additional check of the first and second order terms by studying the angular momentum in field theory

Importance of the alignment of polar π conjugated molecules inside carbon nanotubes in determining second-order non-linear optical properties.

Science.gov (United States)

Yumura, Takashi; Yamamoto, Wataru

2017-09-20

We employed density functional theory (DFT) calculations with dispersion corrections to investigate energetically preferred alignments of certain p,p'-dimethylaminonitrostilbene (DANS) molecules inside an armchair (m,m) carbon nanotube (n × DANS@(m,m)), where the number of inner molecules (n) is no greater than 3. Here, three types of alignments of DANS are considered: a linear alignment in a parallel fashion and stacking alignments in parallel and antiparallel fashions. According to DFT calculations, a threshold tube diameter for containing DANS molecules in linear or stacking alignments was found to be approximately 1.0 nm. Nanotubes with diameters smaller than 1.0 nm result in the selective formation of linearly aligned DANS molecules due to strong confinement effects within the nanotubes. By contrast, larger diameter nanotubes allow DANS molecules to align in a stacking and linear fashion. The type of alignment adopted by the DANS molecules inside a nanotube is responsible for their second-order non-linear optical properties represented by their static hyperpolarizability (β 0 values). In fact, we computed β 0 values of DANS assemblies taken from optimized n × DANS@(m,m) structures, and their values were compared with those of a single DANS molecule. DFT calculations showed that β 0 values of DANS molecules depend on their alignment, which decrease in the following order: linear alignment > parallel stacking alignment > antiparallel stacking alignment. In particular, a linear alignment has a β 0 value more significant than that of the same number of isolated molecules. Therefore, the linear alignment of DANS molecules, which is only allowed inside smaller diameter nanotubes, can strongly enhance their second-order non-linear optical properties. Since the nanotube confinement determines the alignment of DANS molecules, a restricted nanospace can be utilized to control their second-order non-linear optical properties. These DFT findings can assist in the

Asymptotic behavior and stability of second order neutral delay differential equations

NARCIS (Netherlands)

Chen, G.L.; van Gaans, O.W.; Verduyn Lunel, Sjoerd

2014-01-01

We study the asymptotic behavior of a class of second order neutral delay differential equations by both a spectral projection method and an ordinary differential equation method approach. We discuss the relation of these two methods and illustrate some features using examples. Furthermore, a fixed

Composition between mecd and runge-Kutta algorithm method for large system of second order differential equations

International Nuclear Information System (INIS)

Supriyono; Miyoshi, T.

1997-01-01

NECD Method and runge-Kutta method for large system of second order ordinary differential equations in comparing algorithm. The paper introduce a extrapolation method used for solving the large system of second order ordinary differential equation. We call this method the modified extrapolated central difference (MECD) method. for the accuracy and efficiency MECD method. we compare the method with 4-th order runge-Kutta method. The comparison results show that, this method has almost the same accuracy as the 4-th order runge-Kutta method, but the computation time is about half of runge-Kutta. The MECD was declare by the author and Tetsuhiko Miyoshi of the Dept. Applied Science Yamaguchi University Japan

Linearization: Geometric, Complex, and Conditional

Directory of Open Access Journals (Sweden)

Asghar Qadir

2012-01-01

Full Text Available Lie symmetry analysis provides a systematic method of obtaining exact solutions of nonlinear (systems of differential equations, whether partial or ordinary. Of special interest is the procedure that Lie developed to transform scalar nonlinear second-order ordinary differential equations to linear form. Not much work was done in this direction to start with, but recently there have been various developments. Here, first the original work of Lie (and the early developments on it, and then more recent developments based on geometry and complex analysis, apart from Lie’s own method of algebra (namely, Lie group theory, are reviewed. It is relevant to mention that much of the work is not linearization but uses the base of linearization.

Constructive Development of the Solutions of Linear Equations in Introductory Ordinary Differential Equations

Science.gov (United States)

Mallet, D. G.; McCue, S. W.

2009-01-01

The solution of linear ordinary differential equations (ODEs) is commonly taught in first-year undergraduate mathematics classrooms, but the understanding of the concept of a solution is not always grasped by students until much later. Recognizing what it is to be a solution of a linear ODE and how to postulate such solutions, without resorting to…

Path integral solution of linear second order partial differential equations I: the general construction

International Nuclear Information System (INIS)

LaChapelle, J.

2004-01-01

A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schroedinger/diffusion equation in unbounded space. The construction is based on a framework for functional integration introduced by Cartier/DeWitt-Morette

A New Factorisation of a General Second Order Differential Equation

Science.gov (United States)

Clegg, Janet

2006-01-01

A factorisation of a general second order ordinary differential equation is introduced from which the full solution to the equation can be obtained by performing two integrations. The method is compared with traditional methods for solving these type of equations. It is shown how the Green's function can be derived directly from the factorisation…

Relative boundedness and compactness theory for second-order differential operators

Directory of Open Access Journals (Sweden)

Don B. Hinton

1997-01-01

Full Text Available The problem considered is to give necessary and sufficient conditions for perturbations of a second-order ordinary differential operator to be either relatively bounded or relatively compact. Such conditions are found for three classes of operators. The conditions are expressed in terms of integral averages of the coefficients of the perturbing operator.

LSODE, 1. Order Stiff or Non-Stiff Ordinary Differential Equations System Initial Value Problems

International Nuclear Information System (INIS)

Hindmarsh, A.C.; Petzold, L.R.

2005-01-01

1 - Description of program or function: LSODE (Livermore Solver for Ordinary Differential Equations) solves stiff and non-stiff systems of the form dy/dt = f. In the stiff case, it treats the Jacobian matrix df/dy as either a dense (full) or a banded matrix, and as either user-supplied or internally approximated by difference quotients. It uses Adams methods (predictor-corrector) in the non-stiff case, and Backward Differentiation Formula (BDF) methods (the Gear methods) in the stiff case. The linear systems that arise are solved by direct methods (LU factor/solve). The LSODE source is commented extensively to facilitate modification. Both a single-precision version and a double-precision version are available. 2 - Methods: It is assumed that the ODEs are given explicitly, so that the system can be written in the form dy/dt = f(t,y), where y is the vector of dependent variables, and t is the independent variable. LSODE contains two variable-order, variable- step (with interpolatory step-changing) integration methods. The first is the implicit Adams or non-stiff method, of orders one through twelve. The second is the backward differentiation or stiff method (or BDF method, or Gear's method), of orders one through five. 3 - Restrictions on the complexity of the problem: The differential equations must be given in explicit form, i.e., dy/dt = f(y,t). Problems with intermittent high-speed transients may cause inefficient or unstable performance

On periodic bounded and unbounded solutions of second order nonlinear ordinary differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander

2017-01-01

Roč. 24, č. 2 (2017), s. 241-263 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : nonlinear ordinary differential equations * periodic boundary value problem * solvability Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2017-0009/gmj-2017-0009. xml

On periodic bounded and unbounded solutions of second order nonlinear ordinary differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander

2017-01-01

Roč. 24, č. 2 (2017), s. 241-263 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : nonlinear ordinary differential equations * periodic boundary value problem * solvability Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2017-0009/gmj-2017-0009.xml

Algorithmic Verification of Linearizability for Ordinary Differential Equations

KAUST Repository

Lyakhov, Dmitry A.

2017-07-19

For a nonlinear ordinary differential equation solved with respect to the highest order derivative and rational in the other derivatives and in the independent variable, we devise two algorithms to check if the equation can be reduced to a linear one by a point transformation of the dependent and independent variables. The first algorithm is based on a construction of the Lie point symmetry algebra and on the computation of its derived algebra. The second algorithm exploits the differential Thomas decomposition and allows not only to test the linearizability, but also to generate a system of nonlinear partial differential equations that determines the point transformation and the coefficients of the linearized equation. The implementation of both algorithms is discussed and their application is illustrated using several examples.

Parallels between control PDE's (Partial Differential Equations) and systems of ODE's (Ordinary Differential Equations)

Science.gov (United States)

Hunt, L. R.; Villarreal, Ramiro

1987-01-01

System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.

Second-order processing of four-stroke apparent motion.

Science.gov (United States)

Mather, G; Murdoch, L

1999-05-01

In four-stroke apparent motion displays, pattern elements oscillate between two adjacent positions and synchronously reverse in contrast, but appear to move unidirectionally. For example, if rightward shifts preserve contrast but leftward shifts reverse contrast, consistent rightward motion is seen. In conventional first-order displays, elements reverse in luminance contrast (e.g. light elements become dark, and vice-versa). The resulting perception can be explained by responses in elementary motion detectors turned to spatio-temporal orientation. Second-order motion displays contain texture-defined elements, and there is some evidence that they excite second-order motion detectors that extract spatio-temporal orientation following the application of a non-linear 'texture-grabbing' transform by the visual system. We generated a variety of second-order four-stroke displays, containing texture-contrast reversals instead of luminance contrast reversals, and used their effectiveness as a diagnostic test for the presence of various forms of non-linear transform in the second-order motion system. Displays containing only forward or only reversed phi motion sequences were also tested. Displays defined by variation in luminance, contrast, orientation, and size were effective. Displays defined by variation in motion, dynamism, and stereo were partially or wholly ineffective. Results obtained with contrast-reversing and four-stroke displays indicate that only relatively simple non-linear transforms (involving spatial filtering and rectification) are available during second-order energy-based motion analysis.

Linear Ordinary Differential Equations with Constant Coefficients. Revisiting the Impulsive Response Method Using Factorization

Science.gov (United States)

Camporesi, Roberto

2011-01-01

We present an approach to the impulsive response method for solving linear constant-coefficient ordinary differential equations based on the factorization of the differential operator. The approach is elementary, we only assume a basic knowledge of calculus and linear algebra. In particular, we avoid the use of distribution theory, as well as of…

Hojman's theorem of the third-order ordinary differential equation

International Nuclear Information System (INIS)

Hong-Sheng, Lü; Hong-Bin, Zhang; Shu-Long, Gu

2009-01-01

This paper extends Hojman's conservation law to the third-order differential equation. A new conserved quantity is constructed based on the Lie group of transformation generators of the equations of motion. The generators contain variations of the time and generalized coordinates. Two independent non-trivial conserved quantities of the third-order ordinary differential equation are obtained. A simple example is presented to illustrate the applications of the results. (general)

Second-order numerical methods for multi-term fractional differential equations: Smooth and non-smooth solutions

Science.gov (United States)

Zeng, Fanhai; Zhang, Zhongqiang; Karniadakis, George Em

2017-12-01

Starting with the asymptotic expansion of the error equation of the shifted Gr\\"{u}nwald--Letnikov formula, we derive a new modified weighted shifted Gr\\"{u}nwald--Letnikov (WSGL) formula by introducing appropriate correction terms. We then apply one special case of the modified WSGL formula to solve multi-term fractional ordinary and partial differential equations, and we prove the linear stability and second-order convergence for both smooth and non-smooth solutions. We show theoretically and numerically that numerical solutions up to certain accuracy can be obtained with only a few correction terms. Moreover, the correction terms can be tuned according to the fractional derivative orders without explicitly knowing the analytical solutions. Numerical simulations verify the theoretical results and demonstrate that the new formula leads to better performance compared to other known numerical approximations with similar resolution.

Remarks on second-order quadratic systems in algebras

Directory of Open Access Journals (Sweden)

Art Sagle

2017-10-01

Full Text Available This paper is an addendum to our earlier paper [8], where a systematic study of quadratic systems of second order ordinary differential equations defined in commutative algebras was presented. Here we concentrate on special solutions and energy considerations of some quadratic systems defined in algebras which need not be commutative, however, we shall throughout assume the algebra to be associative. We here also give a positive answer to an open question, concerning periodic motions of such systems, posed in our earlier paper.

Variations in the Solution of Linear First-Order Differential Equations. Classroom Notes

Science.gov (United States)

Seaman, Brian; Osler, Thomas J.

2004-01-01

A special project which can be given to students of ordinary differential equations is described in detail. Students create new differential equations by changing the dependent variable in the familiar linear first-order equation (dv/dx)+p(x)v=q(x) by means of a substitution v=f(y). The student then creates a table of the new equations and…

Algebraic Properties of First Integrals for Scalar Linear Third-Order ODEs of Maximal Symmetry

Directory of Open Access Journals (Sweden)

K. S. Mahomed

2013-01-01

Full Text Available By use of the Lie symmetry group methods we analyze the relationship between the first integrals of the simplest linear third-order ordinary differential equations (ODEs and their point symmetries. It is well known that there are three classes of linear third-order ODEs for maximal cases of point symmetries which are 4, 5, and 7. The simplest scalar linear third-order equation has seven-point symmetries. We obtain the classifying relation between the symmetry and the first integral for the simplest equation. It is shown that the maximal Lie algebra of a first integral for the simplest equation y′′′=0 is unique and four-dimensional. Moreover, we show that the Lie algebra of the simplest linear third-order equation is generated by the symmetries of the two basic integrals. We also obtain counting theorems of the symmetry properties of the first integrals for such linear third-order ODEs. Furthermore, we provide insights into the manner in which one can generate the full Lie algebra of higher-order ODEs of maximal symmetry from two of their basic integrals.

Hypergeometric Series Solution to a Class of Second-Order Boundary Value Problems via Laplace Transform with Applications to Nanofluids

Science.gov (United States)

Ebaid, Abdelhalim; Wazwaz, Abdul-Majid; Alali, Elham; Masaedeh, Basem S.

2017-03-01

Very recently, it was observed that the temperature of nanofluids is finally governed by second-order ordinary differential equations with variable coefficients of exponential orders. Such coefficients were then transformed to polynomials type by using new independent variables. In this paper, a class of second-order ordinary differential equations with variable coefficients of polynomials type has been solved analytically. The analytical solution is expressed in terms of a hypergeometric function with generalized parameters. Moreover, applications of the present results have been applied on some selected nanofluids problems in the literature. The exact solutions in the literature were derived as special cases of our generalized analytical solution.

Finding higher order Darboux polynomials for a family of rational first order ordinary differential equations

Science.gov (United States)

Avellar, J.; Claudino, A. L. G. C.; Duarte, L. G. S.; da Mota, L. A. C. P.

2015-10-01

For the Darbouxian methods we are studying here, in order to solve first order rational ordinary differential equations (1ODEs), the most costly (computationally) step is the finding of the needed Darboux polynomials. This can be so grave that it can render the whole approach unpractical. Hereby we introduce a simple heuristics to speed up this process for a class of 1ODEs.

An algorithm for solving initial value problems of third order ordinary ...

African Journals Online (AJOL)

Abstract. We propose an implicit multi-step method for the solution of initial value problems (IVPs) of third order ordinary differential equations (ODE) which does not require reducing the ODE to first order before solving. The development of the method is based on collocation of the differential system and interpolation of the ...

A first course in ordinary differential equations analytical and numerical methods

CERN Document Server

Hermann, Martin

2014-01-01

This book presents a modern introduction to analytical and numerical techniques for solving ordinary differential equations (ODEs). Contrary to the traditional format—the theorem-and-proof format—the book is focusing on analytical and numerical methods. The book supplies a variety of problems and examples, ranging from the elementary to the advanced level, to introduce and study the mathematics of ODEs. The analytical part of the book deals with solution techniques for scalar first-order and second-order linear ODEs, and systems of linear ODEs—with a special focus on the Laplace transform, operator techniques and power series solutions. In the numerical part, theoretical and practical aspects of Runge-Kutta methods for solving initial-value problems and shooting methods for linear two-point boundary-value problems are considered. The book is intended as a primary text for courses on the theory of ODEs and numerical treatment of ODEs for advanced undergraduate and early graduate students. It is assumed t...

First order linear ordinary differential equations in associative algebras

Directory of Open Access Journals (Sweden)

Gordon Erlebacher

2004-01-01

Full Text Available In this paper, we study the linear differential equation $$ frac{dx}{dt}=sum_{i=1}^n a_i(t x b_i(t + f(t $$ in an associative but non-commutative algebra $mathcal{A}$, where the $b_i(t$ form a set of commuting $mathcal{A}$-valued functions expressed in a time-independent spectral basis consisting of mutually annihilating idempotents and nilpotents. Explicit new closed solutions are derived, and examples are presented to illustrate the theory.

On the multisummability of WKB solutions of certain singularly perturbed linear ordinary differential equations

Directory of Open Access Journals (Sweden)

Yoshitsugu Takei

2015-01-01

Full Text Available Using two concrete examples, we discuss the multisummability of WKB solutions of singularly perturbed linear ordinary differential equations. Integral representations of solutions and a criterion for the multisummability based on the Cauchy-Heine transform play an important role in the proof.

Symmetries of th-Order Approximate Stochastic Ordinary Differential Equations

OpenAIRE

Fredericks, E.; Mahomed, F. M.

2012-01-01

Symmetries of $n$ th-order approximate stochastic ordinary differential equations (SODEs) are studied. The determining equations of these SODEs are derived in an Itô calculus context. These determining equations are not stochastic in nature. SODEs are normally used to model nature (e.g., earthquakes) or for testing the safety and reliability of models in construction engineering when looking at the impact of random perturbations.

Rethinking pedagogy for second-order differential equations: a simplified approach to understanding well-posed problems

Science.gov (United States)

Tisdell, Christopher C.

2017-07-01

Knowing an equation has a unique solution is important from both a modelling and theoretical point of view. For over 70 years, the approach to learning and teaching 'well posedness' of initial value problems (IVPs) for second- and higher-order ordinary differential equations has involved transforming the problem and its analysis to a first-order system of equations. We show that this excursion is unnecessary and present a direct approach regarding second- and higher-order problems that does not require an understanding of systems.

A Four-Stage Fifth-Order Trigonometrically Fitted Semi-Implicit Hybrid Method for Solving Second-Order Delay Differential Equations

Directory of Open Access Journals (Sweden)

Sufia Zulfa Ahmad

2016-01-01

Full Text Available We derived a two-step, four-stage, and fifth-order semi-implicit hybrid method which can be used for solving special second-order ordinary differential equations. The method is then trigonometrically fitted so that it is suitable for solving problems which are oscillatory in nature. The methods are then used for solving oscillatory delay differential equations. Numerical results clearly show the efficiency of the new method when compared to the existing explicit and implicit methods in the scientific literature.

On the boundedness and integration of non-oscillatory solutions of certain linear differential equations of second order.

Science.gov (United States)

Tunç, Cemil; Tunç, Osman

2016-01-01

In this paper, certain system of linear homogeneous differential equations of second-order is considered. By using integral inequalities, some new criteria for bounded and [Formula: see text]-solutions, upper bounds for values of improper integrals of the solutions and their derivatives are established to the considered system. The obtained results in this paper are considered as extension to the results obtained by Kroopnick (2014) [1]. An example is given to illustrate the obtained results.

A block Krylov subspace time-exact solution method for linear ordinary differential equation systems

NARCIS (Netherlands)

Bochev, Mikhail A.

2013-01-01

We propose a time-exact Krylov-subspace-based method for solving linear ordinary differential equation systems of the form $y'=-Ay+g(t)$ and $y"=-Ay+g(t)$, where $y(t)$ is the unknown function. The method consists of two stages. The first stage is an accurate piecewise polynomial approximation of

On Sequences of Numbers and Polynomials Defined by Linear Recurrence Relations of Order 2

Directory of Open Access Journals (Sweden)

Tian-Xiao He

2009-01-01

Full Text Available Here we present a new method to construct the explicit formula of a sequence of numbers and polynomials generated by a linear recurrence relation of order 2. The applications of the method to the Fibonacci and Lucas numbers, Chebyshev polynomials, the generalized Gegenbauer-Humbert polynomials are also discussed. The derived idea provides a general method to construct identities of number or polynomial sequences defined by linear recurrence relations. The applications using the method to solve some algebraic and ordinary differential equations are presented.

Nontrivial Periodic Solutions for Nonlinear Second-Order Difference Equations

Directory of Open Access Journals (Sweden)

Tieshan He

2011-01-01

Full Text Available This paper is concerned with the existence of nontrivial periodic solutions and positive periodic solutions to a nonlinear second-order difference equation. Under some conditions concerning the first positive eigenvalue of the linear equation corresponding to the nonlinear second-order equation, we establish the existence results by using the topological degree and fixed point index theories.

Oscillation of two-dimensional linear second-order differential systems

International Nuclear Information System (INIS)

Kwong, M.K.; Kaper, H.G.

1985-01-01

This article is concerned with the oscillatory behavior at infinity of the solution y: [a, ∞) → R 2 of a system of two second-order differential equations, y''(t) + Q(t) y(t) = 0, t epsilon[a, ∞); Q is a continuous matrix-valued function on [a, ∞) whose values are real symmetric matrices of order 2. It is shown that the solution is oscillatory at infinity if the largest eigenvalue of the matrix integral/sub a//sup t/ Q(s) ds tends to infinity as t → ∞. This proves a conjecture of D. Hinton and R.T. Lewis for the two-dimensional case. Furthermore, it is shown that considerably weaker forms of the condition still suffice for oscillatory behavior at infinity. 7 references

Oscillation theory for second order dynamic equations

CERN Document Server

Agarwal, Ravi P; O''Regan, Donal

2003-01-01

The qualitative theory of dynamic equations is a rapidly developing area of research. In the last 50 years, the Oscillation Theory of ordinary, functional, neutral, partial and impulsive differential equations, and their discrete versions, has inspired many scholars. Hundreds of research papers have been published in every major mathematical journal. Many books deal exclusively with the oscillation of solutions of differential equations, but most of these books appeal only to researchers who already know the subject. In an effort to bring Oscillation Theory to a new and broader audience, the authors present a compact, but thorough, understanding of Oscillation Theory for second order differential equations. They include several examples throughout the text not only to illustrate the theory, but also to provide new direction.

A second-order virtual node algorithm for nearly incompressible linear elasticity in irregular domains

Science.gov (United States)

Zhu, Yongning; Wang, Yuting; Hellrung, Jeffrey; Cantarero, Alejandro; Sifakis, Eftychios; Teran, Joseph M.

2012-08-01

We present a cut cell method in R2 for enforcing Dirichlet and Neumann boundary conditions with nearly incompressible linear elastic materials in irregular domains. Virtual nodes on cut uniform grid cells are used to provide geometric flexibility in the domain boundary shape without sacrificing accuracy. We use a mixed formulation utilizing a MAC-type staggered grid with piecewise bilinear displacements centered at cell faces and piecewise constant pressures at cell centers. These discretization choices provide the necessary stability in the incompressible limit and the necessary accuracy in cut cells. Numerical experiments suggest second order accuracy in L∞. We target high-resolution problems and present a class of geometric multigrid methods for solving the discrete equations for displacements and pressures that achieves nearly optimal convergence rates independent of grid resolution.

SIVA/DIVA- INITIAL VALUE ORDINARY DIFFERENTIAL EQUATION SOLUTION VIA A VARIABLE ORDER ADAMS METHOD

Science.gov (United States)

Krogh, F. T.

1994-01-01

The SIVA/DIVA package is a collection of subroutines for the solution of ordinary differential equations. There are versions for single precision and double precision arithmetic. These solutions are applicable to stiff or nonstiff differential equations of first or second order. SIVA/DIVA requires fewer evaluations of derivatives than other variable order Adams predictor-corrector methods. There is an option for the direct integration of second order equations which can make integration of trajectory problems significantly more efficient. Other capabilities of SIVA/DIVA include: monitoring a user supplied function which can be separate from the derivative; dynamically controlling the step size; displaying or not displaying output at initial, final, and step size change points; saving the estimated local error; and reverse communication where subroutines return to the user for output or computation of derivatives instead of automatically performing calculations. The user must supply SIVA/DIVA with: 1) the number of equations; 2) initial values for the dependent and independent variables, integration stepsize, error tolerance, etc.; and 3) the driver program and operational parameters necessary for subroutine execution. SIVA/DIVA contains an extensive diagnostic message library should errors occur during execution. SIVA/DIVA is written in FORTRAN 77 for batch execution and is machine independent. It has a central memory requirement of approximately 120K of 8 bit bytes. This program was developed in 1983 and last updated in 1987.

Optimal explicit strong stability preserving Runge–Kutta methods with high linear order and optimal nonlinear order

KAUST Repository

Gottlieb, Sigal

2015-04-10

High order spatial discretizations with monotonicity properties are often desirable for the solution of hyperbolic PDEs. These methods can advantageously be coupled with high order strong stability preserving time discretizations. The search for high order strong stability time-stepping methods with large allowable strong stability coefficient has been an active area of research over the last two decades. This research has shown that explicit SSP Runge-Kutta methods exist only up to fourth order. However, if we restrict ourselves to solving only linear autonomous problems, the order conditions simplify and this order barrier is lifted: explicit SSP Runge-Kutta methods of any linear order exist. These methods reduce to second order when applied to nonlinear problems. In the current work we aim to find explicit SSP Runge-Kutta methods with large allowable time-step, that feature high linear order and simultaneously have the optimal fourth order nonlinear order. These methods have strong stability coefficients that approach those of the linear methods as the number of stages and the linear order is increased. This work shows that when a high linear order method is desired, it may still be worthwhile to use methods with higher nonlinear order.

Rethinking Pedagogy for Second-Order Differential Equations: A Simplified Approach to Understanding Well-Posed Problems

Science.gov (United States)

Tisdell, Christopher C.

2017-01-01

Knowing an equation has a unique solution is important from both a modelling and theoretical point of view. For over 70 years, the approach to learning and teaching "well posedness" of initial value problems (IVPs) for second- and higher-order ordinary differential equations has involved transforming the problem and its analysis to a…

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

Science.gov (United States)

Park, K. C.; Belvin, W. Keith

1990-01-01

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Second-order splitting schemes for a class of reactive systems

International Nuclear Information System (INIS)

Ren Zhuyin; Pope, Stephen B.

2008-01-01

We consider the numerical time integration of a class of reaction-transport systems that are described by a set of ordinary differential equations for primary variables. In the governing equations, the terms involved may require the knowledge of secondary variables, which are functions of the primary variables. Specifically, we consider the case where, given the primary variables, the evaluation of the secondary variables is computationally expensive. To solve this class of reaction-transport equations, we develop and demonstrate several computationally efficient splitting schemes, wherein the portions of the governing equations containing chemical reaction terms are separated from those parts containing the transport terms. A computationally efficient solution to the transport sub-step is achieved through the use of linearization or predictor-corrector methods. The splitting schemes are applied to the reactive flow in a continuously stirred tank reactor (CSTR) with the Davis-Skodjie reaction model, to the CO+H 2 oxidation in a CSTR with detailed chemical kinetics, and to a reaction-diffusion system with an extension of the Oregonator model of the Belousov-Zhabotinsky reaction. As demonstrated in the test problems, the proposed splitting schemes, which yield efficient solutions to the transport sub-step, achieve second-order accuracy in time

Symmetries, Integrals and Solutions of Ordinary Differential ...

Indian Academy of Sciences (India)

Second-and third-order scalar ordinary differential equations of maximal symmetry in the traditional sense of point, respectively contact, symmetry are examined for the mappings they produce in solutions and fundamental first integrals. The properties of the `exceptional symmetries', i.e. those not considered to be generic to ...

Implementing Families of Implicit Chebyshev Methods with Exact Coefficients for the Numerical Integration of First- and Second-Order Differential Equations

National Research Council Canada - National Science Library

Mitchell, Jason

2002-01-01

A method is presented for the generation of exact numerical coefficients found in two families of implicit Chebyshev methods for the numerical integration of first- and second-order ordinary differential equations...

Algorithm for solving the linear Cauchy problem for large systems of ordinary differential equations with the use of parallel computations

Energy Technology Data Exchange (ETDEWEB)

Moryakov, A. V., E-mail: sailor@orc.ru [National Research Centre Kurchatov Institute (Russian Federation)

2016-12-15

An algorithm for solving the linear Cauchy problem for large systems of ordinary differential equations is presented. The algorithm for systems of first-order differential equations is implemented in the EDELWEISS code with the possibility of parallel computations on supercomputers employing the MPI (Message Passing Interface) standard for the data exchange between parallel processes. The solution is represented by a series of orthogonal polynomials on the interval [0, 1]. The algorithm is characterized by simplicity and the possibility to solve nonlinear problems with a correction of the operator in accordance with the solution obtained in the previous iterative process.

A nanolens-type enhancement in the linear and second harmonic response of a metallic dimer

International Nuclear Information System (INIS)

Pustovit, Vitaliy; Biswas, Sushmita; Vaia, Richard; Urbas, Augustine

2014-01-01

In this paper we explore the linear and second-order nonlinear response of gold nanoparticle pairs (dimers). Despite that even-order nonlinear processes are forbidden in bulk centrosymmetric media like metals, second order nonlinear response exhibits a high degree of sensitivity for spherical nanoparticles where inversion symmetry is broken at the surface. Recent experiments demonstrate significant dependence of linear response and second-harmonic surface nonlinear response arising from the local fundamental field distribution in a dimer configuration. Our calculations are carried out taking into account high order multipolar interactions between metal nanoparticles, and demonstrate that linear and nonlinear optical responses of the dimer exhibit periodic behavior dependent on the separation distance between nanoparticles. This response increases for dimers with a large difference between particle sizes. (paper)

Multi-octave analog photonic link with improved second- and third-order SFDRs

Science.gov (United States)

Tan, Qinggui; Gao, Yongsheng; Fan, Yangyu; He, You

2018-03-01

The second- and third-order spurious free dynamic ranges (SFDRs) are two key performance indicators for a multi-octave analogy photonic link (APL). The linearization methods for either second- or third-order intermodulation distortion (IMD2 or IMD3) have been intensively studied, but the simultaneous suppression for the both were merely reported. In this paper, we propose an APL with improved second- and third-order SFDRs for multi-octave applications based on two parallel DPMZM-based sub-APLs. The IMD3 in each sub-APL is suppressed by properly biasing the DPMZM, and the IMD2 is suppressed by balanced detecting the two sub-APLs. The experiment demonstrates significant suppression ratios for both the IMD2 and IMD3 after linearization in the proposed link, and the measured second- and third-order SFDRs with the operating frequency from 6 to 40 GHz are above 91 dB ṡHz 1 / 2 and 116 dB ṡHz 2 / 3, respectively.

FORCED OSCILLATIONS OF SECOND ORDER SUPER-LINEAR DIFFERENTIAL EQUATION WITH IMPULSES

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

At first,by means of Kartsatos technique,we reduce the impulsive differential equation to a second order nonlinear impulsive homogeneous equation.We find some suitable impulse functions such that all the solutions to the equation are oscillatory.Several criteria on the oscillations of solutions are given.At last,we give an example to demonstrate our results.

An A Posteriori Error Analysis of Mixed Finite Element Galerkin Approximations to Second Order Linear Parabolic Problems

KAUST Repository

Memon, Sajid; Nataraj, Neela; Pani, Amiya Kumar

2012-01-01

In this article, a posteriori error estimates are derived for mixed finite element Galerkin approximations to second order linear parabolic initial and boundary value problems. Using mixed elliptic reconstructions, a posteriori error estimates in L∞(L2)- and L2(L2)-norms for the solution as well as its flux are proved for the semidiscrete scheme. Finally, based on a backward Euler method, a completely discrete scheme is analyzed and a posteriori error bounds are derived, which improves upon earlier results on a posteriori estimates of mixed finite element approximations to parabolic problems. Results of numerical experiments verifying the efficiency of the estimators have also been provided. © 2012 Society for Industrial and Applied Mathematics.

Consensus Algorithms for Networks of Systems with Second- and Higher-Order Dynamics

Science.gov (United States)

Fruhnert, Michael

This thesis considers homogeneous networks of linear systems. We consider linear feedback controllers and require that the directed graph associated with the network contains a spanning tree and systems are stabilizable. We show that, in continuous-time, consensus with a guaranteed rate of convergence can always be achieved using linear state feedback. For networks of continuous-time second-order systems, we provide a new and simple derivation of the conditions for a second-order polynomials with complex coefficients to be Hurwitz. We apply this result to obtain necessary and sufficient conditions to achieve consensus with networks whose graph Laplacian matrix may have complex eigenvalues. Based on the conditions found, methods to compute feedback gains are proposed. We show that gains can be chosen such that consensus is achieved robustly over a variety of communication structures and system dynamics. We also consider the use of static output feedback. For networks of discrete-time second-order systems, we provide a new and simple derivation of the conditions for a second-order polynomials with complex coefficients to be Schur. We apply this result to obtain necessary and sufficient conditions to achieve consensus with networks whose graph Laplacian matrix may have complex eigenvalues. We show that consensus can always be achieved for marginally stable systems and discretized systems. Simple conditions for consensus achieving controllers are obtained when the Laplacian eigenvalues are all real. For networks of continuous-time time-variant higher-order systems, we show that uniform consensus can always be achieved if systems are quadratically stabilizable. In this case, we provide a simple condition to obtain a linear feedback control. For networks of discrete-time higher-order systems, we show that constant gains can be chosen such that consensus is achieved for a variety of network topologies. First, we develop simple results for networks of time

Modeling of second order space charge driven coherent sum and difference instabilities

Directory of Open Access Journals (Sweden)

Yao-Shuo Yuan

2017-10-01

Full Text Available Second order coherent oscillation modes in intense particle beams play an important role for beam stability in linear or circular accelerators. In addition to the well-known second order even envelope modes and their instability, coupled even envelope modes and odd (skew modes have recently been shown in [Phys. Plasmas 23, 090705 (2016PHPAEN1070-664X10.1063/1.4963851] to lead to parametric instabilities in periodic focusing lattices with sufficiently different tunes. While this work was partly using the usual envelope equations, partly also particle-in-cell (PIC simulation, we revisit these modes here and show that the complete set of second order even and odd mode phenomena can be obtained in a unifying approach by using a single set of linearized rms moment equations based on “Chernin’s equations.” This has the advantage that accurate information on growth rates can be obtained and gathered in a “tune diagram.” In periodic focusing we retrieve the parametric sum instabilities of coupled even and of odd modes. The stop bands obtained from these equations are compared with results from PIC simulations for waterbag beams and found to show very good agreement. The “tilting instability” obtained in constant focusing confirms the equivalence of this method with the linearized Vlasov-Poisson system evaluated in second order.

ACCURATE ESTIMATES OF CHARACTERISTIC EXPONENTS FOR SECOND ORDER DIFFERENTIAL EQUATION

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

In this paper, a second order linear differential equation is considered, and an accurate estimate method of characteristic exponent for it is presented. Finally, we give some examples to verify the feasibility of our result.

Second-order differential-delay equation to describe a hybrid bistable device

Science.gov (United States)

Vallee, R.; Dubois, P.; Cote, M.; Delisle, C.

1987-08-01

The problem of a dynamical system with delayed feedback, a hybrid bistable device, characterized by n response times and described by an nth-order differential-delay equation (DDE) is discussed. Starting from a linear-stability analysis of the DDE, the effects of the second-order differential terms on the position of the first bifurcation and on the frequency of the resulting self-oscillation are shown. The effects of the third-order differential terms on the first bifurcation are also considered. Experimental results are shown to support the linear analysis.

On the Linear Stability of the Fifth-Order WENO Discretization

KAUST Repository

Motamed, Mohammad

2010-10-03

We study the linear stability of the fifth-order Weighted Essentially Non-Oscillatory spatial discretization (WENO5) combined with explicit time stepping applied to the one-dimensional advection equation. We show that it is not necessary for the stability domain of the time integrator to include a part of the imaginary axis. In particular, we show that the combination of WENO5 with either the forward Euler method or a two-stage, second-order Runge-Kutta method is linearly stable provided very small time step-sizes are taken. We also consider fifth-order multistep time discretizations whose stability domains do not include the imaginary axis. These are found to be linearly stable with moderate time steps when combined with WENO5. In particular, the fifth-order extrapolated BDF scheme gave superior results in practice to high-order Runge-Kutta methods whose stability domain includes the imaginary axis. Numerical tests are presented which confirm the analysis. © Springer Science+Business Media, LLC 2010.

On method of solving third-order ordinary differential equations directly using Bernstein polynomials

Science.gov (United States)

Khataybeh, S. N.; Hashim, I.

2018-04-01

In this paper, we propose for the first time a method based on Bernstein polynomials for solving directly a class of third-order ordinary differential equations (ODEs). This method gives a numerical solution by converting the equation into a system of algebraic equations which is solved directly. Some numerical examples are given to show the applicability of the method.

The solutions of second-order linear differential systems with constant delays

Science.gov (United States)

Diblík, Josef; Svoboda, Zdeněk

2017-07-01

The representations of solutions to initial problems for non-homogenous n-dimensional second-order differential equations with delays x″(t )-2 A x'(t -τ )+(A2+B2)x (t -2 τ )=f (t ) by means of special matrix delayed functions are derived. Square matrices A and B are commuting and τ > 0. Derived representations use what is called a delayed exponential of a matrix and results generalize some of known results previously derived for homogenous systems.

High Order A-stable Continuous General Linear Methods for Solution of Systems of Initial Value Problems in ODEs

Directory of Open Access Journals (Sweden)

Dauda GuliburYAKUBU

2012-12-01

Full Text Available Accurate solutions to initial value systems of ordinary differential equations may be approximated efficiently by Runge-Kutta methods or linear multistep methods. Each of these has limitations of one sort or another. In this paper we consider, as a middle ground, the derivation of continuous general linear methods for solution of stiff systems of initial value problems in ordinary differential equations. These methods are designed to combine the advantages of both Runge-Kutta and linear multistep methods. Particularly, methods possessing the property of A-stability are identified as promising methods within this large class of general linear methods. We show that the continuous general linear methods are self-starting and have more ability to solve the stiff systems of ordinary differential equations, than the discrete ones. The initial value systems of ordinary differential equations are solved, for instance, without looking for any other method to start the integration process. This desirable feature of the proposed approach leads to obtaining very high accuracy of the solution of the given problem. Illustrative examples are given to demonstrate the novelty and reliability of the methods.

Magnetohydrodynamics (MHD flow of a tangent hyperbolic fluid with nanoparticles past a stretching sheet with second order slip and convective boundary condition

Directory of Open Access Journals (Sweden)

Wubshet Ibrahim

Full Text Available This article presents the effect of thermal radiation on magnetohydrodynamic flow of tangent hyperbolic fluid with nanoparticle past an enlarging sheet with second order slip and convective boundary condition. Condition of zero normal flux of nanoparticles at the wall is used for the concentration boundary condition, which is the current topic that have yet to be studied extensively. The solution for the velocity, temperature and nanoparticle concentration is governed by parameters viz. power-law index (n, Weissenberg number We, Biot number Bi, Prandtl number Pr, velocity slip parameters δ and γ, Lewis number Le, Brownian motion parameter Nb and the thermophoresis parameter Nt. Similarity transformation is used to metamorphosed the governing non-linear boundary-value problem into coupled higher order non-linear ordinary differential equation. The succeeding equations were numerically solved using the function bvp4c from the matlab for different values of emerging parameters. Numerical results are deliberated through graphs and tables for velocity, temperature, concentration, the skin friction coefficient and local Nusselt number. The results designate that the skin friction coefficient Cf deplete as the values of Weissenberg number We, slip parameters γ and δ upturn and it rises as the values of power-law index n increase. The local Nusselt number -θ′(0 decreases as slip parameters γ and δ, radiation parameter Nr, Weissenberg number We, thermophoresis parameter Nt and power-law index n increase. However, the local Nusselt number increases as the Biot number Bi increase. Keywords: Tangent hyperbolic fluid, Second order slip flow, MHD, Convective boundary condition, Radiation effect, Passive control of nanoparticles

Customized Steady-State Constraints for Parameter Estimation in Non-Linear Ordinary Differential Equation Models.

Science.gov (United States)

Rosenblatt, Marcus; Timmer, Jens; Kaschek, Daniel

2016-01-01

Ordinary differential equation models have become a wide-spread approach to analyze dynamical systems and understand underlying mechanisms. Model parameters are often unknown and have to be estimated from experimental data, e.g., by maximum-likelihood estimation. In particular, models of biological systems contain a large number of parameters. To reduce the dimensionality of the parameter space, steady-state information is incorporated in the parameter estimation process. For non-linear models, analytical steady-state calculation typically leads to higher-order polynomial equations for which no closed-form solutions can be obtained. This can be circumvented by solving the steady-state equations for kinetic parameters, which results in a linear equation system with comparatively simple solutions. At the same time multiplicity of steady-state solutions is avoided, which otherwise is problematic for optimization. When solved for kinetic parameters, however, steady-state constraints tend to become negative for particular model specifications, thus, generating new types of optimization problems. Here, we present an algorithm based on graph theory that derives non-negative, analytical steady-state expressions by stepwise removal of cyclic dependencies between dynamical variables. The algorithm avoids multiple steady-state solutions by construction. We show that our method is applicable to most common classes of biochemical reaction networks containing inhibition terms, mass-action and Hill-type kinetic equations. Comparing the performance of parameter estimation for different analytical and numerical methods of incorporating steady-state information, we show that our approach is especially well-tailored to guarantee a high success rate of optimization.

A course in ordinary differential equations

CERN Document Server

Swift, Randall J

2014-01-01

Praise for the First Edition:"A Course in Ordinary Differential Equations deserves to be on the MAA's Basic Library List … the book with its layout, is very student friendly-it is easy to read and understand; every chapter and explanations flow smoothly and coherently … the reviewer would recommend this book highly for undergraduate introductory differential equation courses." -Srabasti Dutta, College of Saint Elizabeth, MAA Online, July 2008"An important feature is that the exposition is richly accompanied by computer algebra code (equally distributed between MATLAB, Mathematica, and Maple). The major part of the book is devoted to classical theory (both for systems and higher order equations). The necessary material from linear algebra is also covered. More advanced topics include numerical methods, stability of equilibria, bifurcations, Laplace transforms, and the power series method."-EMS Newsletter, June 2007"This is a delightful textbook for a standard one-semester undergraduate course in ordinary d...

Exploring inductive linearization for pharmacokinetic-pharmacodynamic systems of nonlinear ordinary differential equations.

Science.gov (United States)

Hasegawa, Chihiro; Duffull, Stephen B

2018-02-01

Pharmacokinetic-pharmacodynamic systems are often expressed with nonlinear ordinary differential equations (ODEs). While there are numerous methods to solve such ODEs these methods generally rely on time-stepping solutions (e.g. Runge-Kutta) which need to be matched to the characteristics of the problem at hand. The primary aim of this study was to explore the performance of an inductive approximation which iteratively converts nonlinear ODEs to linear time-varying systems which can then be solved algebraically or numerically. The inductive approximation is applied to three examples, a simple nonlinear pharmacokinetic model with Michaelis-Menten elimination (E1), an integrated glucose-insulin model and an HIV viral load model with recursive feedback systems (E2 and E3, respectively). The secondary aim of this study was to explore the potential advantages of analytically solving linearized ODEs with two examples, again E3 with stiff differential equations and a turnover model of luteinizing hormone with a surge function (E4). The inductive linearization coupled with a matrix exponential solution provided accurate predictions for all examples with comparable solution time to the matched time-stepping solutions for nonlinear ODEs. The time-stepping solutions however did not perform well for E4, particularly when the surge was approximated by a square wave. In circumstances when either a linear ODE is particularly desirable or the uncertainty in matching the integrator to the ODE system is of potential risk, then the inductive approximation method coupled with an analytical integration method would be an appropriate alternative.

Pointwise second-order necessary optimality conditions and second-order sensitivity relations in optimal control

Science.gov (United States)

Frankowska, Hélène; Hoehener, Daniel

2017-06-01

This paper is devoted to pointwise second-order necessary optimality conditions for the Mayer problem arising in optimal control theory. We first show that with every optimal trajectory it is possible to associate a solution p (ṡ) of the adjoint system (as in the Pontryagin maximum principle) and a matrix solution W (ṡ) of an adjoint matrix differential equation that satisfy a second-order transversality condition and a second-order maximality condition. These conditions seem to be a natural second-order extension of the maximum principle. We then prove a Jacobson like necessary optimality condition for general control systems and measurable optimal controls that may be only ;partially singular; and may take values on the boundary of control constraints. Finally we investigate the second-order sensitivity relations along optimal trajectories involving both p (ṡ) and W (ṡ).

Approximating second-order vector differential operators on distorted meshes in two space dimensions

International Nuclear Information System (INIS)

Hermeline, F.

2008-01-01

A new finite volume method is presented for approximating second-order vector differential operators in two space dimensions. This method allows distorted triangle or quadrilateral meshes to be used without the numerical results being too much altered. The matrices that need to be inverted are symmetric positive definite therefore, the most powerful linear solvers can be applied. The method has been tested on a few second-order vector partial differential equations coming from elasticity and fluids mechanics areas. These numerical experiments show that it is second-order accurate and locking-free. (authors)

A canonical form of the equation of motion of linear dynamical systems

Science.gov (United States)

Kawano, Daniel T.; Salsa, Rubens Goncalves; Ma, Fai; Morzfeld, Matthias

2018-03-01

The equation of motion of a discrete linear system has the form of a second-order ordinary differential equation with three real and square coefficient matrices. It is shown that, for almost all linear systems, such an equation can always be converted by an invertible transformation into a canonical form specified by two diagonal coefficient matrices associated with the generalized acceleration and displacement. This canonical form of the equation of motion is unique up to an equivalence class for non-defective systems. As an important by-product, a damped linear system that possesses three symmetric and positive definite coefficients can always be recast as an undamped and decoupled system.

An Online Method for Interpolating Linear Parametric Reduced-Order Models

KAUST Repository

Amsallem, David; Farhat, Charbel

2011-01-01

A two-step online method is proposed for interpolating projection-based linear parametric reduced-order models (ROMs) in order to construct a new ROM for a new set of parameter values. The first step of this method transforms each precomputed ROM into a consistent set of generalized coordinates. The second step interpolates the associated linear operators on their appropriate matrix manifold. Real-time performance is achieved by precomputing inner products between the reduced-order bases underlying the precomputed ROMs. The proposed method is illustrated by applications in mechanical and aeronautical engineering. In particular, its robustness is demonstrated by its ability to handle the case where the sampled parameter set values exhibit a mode veering phenomenon. © 2011 Society for Industrial and Applied Mathematics.

The shooting method and multiple solutions of two/multi-point BVPs of second-order ODE

Directory of Open Access Journals (Sweden)

Man Kam Kwong

2006-06-01

Full Text Available Within the last decade, there has been growing interest in the study of multiple solutions of two- and multi-point boundary value problems of nonlinear ordinary differential equations as fixed points of a cone mapping. Undeniably many good results have emerged. The purpose of this paper is to point out that, in the special case of second-order equations, the shooting method can be an effective tool, sometimes yielding better results than those obtainable via fixed point techniques.

The invariance of second-order functionals revisited

International Nuclear Information System (INIS)

Battezzati, M.

1984-01-01

In this paper some invariance properties of certain homogeneous functional forms of perturbative second-order energies with respect to transformations on the arguments are briefly considered. It has been shown that, if this energy is regarded as an Hamiltonian governing the time evolution of the arguments, which are the components of the first-order perturbed functions, the x and y couples play naturally the role of canonically conjugated co-ordinates and momenta. A search has been made for those linear transformations on these functions which preserve the above duality or reciprocity relations. It has been found that certain canonical transformations are of this type. In particular, the spinorial covariant-contravariant transformations for rotations in four-dimensional space-time

Second Order Ideal-Ward Continuity

Directory of Open Access Journals (Sweden)

Bipan Hazarika

2014-01-01

Full Text Available The main aim of the paper is to introduce a concept of second order ideal-ward continuity in the sense that a function f is second order ideal-ward continuous if I-limn→∞Δ2f(xn=0 whenever I-limn→∞Δ2xn=0 and a concept of second order ideal-ward compactness in the sense that a subset E of R is second order ideal-ward compact if any sequence x=(xn of points in E has a subsequence z=(zk=(xnk of the sequence x such that I-limk→∞Δ2zk=0 where Δ2zk=zk+2-2zk+1+zk. We investigate the impact of changing the definition of convergence of sequences on the structure of ideal-ward continuity in the sense of second order ideal-ward continuity and compactness of sets in the sense of second order ideal-ward compactness and prove related theorems.

Searching fundamental information in ordinary differential equations. Nondimensionalization technique.

Science.gov (United States)

Sánchez Pérez, J F; Conesa, M; Alhama, I; Alhama, F; Cánovas, M

2017-01-01

Classical dimensional analysis and nondimensionalization are assumed to be two similar approaches in the search for dimensionless groups. Both techniques, simplify the study of many problems. The first approach does not need to know the mathematical model, being sufficient a deep understanding of the physical phenomenon involved, while the second one begins with the governing equations and reduces them to their dimensionless form by simple mathematical manipulations. In this work, a formal protocol is proposed for applying the nondimensionalization process to ordinary differential equations, linear or not, leading to dimensionless normalized equations from which the resulting dimensionless groups have two inherent properties: In one hand, they are physically interpreted as balances between counteracting quantities in the problem, and on the other hand, they are of the order of magnitude unity. The solutions provided by nondimensionalization are more precise in every case than those from dimensional analysis, as it is illustrated by the applications studied in this work.

Extending the Constant Coefficient Solution Technique to Variable Coefficient Ordinary Differential Equations

Science.gov (United States)

Mohammed, Ahmed; Zeleke, Aklilu

2015-01-01

We introduce a class of second-order ordinary differential equations (ODEs) with variable coefficients whose closed-form solutions can be obtained by the same method used to solve ODEs with constant coefficients. General solutions for the homogeneous case are discussed.

Time-averaged second-order pressure and velocity measurements in a pressurized oscillating flow prime mover

Energy Technology Data Exchange (ETDEWEB)

Paridaens, Richard [DynFluid, Arts et Metiers, 151 boulevard de l' Hopital, Paris (France); Kouidri, Smaine [LIMSI-CNRS, Orsay Cedex (France)

2016-11-15

Nonlinear phenomena in oscillating flow devices cause the appearance of a relatively minor secondary flow known as acoustic streaming, which is superimposed on the primary oscillating flow. Knowledge of control parameters, such as the time-averaged second-order velocity and pressure, would elucidate the non-linear phenomena responsible for this part of the decrease in the system's energetic efficiency. This paper focuses on the characterization of a travelling wave oscillating flow engine by measuring the time-averaged second order pressure and velocity. Laser Doppler velocimetry technique was used to measure the time-averaged second-order velocity. As streaming is a second-order phenomenon, its measurement requires specific settings especially in a pressurized device. Difficulties in obtaining the proper settings are highlighted in this study. The experiments were performed for mean pressures varying from 10 bars to 22 bars. Non-linear effect does not constantly increase with pressure.

Time-course window estimator for ordinary differential equations linear in the parameters

NARCIS (Netherlands)

Vujacic, Ivan; Dattner, Itai; Gonzalez, Javier; Wit, Ernst

In many applications obtaining ordinary differential equation descriptions of dynamic processes is scientifically important. In both, Bayesian and likelihood approaches for estimating parameters of ordinary differential equations, the speed and the convergence of the estimation procedure may

Valuing a gas-fired power plant: A comparison of ordinary linear models, regime-switching approaches, and models with stochastic volatility

International Nuclear Information System (INIS)

Heydari, Somayeh; Siddiqui, Afzal

2010-01-01

Energy prices are often highly volatile with unexpected spikes. Capturing these sudden spikes may lead to more informed decision-making in energy investments, such as valuing gas-fired power plants, than ignoring them. In this paper, non-linear regime-switching models and models with mean-reverting stochastic volatility are compared with ordinary linear models. The study is performed using UK electricity and natural gas daily spot prices and suggests that with the aim of valuing a gas-fired power plant with and without operational flexibility, non-linear models with stochastic volatility, specifically for logarithms of electricity prices, provide better out-of-sample forecasts than both linear models and regime-switching models.

A novel algebraic procedure for solving non-linear evolution equations of higher order

International Nuclear Information System (INIS)

Huber, Alfred

2007-01-01

We report here a systematic approach that can easily be used for solving non-linear partial differential equations (nPDE), especially of higher order. We restrict the analysis to the so called evolution equations describing any wave propagation. The proposed new algebraic approach leads us to traveling wave solutions and moreover, new class of solution can be obtained. The crucial step of our method is the basic assumption that the solutions satisfy an ordinary differential equation (ODE) of first order that can be easily integrated. The validity and reliability of the method is tested by its application to some non-linear evolution equations. The important aspect of this paper however is the fact that we are able to calculate distinctive class of solutions which cannot be found in the current literature. In other words, using this new algebraic method the solution manifold is augmented to new class of solution functions. Simultaneously we would like to stress the necessity of such sophisticated methods since a general theory of nPDE does not exist. Otherwise, for practical use the algebraic construction of new class of solutions is of fundamental interest

Introduction to ordinary differential equations

CERN Document Server

Rabenstein, Albert L

1966-01-01

Introduction to Ordinary Differential Equations is a 12-chapter text that describes useful elementary methods of finding solutions using ordinary differential equations. This book starts with an introduction to the properties and complex variable of linear differential equations. Considerable chapters covered topics that are of particular interest in applications, including Laplace transforms, eigenvalue problems, special functions, Fourier series, and boundary-value problems of mathematical physics. Other chapters are devoted to some topics that are not directly concerned with finding solutio

Wetting transitions: First order or second order

International Nuclear Information System (INIS)

Teletzke, G.F.; Scriven, L.E.; Davis, H.T.

1982-01-01

A generalization of Sullivan's recently proposed theory of the equilibrium contact angle, the angle at which a fluid interface meets a solid surface, is investigated. The generalized theory admits either a first-order or second-order transition from a nonzero contact angle to perfect wetting as a critical point is approached, in contrast to Sullivan's original theory, which predicts only a second-order transition. The predictions of this computationally convenient theory are in qualitative agreement with a more rigorous theory to be presented in a future publication

General solutions of second-order linear difference equations of Euler type

Directory of Open Access Journals (Sweden)

Akane Hongyo

2017-01-01

Full Text Available The purpose of this paper is to give general solutions of linear difference equations which are related to the Euler-Cauchy differential equation \\(y^{\\prime\\prime}+(\\lambda/t^2y=0\\ or more general linear differential equations. We also show that the asymptotic behavior of solutions of the linear difference equations are similar to solutions of the linear differential equations.

Nonparametric Second-Order Theory of Error Propagation on Motion Groups.

Science.gov (United States)

Wang, Yunfeng; Chirikjian, Gregory S

2008-01-01

Error propagation on the Euclidean motion group arises in a number of areas such as in dead reckoning errors in mobile robot navigation and joint errors that accumulate from the base to the distal end of kinematic chains such as manipulators and biological macromolecules. We address error propagation in rigid-body poses in a coordinate-free way. In this paper we show how errors propagated by convolution on the Euclidean motion group, SE(3), can be approximated to second order using the theory of Lie algebras and Lie groups. We then show how errors that are small (but not so small that linearization is valid) can be propagated by a recursive formula derived here. This formula takes into account errors to second-order, whereas prior efforts only considered the first-order case. Our formulation is nonparametric in the sense that it will work for probability density functions of any form (not only Gaussians). Numerical tests demonstrate the accuracy of this second-order theory in the context of a manipulator arm and a flexible needle with bevel tip.

An implicit second order numerical method for two-fluid models

International Nuclear Information System (INIS)

Toumi, I.

1995-01-01

We present an implicit upwind numerical method for a six equation two-fluid model based on a linearized Riemann solver. The construction of this approximate Riemann solver uses an extension of Roe's scheme. Extension to second order accurate method is achieved using a piecewise linear approximation of the solution and a slope limiter method. For advancing in time, a linearized implicit integrating step is used. In practice this new numerical method has proved to be stable and capable of generating accurate non-oscillating solutions for two-phase flow calculations. The scheme was applied both to shock tube problems and to standard tests for two-fluid codes. (author)

On Solvability Theorems of Second-Order Ordinary Differential Equations with Delay

Directory of Open Access Journals (Sweden)

Nai-Sher Yeh

2018-01-01

Full Text Available For each x0∈[0,2π and k∈N, we obtain some existence theorems of periodic solutions to the two-point boundary value problem u′′(x+k2u(x-x0+g(x,u(x-x0=h(x in (0,2π with u(0-u(2π=u′(0-u′(2π=0 when g:(0,2π×R→R is a Caratheodory function which grows linearly in u as u→∞, and h∈L1(0,2π may satisfy a generalized Landesman-Lazer condition (1+sign(β∫02πh(xv(xdx0gβ+(xvx1-βdx+∫v(x<0gβ-(xvx1-βdx for all v∈N(L\\{0}. Here N(L denotes the subspace of L1(0,2π spanned by sin⁡kx and cos⁡kx, -1<β≤0, gβ+(x=lim  infu→∞(gx,uu/u1-β, and gβ-(x=lim  infu→-∞(gx,uu/u1-β.

Self-Tuning Control of Linear Systems Followed by Deadzones

Directory of Open Access Journals (Sweden)

K. Kazlauskas

2014-02-01

Full Text Available The aim of the present paper is to increase the efficiency of self-tuning generalized minimum variance (GMV control of linear time-invariant (LTI systems followed by deadzone nonlinearities. An approach, based on reordering of observations to be processed for the reconstruction of an unknown internal signal that acts between LTI system and a static nonlinear block of the closed-loop Wiener system, has been developed. The results of GMV self-tuning control of the second order LTI system with an ordinary deadzone are given.

Second-order nonlinearity induced transparency.

Science.gov (United States)

Zhou, Y H; Zhang, S S; Shen, H Z; Yi, X X

2017-04-01

In analogy to electromagnetically induced transparency, optomechanically induced transparency was proposed recently in [Science330, 1520 (2010)SCIEAS0036-807510.1126/science.1195596]. In this Letter, we demonstrate another form of induced transparency enabled by second-order nonlinearity. A practical application of the second-order nonlinearity induced transparency is to measure the second-order nonlinear coefficient. Our scheme might find applications in quantum optics and quantum information processing.

Existence of infinitely many periodic solutions for second-order nonautonomous Hamiltonian systems

Directory of Open Access Journals (Sweden)

Wen Guan

2015-04-01

Full Text Available By using minimax methods and critical point theory, we obtain infinitely many periodic solutions for a second-order nonautonomous Hamiltonian systems, when the gradient of potential energy does not exceed linear growth.

First integrals and parametric solutions of third-order ODEs admitting {\\mathfrak{sl}(2, {R})}

Science.gov (United States)

Ruiz, A.; Muriel, C.

2017-05-01

A complete set of first integrals for any third-order ordinary differential equation admitting a Lie symmetry algebra isomorphic to sl(2, {R}) is explicitly computed. These first integrals are derived from two linearly independent solutions of a linear second-order ODE, without additional integration. The general solution in parametric form can be obtained by using the computed first integrals. The study includes a parallel analysis of the four inequivalent realizations of sl(2, {R}) , and it is applied to several particular examples. These include the generalized Chazy equation, as well as an example of an equation which admits the most complicated of the four inequivalent realizations.

A new implementation of the second-order polarization propagator approximation (SOPPA)

DEFF Research Database (Denmark)

Packer, Martin J.; Dalskov, Erik K.; Enevoldsen, Thomas

1996-01-01

We present a new implementation of the second-order polarization propagator approximation (SOPPA) using a direct linear transformation approach, in which the SOPPA equations are solved iteratively. This approach has two important advantages over its predecessors. First, the direct linear...... and triplet transitions for benzene and naphthalene. The results compare well with experiment and CASPT2 values, calculated with identical basis sets and molecular geometries. This indicates that SOPPA can provide reliable values for excitation energies and response properties for relatively large molecular...

Nonoscillation criteria for half-linear second order difference equations

Czech Academy of Sciences Publication Activity Database

Došlý, Ondřej; Řehák, Pavel

2001-01-01

Roč. 42, - (2001), s. 453-464 ISSN 0898-1221 R&D Projects: GA ČR GA201/98/0677; GA ČR GA201/99/0295 Keywords : half-linear difference equation%nonoscillation criteria%variational principle Subject RIV: BA - General Mathematics Impact factor: 0.383, year: 2001

Green's matrix for a second-order self-adjoint matrix differential operator

International Nuclear Information System (INIS)

Sisman, Tahsin Cagri; Tekin, Bayram

2010-01-01

A systematic construction of the Green's matrix for a second-order self-adjoint matrix differential operator from the linearly independent solutions of the corresponding homogeneous differential equation set is carried out. We follow the general approach of extracting the Green's matrix from the Green's matrix of the corresponding first-order system. This construction is required in the cases where the differential equation set cannot be turned to an algebraic equation set via transform techniques.

Symmetries and integrability of a fourth-order Euler-Bernoulli beam equation

International Nuclear Information System (INIS)

Bokhari, Ashfaque H.; Zaman, F. D.; Mahomed, F. M.

2010-01-01

The complete symmetry group classification of the fourth-order Euler-Bernoulli ordinary differential equation, where the elastic modulus and the area moment of inertia are constants and the applied load is a function of the normal displacement, is obtained. We perform the Lie and Noether symmetry analysis of this problem. In the Lie analysis, the principal Lie algebra which is one dimensional extends in four cases, viz. the linear, exponential, general power law, and a negative fractional power law. It is further shown that two cases arise in the Noether classification with respect to the standard Lagrangian. That is, the linear case for which the Noether algebra dimension is one less than the Lie algebra dimension as well as the negative fractional power law. In the latter case the Noether algebra is three dimensional and is isomorphic to the Lie algebra which is sl(2,R). This exceptional case, although admitting the nonsolvable algebra sl(2,R), remarkably allows for a two-parameter family of exact solutions via the Noether integrals. The Lie reduction gives a second-order ordinary differential equation which has nonlocal symmetry.

Second-Order Science of Interdisciplinary Research

DEFF Research Database (Denmark)

Alrøe, Hugo Fjelsted; Noe, Egon

2014-01-01

require and challenge interdisciplinarity. Problem: The conventional methods of interdisciplinary research fall short in the case of wicked problems because they remain first-order science. Our aim is to present workable methods and research designs for doing second-order science in domains where...... there are many different scientific knowledges on any complex problem. Method: We synthesize and elaborate a framework for second-order science in interdisciplinary research based on a number of earlier publications, experiences from large interdisciplinary research projects, and a perspectivist theory...... of science. Results: The second-order polyocular framework for interdisciplinary research is characterized by five principles. Second-order science of interdisciplinary research must: 1. draw on the observations of first-order perspectives, 2. address a shared dynamical object, 3. establish a shared problem...

On the asymptotic expansions of solutions of an nth order linear differential equation with power coefficients

International Nuclear Information System (INIS)

Paris, R.B.; Wood, A.D.

1984-11-01

The asymptotic expansions of solutions of a class of linear ordinary differential equations of arbitrary order n, containing a factor zsup(m) multiplying the lower order derivatives, are investigated for large values of z in the complex plane. Four classes of solutions are considered which exhibit the following behaviour as /z/ → infinity in certain sectors: (i) solutions whose behaviour is either exponentially large or algebraic (involving p ( < n) algebraic expansions), (ii) solutions which are exponentially small (iii) solutions with a single algebraic expansion and (iv) solutions which are even and odd functions of z whenever n+m is even. The asymptotic expansions of these solutions in a full neigbourhood of the point at infinity are obtained by means of the theory of the solutions in the case m=O developed in a previous paper

Second-Order Systems of ODEs Admitting Three-Dimensional Lie Algebras and Integrability

Directory of Open Access Journals (Sweden)

Muhammad Ayub

2013-01-01

the case of k≥3. We discuss the singular invariant representations of canonical forms for systems of two second-order ODEs admitting three-dimensional Lie algebras. Furthermore, we give an integration procedure for canonical forms for systems of two second-order ODEs admitting three-dimensional Lie algebras which comprises of two approaches, namely, division into four types I, II, III, and IV and that of integrability of the invariant representations. We prove that if a system of two second-order ODEs has a three-dimensional solvable Lie algebra, then, its general solution can be obtained from a partially linear, partially coupled or reduced invariantly represented system of equations. A natural extension of this result is provided for a system of two kth-order (k≥3 ODEs. We present illustrative examples of familiar integrable physical systems which admit three-dimensional Lie algebras such as the classical Kepler problem and the generalized Ermakov systems that give rise to closed trajectories.

Ordinary differential equations

CERN Document Server

Greenberg, Michael D

2014-01-01

Features a balance between theory, proofs, and examples and provides applications across diverse fields of study Ordinary Differential Equations presents a thorough discussion of first-order differential equations and progresses to equations of higher order. The book transitions smoothly from first-order to higher-order equations, allowing readers to develop a complete understanding of the related theory. Featuring diverse and interesting applications from engineering, bioengineering, ecology, and biology, the book anticipates potential difficulties in understanding the various solution steps

Uncertainty of pesticide residue concentration determined from ordinary and weighted linear regression curve.

Science.gov (United States)

Yolci Omeroglu, Perihan; Ambrus, Árpad; Boyacioglu, Dilek

2018-03-28

Determination of pesticide residues is based on calibration curves constructed for each batch of analysis. Calibration standard solutions are prepared from a known amount of reference material at different concentration levels covering the concentration range of the analyte in the analysed samples. In the scope of this study, the applicability of both ordinary linear and weighted linear regression (OLR and WLR) for pesticide residue analysis was investigated. We used 782 multipoint calibration curves obtained for 72 different analytical batches with high-pressure liquid chromatography equipped with an ultraviolet detector, and gas chromatography with electron capture, nitrogen phosphorus or mass spectrophotometer detectors. Quality criteria of the linear curves including regression coefficient, standard deviation of relative residuals and deviation of back calculated concentrations were calculated both for WLR and OLR methods. Moreover, the relative uncertainty of the predicted analyte concentration was estimated for both methods. It was concluded that calibration curve based on WLR complies with all the quality criteria set by international guidelines compared to those calculated with OLR. It means that all the data fit well with WLR for pesticide residue analysis. It was estimated that, regardless of the actual concentration range of the calibration, relative uncertainty at the lowest calibrated level ranged between 0.3% and 113.7% for OLR and between 0.2% and 22.1% for WLR. At or above 1/3 of the calibrated range, uncertainty of calibration curve ranged between 0.1% and 16.3% for OLR and 0% and 12.2% for WLR, and therefore, the two methods gave comparable results.

Stream Kriging: Incremental and recursive ordinary Kriging over spatiotemporal data streams

Science.gov (United States)

Zhong, Xu; Kealy, Allison; Duckham, Matt

2016-05-01

Ordinary Kriging is widely used for geospatial interpolation and estimation. Due to the O (n3) time complexity of solving the system of linear equations, ordinary Kriging for a large set of source points is computationally intensive. Conducting real-time Kriging interpolation over continuously varying spatiotemporal data streams can therefore be especially challenging. This paper develops and tests two new strategies for improving the performance of an ordinary Kriging interpolator adapted to a stream-processing environment. These strategies rely on the expectation that, over time, source data points will frequently refer to the same spatial locations (for example, where static sensor nodes are generating repeated observations of a dynamic field). First, an incremental strategy improves efficiency in cases where a relatively small proportion of previously processed spatial locations are absent from the source points at any given iteration. Second, a recursive strategy improves efficiency in cases where there is substantial set overlap between the sets of spatial locations of source points at the current and previous iterations. These two strategies are evaluated in terms of their computational efficiency in comparison to ordinary Kriging algorithm. The results show that these two strategies can reduce the time taken to perform the interpolation by up to 90%, and approach average-case time complexity of O (n2) when most but not all source points refer to the same locations over time. By combining the approaches developed in this paper with existing heuristic ordinary Kriging algorithms, the conclusions indicate how further efficiency gains could potentially be accrued. The work ultimately contributes to the development of online ordinary Kriging interpolation algorithms, capable of real-time spatial interpolation with large streaming data sets.

Access is mainly a second-order process: SDT models whether phenomenally (first-order) conscious states are accessed by reflectively (second-order) conscious processes.

Science.gov (United States)

Snodgrass, Michael; Kalaida, Natasha; Winer, E Samuel

2009-06-01

Access can either be first-order or second-order. First order access concerns whether contents achieve representation in phenomenal consciousness at all; second-order access concerns whether phenomenally conscious contents are selected for metacognitive, higher order processing by reflective consciousness. When the optional and flexible nature of second-order access is kept in mind, there remain strong reasons to believe that exclusion failure can indeed isolate phenomenally conscious stimuli that are not so accessed. Irvine's [Irvine, E. (2009). Signal detection theory, the exclusion failure paradigm and weak consciousness-Evidence for the access/phenomenal distinction? Consciousness and Cognition.] partial access argument fails because exclusion failure is indeed due to lack of second-order access, not insufficient phenomenally conscious information. Further, the enable account conforms with both qualitative differences and subjective report, and is simpler than the endow account. Finally, although first-order access may be a distinct and important process, second-order access arguably reflects the core meaning of access generally.

A Second Look at Second-Order Belief Attribution in Autism.

Science.gov (United States)

Tager-Flusberg, Helen; Sullivan, Kate

1994-01-01

Twelve students with autism and 12 with mental retardation, who had passed a first-order test of false belief, were given a second-order reasoning task. No intergroup performance differences were seen. Findings suggest that the difficulty for both groups with the second-order task lies in information processing demands rather than conceptual…

Second-order gauge-invariant perturbations during inflation

International Nuclear Information System (INIS)

Finelli, F.; Marozzi, G.; Vacca, G. P.; Venturi, G.

2006-01-01

The evolution of gauge invariant second-order scalar perturbations in a general single field inflationary scenario are presented. Different second-order gauge-invariant expressions for the curvature are considered. We evaluate perturbatively one of these second order curvature fluctuations and a second-order gauge-invariant scalar field fluctuation during the slow-roll stage of a massive chaotic inflationary scenario, taking into account the deviation from a pure de Sitter evolution and considering only the contribution of super-Hubble perturbations in mode-mode coupling. The spectra resulting from their contribution to the second order quantum correlation function are nearly scale-invariant, with additional logarithmic corrections with respect to the first order spectrum. For all scales of interest the amplitude of these spectra depends on the total number of e-folds. We find, on comparing first and second order perturbation results, an upper limit to the total number of e-folds beyond which the two orders are comparable

Spatial Processes in Linear Ordering

Science.gov (United States)

von Hecker, Ulrich; Klauer, Karl Christoph; Wolf, Lukas; Fazilat-Pour, Masoud

2016-01-01

Memory performance in linear order reasoning tasks (A > B, B > C, C > D, etc.) shows quicker, and more accurate responses to queries on wider (AD) than narrower (AB) pairs on a hypothetical linear mental model (A -- B -- C -- D). While indicative of an analogue representation, research so far did not provide positive evidence for spatial…

Relaxation approximations to second-order traffic flow models by high-resolution schemes

International Nuclear Information System (INIS)

Nikolos, I.K.; Delis, A.I.; Papageorgiou, M.

2015-01-01

A relaxation-type approximation of second-order non-equilibrium traffic models, written in conservation or balance law form, is considered. Using the relaxation approximation, the nonlinear equations are transformed to a semi-linear diagonilizable problem with linear characteristic variables and stiff source terms with the attractive feature that neither Riemann solvers nor characteristic decompositions are in need. In particular, it is only necessary to provide the flux and source term functions and an estimate of the characteristic speeds. To discretize the resulting relaxation system, high-resolution reconstructions in space are considered. Emphasis is given on a fifth-order WENO scheme and its performance. The computations reported demonstrate the simplicity and versatility of relaxation schemes as numerical solvers

Estimation of Multiple Point Sources for Linear Fractional Order Systems Using Modulating Functions

KAUST Repository

Belkhatir, Zehor

2017-06-28

This paper proposes an estimation algorithm for the characterization of multiple point inputs for linear fractional order systems. First, using polynomial modulating functions method and a suitable change of variables the problem of estimating the locations and the amplitudes of a multi-pointwise input is decoupled into two algebraic systems of equations. The first system is nonlinear and solves for the time locations iteratively, whereas the second system is linear and solves for the input’s amplitudes. Second, closed form formulas for both the time location and the amplitude are provided in the particular case of single point input. Finally, numerical examples are given to illustrate the performance of the proposed technique in both noise-free and noisy cases. The joint estimation of pointwise input and fractional differentiation orders is also presented. Furthermore, a discussion on the performance of the proposed algorithm is provided.

The effects of second-order hydrodynamics on a semisubmersible floating offshore wind turbine

International Nuclear Information System (INIS)

Bayati, I; Jonkman, J; Robertson, A; Platt, A

2014-01-01

The objective of this paper is to assess the second-order hydrodynamic effects on a semisubmersible floating offshore wind turbine. Second-order hydrodynamics induce loads and motions at the sum- and difference-frequencies of the incident waves. These effects have often been ignored in offshore wind analysis, under the assumption that they are significantly smaller than first-order effects. The sum- and difference-frequency loads can, however, excite eigenfrequencies of a floating system, leading to large oscillations that strain the mooring system or vibrations that cause fatigue damage to the structure. Observations of supposed second-order responses in wave-tank tests performed by the DeepCwind consortium at the Maritime Research Institute Netherlands (MARIN) offshore basin suggest that these effects might be more important than originally expected. These observations inspired interest in investigating how second-order excitation affects floating offshore wind turbines and whether second-order hydrodynamics should be included in offshore wind simulation tools like FAST. In this work, the effects of second-order hydrodynamics on a floating semisubmersible offshore wind turbine are investigated. Because FAST is currently unable to account for second-order effects, a method to assess these effects was applied in which linearized properties of the floating wind system derived from FAST (including the 6x6 mass and stiffness matrices) are used by WAMIT to solve the first- and second-order hydrodynamics problems in the frequency domain. The method was applied to the Offshore Code Comparison Collaboration Continuation OC4-DeepCwind semisubmersible platform, supporting the National Renewable Energy Laboratory's 5-MW baseline wind turbine. In this paper, the loads and response of the system caused by the second-order hydrodynamics are analysed and compared to the first-order hydrodynamic loads and induced motions in the frequency domain. Further, the second-order

Mixed H2/Hinfinity output-feedback control of second-order neutral systems with time-varying state and input delays.

Science.gov (United States)

Karimi, Hamid Reza; Gao, Huijun

2008-07-01

A mixed H2/Hinfinity output-feedback control design methodology is presented in this paper for second-order neutral linear systems with time-varying state and input delays. Delay-dependent sufficient conditions for the design of a desired control are given in terms of linear matrix inequalities (LMIs). A controller, which guarantees asymptotic stability and a mixed H2/Hinfinity performance for the closed-loop system of the second-order neutral linear system, is then developed directly instead of coupling the model to a first-order neutral system. A Lyapunov-Krasovskii method underlies the LMI-based mixed H2/Hinfinity output-feedback control design using some free weighting matrices. The simulation results illustrate the effectiveness of the proposed methodology.

Spectrum of Discrete Second-Order Difference Operator with Sign-Changing Weight and Its Applications

Directory of Open Access Journals (Sweden)

Ruyun Ma

2014-01-01

Full Text Available Let T>1 be an integer, and let=1,2,…,T. We discuss the spectrum of discrete linear second-order eigenvalue problems Δ2ut-1+λmtut=0, t∈,  u0=uT+1=0, where λ≠0 is a parameter, m:→ℝ changes sign and mt≠0 on . At last, as an application of this spectrum result, we show the existence of sign-changing solutions of discrete nonlinear second-order problems by using bifurcate technique.

Source of second order chromaticity in RHIC

International Nuclear Information System (INIS)

Luo, Y.; Gu, X.; Fischer, W.; Trbojevic, D.

2011-01-01

In this note we will answer the following questions: (1) what is the source of second order chromaticities in RHIC? (2) what is the dependence of second order chromaticity on the on-momentum β-beat? (3) what is the dependence of second order chromaticity on β* at IP6 and IP8? To answer these questions, we use the perturbation theory to numerically calculate the contributions of each quadrupole and sextupole to the first, second, and third order chromaticities.

An Analysis of Second-Order Autoshaping

Science.gov (United States)

Ward-Robinson, Jasper

2004-01-01

Three mechanisms can explain second-order conditioning: (1) The second-order conditioned stimulus (CS2) could activate a representation of the first-order conditioned stimulus (CS1), thereby provoking the conditioned response (CR); The CS2 could enter into an excitatory association with either (2) the representation governing the CR, or (3) with a…

Iterative oscillation results for second-order differential equations with advanced argument

Directory of Open Access Journals (Sweden)

Irena Jadlovska

2017-07-01

Full Text Available This article concerns the oscillation of solutions to a linear second-order differential equation with advanced argument. Sufficient oscillation conditions involving limit inferior are given which essentially improve known results. We base our technique on the iterative construction of solution estimates and some of the recent ideas developed for first-order advanced differential equations. We demonstrate the advantage of our results on Euler-type advanced equation. Using MATLAB software, a comparison of the effectiveness of newly obtained criteria as well as the necessary iteration length in particular cases are discussed.

Pap-smear Classification Using Efficient Second Order Neural Network Training Algorithms

DEFF Research Database (Denmark)

Ampazis, Nikolaos; Dounias, George; Jantzen, Jan

2004-01-01

In this paper we make use of two highly efficient second order neural network training algorithms, namely the LMAM (Levenberg-Marquardt with Adaptive Momentum) and OLMAM (Optimized Levenberg-Marquardt with Adaptive Momentum), for the construction of an efficient pap-smear test classifier. The alg......In this paper we make use of two highly efficient second order neural network training algorithms, namely the LMAM (Levenberg-Marquardt with Adaptive Momentum) and OLMAM (Optimized Levenberg-Marquardt with Adaptive Momentum), for the construction of an efficient pap-smear test classifier....... The algorithms are methodologically similar, and are based on iterations of the form employed in the Levenberg-Marquardt (LM) method for non-linear least squares problems with the inclusion of an additional adaptive momentum term arising from the formulation of the training task as a constrained optimization...

Existence of solutions to second-order nonlinear coupled systems with nonlinear coupled boundary conditions

Directory of Open Access Journals (Sweden)

Imran Talib

2015-12-01

Full Text Available In this article, study the existence of solutions for the second-order nonlinear coupled system of ordinary differential equations $$\\displaylines{ u''(t=f(t,v(t,\\quad t\\in [0,1],\\cr v''(t=g(t,u(t,\\quad t\\in [0,1], }$$ with nonlinear coupled boundary conditions $$\\displaylines{ \\phi(u(0,v(0,u(1,v(1,u'(0,v'(0=(0,0, \\cr \\psi(u(0,v(0,u(1,v(1,u'(1,v'(1=(0,0, }$$ where $f,g:[0,1]\\times \\mathbb{R}\\to \\mathbb{R}$ and $\\phi,\\psi:\\mathbb{R}^6\\to \\mathbb{R}^2$ are continuous functions. Our main tools are coupled lower and upper solutions, Arzela-Ascoli theorem, and Schauder's fixed point theorem.

On the existence of eigenmodes of linear quasi-periodic differential equations and their relation to the MHD continuum

International Nuclear Information System (INIS)

Salat, A.

1981-12-01

The existence of quasi-periodic eigensolutions of a linear second order ordinary differential equation with quasi-periodic coefficient f(ω 1 t,ω 2 t) is investigated numerically and graphically. For sufficiently incommensurate frequencies ω 1 , ω 2 a doubly indexed infinite sequence of eigenvalues and eigenmodes is obtained. The equation considered is a model for the magneto-hydrodynamic 'continuum' in general toroidal geometry. The result suggests that continuum modes exist at least on sufficiently irrational magnetic surfaces. (orig.)

Generalized Second-Order Parametric Optimality Conditions in Semiinfinite Discrete Minmax Fractional Programming and Second-Order Univexity

Directory of Open Access Journals (Sweden)

Ram Verma

2016-02-01

Full Text Available This paper deals with mainly establishing numerous sets of generalized second order paramertic sufficient optimality conditions for a semiinfinite discrete minmax fractional programming problem, while the results on semiinfinite discrete minmax fractional programming problem achieved based on some partitioning schemes under various types of generalized second order univexity assumptions. 

On the Robustness of Hysteretic Second-Order Systems with PID : iISS approach

NARCIS (Netherlands)

Ouyang, Ruiyue; Jayawardhana, Bayu; Andrieu, Vincent

2012-01-01

In this paper, we study the robustness property of a second-order linear plant controlled by a proportional, integral and derivative (PID) controller with a hysteretic actuator. The hysteretic actuator is modeled by a Duhem model that exhibits clockwise (CW) input-output (I/O) dynamics (such as the

Second order approximation for optical polaron in the strong coupling case

International Nuclear Information System (INIS)

Bogolubov, N.N. Jr.

1993-11-01

Here we propose a method of construction second order approximation for ground state energy for class of model Hamiltonian with linear type interaction on Bose operators in strong coupling case. For the application of the above method we have considered polaron model and propose construction set of nonlinear differential equations for definition ground state energy in strong coupling case. We have considered also radial symmetry case. (author). 10 refs

First-order systems of linear partial differential equations: normal forms, canonical systems, transform methods

Directory of Open Access Journals (Sweden)

Heinz Toparkus

2014-04-01

Full Text Available In this paper we consider first-order systems with constant coefficients for two real-valued functions of two real variables. This is both a problem in itself, as well as an alternative view of the classical linear partial differential equations of second order with constant coefficients. The classification of the systems is done using elementary methods of linear algebra. Each type presents its special canonical form in the associated characteristic coordinate system. Then you can formulate initial value problems in appropriate basic areas, and you can try to achieve a solution of these problems by means of transform methods.

Analysis of second order harmonic distortion due to transmitter non-linearity and chromatic and modal dispersion of optical OFDM SSB modulated signals in SMF-MMF fiber links

Science.gov (United States)

Patel, Dhananjay; Singh, Vinay Kumar; Dalal, U. D.

2017-01-01

Single mode fibers (SMF) are typically used in Wide Area Networks (WAN), Metropolitan Area Networks (MAN) and also find applications in Radio over Fiber (RoF) architectures supporting data transmission in Fiber to the Home (FTTH), Remote Antenna Units (RAUs), in-building networks etc. Multi-mode fibers (MMFs) with low cost, ease of installation and low maintenance are predominantly (85-90%) deployed in-building networks providing data access in local area networks (LANs). The transmission of millimeter wave signals through the SMF in WAN and MAN, along with the reuse of MMF in-building networks will not levy fiber reinstallation cost. The transmission of the millimeter waves experiences signal impairments due to the transmitter non-linearity and modal dispersion of the MMF. The MMF exhibiting large modal dispersion limits the bandwidth-length product of the fiber. The second and higher-order harmonics present in the optical signal fall within the system bandwidth. This causes degradation in the received signal and an unwanted radiation of power at the RAU. The power of these harmonics is proportional to the non-linearity of the transmitter and the modal dispersion of the MMF and should be maintained below the standard values as per the international norms. In this paper, a mathematical model is developed for Second-order Harmonic Distortion (HD2) generated due to non-linearity of the transmitter and chromatic-modal dispersion of the SMF-MMF optic link. This is also verified using a software simulation. The model consists of a Mach Zehnder Modulator (MZM) that generates two m-QAM OFDM Single Sideband (SSB) signals based on phase shift of the hybrid coupler (90° and 120°). Our results show that the SSB signal with 120° hybrid coupler has suppresses the higher-order harmonics and makes the system more robust against the HD2 in the SMF-MMF optic link.

Class of unconditionally stable second-order implicit schemes for hyperbolic and parabolic equations

International Nuclear Information System (INIS)

Lui, H.C.

The linearized Burgers equation is considered as a model u/sub t/ tau/sub x/ = bu/sub xx/, where the subscripts t and x denote the derivatives of the function u with respect to time t and space x; a and b are constants (b greater than or equal to 0). Numerical schemes for solving the equation are described that are second-order accurate, unconditionally stable, and dissipative of higher order. (U.S.)

Second order finite volume scheme for Maxwell's equations with discontinuous electromagnetic properties on unstructured meshes

Energy Technology Data Exchange (ETDEWEB)

Ismagilov, Timur Z., E-mail: ismagilov@academ.org

2015-02-01

This paper presents a second order finite volume scheme for numerical solution of Maxwell's equations with discontinuous dielectric permittivity and magnetic permeability on unstructured meshes. The scheme is based on Godunov scheme and employs approaches of Van Leer and Lax–Wendroff to increase the order of approximation. To keep the second order of approximation near dielectric permittivity and magnetic permeability discontinuities a novel technique for gradient calculation and limitation is applied near discontinuities. Results of test computations for problems with linear and curvilinear discontinuities confirm second order of approximation. The scheme was applied to modelling propagation of electromagnetic waves inside photonic crystal waveguides with a bend.

Structural changes of small amplitude kinetic Alfvén solitary waves due to second-order corrections

International Nuclear Information System (INIS)

Choi, Cheong R.

2015-01-01

The structural changes of kinetic Alfvén solitary waves (KASWs) due to higher-order terms are investigated. While the first-order differential equation for KASWs provides the dispersion relation for kinetic Alfvén waves, the second-order differential equation describes the structural changes of the solitary waves due to higher-order nonlinearity. The reductive perturbation method is used to obtain the second-order and third-order partial differential equations; then, Kodama and Taniuti's technique [J. Phys. Soc. Jpn. 45, 298 (1978)] is applied in order to remove the secularities in the third-order differential equations and derive a linear second-order inhomogeneous differential equation. The solution to this new second-order equation indicates that, as the amplitude increases, the hump-type Korteweg-de Vries solution is concentrated more around the center position of the soliton and that dip-type structures form near the two edges of the soliton. This result has a close relationship with the interpretation of the complex KASW structures observed in space with satellites

A unified model for transfer alignment at random misalignment angles based on second-order EKF

International Nuclear Information System (INIS)

Cui, Xiao; Qin, Yongyuan; Yan, Gongmin; Liu, Zhenbo; Mei, Chunbo

2017-01-01

In the transfer alignment process of inertial navigation systems (INSs), the conventional linear error model based on the small misalignment angle assumption cannot be applied to large misalignment situations. Furthermore, the nonlinear model based on the large misalignment angle suffers from redundant computation with nonlinear filters. This paper presents a unified model for transfer alignment suitable for arbitrary misalignment angles. The alignment problem is transformed into an estimation of the relative attitude between the master INS (MINS) and the slave INS (SINS), by decomposing the attitude matrix of the latter. Based on the Rodriguez parameters, a unified alignment model in the inertial frame with the linear state-space equation and a second order nonlinear measurement equation are established, without making any assumptions about the misalignment angles. Furthermore, we employ the Taylor series expansions on the second-order nonlinear measurement equation to implement the second-order extended Kalman filter (EKF2). Monte-Carlo simulations demonstrate that the initial alignment can be fulfilled within 10 s, with higher accuracy and much smaller computational cost compared with the traditional unscented Kalman filter (UKF) at large misalignment angles. (paper)

A unified model for transfer alignment at random misalignment angles based on second-order EKF

Science.gov (United States)

Cui, Xiao; Mei, Chunbo; Qin, Yongyuan; Yan, Gongmin; Liu, Zhenbo

2017-04-01

In the transfer alignment process of inertial navigation systems (INSs), the conventional linear error model based on the small misalignment angle assumption cannot be applied to large misalignment situations. Furthermore, the nonlinear model based on the large misalignment angle suffers from redundant computation with nonlinear filters. This paper presents a unified model for transfer alignment suitable for arbitrary misalignment angles. The alignment problem is transformed into an estimation of the relative attitude between the master INS (MINS) and the slave INS (SINS), by decomposing the attitude matrix of the latter. Based on the Rodriguez parameters, a unified alignment model in the inertial frame with the linear state-space equation and a second order nonlinear measurement equation are established, without making any assumptions about the misalignment angles. Furthermore, we employ the Taylor series expansions on the second-order nonlinear measurement equation to implement the second-order extended Kalman filter (EKF2). Monte-Carlo simulations demonstrate that the initial alignment can be fulfilled within 10 s, with higher accuracy and much smaller computational cost compared with the traditional unscented Kalman filter (UKF) at large misalignment angles.

Systemic Design for Second-Order Effects

Directory of Open Access Journals (Sweden)

Evan Barba

2017-04-01

Full Text Available Second-order effects refer to changes within a system that are the result of changes made somewhere else in the system (the first-order effects. Second-order effects can occur at different spatial, temporal, or organizational scales from the original interventions, and are difficult to control. Some organizational theorists suggest that careful management of feedback processes can facilitate controlled change from one organizational configuration to another. Recognizing that skill in managing feedback processes is a core competency of design suggests that design skills are potentially useful tools in achieving organizational change. This paper describes a case study in which a co-design methodology was used to control the second-order effects resulting from a classroom intervention to create organizational change. This approach is then theorized as the Instigator Systems approach.

Range-separated time-dependent density-functional theory with a frequency-dependent second-order Bethe-Salpeter correlation kernel

Energy Technology Data Exchange (ETDEWEB)

Rebolini, Elisa, E-mail: elisa.rebolini@kjemi.uio.no; Toulouse, Julien, E-mail: julien.toulouse@upmc.fr [Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Univ Paris 06, CNRS, 4 place Jussieu, F-75005 Paris (France)

2016-03-07

We present a range-separated linear-response time-dependent density-functional theory (TDDFT) which combines a density-functional approximation for the short-range response kernel and a frequency-dependent second-order Bethe-Salpeter approximation for the long-range response kernel. This approach goes beyond the adiabatic approximation usually used in linear-response TDDFT and aims at improving the accuracy of calculations of electronic excitation energies of molecular systems. A detailed derivation of the frequency-dependent second-order Bethe-Salpeter correlation kernel is given using many-body Green-function theory. Preliminary tests of this range-separated TDDFT method are presented for the calculation of excitation energies of the He and Be atoms and small molecules (H{sub 2}, N{sub 2}, CO{sub 2}, H{sub 2}CO, and C{sub 2}H{sub 4}). The results suggest that the addition of the long-range second-order Bethe-Salpeter correlation kernel overall slightly improves the excitation energies.

On the second-order temperature jump coefficient of a dilute gas

Science.gov (United States)

Radtke, Gregg A.; Hadjiconstantinou, N. G.; Takata, S.; Aoki, K.

2012-09-01

We use LVDSMC simulations to calculate the second-order temperature jump coefficient for a dilute gas whose temperature is governed by the Poisson equation with a constant forcing term. Both the hard sphere gas and the BGK model of the Boltzmann equation are considered. Our results show that the temperature jump coefficient is different from the well known linear and steady case where the temperature is governed by the homogeneous heat conduction (Laplace) equation.

First- and second-order processing in transient stereopsis.

Science.gov (United States)

Edwards, M; Pope, D R; Schor, C M

2000-01-01

Large-field stimuli were used to investigate the interaction of first- and second-order pathways in transient-stereo processing. Stimuli consisted of sinewave modulations in either the mean luminance (first-order stimulus) or the contrast (second-order stimulus) of a dynamic-random-dot field. The main results of the present study are that: (1) Depth could be extracted with both the first-order and second-order stimuli; (2) Depth could be extracted from dichoptically mixed first- and second-order stimuli, however, the same stimuli, when presented as a motion sequence, did not result in a motion percept. Based upon these findings we conclude that the transient-stereo system processes both first- and second-order signals, and that these two signals are pooled prior to the extraction of transient depth. This finding of interaction between first- and second-order stereoscopic processing is different from the independence that has been found with the motion system.

Motion induced second order temperature and y-type anisotropies after the subtraction of linear dipole in the CMB maps

International Nuclear Information System (INIS)

Sunyaev, Rashid A.; Khatri, Rishi

2013-01-01

y-type spectral distortions of the cosmic microwave background allow us to detect clusters and groups of galaxies, filaments of hot gas and the non-uniformities in the warm hot intergalactic medium. Several CMB experiments (on small areas of sky) and theoretical groups (for full sky) have recently published y-type distortion maps. We propose to search for two artificial hot spots in such y-type maps resulting from the incomplete subtraction of the effect of the motion induced dipole on the cosmic microwave background sky. This dipole introduces, at second order, additional temperature and y-distortion anisotropy on the sky of amplitude few μK which could potentially be measured by Planck HFI and Pixie experiments and can be used as a source of cross channel calibration by CMB experiments. This y-type distortion is present in every pixel and is not the result of averaging the whole sky. This distortion, calculated exactly from the known linear dipole, can be subtracted from the final y-type maps, if desired

Modulation masking produced by second-order modulators

DEFF Research Database (Denmark)

Füllgrabe, Christian; Moore, Brian C.J.; Demany, Laurent

2005-01-01

Recent studies suggest that an auditory nonlinearity converts second-order sinusoidal amplitude modulation (SAM) (i.e., modulation of SAM depth) into a first-order SAM component, which contributes to the perception of second-order SAM. However, conversion may also occur in other ways such as coch...

A second order discontinuous Galerkin fast sweeping method for Eikonal equations

Science.gov (United States)

Li, Fengyan; Shu, Chi-Wang; Zhang, Yong-Tao; Zhao, Hongkai

2008-09-01

In this paper, we construct a second order fast sweeping method with a discontinuous Galerkin (DG) local solver for computing viscosity solutions of a class of static Hamilton-Jacobi equations, namely the Eikonal equations. Our piecewise linear DG local solver is built on a DG method developed recently [Y. Cheng, C.-W. Shu, A discontinuous Galerkin finite element method for directly solving the Hamilton-Jacobi equations, Journal of Computational Physics 223 (2007) 398-415] for the time-dependent Hamilton-Jacobi equations. The causality property of Eikonal equations is incorporated into the design of this solver. The resulting local nonlinear system in the Gauss-Seidel iterations is a simple quadratic system and can be solved explicitly. The compactness of the DG method and the fast sweeping strategy lead to fast convergence of the new scheme for Eikonal equations. Extensive numerical examples verify efficiency, convergence and second order accuracy of the proposed method.

Local linearization methods for the numerical integration of ordinary differential equations: An overview

International Nuclear Information System (INIS)

Jimenez, J.C.

2009-06-01

Local Linearization (LL) methods conform a class of one-step explicit integrators for ODEs derived from the following primary and common strategy: the vector field of the differential equation is locally (piecewise) approximated through a first-order Taylor expansion at each time step, thus obtaining successive linear equations that are explicitly integrated. Hereafter, the LL approach may include some additional strategies to improve that basic affine approximation. Theoretical and practical results have shown that the LL integrators have a number of convenient properties. These include arbitrary order of convergence, A-stability, linearization preserving, regularity under quite general conditions, preservation of the dynamics of the exact solution around hyperbolic equilibrium points and periodic orbits, integration of stiff and high-dimensional equations, low computational cost, and others. In this paper, a review of the LL methods and their properties is presented. (author)

A Second-Order Maximum Principle Preserving Lagrange Finite Element Technique for Nonlinear Scalar Conservation Equations

KAUST Repository

Guermond, Jean-Luc; Nazarov, Murtazo; Popov, Bojan; Yang, Yong

2014-01-01

© 2014 Society for Industrial and Applied Mathematics. This paper proposes an explicit, (at least) second-order, maximum principle satisfying, Lagrange finite element method for solving nonlinear scalar conservation equations. The technique is based on a new viscous bilinear form introduced in Guermond and Nazarov [Comput. Methods Appl. Mech. Engrg., 272 (2014), pp. 198-213], a high-order entropy viscosity method, and the Boris-Book-Zalesak flux correction technique. The algorithm works for arbitrary meshes in any space dimension and for all Lipschitz fluxes. The formal second-order accuracy of the method and its convergence properties are tested on a series of linear and nonlinear benchmark problems.

On moduli for which certain second-order linear reczrrebces contain a complete system of residues modulo m

Czech Academy of Sciences Publication Activity Database

Somer, L.; Křížek, Michal

2017-01-01

Roč. 55, č. 3 (2017), s. 209-228 ISSN 0015-0517 Institutional support: RVO:67985840 Keywords : Lucas sequence * second-order * Fibonacci sequence Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics

Second-order perturbations of cosmological fluids: Relativistic effects of pressure, multicomponent, curvature, and rotation

International Nuclear Information System (INIS)

Hwang, Jai-chan; Noh, Hyerim

2007-01-01

We present general relativistic correction terms appearing in Newton's gravity to the second-order perturbations of cosmological fluids. In our previous work we have shown that to the second-order perturbations, the density and velocity perturbation equations of general relativistic zero-pressure, irrotational, single-component fluid in a spatially flat background coincide exactly with the ones known in Newton's theory without using the gravitational potential. We also have shown the effect of gravitational waves to the second order, and pure general relativistic correction terms appearing in the third-order perturbations. Here, we present results of second-order perturbations relaxing all the assumptions made in our previous works. We derive the general relativistic correction terms arising due to (i) pressure, (ii) multicomponent, (iii) background spatial curvature, and (iv) rotation. In the case of multicomponent zero-pressure, irrotational fluids under the flat background, we effectively do not have relativistic correction terms, thus the relativistic equations expressed in terms of density and velocity perturbations again coincide with the Newtonian ones. In the other three cases we generally have pure general relativistic correction terms. In the case of pressure, the relativistic corrections appear even in the level of background and linear perturbation equations. In the presence of background spatial curvature, or rotation, pure relativistic correction terms directly appear in the Newtonian equations of motion of density and velocity perturbations to the second order; to the linear order, without using the gravitational potential (or metric perturbations), we have relativistic/Newtonian correspondences for density and velocity perturbations of a single-component fluid including the rotation even in the presence of background spatial curvature. In the small-scale limit (far inside the horizon), to the second-order, relativistic equations of density and

Integration of differential equations by the pseudo-linear (PL) approximation

International Nuclear Information System (INIS)

Bonalumi, Riccardo A.

1998-01-01

A new method of integrating differential equations was originated with the technique of approximately calculating the integrals called the pseudo-linear (PL) procedure: this method is A-stable. This article contains the following examples: 1st order ordinary differential equations (ODEs), 2nd order linear ODEs, stiff system of ODEs (neutron kinetics), one-dimensional parabolic (diffusion) partial differential equations. In this latter case, this PL method coincides with the Crank-Nicholson method

The non-ordinary Regge behavior of the K{sup *}{sub 0}(800) or κ-meson versus the ordinary K{sup *}{sub 0}(1430)

Energy Technology Data Exchange (ETDEWEB)

Pelaez, J.R.; Rodas, A. [Universidad Complutense de Madrid, Departamento de Fisica Teorica II and UPARCOS, Madrid (Spain)

2017-06-15

The Regge trajectory of an elastic resonance can be calculated from dispersion theory, instead of fitted phenomenologically, using only its pole parameters as input. This also provides a correct treatment of resonance widths in Regge trajectories, essential for very wide resonances. In this work we first calculate the K{sup *}{sub 0}(1430) Regge trajectory, finding the ordinary almost real and linear behavior, typical of q anti q resonances. In contrast, for the K{sup *}{sub 0}(800) meson, the resulting Regge trajectory is non-linear and has a much smaller slope than ordinary resonances, being remarkably similar to that of the f{sub 0}(500) or σ meson. The slope of these unusual Regge trajectories seems to scale with the meson masses rather than with scales typical of quark degrees of freedom. We also calculate the range of the interaction responsible for the formation of these resonances. Our results strongly support a non-ordinary, predominantly meson-meson-like, interpretation for the lightest strange and non-strange resonances. (orig.)

From ordinary to partial differential equations

CERN Document Server

Esposito, Giampiero

2017-01-01

This book is addressed to mathematics and physics students who want to develop an interdisciplinary view of mathematics, from the age of Riemann, Poincaré and Darboux to basic tools of modern mathematics. It enables them to acquire the sensibility necessary for the formulation and solution of difficult problems, with an emphasis on concepts, rigour and creativity. It consists of eight self-contained parts: ordinary differential equations; linear elliptic equations; calculus of variations; linear and non-linear hyperbolic equations; parabolic equations; Fuchsian functions and non-linear equations; the functional equations of number theory; pseudo-differential operators and pseudo-differential equations. The author leads readers through the original papers and introduces new concepts, with a selection of topics and examples that are of high pedagogical value.

Recoil velocity at second post-Newtonian order for spinning black hole binaries

International Nuclear Information System (INIS)

Racine, Etienne; Buonanno, Alessandra; Kidder, Larry

2009-01-01

We compute the flux of linear momentum carried by gravitational waves emitted from spinning binary black holes at second post-Newtonian (2PN) order for generic orbits. In particular we provide explicit expressions of three new types of terms, namely, next-to-leading order spin-orbit terms at 1.5 post-Newtonian (1.5PN) order, spin-orbit tail terms at 2PN order, and spin-spin terms at 2PN order. Restricting ourselves to quasicircular orbits, we integrate the linear-momentum flux over time to obtain the recoil velocity as function of orbital frequency. We find that in the so-called superkick configuration the higher-order spin corrections can increase the recoil velocity up to a factor ∼3 with respect to the leading-order PN prediction. Whereas the recoil velocity computed in PN theory within the adiabatic approximation can accurately describe the early inspiral phase, we find that its fast increase during the late inspiral and plunge, and the arbitrariness in determining until when it should be trusted, makes the PN predictions for the total recoil not very accurate and robust. Nevertheless, the linear-momentum flux at higher PN orders can be employed to build more reliable resummed expressions aimed at capturing the nonperturbative effects until merger. Furthermore, we provide expressions valid for generic orbits, and accurate at 2PN order, for the energy and angular momentum carried by gravitational waves emitted from spinning binary black holes. Specializing to quasicircular orbits we compute the spin-spin terms at 2PN order in the expression for the evolution of the orbital frequency and found agreement with Mikoczi, Vasuth, and Gergely. We also verified that in the limit of extreme mass ratio our expressions for the energy and angular momentum fluxes match the ones of Tagoshi, Shibata, Tanaka, and Sasaki obtained in the context of black hole perturbation theory.

PSsolver: A Maple implementation to solve first order ordinary differential equations with Liouvillian solutions

Science.gov (United States)

Avellar, J.; Duarte, L. G. S.; da Mota, L. A. C. P.

2012-10-01

of revisions: • As time went by, many commands in Maple were deprecated. So, in order to make the program able to run with the newer versions, we have checked and changed some of those. For instance, the command sum had changed, and some program lines were substituted so that the package works properly. • In the old version we must supply the degree of the Darboux polynomials we want to determine. In the present version the user can set the degree by typing Deg = number in the command call (e.g., PSsolve(ode, Deg =3); telling the command PSsolve that it must use Darboux polynomials of degree up to three). If the user does not specify the degree, the routines use, as default, the degree 1. Restrictions: If the integrating factor for the FOODE under consideration has factors of high degree in the dependent and independent variables and in the elementary functions appearing in the FOODE, the package may spend a long time finding the solution. Also, when dealing with FOODEs containing elementary functions, it is essential that the algebraic dependency between them is recognized. If that does not happen, our program can miss some solutions. Unusual features: Our implementation of the Prelle-Singer approach not only solves FOODEs, but can also be used as a research tool that allows the user to follow all the steps of the procedure. For example, the Darboux polynomials (eigenpolynomials) of the D-operator associated with a FOODE (see Section 4) can be calculated. In addition, our package is successful in solving FOODEs that were not solved by some of the most commonly available solvers. Finally, our package implements a theoretical extension (for details, see [1,2]) to the original Prelle-Singer approach that enhances its scope, allowing it to tackle some FOODEs whose solutions involve non-elementary Liouvillian functions. Running time: This depends strongly on the FOODE, but usually under 2 seconds when running our 'arena' test file: The non linear FOODEs presented in

A Unified Introduction to Ordinary Differential Equations

Science.gov (United States)

Lutzer, Carl V.

2006-01-01

This article describes how a presentation from the point of view of differential operators can be used to (partially) unify the myriad techniques in an introductory course in ordinary differential equations by providing students with a powerful, flexible paradigm that extends into (or from) linear algebra. (Contains 1 footnote.)

Geometrical foundations of continuum mechanics an application to first- and second-order elasticity and elasto-plasticity

CERN Document Server

Steinmann, Paul

2015-01-01

This book illustrates the deep roots of the geometrically nonlinear kinematics of generalized continuum mechanics in differential geometry. Besides applications to first- order elasticity and elasto-plasticity an appreciation thereof is particularly illuminating for generalized models of continuum mechanics such as second-order (gradient-type) elasticity and elasto-plasticity.   After a motivation that arises from considering geometrically linear first- and second- order crystal plasticity in Part I several concepts from differential geometry, relevant for what follows, such as connection, parallel transport, torsion, curvature, and metric for holonomic and anholonomic coordinate transformations are reiterated in Part II. Then, in Part III, the kinematics of geometrically nonlinear continuum mechanics are considered. There various concepts of differential geometry, in particular aspects related to compatibility, are generically applied to the kinematics of first- and second- order geometrically nonlinear con...

The linear ordering problem: an algorithm for the optimal solution ...

African Journals Online (AJOL)

In this paper we describe and implement an algorithm for the exact solution of the Linear Ordering problem. Linear Ordering is the problem of finding a linear order of the nodes of a graph such that the sum of the weights which are consistent with this order is as large as possible. It is an NP - Hard combinatorial optimisation ...

Recursive belief manipulation and second-order false-beliefs

DEFF Research Database (Denmark)

Braüner, Torben; Blackburn, Patrick Rowan; Polyanskaya, Irina

2016-01-01

it indicate that a more fundamental *conceptual change* has taken place? In this paper we extend Braüner's hybrid-logical analysis of first-order false-belief tasks to the second-order case, and argue that our analysis supports a version of the conceptual change position.......The literature on first-order false-belief is extensive, but less is known about the second-order case. The ability to handle second-order false-beliefs correctly seems to mark a cognitively significant step, but what is its status? Is it an example of *complexity only* development, or does...

Axisymmetric flow and heat transfer to modified second grade fluid over a radially stretching sheet

Directory of Open Access Journals (Sweden)

Masood Khan

Full Text Available In the present work, an analysis is made to the two-dimensional axisymmetric flow and heat transfer of a modified second grade fluid over an isothermal non-linear radially stretching sheet. The momentum and energy equations are modelled and the boundary layer equations are derived. The governing equations for velocity and temperature are turned down into a system of ordinary differential equations by invoking appropriate transformations which are then solved numerically via fourth and fifth order Runge-Kutta Fehlberg method. Moreover, the influence of the pertinent parameters namely the generalized second grade parameter, stretching parameter, the power-law index and the generalized Prandtl number is graphically portrayed. It is inferred that the generalized second grade parameter uplifted the momentum boundary layer while lessened the thermal boundary layer. Furthermore, the impact of stretching parameter is more pronounced for the second grade fluid (m = 0 in contrast with the power-law fluid (k = 0. For some special cases, comparisons are made with previously reported results and an excellent agreement is established. Keywords: Modified second grade fluid, Axisymmetric flow, Heat transfer, Non-linear stretching sheet

Calculating Second-Order Effects in MOSFET's

Science.gov (United States)

Benumof, Reuben; Zoutendyk, John A.; Coss, James R.

1990-01-01

Collection of mathematical models includes second-order effects in n-channel, enhancement-mode, metal-oxide-semiconductor field-effect transistors (MOSFET's). When dimensions of circuit elements relatively large, effects neglected safely. However, as very-large-scale integration of microelectronic circuits leads to MOSFET's shorter or narrower than 2 micrometer, effects become significant in design and operation. Such computer programs as widely-used "Simulation Program With Integrated Circuit Emphasis, Version 2" (SPICE 2) include many of these effects. In second-order models of n-channel, enhancement-mode MOSFET, first-order gate-depletion region diminished by triangular-cross-section deletions on end and augmented by circular-wedge-cross-section bulges on sides.

Universal formats for nonlinear ordinary differential systems

International Nuclear Information System (INIS)

Kerner, E.H.

1981-01-01

It is shown that very general nonlinear ordinary differential systems (embracing all that arise in practice) may, first, be brought down to polynomial systems (where the nonlinearities occur only as polynomials in the dependent variables) by introducing suitable new variables into the original system; second, that polynomial systems are reducible to ''Riccati systems,'' where the nonlinearities are quadratic at most; third, that Riccati systems may be brought to elemental universal formats containing purely quadratic terms with simple arrays of coefficients that are all zero or unity. The elemental systems have representations as novel types of matrix Riccati equations. Different starting systems and their associated Riccati systems differ from one another, at the final elemental level, in order and in initial data, but not in format

Generation and Identification of Ordinary Differential Equations of Maximal Symmetry Algebra

Directory of Open Access Journals (Sweden)

J. C. Ndogmo

2016-01-01

Full Text Available An effective method for generating linear ordinary differential equations of maximal symmetry in their most general form is found, and an explicit expression for the point transformation reducing the equation to its canonical form is obtained. New expressions for the general solution are also found, as well as several identification and other results and a direct proof of the fact that a linear ordinary differential equation is iterative if and only if it is reducible to the canonical form by a point transformation. New classes of solvable equations parameterized by an arbitrary function are also found, together with simple algebraic expressions for the corresponding general solution.

On sign constant solutions of certain boundary value problems for second-order functional differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander; Vodstrčil, Petr

2005-01-01

Roč. 84, č. 2 (2005), s. 197-209 ISSN 0003-6811 Institutional research plan: CEZ:AV0Z10190503 Keywords : second order linear functional differential equations * nonnegative solution * two-point boundary value problem Subject RIV: BA - General Mathematics http://www.tandfonline.com/doi/full/10.1080/00036810410001724427

Second-Order Risk Constraints in Decision Analysis

Directory of Open Access Journals (Sweden)

Love Ekenberg

2014-01-01

Full Text Available Recently, representations and methods aimed at analysing decision problems where probabilities and values (utilities are associated with distributions over them (second-order representations have been suggested. In this paper we present an approach to how imprecise information can be modelled by means of second-order distributions and how a risk evaluation process can be elaborated by integrating procedures for numerically imprecise probabilities and utilities. We discuss some shortcomings of the use of the principle of maximising the expected utility and of utility theory in general, and offer remedies by the introduction of supplementary decision rules based on a concept of risk constraints taking advantage of second-order distributions.

Scintillation camera with second order resolution

International Nuclear Information System (INIS)

Muehllehner, G.

1976-01-01

A scintillation camera for use in radioisotope imaging to determine the concentration of radionuclides in a two-dimensional area is described in which means is provided for second order positional resolution. The phototubes, which normally provide only a single order of resolution, are modified to provide second order positional resolution of radiation within an object positioned for viewing by the scintillation camera. The phototubes are modified in that multiple anodes are provided to receive signals from the photocathode in a manner such that each anode is particularly responsive to photoemissions from a limited portion of the photocathode. Resolution of radioactive events appearing as an output of this scintillation camera is thereby improved

Scintillation camera with second order resolution

International Nuclear Information System (INIS)

1975-01-01

A scintillation camera is described for use in radioisotope imaging to determine the concentration of radionuclides in a two-dimensional area in which means is provided for second-order positional resolution. The phototubes which normally provide only a single order of resolution, are modified to provide second-order positional resolution of radiation within an object positioned for viewing by the scintillation camera. The phototubes are modified in that multiple anodes are provided to receive signals from the photocathode in a manner such that each anode is particularly responsive to photoemissions from a limited portion of the photocathode. Resolution of radioactive events appearing as an output of this scintillation camera is thereby improved

Periodic oscillations in linear continuous media coupled with nonlinear discrete systems

International Nuclear Information System (INIS)

Lupini, R.

1998-01-01

A general derivation of partial differential equations with boundary conditions in the form of ordinary differential equations is obtained using the principle of stationary action for a Lagrangian function composed of continuous plus discrete parts in interaction across the boundaries of a 1-dimensional medium. This approach leads directly to the theorem of energy conservation. For linear continuous medium, homogeneous Dirichlet condition at one boundary, and nonlinear oscillator at the other boundary, the entire differential problem reduces to a nonlinear differential-difference equation of neutral type and of the second order. The lag parameter is τ = l/c, where c is the phase speed, l the length of the continuum. The Author investigate the problem of the occurrence of periodic solutions of period integer multiple of the lag (super harmonic solutions) in the case of zero inertia of the boundary system. The problem for such oscillations is shown to reduce to systems of ordinary differential equations with matching conditions in a phase space of lower dimensionality: Phase-plane techniques are used to determine solutions of period 4τ, 8τ and 6τ

Investigating local network interactions underlying first- and second-order processing.

Science.gov (United States)

Ellemberg, Dave; Allen, Harriet A; Hess, Robert F

2004-01-01

We compared the spatial lateral interactions for first-order cues to those for second-order cues, and investigated spatial interactions between these two types of cues. We measured the apparent modulation depth of a target Gabor at fixation, in the presence and the absence of horizontally flanking Gabors. The Gabors' gratings were either added to (first-order) or multiplied with (second-order) binary 2-D noise. Apparent "contrast" or modulation depth (i.e., the perceived difference between the high and low luminance regions for the first-order stimulus, or between the high and low contrast regions for the second-order stimulus) was measured with a modulation depth-matching paradigm. For each observer, the first- and second-order Gabors were equated for apparent modulation depth without the flankers. Our results indicate that at the smallest inter-element spacing, the perceived reduction in modulation depth is significantly smaller for the second-order than for the first-order stimuli. Further, lateral interactions operate over shorter distances and the spatial frequency and orientation tuning of the suppression effect are broader for second- than first-order stimuli. Finally, first- and second-order information interact in an asymmetrical fashion; second-order flankers do not reduce the apparent modulation depth of the first-order target, whilst first-order flankers reduce the apparent modulation depth of the second-order target.

Linear Matrix Inequality Based Fuzzy Synchronization for Fractional Order Chaos

Directory of Open Access Journals (Sweden)

Bin Wang

2015-01-01

Full Text Available This paper investigates fuzzy synchronization for fractional order chaos via linear matrix inequality. Based on generalized Takagi-Sugeno fuzzy model, one efficient stability condition for fractional order chaos synchronization or antisynchronization is given. The fractional order stability condition is transformed into a set of linear matrix inequalities and the rigorous proof details are presented. Furthermore, through fractional order linear time-invariant (LTI interval theory, the approach is developed for fractional order chaos synchronization regardless of the system with uncertain parameters. Three typical examples, including synchronization between an integer order three-dimensional (3D chaos and a fractional order 3D chaos, anti-synchronization of two fractional order hyperchaos, and the synchronization between an integer order 3D chaos and a fractional order 4D chaos, are employed to verify the theoretical results.

Second order logic, set theory and foundations of mathematics

NARCIS (Netherlands)

Väänänen, J.A.; Dybjer, P; Lindström, S; Palmgren, E; Sundholm, G

2012-01-01

The question, whether second order logic is a better foundation for mathematics than set theory, is addressed. The main difference between second order logic and set theory is that set theory builds up a transfinite cumulative hierarchy while second order logic stays within one application of the

Decomposition of a symmetric second-order tensor

Science.gov (United States)

Heras, José A.

2018-05-01

In the three-dimensional space there are different definitions for the dot and cross products of a vector with a second-order tensor. In this paper we show how these products can uniquely be defined for the case of symmetric tensors. We then decompose a symmetric second-order tensor into its ‘dot’ part, which involves the dot product, and the ‘cross’ part, which involves the cross product. For some physical applications, this decomposition can be interpreted as one in which the dot part identifies with the ‘parallel’ part of the tensor and the cross part identifies with the ‘perpendicular’ part. This decomposition of a symmetric second-order tensor may be suitable for undergraduate courses of vector calculus, mechanics and electrodynamics.

The second-order decomposition model of nonlinear irregular waves

DEFF Research Database (Denmark)

Yang, Zhi Wen; Bingham, Harry B.; Li, Jin Xuan

2013-01-01

into the first- and the second-order super-harmonic as well as the second-order sub-harmonic components by transferring them into an identical Fourier frequency-space and using a Newton-Raphson iteration method. In order to evaluate the present model, a variety of monochromatic waves and the second...

Investigation of second-order hyperpolarizability of some organic compounds

Science.gov (United States)

Tajalli, H.; Zirak, P.; Ahmadi, S.

2003-04-01

In this work, we have measured the second order hyperpolarizability of some organic materials with (EFISH) method and also calculated the second order hyperpolarizability of 13 organic compound with Mopac6 software and investigated the different factors that affect the amount of second order hyperpolarizability and ways to increase it.

Solution of linear ordinary differential equations by means of the method of variation of arbitrary constants

DEFF Research Database (Denmark)

Mejlbro, Leif

1997-01-01

An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians.......An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians....

Shooting method for third order simultaneous ordinary differential equations with application to magnetohydrodynamic boundary layer

International Nuclear Information System (INIS)

Srivastava, A.C.; Hazarika, G.C.

1990-01-01

An algorithm based on the shooting method has been developed for the solution of a two-point boundary value problem consisting of a system of third order simultaneous ordinary differential equations. The Falkner-Skan equations for electrically conducting viscous fluid with applied magnetic field has been solved by using this algorithm for various values of the wedge angle and magnetic parameters. The shooting method seems to be well convergent for a system as the results are in good agreement with those obtained by other methods. It is observed that both viscous boundary layer and magnetic boundary layer decrease while velocity as well as magnetic field increase with the increase of the wedge angle. (author). 6 tabs., 7 refs

Method of construction of the Riemann function for a second-order hyperbolic equation

Science.gov (United States)

Aksenov, A. V.

2017-12-01

A linear hyperbolic equation of the second order in two independent variables is considered. The Riemann function of the adjoint equation is shown to be invariant with respect to the fundamental solutions transformation group. Symmetries and symmetries of fundamental solutions of the Euler-Poisson-Darboux equation are found. The Riemann function is constructed with the aid of fundamental solutions symmetries. Examples of the application of the algorithm for constructing Riemann function are given.

Second-Order Footsteps Illusions

Directory of Open Access Journals (Sweden)

Akiyoshi Kitaoka

2015-12-01

Full Text Available In the “footsteps illusion”, light and dark squares travel at constant speed across black and white stripes. The squares appear to move faster and slower as their contrast against the stripes varies. We now demonstrate some second-order footsteps illusions, in which all edges are defined by colors or textures—even though luminance-based neural motion detectors are blind to such edges.

Solving polynomial differential equations by transforming them to linear functional-differential equations

OpenAIRE

Nahay, John Michael

2008-01-01

We present a new approach to solving polynomial ordinary differential equations by transforming them to linear functional equations and then solving the linear functional equations. We will focus most of our attention upon the first-order Abel differential equation with two nonlinear terms in order to demonstrate in as much detail as possible the computations necessary for a complete solution. We mention in our section on further developments that the basic transformation idea can be generali...

Second order pedagogy as an example of second order cybernetics

Directory of Open Access Journals (Sweden)

Anne B. Reinertsen

2012-07-01

Full Text Available This article is about seeing/creating/trying out an idea of pedagogy and pedagogical/ educational research in/as/with self-reflexive, circular and diffractive perspectives and about using second order cybernetics as thinking tool. It is a move away from traditional hypothesis driven activities and a move towards data driven pedagogies and research: Teachers, teacher researchers and researchers simultaneously producing and theorizing our practices and ourselves. Deleuzian becomings- eventually becomings with data - theory - theodata is pivotal. It is a move towards a Derridean bricolage. A different science of pedagogy operating as a circular science of self-reflexivity and diffraction in search of quality again and again and again: Theopractical becomings and inspiractionresearch.

Weak value amplification via second-order correlated technique

International Nuclear Information System (INIS)

Cui Ting; Huang Jing-Zheng; Zeng Gui-Hua; Liu Xiang

2016-01-01

We propose a new framework combining weak measurement and second-order correlated technique. The theoretical analysis shows that weak value amplification (WVA) experiment can also be implemented by a second-order correlated system. We then build two-dimensional second-order correlated function patterns for achieving higher amplification factor and discuss the signal-to-noise ratio influence. Several advantages can be obtained by our proposal. For instance, detectors with high resolution are not necessary. Moreover, detectors with low saturation intensity are available in WVA setup. Finally, type-one technical noise can be effectively suppressed. (paper)

An Ordering Linear Unification Algorithm

Institute of Scientific and Technical Information of China (English)

胡运发

1989-01-01

In this paper,we present an ordering linear unification algorithm(OLU).A new idea on substituteion of the binding terms is introduced to the algorithm,which is able to overcome some drawbacks of other algorithms,e.g.,MM algorithm[1],RG1 and RG2 algorithms[2],Particularly,if we use the directed eyclie graphs,the algoritm needs not check the binding order,then the OLU algorithm can also be aplied to the infinite tree data struceture,and a higher efficiency can be expected.The paper focuses upon the discussion of OLU algorithm and a partial order structure with respect to the unification algorithm.This algorithm has been implemented in the GKD-PROLOG/VAX 780 interpreting system.Experimental results have shown that the algorithm is very simple and efficient.

Second-order nonlinear optical metamaterials: ABC-type nanolaminates

International Nuclear Information System (INIS)

Alloatti, L.; Kieninger, C.; Lauermann, M.; Köhnle, K.; Froelich, A.; Wegener, M.; Frenzel, T.; Freude, W.; Leuthold, J.; Koos, C.

2015-01-01

We demonstrate a concept for second-order nonlinear metamaterials that can be obtained from non-metallic centrosymmetric constituents with inherently low optical absorption. The concept is based on iterative atomic-layer deposition of three different materials, A = Al 2 O 3 , B = TiO 2 , and C = HfO 2 . The centrosymmetry of the resulting ABC stack is broken since the ABC and the inverted CBA sequences are not equivalent—a necessary condition for non-zero second-order nonlinearity. In our experiments, we find that the bulk second-order nonlinear susceptibility depends on the density of interfaces, leading to a nonlinear susceptibility of 0.26 pm/V at a wavelength of 800 nm. ABC-type nanolaminates can be deposited on virtually any substrate and offer a promising route towards engineering of second-order optical nonlinearities at both infrared and visible wavelengths

Second-order optical effects in several pyrazolo-quinoline derivatives

International Nuclear Information System (INIS)

Makowska-Janusik, M.; Gondek, E.; Kityk, I.V.; WisIa, J.; Sanetra, J.; Danel, A.

2004-01-01

Using optical poling of several pyazolo-quinoline (PAQ) derivatives we have found an existence of sufficiently high second order optical susceptibility at wavelength 1.76 μm varying in the range 0.9-2.8 pm/V. The performed quantum chemical simulations of the UV-absorption for isolated, solvated and incorporated into the polymethacrylate (PMMA) polymer films have shown that the PM3 method is the best among the semi-empirical ones to simulate the optical properties. The calculations of the hyperpolarizabilites have shown a good correlation with experimentally measured susceptibilities obtained from the optical poling. We have found that experimental susceptibility depends on linear molecular polarizability and photoinducing changes of the molecular dipole moment. It is clearly seen for the PAQ4-PAQ6 molecules possessing halogen atoms with relatively large polarizabilities

Second-order optical effects in several pyrazolo-quinoline derivatives

Energy Technology Data Exchange (ETDEWEB)

Makowska-Janusik, M. [Solid State Department, Institute of Physics, WSP Czestochowa, Al. Armii Krajowej 13/15, Czestochowa PL42201 (Poland); Gondek, E. [Institute of Physics, Cracow University of Technology, ul. Podchorazych 1, 30-084 (Poland); Kityk, I.V. [Department of Biology and Biophysics, Technical University of Czestochowa, Al. Armii Krajowej 36, Czestochowa PL-42210 (Poland)]. E-mail: i.kityk@wsp.czest.pl; WisIa, J. [Departament of Chemistry, Hugon Kollataj Agricultural University, Al. Mickiewicza 24/28, 30-059 Cracow (Poland); Sanetra, J. [Institute of Physics, Cracow University of Technology, ul. Podchorazych 1, 30-084 (Poland); Danel, A. [Department of Chemistry, Hugon Kollataj Agricultural University, Al. Mickiewicza 24/28, 30-059 Cracow (Poland)

2004-11-15

Using optical poling of several pyazolo-quinoline (PAQ) derivatives we have found an existence of sufficiently high second order optical susceptibility at wavelength 1.76 {mu}m varying in the range 0.9-2.8 pm/V. The performed quantum chemical simulations of the UV-absorption for isolated, solvated and incorporated into the polymethacrylate (PMMA) polymer films have shown that the PM3 method is the best among the semi-empirical ones to simulate the optical properties. The calculations of the hyperpolarizabilites have shown a good correlation with experimentally measured susceptibilities obtained from the optical poling. We have found that experimental susceptibility depends on linear molecular polarizability and photoinducing changes of the molecular dipole moment. It is clearly seen for the PAQ4-PAQ6 molecules possessing halogen atoms with relatively large polarizabilities.

Second-order optical effects in several pyrazolo-quinoline derivatives

Science.gov (United States)

Makowska-Janusik, M.; Gondek, E.; Kityk, I. V.; Wisła, J.; Sanetra, J.; Danel, A.

2004-11-01

Using optical poling of several pyazolo-quinoline (PAQ) derivatives we have found an existence of sufficiently high second order optical susceptibility at wavelength 1.76 μm varying in the range 0.9-2.8 pm/V. The performed quantum chemical simulations of the UV-absorption for isolated, solvated and incorporated into the polymethacrylate (PMMA) polymer films have shown that the PM3 method is the best among the semi-empirical ones to simulate the optical properties. The calculations of the hyperpolarizabilites have shown a good correlation with experimentally measured susceptibilities obtained from the optical poling. We have found that experimental susceptibility depends on linear molecular polarizability and photoinducing changes of the molecular dipole moment. It is clearly seen for the PAQ4-PAQ6 molecules possessing halogen atoms with relatively large polarizabilities.

Kubo Formulas for Second-Order Hydrodynamic Coefficients

International Nuclear Information System (INIS)

Moore, Guy D.; Sohrabi, Kiyoumars A.

2011-01-01

At second order in gradients, conformal relativistic hydrodynamics depends on the viscosity η and on five additional ''second-order'' hydrodynamical coefficients τ Π , κ, λ 1 , λ 2 , and λ 3 . We derive Kubo relations for these coefficients, relating them to equilibrium, fully retarded three-point correlation functions of the stress tensor. We show that the coefficient λ 3 can be evaluated directly by Euclidean means and does not in general vanish.

Controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with delay and Poisson jumps

Directory of Open Access Journals (Sweden)

Diem Dang Huan

2015-12-01

Full Text Available The current paper is concerned with the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps in Hilbert spaces. Using the theory of a strongly continuous cosine family of bounded linear operators, stochastic analysis theory and with the help of the Banach fixed point theorem, we derive a new set of sufficient conditions for the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps. Finally, an application to the stochastic nonlinear wave equation with infinite delay and Poisson jumps is given.

Second order nonlinear optical properties of zinc oxide films deposited by low temperature dual ion beam sputtering

International Nuclear Information System (INIS)

Larciprete, M.C.; Passeri, D.; Michelotti, F.; Paoloni, S.; Sibilia, C.; Bertolotti, M.; Belardini, A.; Sarto, F.; Somma, F.; Lo Mastro, S.

2005-01-01

We investigated second order optical nonlinearity of zinc oxide thin films, grown on glass substrates by the dual ion beam sputtering technique under different deposition conditions. Linear optical characterization of the films was carried out by spectrophotometric optical transmittance and reflectance measurements, giving the complex refractive index dispersion. Resistivity of the films was determined using the four-point probe sheet resistance method. Second harmonic generation measurements were performed by means of the Maker fringes technique where the fundamental beam was originated by nanosecond laser at λ=1064 nm. We found a relatively high nonlinear optical response, and evidence of a dependence of the nonlinear coefficient on the deposition parameters for each sample. Moreover, the crystalline properties of the films were investigated by x-ray diffraction measurements and correlation with second order nonlinearity were analyzed. Finally, we investigated the influence of the oxygen flow rate during the deposition process on both the second order nonlinearity and the structural properties of the samples

Third-order ordinary differential equations Y”' = f(x, y, y'', y′”) with ...

African Journals Online (AJOL)

dimensional symmetry algebra. Mathematics Subject Classication (2010): 34A05, 34A25, 53A55, 76M60. Key words: Linearization, third order ODEs, point transformation, contact transformation, Lie symmetries, relative differential invariants.

Linearly Ordered Attribute Grammar Scheduling Using SAT-Solving

NARCIS (Netherlands)

Bransen, Jeroen; van Binsbergen, L.Thomas; Claessen, Koen; Dijkstra, Atze

2015-01-01

Many computations over trees can be specified using attribute grammars. Compilers for attribute grammars need to find an evaluation order (or schedule) in order to generate efficient code. For the class of linearly ordered attribute grammars such a schedule can be found statically, but this problem

Estimation of ibuprofen and famotidine in tablets by second order derivative spectrophotometery method

Directory of Open Access Journals (Sweden)

Dimal A. Shah

2017-02-01

Full Text Available A simple and accurate method for the analysis of ibuprofen (IBU and famotidine (FAM in their combined dosage form was developed using second order derivative spectrophotometery. IBU and FAM were quantified using second derivative responses at 272.8 nm and 290 nm in the spectra of their solutions in methanol. The calibration curves were linear in the concentration range of 100–600 μg/mL for IBU and 5–25 μg/mL for FAM. The method was validated and found to be accurate and precise. Developed method was successfully applied for the estimation of IBU and FAM in their combined dosage form.

Second-order impartiality and public sphere

Directory of Open Access Journals (Sweden)

Sládeček Michal

2016-01-01

Full Text Available In the first part of the text the distinction between first- and second-order impartiality, along with Brian Barry’s thorough elaboration of their characteristics and the differences between them, is examined. While the former impartiality is related to non-favoring fellow-persons in everyday occasions, the latter is manifested in the institutional structure of society and its political and public morality. In the second part of the article, the concept of public impartiality is introduced through analysis of two examples. In the first example, a Caledonian Club with its exclusive membership is considered as a form of association which is partial, but nevertheless morally acceptable. In the second example, the so-called Heinz dilemma has been reconsidered and the author points to some flaws in Barry’s interpretation, arguing that Heinz’s right of giving advantage to his wife’s life over property rights can be recognized through mitigating circum-stances, and this partiality can be appreciated in the public sphere. Thus, public impartiality imposes limits to the restrictiveness and rigidity of political impartiality implied in second-order morality. [Projekat Ministarstva nauke Republike Srbije, br. 179049

Generalized principal resonance in oscillatory systems of second order; Resonancia principal generalizada en sistemas oscilatorios de segundo orden

Energy Technology Data Exchange (ETDEWEB)

Munoz Aguirre, E. [Universidad Autonoma de Puebla, Oaxaca (Mexico); Alexandrov, V. V. [Benemerita Universidad Autonoma de Puebla, Puebla (Mexico)

2001-02-01

This paper will describe the generalized principal resonance of systems as described by the second order of ordinary differential equations and proved by the Pontriaguin maximal principle to coincide with the lengthened solution of an external problem of the same system. The results are verified in special cases of general resonance and parametric resonance for the Mathieu equation. [Spanish] En el presente articulo se estudia la resonancia principal generalizada para sistemas descritos por ecuaciones diferenciales ordinarias de segundo orden y se demuestra con ayuda del principio del maximo de Pontriaguin, la coincidencia de esta con la solucion prolongada de un problema extremal para el mismo sistema. Ademas se verifican estos resultados en los casos particulares de resonancia general y resonancia parametrica para la ecuacion de Mathieu.

Hierarchy among Automata on Linear Orderings

OpenAIRE

Bruyère , Véronique; Carton , Olivier

2005-01-01

In a preceding paper, automata and rational expressions have been introduced for words indexed by linear orderings, together with a Kleene-like theorem. We here pursue this work by proposing a hierarchy among the rational sets. Each class of the hierarchy is defined by a subset of the rational operations that can be used. We then characterize any class by an appropriate class of automata, leading to a Kleene theorem inside the class. A characterization by particular classes of orderings is al...

Fluid-structure interactions in one-dimensional linear cases

International Nuclear Information System (INIS)

Schumann, U.

1979-01-01

The interaction of pressure waves in a pipe with an elastic endwall (piston) is analyzed using a linear ('acoustic') model. Two transient and two periodic cases are investigated. In the transient cases the motions are initiated by either a sudden pressure drop at the opeen end (breaking membrane) or by a sudden release of the piston from a non-equilibrium position ('snapback'); in the latter case the other end of the pipe is closed. In the periodic cases harmonic oscillations of the piston and the fluid are investigated with the other end of the pipe being either closed or open (kept at constant pressure). The problem is characterized by three non-dimensional numbers (e.g.: Mach-, Strouhal-, and an interaction-number). The solution of the wave equation for the pressure accounting for the coupling to the structure can be reduced analytically to the problem of integrating one ordinary differential equation of second order in time. This differential equation in turn can be integrated analytically at least for a certain time period. At later times this ordinary differential equation is integrated numerically. For the periodic cases eigenvalue-problems arise with an infinite number of solutions. The first few eigensolutions are given. (orig./RW) [de

Abnormal Waves Modelled as Second-order Conditional Waves

DEFF Research Database (Denmark)

Jensen, Jørgen Juncher

2005-01-01

The paper presents results for the expected second order short-crested wave conditional of a given wave crest at a specific point in time and space. The analysis is based on the second order Sharma and Dean shallow water wave theory. Numerical results showing the importance of the spectral densit...

'Second' Ehrenfest equation for second order phase transition under hydrostatic pressure

Science.gov (United States)

Moin, Ph. B.

2018-02-01

It is shown that the fundamental conditions for the second-order phase transitions ? and ?, from which the two Ehrenfest equations follow (the 'usual' and the 'second' ones), are realised only at zero hydrostatic pressure (?). At ? the volume jump ΔV at the transition is proportional to the pressure and to the jump of the compressibility ΔζV, whereas the entropy jump ΔS is proportional to the pressure and to the jump of the thermal expansion coefficient ΔαV. This means that at non-zero hydrostatic pressure the phase transition is of the first order and is described by the Clausius-Clapeyron equation. At small pressure this equation coincides with the 'second' Ehrenfest equation ?. At high P, the Clausius-Clapeyron equation describes qualitatively the caused by the crystal compression positive curvature of the ? dependence.

Second-Order Conditioning in "Drosophila"

Science.gov (United States)

Tabone, Christopher J.; de Belle, J. Steven

2011-01-01

Associative conditioning in "Drosophila melanogaster" has been well documented for several decades. However, most studies report only simple associations of conditioned stimuli (CS, e.g., odor) with unconditioned stimuli (US, e.g., electric shock) to measure learning or establish memory. Here we describe a straightforward second-order conditioning…

First-principles prediction of optical second-order harmonic generation in the endohedral N-C{sub 60} compound

Energy Technology Data Exchange (ETDEWEB)

Zhang, G. P.; Strubbe, David A.; Louie, Steven G.; George, Thomas F. [Department of Physics, Indiana State University, Terre Haute, Indiana 47809 (United States); Department of Physics, University of California, Berkeley, California 94720 (United States) and Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720 (United States); Department of Chemistry and Biochemistry and Department of Physics and Astronomy, Office of the Chancellor and Center for Nanoscience, University of Missouri-St. Louis, St. Louis, Missouri 63121 (United States)

2011-08-15

Non-linear-optical properties in C{sub 60} have attracted enormous attention for over two decades. The endohedral complex N-C{sub 60}, with its remarkable thermal stability and spin-quartet ground state, is a candidate for future room-temperature quantum computing, but there has been no investigation of its non-linear-optical properties. Here, a first-principles calculation shows that N-C{sub 60} is a promising material for nanoscale and ultrafast modulations. Excitation by a pump laser pulse of the nitrogen-atom vibration inside the C{sub 60} cage transiently breaks inversion symmetry and can enable second-harmonic generation (SHG) from a probe pulse. Unlike the SHG observed in C{sub 60} thin films, this harmonic signal is switched on and off periodically every 345 fs. For an fcc crystal of N-C{sub 60}, the second-order susceptibility {chi}{sup (2)} is on the order of 10{sup -8} esu, similar to commercially used nonlinear materials.

Second order sliding mode control for a quadrotor UAV.

Science.gov (United States)

Zheng, En-Hui; Xiong, Jing-Jing; Luo, Ji-Liang

2014-07-01

A method based on second order sliding mode control (2-SMC) is proposed to design controllers for a small quadrotor UAV. For the switching sliding manifold design, the selection of the coefficients of the switching sliding manifold is in general a sophisticated issue because the coefficients are nonlinear. In this work, in order to perform the position and attitude tracking control of the quadrotor perfectly, the dynamical model of the quadrotor is divided into two subsystems, i.e., a fully actuated subsystem and an underactuated subsystem. For the former, a sliding manifold is defined by combining the position and velocity tracking errors of one state variable, i.e., the sliding manifold has two coefficients. For the latter, a sliding manifold is constructed via a linear combination of position and velocity tracking errors of two state variables, i.e., the sliding manifold has four coefficients. In order to further obtain the nonlinear coefficients of the sliding manifold, Hurwitz stability analysis is used to the solving process. In addition, the flight controllers are derived by using Lyapunov theory, which guarantees that all system state trajectories reach and stay on the sliding surfaces. Extensive simulation results are given to illustrate the effectiveness of the proposed control method. Copyright © 2014 ISA. Published by Elsevier Ltd. All rights reserved.

First and second order operator splitting methods for the phase field crystal equation

International Nuclear Information System (INIS)

Lee, Hyun Geun; Shin, Jaemin; Lee, June-Yub

2015-01-01

In this paper, we present operator splitting methods for solving the phase field crystal equation which is a model for the microstructural evolution of two-phase systems on atomic length and diffusive time scales. A core idea of the methods is to decompose the original equation into linear and nonlinear subequations, in which the linear subequation has a closed-form solution in the Fourier space. We apply a nonlinear Newton-type iterative method to solve the nonlinear subequation at the implicit time level and thus a considerably large time step can be used. By combining these subequations, we achieve the first- and second-order accuracy in time. We present numerical experiments to show the accuracy and efficiency of the proposed methods

A second-order class-D audio amplifier

OpenAIRE

Cox, Stephen M.; Tan, M.T.; Yu, J.

2011-01-01

Class-D audio amplifiers are particularly efficient, and this efficiency has led to their ubiquity in a wide range of modern electronic appliances. Their output takes the form of a high-frequency square wave whose duty cycle (ratio of on-time to off-time) is modulated at low frequency according to the audio signal. A mathematical model is developed here for a second-order class-D amplifier design (i.e., containing one second-order integrator) with negative feedback. We derive exact expression...

Second-order wave diffraction by a circular cylinder using scaled boundary finite element method

International Nuclear Information System (INIS)

Song, H; Tao, L

2010-01-01

The scaled boundary finite element method (SBFEM) has achieved remarkable success in structural mechanics and fluid mechanics, combing the advantage of both FEM and BEM. Most of the previous works focus on linear problems, in which superposition principle is applicable. However, many physical problems in the real world are nonlinear and are described by nonlinear equations, challenging the application of the existing SBFEM model. A popular idea to solve a nonlinear problem is decomposing the nonlinear equation to a number of linear equations, and then solves them individually. In this paper, second-order wave diffraction by a circular cylinder is solved by SBFEM. By splitting the forcing term into two parts, the physical problem is described as two second-order boundary-value problems with different asymptotic behaviour at infinity. Expressing the velocity potentials as a series of depth-eigenfunctions, both of the 3D boundary-value problems are decomposed to a number of 2D boundary-value sub-problems, which are solved semi-analytically by SBFEM. Only the cylinder boundary is discretised with 1D curved finite-elements on the circumference of the cylinder, while the radial differential equation is solved completely analytically. The method can be extended to solve more complex wave-structure interaction problems resulting in direct engineering applications.

Testing for one Generalized Linear Single Order Parameter

DEFF Research Database (Denmark)

Ellegaard, Niels Langager; Christensen, Tage Emil; Dyre, Jeppe

We examine a linear single order parameter model for thermoviscoelastic relaxation in viscous liquids, allowing for a distribution of relaxation times. In this model the relaxation of volume and entalpy is completely described by the relaxation of one internal order parameter. In contrast to prior...... work the order parameter may be chosen to have a non-exponential relaxation. The model predictions contradict the general consensus of the properties of viscous liquids in two ways: (i) The model predicts that following a linear isobaric temperature step, the normalized volume and entalpy relaxation...... responses or extrapolate from measurements of a glassy state away from equilibrium. Starting from a master equation description of inherent dynamics, we calculate the complex thermodynamic response functions. We device a way of testing for the generalized single order parameter model by measuring 3 complex...

First order and second order fermi acceleration of energetic charged particles by shock waves

International Nuclear Information System (INIS)

Webb, G.M.

1983-01-01

Steady state solutions of the cosmic ray transport equation describing first order Fermi acceleration of energetic charged particles at a plane shock (without losses) and second order Fermi acceleration in the downstream region of the shock are derived. The solutions for the isotropic part of the phase space distribution function are expressible as eigenfunction expansions, being superpositions of series of power law momentum spectra, with the power law indices being the roots of an eigenvalue equation. The above exact analytic solutions are for the case where the spatial diffusion coefficient kappa is independent of momentum. The solutions in general depend on the shock compression ratio, the modulation parameters V 1 L/kappa 1 , V 2 L/kappa 2 (V is the plasma velocity, kappa is the energetic particle diffusion coefficient, and L a characteristic length over which second order Fermi acceleration is effective) in the upstream and downstream regions of the shock, respectively, and also on a further dimensionless parameter, zeta, characterizing second order Fermi acceleration. In the limit as zeta→0 (no second order Fermi acceleration) the power law momentum spectrum characteristic of first order Fermi acceleration (depending only on the shock compression ratio) obtained previously is recovered. Perturbation solutions for the case where second order Fermi effects are small, and for realistic diffusion coefficients (kappainfinityp/sup a/, a>0, p = particle momentum), applicable at high momenta, are also obtained

Efficient education policy: A second-order elasticity rule

OpenAIRE

Richter, Wolfram F.

2010-01-01

Assuming a two-period model with endogenous choices of labour, education, and saving, efficient education policy is characterized for a Ramsey-like scenario in which the government is constrained to use linear instruments. It is shown that education should be effectively subsidized if, and only if, the elasticity of the earnings function is increasing in education. The strength of second-best subsidization increases in the elasticity of the elasticity of the earnings function. This second-ord...

Skyrme interaction to second order in nuclear matter

Science.gov (United States)

Kaiser, N.

2015-09-01

Based on the phenomenological Skyrme interaction various density-dependent nuclear matter quantities are calculated up to second order in many-body perturbation theory. The spin-orbit term as well as two tensor terms contribute at second order to the energy per particle. The simultaneous calculation of the isotropic Fermi-liquid parameters provides a rigorous check through the validity of the Landau relations. It is found that published results for these second order contributions are incorrect in most cases. In particular, interference terms between s-wave and p-wave components of the interaction can contribute only to (isospin or spin) asymmetry energies. Even with nine adjustable parameters, one does not obtain a good description of the empirical nuclear matter saturation curve in the low density region 0\\lt ρ \\lt 2{ρ }0. The reason for this feature is the too strong density-dependence {ρ }8/3 of several second-order contributions. The inclusion of the density-dependent term \\frac{1}{6}{t}3{ρ }1/6 is therefore indispensable for a realistic description of nuclear matter in the Skyrme framework.

Inverse problems in ordinary differential equations and applications

CERN Document Server

Llibre, Jaume

2016-01-01

This book is dedicated to study the inverse problem of ordinary differential equations, that is it focuses in finding all ordinary differential equations that satisfy a given set of properties. The Nambu bracket is the central tool in developing this approach. The authors start characterizing the ordinary differential equations in R^N which have a given set of partial integrals or first integrals. The results obtained are applied first to planar polynomial differential systems with a given set of such integrals, second to solve the 16th Hilbert problem restricted to generic algebraic limit cycles, third for solving the inverse problem for constrained Lagrangian and Hamiltonian mechanical systems, fourth for studying the integrability of a constrained rigid body. Finally the authors conclude with an analysis on nonholonomic mechanics, a generalization of the Hamiltonian principle, and the statement an solution of the inverse problem in vakonomic mechanics.

On holographic entanglement entropy with second order excitations

Science.gov (United States)

He, Song; Sun, Jia-Rui; Zhang, Hai-Qing

2018-03-01

We study the low-energy corrections to the holographic entanglement entropy (HEE) in the boundary CFT by perturbing the bulk geometry up to second order excitations. Focusing on the case that the boundary subsystem is a strip, we show that the area of the bulk minimal surface can be expanded in terms of the conserved charges, such as mass, angular momentum and electric charge of the AdS black brane. We also calculate the variation of the energy in the subsystem and verify the validity of the first law-like relation of thermodynamics at second order. Moreover, the HEE is naturally bounded at second order perturbations if the cosmic censorship conjecture for the dual black hole still holds.

First- and second-order charged particle optics

International Nuclear Information System (INIS)

Brown, K.L.; Servranckx, R.V.

1984-07-01

Since the invention of the alternating gradient principle there has been a rapid evolution of the mathematics and physics techniques applicable to charged particle optics. In this publication we derive a differential equation and a matrix algebra formalism valid to second-order to present the basic principles governing the design of charged particle beam transport systems. A notation first introduced by John Streib is used to convey the essential principles dictating the design of such beam transport systems. For example the momentum dispersion, the momentum resolution, and all second-order aberrations are expressed as simple integrals of the first-order trajectories (matrix elements) and of the magnetic field parameters (multipole components) characterizing the system. 16 references, 30 figures

Second-Order Learning Methods for a Multilayer Perceptron

International Nuclear Information System (INIS)

Ivanov, V.V.; Purehvdorzh, B.; Puzynin, I.V.

1994-01-01

First- and second-order learning methods for feed-forward multilayer neural networks are studied. Newton-type and quasi-Newton algorithms are considered and compared with commonly used back-propagation algorithm. It is shown that, although second-order algorithms require enhanced computer facilities, they provide better convergence and simplicity in usage. 13 refs., 2 figs., 2 tabs

A textbook on ordinary differential equations

CERN Document Server

Ahmad, Shair

2015-01-01

This book offers readers a primer on the theory and applications of Ordinary Differential Equations. The style used is simple, yet thorough and rigorous. Each chapter ends with a broad set of exercises that range from the routine to the more challenging and thought-provoking. Solutions to selected exercises can be found at the end of the book. The book contains many interesting examples on topics such as electric circuits, the pendulum equation, the logistic equation, the Lotka-Volterra system, the Laplace Transform, etc., which introduce students to a number of interesting aspects of the theory and applications. The work is mainly intended for students of Mathematics, Physics, Engineering, Computer Science and other areas of the natural and social sciences that use ordinary differential equations, and who have a firm grasp of Calculus and a minimal understanding of the basic concepts used in Linear Algebra. It also studies a few more advanced topics, such as Stability Theory and Boundary Value Problems, whic...

Comparison of third-order plasma wave echoes with ballistic second-order plasma wave echoes

International Nuclear Information System (INIS)

Leppert, H.D.; Schuelter, H.; Wiesemann, K.

1982-01-01

The apparent dispersion of third-order plasma wave echoes observed in a high frequency plasma is compared with that of simultaneously observed ballistic second-order echoes. Amplitude and wavelength of third-order echoes are found to be always smaller than those of second-order echoes, however, the dispersion curves of both types of echoes are very similar. These observations are in qualitative agreement with calculations of special ballistic third-order echoes. The ballistic nature of the observed third-order echoes may, therefore, be concluded from these measurements. (author)

Linear Parameter Varying Versus Linear Time Invariant Reduced Order Controller Design of Turboprop Aircraft Dynamics

Directory of Open Access Journals (Sweden)

Widowati

2012-07-01

Full Text Available The applicability of parameter varying reduced order controllers to aircraft model is proposed. The generalization of the balanced singular perturbation method of linear time invariant (LTI system is used to reduce the order of linear parameter varying (LPV system. Based on the reduced order model the low-order LPV controller is designed by using synthesis technique. The performance of the reduced order controller is examined by applying it to lateral-directional control of aircraft model having 20th order. Furthermore, the time responses of the closed loop system with reduced order LPV controllers and reduced order LTI controller is compared. From the simulation results, the 8th order LPV controller can maintain stability and to provide the same level of closed-loop systems performance as the full-order LPV controller. It is different with the reduced-order LTI controller that cannot maintain stability and performance for all allowable parameter trajectories.

Second order time evolution of the multigroup diffusion and P1 equations for radiation transport

International Nuclear Information System (INIS)

Olson, Gordon L.

2011-01-01

Highlights: → An existing multigroup transport algorithm is extended to be second-order in time. → A new algorithm is presented that does not require a grey acceleration solution. → The two algorithms are tested with 2D, multi-material problems. → The two algorithms have comparable computational requirements. - Abstract: An existing solution method for solving the multigroup radiation equations, linear multifrequency-grey acceleration, is here extended to be second order in time. This method works for simple diffusion and for flux-limited diffusion, with or without material conduction. A new method is developed that does not require the solution of an averaged grey transport equation. It is effective solving both the diffusion and P 1 forms of the transport equation. Two dimensional, multi-material test problems are used to compare the solution methods.

A second order anti-diffusive Lagrange-remap scheme for two-component flows

Directory of Open Access Journals (Sweden)

Lagoutière Frédéric

2011-11-01

Full Text Available We build a non-dissipative second order algorithm for the approximate resolution of the one-dimensional Euler system of compressible gas dynamics with two components. The considered model was proposed in [1]. The algorithm is based on [8] which deals with a non-dissipative first order resolution in Lagrange-remap formalism. In the present paper we describe, in the same framework, an algorithm that is second order accurate in time and space, and that preserves sharp interfaces. Numerical results reported at the end of the paper are very encouraging, showing the interest of the second order accuracy for genuinely non-linear waves. Nous construisons un algorithme d’ordre deux et non dissipatif pour la résolution approchée des équations d’Euler de la dynamique des gaz compressibles à deux constituants en dimension un. Le modèle que nous considérons est celui à cinq équations proposé et analysé dans [1]. L’algorithme est basé sur [8] qui proposait une résolution approchée à l’ordre un et non dissipative au moyen d’un splitting de type Lagrange-projection. Dans le présent article, nous décrivons, dans le même formalisme, un algorithme d’ordre deux en temps et en espace, qui préserve des interfaces « parfaites » entre les constituants. Les résultats numériques rapportés à la fin de l’article sont très encourageants ; ils montrent clairement les avantages d’un schéma d’ordre deux pour les ondes vraiment non linéaires.

The known unknowns: neural representation of second-order uncertainty, and ambiguity

Science.gov (United States)

Bach, Dominik R.; Hulme, Oliver; Penny, William D.; Dolan, Raymond J.

2011-01-01

Predictions provided by action-outcome probabilities entail a degree of (first-order) uncertainty. However, these probabilities themselves can be imprecise and embody second-order uncertainty. Tracking second-order uncertainty is important for optimal decision making and reinforcement learning. Previous functional magnetic resonance imaging investigations of second-order uncertainty in humans have drawn on an economic concept of ambiguity, where action-outcome associations in a gamble are either known (unambiguous) or completely unknown (ambiguous). Here, we relaxed the constraints associated with a purely categorical concept of ambiguity and varied the second-order uncertainty of gambles continuously, quantified as entropy over second-order probabilities. We show that second-order uncertainty influences decisions in a pessimistic way by biasing second-order probabilities, and that second-order uncertainty is negatively correlated with posterior cingulate cortex activity. The category of ambiguous (compared to non-ambiguous) gambles also biased choice in a similar direction, but was associated with distinct activation of a posterior parietal cortical area; an activation that we show reflects a different computational mechanism. Our findings indicate that behavioural and neural responses to second-order uncertainty are distinct from those associated with ambiguity and may call for a reappraisal of previous data. PMID:21451019

On holographic entanglement entropy with second order excitations

Directory of Open Access Journals (Sweden)

Song He

2018-03-01

Full Text Available We study the low-energy corrections to the holographic entanglement entropy (HEE in the boundary CFT by perturbing the bulk geometry up to second order excitations. Focusing on the case that the boundary subsystem is a strip, we show that the area of the bulk minimal surface can be expanded in terms of the conserved charges, such as mass, angular momentum and electric charge of the AdS black brane. We also calculate the variation of the energy in the subsystem and verify the validity of the first law-like relation of thermodynamics at second order. Moreover, the HEE is naturally bounded at second order perturbations if the cosmic censorship conjecture for the dual black hole still holds.

Parasupersymmetry and N-fold supersymmetry in quantum many-body systems. I: General formalism and second order

International Nuclear Information System (INIS)

Tanaka, Toshiaki

2007-01-01

We propose an elegant formulation of parafermionic algebra and parasupersymmetry of arbitrary order in quantum many-body systems without recourse to any specific matrix representation of parafermionic operators and any kind of deformed algebra. Within our formulation, we show generically that every parasupersymmetric quantum system of order p consists of N-fold supersymmetric pairs with N≤p and thus has weak quasi-solvability and isospectral property. We also propose a new type of non-linear supersymmetries, called quasi-parasupersymmetry, which is less restrictive than parasupersymmetry and is different from N-fold supersymmetry even in one-body systems though the conserved charges are represented by higher-order linear differential operators. To illustrate how our formulation works, we construct second-order parafermionic algebra and three simple examples of parasupersymmetric quantum systems of order 2, one is essentially equivalent to the one-body Rubakov-Spiridonov type and the others are two-body systems in which two supersymmetries are folded. In particular, we show that the first model admits a generalized 2-fold superalgebra

Linear matrix differential equations of higher-order and applications

Directory of Open Access Journals (Sweden)

Mustapha Rachidi

2008-07-01

Full Text Available In this article, we study linear differential equations of higher-order whose coefficients are square matrices. The combinatorial method for computing the matrix powers and exponential is adopted. New formulas representing auxiliary results are obtained. This allows us to prove properties of a large class of linear matrix differential equations of higher-order, in particular results of Apostol and Kolodner are recovered. Also illustrative examples and applications are presented.

Ordinary differential equations

CERN Document Server

Miller, Richard K

1982-01-01

Ordinary Differential Equations is an outgrowth of courses taught for a number of years at Iowa State University in the mathematics and the electrical engineering departments. It is intended as a text for a first graduate course in differential equations for students in mathematics, engineering, and the sciences. Although differential equations is an old, traditional, and well-established subject, the diverse backgrounds and interests of the students in a typical modern-day course cause problems in the selection and method of presentation of material. In order to compensate for this diversity,

Conformal symmetry and non-relativistic second-order fluid dynamics

International Nuclear Information System (INIS)

Chao Jingyi; Schäfer, Thomas

2012-01-01

We study the constraints imposed by conformal symmetry on the equations of fluid dynamics at second order in the gradients of the hydrodynamic variables. At zeroth order, conformal symmetry implies a constraint on the equation of state, E 0 =2/3 P, where E 0 is the energy density and P is the pressure. At first order, conformal symmetry implies that the bulk viscosity must vanish. We show that at second order, conformal invariance requires that two-derivative terms in the stress tensor must be traceless, and that it determines the relaxation of dissipative stresses to the Navier–Stokes form. We verify these results by solving the Boltzmann equation at second order in the gradient expansion. We find that only a subset of the terms allowed by conformal symmetry appear. - Highlights: ► We derive conformal constraints for the stress tensor of a scale invariant fluid. ► We determine the relaxation time in kinetic theory. ► We compute the rate of entropy production in second-order fluid dynamics.

Solutions manual to accompany Ordinary differential equations

CERN Document Server

Greenberg, Michael D

2014-01-01

Features a balance between theory, proofs, and examples and provides applications across diverse fields of study Ordinary Differential Equations presents a thorough discussion of first-order differential equations and progresses to equations of higher order. The book transitions smoothly from first-order to higher-order equations, allowing readers to develop a complete understanding of the related theory. Featuring diverse and interesting applications from engineering, bioengineering, ecology, and biology, the book anticipates potential difficulties in understanding the various solution steps

Experimental verification of a self-consistent theory of the first-, second-, and third-order (non)linear optical response

International Nuclear Information System (INIS)

Perez-Moreno, Javier; Hung, Sheng-Ting; Kuzyk, Mark G.; Zhou, Juefei; Ramini, Shiva K.; Clays, Koen

2011-01-01

We show that a combination of linear absorption spectroscopy, hyper-Rayleigh scattering, and a theoretical analysis using sum rules to reduce the size of the parameter space leads to a prediction of the imaginary part of the second hyperpolarizability of the dye AF-455 that agrees with the experimental data gathered through two-photon absorption spectroscopy. Our procedure, which demands self-consistency between several measurement techniques and does not use adjustable parameters, provides a means for determining transition moments between the dominant excited states based strictly on experimental characterization. This is made possible by our new approach that uses sum rules and molecular symmetry to rigorously reduce the number of required physical quantities.

Second-Order Assortative Mixing in Social Networks

DEFF Research Database (Denmark)

Zhou, Shi; Cox, Ingemar; Hansen, Lars Kai

2017-01-01

In a social network, the number of links of a node, or node degree, is often assumed as a proxy for the node’s importance or prominence within the network. It is known that social networks exhibit the (first-order) assortative mixing, i.e. if two nodes are connected, they tend to have similar node...... degrees, suggesting that people tend to mix with those of comparable prominence. In this paper, we report the second-order assortative mixing in social networks. If two nodes are connected, we measure the degree correlation between their most prominent neighbours, rather than between the two nodes...... themselves. We observe very strong second-order assortative mixing in social networks, often significantly stronger than the first-order assortative mixing. This suggests that if two people interact in a social network, then the importance of the most prominent person each knows is very likely to be the same...

Numerical method for solving linear Fredholm fuzzy integral equations of the second kind

Energy Technology Data Exchange (ETDEWEB)

Abbasbandy, S. [Department of Mathematics, Imam Khomeini International University, P.O. Box 288, Ghazvin 34194 (Iran, Islamic Republic of)]. E-mail: saeid@abbasbandy.com; Babolian, E. [Faculty of Mathematical Sciences and Computer Engineering, Teacher Training University, Tehran 15618 (Iran, Islamic Republic of); Alavi, M. [Department of Mathematics, Arak Branch, Islamic Azad University, Arak 38135 (Iran, Islamic Republic of)

2007-01-15

In this paper we use parametric form of fuzzy number and convert a linear fuzzy Fredholm integral equation to two linear system of integral equation of the second kind in crisp case. We can use one of the numerical method such as Nystrom and find the approximation solution of the system and hence obtain an approximation for fuzzy solution of the linear fuzzy Fredholm integral equations of the second kind. The proposed method is illustrated by solving some numerical examples.

SECOND-ORDER SOLUTIONS OF COSMOLOGICAL PERTURBATION IN THE MATTER-DOMINATED ERA

International Nuclear Information System (INIS)

Hwang, Jai-chan; Noh, Hyerim; Gong, Jinn-Ouk

2012-01-01

We present the growing mode solutions of cosmological perturbations to the second order in the matter-dominated era. We also present several gauge-invariant combinations of perturbation variables to the second order in the most general fluid context. Based on these solutions, we study the Newtonian correspondence of relativistic perturbations to the second order. In addition to the previously known exact relativistic/Newtonian correspondence of density and velocity perturbations to the second order in the comoving gauge, here we show that in the sub-horizon limit we have the correspondences for density, velocity, and potential perturbations in the zero-shear gauge and in the uniform-expansion gauge to the second order. Density perturbation in the uniform-curvature gauge also shows the correspondence to the second order in the sub-horizon scale. We also identify the relativistic gravitational potential that shows exact correspondence to the Newtonian one to the second order.

Nonlinear ordinary differential equations analytical approximation and numerical methods

CERN Document Server

Hermann, Martin

2016-01-01

The book discusses the solutions to nonlinear ordinary differential equations (ODEs) using analytical and numerical approximation methods. Recently, analytical approximation methods have been largely used in solving linear and nonlinear lower-order ODEs. It also discusses using these methods to solve some strong nonlinear ODEs. There are two chapters devoted to solving nonlinear ODEs using numerical methods, as in practice high-dimensional systems of nonlinear ODEs that cannot be solved by analytical approximate methods are common. Moreover, it studies analytical and numerical techniques for the treatment of parameter-depending ODEs. The book explains various methods for solving nonlinear-oscillator and structural-system problems, including the energy balance method, harmonic balance method, amplitude frequency formulation, variational iteration method, homotopy perturbation method, iteration perturbation method, homotopy analysis method, simple and multiple shooting method, and the nonlinear stabilized march...

Estimates of solutions of certain classes of second-order differential equations in a Hilbert space

International Nuclear Information System (INIS)

Artamonov, N V

2003-01-01

Linear second-order differential equations of the form u''(t)+(B+iD)u'(t)+(T+iS)u(t)=0 in a Hilbert space are studied. Under certain conditions on the (generally speaking, unbounded) operators T, S, B and D the correct solubility of the equation in the 'energy' space is proved and best possible (in the general case) estimates of the solutions on the half-axis are obtained

Second-Order Conformally Equivariant Quantization in Dimension 1|2

Directory of Open Access Journals (Sweden)

Najla Mellouli

2009-12-01

Full Text Available This paper is the next step of an ambitious program to develop conformally equivariant quantization on supermanifolds. This problem was considered so far in (superdimensions 1 and 1|1. We will show that the case of several odd variables is much more difficult. We consider the supercircle S^{1|2} equipped with the standard contact structure. The conformal Lie superalgebra K(2 of contact vector fields on S^{1|2} contains the Lie superalgebra osp(2|2. We study the spaces of linear differential operators on the spaces of weighted densities as modules over osp(2|2. We prove that, in the non-resonant case, the spaces of second order differential operators are isomorphic to the corresponding spaces of symbols as osp(2|2-modules. We also prove that the conformal equivariant quantization map is unique and calculate its explicit formula.

Normal-mode-based analysis of electron plasma waves with second-order Hermitian formalism

Science.gov (United States)

Ramos, J. J.; White, R. L.

2018-03-01

The classic problem of the dynamic evolution and Landau damping of linear Langmuir electron waves in a collisionless plasma with Maxwellian background is cast as a second-order, self-adjoint problem with a continuum spectrum of real and positive squared frequencies. The corresponding complete basis of singular normal modes is obtained, along with their orthogonality relation. This yields easily the general expression of the time-reversal-invariant solution for any initial-value problem. Examples are given for specific initial conditions that illustrate different behaviors of the Landau-damped macroscopic moments of the perturbations.

A new multi-domain method based on an analytical control surface for linear and second-order mean drift wave loads on floating bodies

Science.gov (United States)

Liang, Hui; Chen, Xiaobo

2017-10-01

A novel multi-domain method based on an analytical control surface is proposed by combining the use of free-surface Green function and Rankine source function. A cylindrical control surface is introduced to subdivide the fluid domain into external and internal domains. Unlike the traditional domain decomposition strategy or multi-block method, the control surface here is not panelized, on which the velocity potential and normal velocity components are analytically expressed as a series of base functions composed of Laguerre function in vertical coordinate and Fourier series in the circumference. Free-surface Green function is applied in the external domain, and the boundary integral equation is constructed on the control surface in the sense of Galerkin collocation via integrating test functions orthogonal to base functions over the control surface. The external solution gives rise to the so-called Dirichlet-to-Neumann [DN2] and Neumann-to-Dirichlet [ND2] relations on the control surface. Irregular frequencies, which are only dependent on the radius of the control surface, are present in the external solution, and they are removed by extending the boundary integral equation to the interior free surface (circular disc) on which the null normal derivative of potential is imposed, and the dipole distribution is expressed as Fourier-Bessel expansion on the disc. In the internal domain, where the Rankine source function is adopted, new boundary integral equations are formulated. The point collocation is imposed over the body surface and free surface, while the collocation of the Galerkin type is applied on the control surface. The present method is valid in the computation of both linear and second-order mean drift wave loads. Furthermore, the second-order mean drift force based on the middle-field formulation can be calculated analytically by using the coefficients of the Fourier-Laguerre expansion.

Uncertainty evaluation for ordinary least-square fitting with arbitrary order polynomial in joule balance method

International Nuclear Information System (INIS)

You, Qiang; Xu, JinXin; Wang, Gang; Zhang, Zhonghua

2016-01-01

The ordinary least-square fitting with polynomial is used in both the dynamic phase of the watt balance method and the weighting phase of joule balance method but few researches have been conducted to evaluate the uncertainty of the fitting data in the electrical balance methods. In this paper, a matrix-calculation method for evaluating the uncertainty of the polynomial fitting data is derived and the properties of this method are studied by simulation. Based on this, another two derived methods are proposed. One is used to find the optimal fitting order for the watt or joule balance methods. Accuracy and effective factors of this method are experimented with simulations. The other is used to evaluate the uncertainty of the integral of the fitting data for joule balance, which is demonstrated with an experiment from the NIM-1 joule balance. (paper)

Lagrange-Noether method for solving second-order differential equations

Institute of Scientific and Technical Information of China (English)

Wu Hui-Bin; Wu Run-Heng

2009-01-01

The purpose of this paper is to provide a new method called the Lagrange-Noether method for solving second-order differential equations. The method is,firstly,to write the second-order differential equations completely or partially in the form of Lagrange equations,and secondly,to obtain the integrals of the equations by using the Noether theory of the Lagrange system. An example is given to illustrate the application of the result.

Second-order sliding mode controller with model reference adaptation for automatic train operation

Science.gov (United States)

Ganesan, M.; Ezhilarasi, D.; Benni, Jijo

2017-11-01

In this paper, a new approach to model reference based adaptive second-order sliding mode control together with adaptive state feedback is presented to control the longitudinal dynamic motion of a high speed train for automatic train operation with the objective of minimal jerk travel by the passengers. The nonlinear dynamic model for the longitudinal motion of the train comprises of a locomotive and coach subsystems is constructed using multiple point-mass model by considering the forces acting on the vehicle. An adaptation scheme using Lyapunov criterion is derived to tune the controller gains by considering a linear, stable reference model that ensures the stability of the system in closed loop. The effectiveness of the controller tracking performance is tested under uncertain passenger load, coupler-draft gear parameters, propulsion resistance coefficients variations and environmental disturbances due to side wind and wet rail conditions. The results demonstrate improved tracking performance of the proposed control scheme with a least jerk under maximum parameter uncertainties when compared to constant gain second-order sliding mode control.

An algebraic fractional order differentiator for a class of signals satisfying a linear differential equation

KAUST Repository

Liu, Da-Yan; Tian, Yang; Boutat, Driss; Laleg-Kirati, Taous-Meriem

2015-01-01

This paper aims at designing a digital fractional order differentiator for a class of signals satisfying a linear differential equation to estimate fractional derivatives with an arbitrary order in noisy case, where the input can be unknown or known with noises. Firstly, an integer order differentiator for the input is constructed using a truncated Jacobi orthogonal series expansion. Then, a new algebraic formula for the Riemann-Liouville derivative is derived, which is enlightened by the algebraic parametric method. Secondly, a digital fractional order differentiator is proposed using a numerical integration method in discrete noisy case. Then, the noise error contribution is analyzed, where an error bound useful for the selection of the design parameter is provided. Finally, numerical examples illustrate the accuracy and the robustness of the proposed fractional order differentiator.

An algebraic fractional order differentiator for a class of signals satisfying a linear differential equation

KAUST Repository

Liu, Da-Yan

2015-04-30

This paper aims at designing a digital fractional order differentiator for a class of signals satisfying a linear differential equation to estimate fractional derivatives with an arbitrary order in noisy case, where the input can be unknown or known with noises. Firstly, an integer order differentiator for the input is constructed using a truncated Jacobi orthogonal series expansion. Then, a new algebraic formula for the Riemann-Liouville derivative is derived, which is enlightened by the algebraic parametric method. Secondly, a digital fractional order differentiator is proposed using a numerical integration method in discrete noisy case. Then, the noise error contribution is analyzed, where an error bound useful for the selection of the design parameter is provided. Finally, numerical examples illustrate the accuracy and the robustness of the proposed fractional order differentiator.

Second-order analysis of structured inhomogeneous spatio-temporal point processes

DEFF Research Database (Denmark)

Møller, Jesper; Ghorbani, Mohammad

Statistical methodology for spatio-temporal point processes is in its infancy. We consider second-order analysis based on pair correlation functions and K-functions for first general inhomogeneous spatio-temporal point processes and second inhomogeneous spatio-temporal Cox processes. Assuming...... spatio-temporal separability of the intensity function, we clarify different meanings of second-order spatio-temporal separability. One is second-order spatio-temporal independence and relates e.g. to log-Gaussian Cox processes with an additive covariance structure of the underlying spatio......-temporal Gaussian process. Another concerns shot-noise Cox processes with a separable spatio-temporal covariance density. We propose diagnostic procedures for checking hypotheses of second-order spatio-temporal separability, which we apply on simulated and real data (the UK 2001 epidemic foot and mouth disease data)....

Second-order generalized perturbation theory for source-driven systems

International Nuclear Information System (INIS)

Greenspan, E.; Gilai, D.; Oblow, E.M.

1978-01-01

A second-order generalized perturbation theory (GPT) for the effect of multiple system variations on a general flux functional in source-driven systems is derived. The derivation is based on a functional Taylor series in which second-order derivatives are retained. The resulting formulation accounts for the nonlinear effect of a given variation accurate to third order in the flux and adjoint perturbations. It also accounts for the effect of interaction between any number of variations. The new formulation is compared with exact perturbation theory as well as with perturbation theory for altered systems. The usefulnes of the second-order GPT formulation is illustrated by applying it to optimization problems. Its applicability to areas of cross-section sensitivity analysis and system design and evaluation is also discussed

Assessment of First- and Second-Order Wave-Excitation Load Models for Cylindrical Substructures: Preprint

Energy Technology Data Exchange (ETDEWEB)

Pereyra, Brandon; Wendt, Fabian; Robertson, Amy; Jonkman, Jason

2017-03-09

The hydrodynamic loads on an offshore wind turbine's support structure present unique engineering challenges for offshore wind. Two typical approaches used for modeling these hydrodynamic loads are potential flow (PF) and strip theory (ST), the latter via Morison's equation. This study examines the first- and second-order wave-excitation surge forces on a fixed cylinder in regular waves computed by the PF and ST approaches to (1) verify their numerical implementations in HydroDyn and (2) understand when the ST approach breaks down. The numerical implementation of PF and ST in HydroDyn, a hydrodynamic time-domain solver implemented as a module in the FAST wind turbine engineering tool, was verified by showing the consistency in the first- and second-order force output between the two methods across a range of wave frequencies. ST is known to be invalid at high frequencies, and this study investigates where the ST solution diverges from the PF solution. Regular waves across a range of frequencies were run in HydroDyn for a monopile substructure. As expected, the solutions for the first-order (linear) wave-excitation loads resulting from these regular waves are similar for PF and ST when the diameter of the cylinder is small compared to the length of the waves (generally when the diameter-to-wavelength ratio is less than 0.2). The same finding applies to the solutions for second-order wave-excitation loads, but for much smaller diameter-to-wavelength ratios (based on wavelengths of first-order waves).

Second order optical nonlinearity in silicon by symmetry breaking

Energy Technology Data Exchange (ETDEWEB)

Cazzanelli, Massimo, E-mail: massimo.cazzanelli@unitn.it [Laboratorio IdEA, Dipartimento di Fisica, Università di Trento, via Sommarive, 14 Povo (Trento) (Italy); Schilling, Joerg, E-mail: joerg.schilling@physik.uni-halle.de [Centre for Innovation Competence SiLi-nano, Martin-Luther-University Halle-Wittenberg, Karl-Freiherr-von-Fritsch Str. 3, 06120 Halle (Germany)

2016-03-15

Although silicon does not possess a dipolar bulk second order nonlinear susceptibility due to its centro-symmetric crystal structure, in recent years several attempts were undertaken to create such a property in silicon. This review presents the different sources of a second order susceptibility (χ{sup (2)}) in silicon and the connected second order nonlinear effects which were investigated up to now. After an introduction, a theoretical overview discusses the second order nonlinearity in general and distinguishes between the dipolar contribution—which is usually dominating in non-centrosymmetric structures—and the quadrupolar contribution, which even exists in centro-symmetric materials. Afterwards, the classic work on second harmonic generation from silicon surfaces in reflection measurements is reviewed. Due to the abrupt symmetry breaking at surfaces and interfaces locally a dipolar second order susceptibility appears, resulting in, e.g., second harmonic generation. Since the bulk contribution is usually small, the study of this second harmonic signal allows a sensitive observation of the surface/interface conditions. The impact of covering films, strain, electric fields, and defect states at the interfaces was already investigated in this way. With the advent of silicon photonics and the search for ever faster electrooptic modulators, the interest turned to the creation of a dipolar bulk χ{sup (2)} in silicon. These efforts have been focussing on several experiments applying an inhomogeneous strain to the silicon lattice to break its centro-symmetry. Recent results suggesting the impact of electric fields which are exerted from fixed charges in adjacent covering layers are also included. After a subsequent summary on “competing” concepts using not Si but Si-related materials, the paper will end with some final conclusions, suggesting possible future research direction in this dynamically developing field.

Existence and controllability results for damped second order impulsive functional differential systems with state-dependent delay

Directory of Open Access Journals (Sweden)

M. Mallika Arjunan

2014-01-01

Full Text Available In this paper, we investigate the existence and controllability of mild solutions for a damped second order impulsive functional differential equation with state-dependent delay in Banach spaces. The results are obtained by using Sadovskii's fixed point theorem combined with the theories of a strongly continuous cosine family of bounded linear operators. Finally, an example is provided to illustrate the main results.

Probabilistic Sophistication, Second Order Stochastic Dominance, and Uncertainty Aversion

OpenAIRE

Simone Cerreia-Vioglio; Fabio Maccheroni; Massimo Marinacci; Luigi Montrucchio

2010-01-01

We study the interplay of probabilistic sophistication, second order stochastic dominance, and uncertainty aversion, three fundamental notions in choice under uncertainty. In particular, our main result, Theorem 2, characterizes uncertainty averse preferences that satisfy second order stochastic dominance, as well as uncertainty averse preferences that are probabilistically sophisticated.

Second-order polarization-mode dispersion in photonic crystal fibers

DEFF Research Database (Denmark)

Larsen, T; Bjarklev, Anders Overgaard; Peterson, A

2003-01-01

We report the first experimental measurements of second-order polarization-mode dispersion in two successive 900 meter pulls of a silica photonic crystal fiber.......We report the first experimental measurements of second-order polarization-mode dispersion in two successive 900 meter pulls of a silica photonic crystal fiber....

Method to render second order beam optics programs symplectic

International Nuclear Information System (INIS)

Douglas, D.; Servranckx, R.V.

1984-10-01

We present evidence that second order matrix-based beam optics programs violate the symplectic condition. A simple method to avoid this difficulty, based on a generating function approach to evaluating transfer maps, is described. A simple example illustrating the non-symplectricity of second order matrix methods, and the effectiveness of our solution to the problem, is provided. We conclude that it is in fact possible to bring second order matrix optics methods to a canonical form. The procedure for doing so has been implemented in the program DIMAT, and could be implemented in programs such as TRANSPORT and TURTLE, making them useful in multiturn applications. 15 refs

Finite-dimensional linear algebra

CERN Document Server

Gockenbach, Mark S

2010-01-01

Some Problems Posed on Vector SpacesLinear equationsBest approximationDiagonalizationSummaryFields and Vector SpacesFields Vector spaces Subspaces Linear combinations and spanning sets Linear independence Basis and dimension Properties of bases Polynomial interpolation and the Lagrange basis Continuous piecewise polynomial functionsLinear OperatorsLinear operatorsMore properties of linear operatorsIsomorphic vector spaces Linear operator equations Existence and uniqueness of solutions The fundamental theorem; inverse operatorsGaussian elimination Newton's method Linear ordinary differential eq

Groups of integral transforms generated by Lie algebras of second-and higher-order differential operators

International Nuclear Information System (INIS)

Steinberg, S.; Wolf, K.B.

1979-01-01

The authors study the construction and action of certain Lie algebras of second- and higher-order differential operators on spaces of solutions of well-known parabolic, hyperbolic and elliptic linear differential equations. The latter include the N-dimensional quadratic quantum Hamiltonian Schroedinger equations, the one-dimensional heat and wave equations and the two-dimensional Helmholtz equation. In one approach, the usual similarity first-order differential operator algebra of the equation is embedded in the larger one, which appears as a quantum-mechanical dynamic algebra. In a second approach, the new algebra is built as the time evolution of a finite-transformation algebra on the initial conditions. In a third approach, the algebra to inhomogeneous similarity algebra is deformed to a noncompact classical one. In every case, we can integrate the algebra to a Lie group of integral transforms acting effectively on the solution space of the differential equation. (author)

A Comparison of Two-Stage Approaches for Fitting Nonlinear Ordinary Differential Equation Models with Mixed Effects.

Science.gov (United States)

Chow, Sy-Miin; Bendezú, Jason J; Cole, Pamela M; Ram, Nilam

2016-01-01

Several approaches exist for estimating the derivatives of observed data for model exploration purposes, including functional data analysis (FDA; Ramsay & Silverman, 2005 ), generalized local linear approximation (GLLA; Boker, Deboeck, Edler, & Peel, 2010 ), and generalized orthogonal local derivative approximation (GOLD; Deboeck, 2010 ). These derivative estimation procedures can be used in a two-stage process to fit mixed effects ordinary differential equation (ODE) models. While the performance and utility of these routines for estimating linear ODEs have been established, they have not yet been evaluated in the context of nonlinear ODEs with mixed effects. We compared properties of the GLLA and GOLD to an FDA-based two-stage approach denoted herein as functional ordinary differential equation with mixed effects (FODEmixed) in a Monte Carlo (MC) study using a nonlinear coupled oscillators model with mixed effects. Simulation results showed that overall, the FODEmixed outperformed both the GLLA and GOLD across all the embedding dimensions considered, but a novel use of a fourth-order GLLA approach combined with very high embedding dimensions yielded estimation results that almost paralleled those from the FODEmixed. We discuss the strengths and limitations of each approach and demonstrate how output from each stage of FODEmixed may be used to inform empirical modeling of young children's self-regulation.

Microscopic cascading of second-order molecular nonlinearity: New design principles for enhancing third-order nonlinearity.

Science.gov (United States)

Baev, Alexander; Autschbach, Jochen; Boyd, Robert W; Prasad, Paras N

2010-04-12

Herein, we develop a phenomenological model for microscopic cascading and substantiate it with ab initio calculations. It is shown that the concept of local microscopic cascading of a second-order nonlinearity can lead to a third-order nonlinearity, without introducing any new loss mechanisms that could limit the usefulness of our approach. This approach provides a new molecular design protocol, in which the current great successes achieved in producing molecules with extremely large second-order nonlinearity can be used in a supra molecular organization in a preferred orientation to generate very large third-order response magnitudes. The results of density functional calculations for a well-known second-order molecule, (para)nitroaniline, show that a head-to-tail dimer configuration exhibits enhanced third-order nonlinearity, in agreement with the phenomenological model which suggests that such an arrangement will produce cascading due to local field effects.

Characterization of Symmetry Properties of First Integrals for Submaximal Linearizable Third-Order ODEs

Directory of Open Access Journals (Sweden)

K. S. Mahomed

2013-01-01

Full Text Available The relationship between first integrals of submaximal linearizable third-order ordinary differential equations (ODEs and their symmetries is investigated. We obtain the classifying relations between the symmetries and the first integral for submaximal cases of linear third-order ODEs. It is known that the maximum Lie algebra of the first integral is achieved for the simplest equation and is four-dimensional. We show that for the other two classes they are not unique. We also obtain counting theorems of the symmetry properties of the first integrals for these classes of linear third-order ODEs. For the 5 symmetry class of linear third-order ODEs, the first integrals can have 0, 1, 2, and 3 symmetries, and for the 4 symmetry class of linear third-order ODEs, they are 0, 1, and 2 symmetries, respectively. In the case of submaximal linear higher-order ODEs, we show that their full Lie algebras can be generated by the subalgebras of certain basic integrals.

Conformal conservation laws for second-order scalar fields

International Nuclear Information System (INIS)

Blakeskee, J.S.; Logan, J.D.

1976-01-01

It is considered an action integral over space-time whose Lagrangian depends upon a scalar field an upon derivatives of the field function up to second order. From invariance identities obtained by the authors in an earlier work it is shown how a new proof of Noether's theorem for this second-order problem follows in the multiple integral case. Finally, conservation laws are written down in the case that the given action integral be invariant under the fifteen-parameter special conformal group

Higher-order Cauchy of the second kind and poly-Cauchy of the second kind mixed type polynomials

OpenAIRE

Kim, Dae San; Kim, Taekyun

2013-01-01

In this paper, we investigate some properties of higher-order Cauchy of the second kind and poly-Cauchy of the second mixed type polynomials with umbral calculus viewpoint. From our investigation, we derive many interesting identities of higher-order Cauchy of the second kind and poly-Cauchy of the second kind mixed type polynomials.

Optimality Conditions in Differentiable Vector Optimization via Second-Order Tangent Sets

International Nuclear Information System (INIS)

Jimenez, Bienvenido; Novo, Vicente

2004-01-01

We provide second-order necessary and sufficient conditions for a point to be an efficient element of a set with respect to a cone in a normed space, so that there is only a small gap between necessary and sufficient conditions. To this aim, we use the common second-order tangent set and the asymptotic second-order cone utilized by Penot. As an application we establish second-order necessary conditions for a point to be a solution of a vector optimization problem with an arbitrary feasible set and a twice Frechet differentiable objective function between two normed spaces. We also establish second-order sufficient conditions when the initial space is finite-dimensional so that there is no gap with necessary conditions. Lagrange multiplier rules are also given

Computation of Optimal Monotonicity Preserving General Linear Methods

KAUST Repository

Ketcheson, David I.

2009-07-01

Monotonicity preserving numerical methods for ordinary differential equations prevent the growth of propagated errors and preserve convex boundedness properties of the solution. We formulate the problem of finding optimal monotonicity preserving general linear methods for linear autonomous equations, and propose an efficient algorithm for its solution. This algorithm reliably finds optimal methods even among classes involving very high order accuracy and that use many steps and/or stages. The optimality of some recently proposed methods is verified, and many more efficient methods are found. We use similar algorithms to find optimal strong stability preserving linear multistep methods of both explicit and implicit type, including methods for hyperbolic PDEs that use downwind-biased operators.

A Linear-Elasticity Solver for Higher-Order Space-Time Mesh Deformation

Science.gov (United States)

Diosady, Laslo T.; Murman, Scott M.

2018-01-01

A linear-elasticity approach is presented for the generation of meshes appropriate for a higher-order space-time discontinuous finite-element method. The equations of linear-elasticity are discretized using a higher-order, spatially-continuous, finite-element method. Given an initial finite-element mesh, and a specified boundary displacement, we solve for the mesh displacements to obtain a higher-order curvilinear mesh. Alternatively, for moving-domain problems we use the linear-elasticity approach to solve for a temporally discontinuous mesh velocity on each time-slab and recover a continuous mesh deformation by integrating the velocity. The applicability of this methodology is presented for several benchmark test cases.

Neural network construction of flow of a viscoelastic fluid of a second order between two eccentric spheres

International Nuclear Information System (INIS)

Elbakry, M.Y.; El-Helly, M.; Elbakry, M.Y.

2010-01-01

Neural networks are widely for solving many scientific linear and non-linear problems. In this work ,we used the artificial neural network (ANN) to simulate and predict the torque and force acting on the outer stationary sphere due to steady state motion of the second order fluid between two eccentric spheres by a rotating inner sphere with an angular velocity Ω. the (ANN) model has been trained based on the experimental data to produce the torque and force at different eccentricities. The experimental and trained torque and force are compared. The designed ANN shows a good match to the experimental data.

Second-order accurate volume-of-fluid algorithms for tracking material interfaces

International Nuclear Information System (INIS)

Pilliod, James Edward; Puckett, Elbridge Gerry

2004-01-01

We introduce two new volume-of-fluid interface reconstruction algorithms and compare the accuracy of these algorithms to four other widely used volume-of-fluid interface reconstruction algorithms. We find that when the interface is smooth (e.g., continuous with two continuous derivatives) the new methods are second-order accurate and the other algorithms are first-order accurate. We propose a design criteria for a volume-of-fluid interface reconstruction algorithm to be second-order accurate. Namely, that it reproduce lines in two space dimensions or planes in three space dimensions exactly. We also introduce a second-order, unsplit, volume-of-fluid advection algorithm that is based on a second-order, finite difference method for scalar conservation laws due to Bell, Dawson and Shubin. We test this advection algorithm by modeling several different interface shapes propagating in two simple incompressible flows and compare the results with the standard second-order, operator-split advection algorithm. Although both methods are second-order accurate when the interface is smooth, we find that the unsplit algorithm exhibits noticeably better resolution in regions where the interface has discontinuous derivatives, such as at corners

Passive longitudinal phase space linearizer

Directory of Open Access Journals (Sweden)

P. Craievich

2010-03-01

Full Text Available We report on the possibility to passively linearize the bunch compression process in electron linacs for the next generation x-ray free electron lasers. This can be done by using the monopole wakefields in a dielectric-lined waveguide. The optimum longitudinal voltage loss over the length of the bunch is calculated in order to compensate both the second-order rf time curvature and the second-order momentum compaction terms. Thus, the longitudinal phase space after the compression process is linearized up to a fourth-order term introduced by the convolution between the bunch and the monopole wake function.

Linear conversion theory on the second harmonic emission from a plasma filament

International Nuclear Information System (INIS)

Tan Weihan; Gu Min

1989-01-01

The linear conversion theory of laser produced plasma filaments is studied. By calculations for the energy flux of the second harmonic emission on the basis of the planar wave-plasma interaction model, it has been found that there exists no 2ω 0 harmonic emission in the direction perpendicular to the incident laser, in contradiction with the experiments. A linear conversion theory is proposed on the second harmonic emission from a plasma filament and discovered the intense 2ω 0 harmonic emission in the direction perpendicular to the incident laser, which is in agreement with the experiments. (author)

Impurity strength and impurity domain modulated frequency-dependent linear and second non-linear response properties of doped quantum dots

Energy Technology Data Exchange (ETDEWEB)

Datta, Nirmal Kumar [Department of Physics, Suri Vidyasagar College, Suri, Birbhum 731 101, West Bengal (India); Ghosh, Manas [Department of Chemistry, Physical Chemistry Section, Visva Bharati University, Santiniketan, Birbhum 731 235, West Bengal (India)

2011-08-15

We explore the pattern of frequency-dependent linear and second non-linear optical responses of repulsive impurity doped quantum dots harmonically confined in two dimensions. The dopant impurity potential chosen assumes a Gaussian form and it is doped into an on-center location. The quantum dot is subject to a periodically oscillating external electric field. For some fixed values of transverse magnetic field strength ({omega}{sub c}) and harmonic confinement potential ({omega}{sub 0}), the influence of impurity strength (V{sub 0}) and impurity domain ({xi}) on the diagonal components of the frequency-dependent linear ({alpha}{sub xx} and {alpha}{sub yy}) and second non-linear ({gamma}{sub xxxx} and {gamma}{sub yyyy}) responses of the dot are computed through a linear variational route. The investigations reveal that the optical responses undergo enhancement with increase in both V{sub 0} and {xi} values. However, in the limitingly small dopant strength regime one observes a drop in the optical responses with increase in V{sub 0}. A time-average rate of energy transfer to the system is often invoked to support the findings. (Copyright copyright 2011 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim)

Application of Hierarchical Linear Models/Linear Mixed-Effects Models in School Effectiveness Research

Science.gov (United States)

Ker, H. W.

2014-01-01

Multilevel data are very common in educational research. Hierarchical linear models/linear mixed-effects models (HLMs/LMEs) are often utilized to analyze multilevel data nowadays. This paper discusses the problems of utilizing ordinary regressions for modeling multilevel educational data, compare the data analytic results from three regression…

Self-triggered rendezvous of gossiping second-order agents

NARCIS (Netherlands)

De Persis, Claudio; Frasca, Paolo; Hendrickx, Julien M.

2013-01-01

A recent paper by some of the authors introduced several self-triggered coordination algorithms for first-order continuous-time systems. The extension of these algorithms to second-order agents is relevant in many practical applications but presents some challenges that are tackled in this

A Three Step Explicit Method for Direct Solution of Third Order ...

African Journals Online (AJOL)

This study produces a three step discrete Linear Multistep Method for Direct solution of third order initial value problems of ordinary differential equations of the form y'''= f(x,y,y',y''). Taylor series expansion technique was adopted in the development of the method. The differential system from the basis polynomial function to ...

FORSIM, Solution of Ordinary or Partial Differential Equation with Initial Conditions

International Nuclear Information System (INIS)

Carver, M.B.

1985-01-01

1 - Description of problem or function: FORSIM is a FORTRAN oriented simulation program which automates the continuous transient solution of systems of ordinary and/or partial differential equations. The user writes his equations in a FORTRAN subroutine, following prescribed rules, and loads this routine along with the executive routines. The executive routines then read in initial data supplied by the user and proceed with the integration. 2 - Method of solution: Partial differential equations are converted to coupled ordinary differential equations by suitable discretization formulae. Integration is done by variable order, variable step-size error controlled algorithms. 3 - Restrictions on the complexity of the problem - Maximum of: 1000 ordinary differential equations

Punishing second-order free riders before first-order free riders: The effect of pool punishment priority on cooperation

OpenAIRE

Ozono, Hiroki; Kamijo, Yoshio; Shimizu, Kazumi

2017-01-01

Second-order free riders, who do not owe punishment cost to first-order free riders in public goods games, lead to low cooperation. Previous studies suggest that for stable cooperation, it is critical to have a pool punishment system with second-order punishment, which gathers resources from group members and punishes second-order free riders as well as first-order free riders. In this study, we focus on the priority of punishment. We hypothesize that the pool punishment system that prioritiz...

Algebraic dynamics solutions and algebraic dynamics algorithm for nonlinear ordinary differential equations

Institute of Scientific and Technical Information of China (English)

WANG; Shunjin; ZHANG; Hua

2006-01-01

The problem of preserving fidelity in numerical computation of nonlinear ordinary differential equations is studied in terms of preserving local differential structure and approximating global integration structure of the dynamical system.The ordinary differential equations are lifted to the corresponding partial differential equations in the framework of algebraic dynamics,and a new algorithm-algebraic dynamics algorithm is proposed based on the exact analytical solutions of the ordinary differential equations by the algebraic dynamics method.In the new algorithm,the time evolution of the ordinary differential system is described locally by the time translation operator and globally by the time evolution operator.The exact analytical piece-like solution of the ordinary differential equations is expressd in terms of Taylor series with a local convergent radius,and its finite order truncation leads to the new numerical algorithm with a controllable precision better than Runge Kutta Algorithm and Symplectic Geometric Algorithm.

Differential effects of exogenous and endogenous attention on second-order texture contrast sensitivity

Science.gov (United States)

Barbot, Antoine; Landy, Michael S.; Carrasco, Marisa

2012-01-01

The visual system can use a rich variety of contours to segment visual scenes into distinct perceptually coherent regions. However, successfully segmenting an image is a computationally expensive process. Previously we have shown that exogenous attention—the more automatic, stimulus-driven component of spatial attention—helps extract contours by enhancing contrast sensitivity for second-order, texture-defined patterns at the attended location, while reducing sensitivity at unattended locations, relative to a neutral condition. Interestingly, the effects of exogenous attention depended on the second-order spatial frequency of the stimulus. At parafoveal locations, attention enhanced second-order contrast sensitivity to relatively high, but not to low second-order spatial frequencies. In the present study we investigated whether endogenous attention—the more voluntary, conceptually-driven component of spatial attention—affects second-order contrast sensitivity, and if so, whether its effects are similar to those of exogenous attention. To that end, we compared the effects of exogenous and endogenous attention on the sensitivity to second-order, orientation-defined, texture patterns of either high or low second-order spatial frequencies. The results show that, like exogenous attention, endogenous attention enhances second-order contrast sensitivity at the attended location and reduces it at unattended locations. However, whereas the effects of exogenous attention are a function of the second-order spatial frequency content, endogenous attention affected second-order contrast sensitivity independent of the second-order spatial frequency content. This finding supports the notion that both exogenous and endogenous attention can affect second-order contrast sensitivity, but that endogenous attention is more flexible, benefitting performance under different conditions. PMID:22895879

Differential effects of exogenous and endogenous attention on second-order texture contrast sensitivity.

Science.gov (United States)

Barbot, Antoine; Landy, Michael S; Carrasco, Marisa

2012-08-15

The visual system can use a rich variety of contours to segment visual scenes into distinct perceptually coherent regions. However, successfully segmenting an image is a computationally expensive process. Previously we have shown that exogenous attention--the more automatic, stimulus-driven component of spatial attention--helps extract contours by enhancing contrast sensitivity for second-order, texture-defined patterns at the attended location, while reducing sensitivity at unattended locations, relative to a neutral condition. Interestingly, the effects of exogenous attention depended on the second-order spatial frequency of the stimulus. At parafoveal locations, attention enhanced second-order contrast sensitivity to relatively high, but not to low second-order spatial frequencies. In the present study we investigated whether endogenous attention-the more voluntary, conceptually-driven component of spatial attention--affects second-order contrast sensitivity, and if so, whether its effects are similar to those of exogenous attention. To that end, we compared the effects of exogenous and endogenous attention on the sensitivity to second-order, orientation-defined, texture patterns of either high or low second-order spatial frequencies. The results show that, like exogenous attention, endogenous attention enhances second-order contrast sensitivity at the attended location and reduces it at unattended locations. However, whereas the effects of exogenous attention are a function of the second-order spatial frequency content, endogenous attention affected second-order contrast sensitivity independent of the second-order spatial frequency content. This finding supports the notion that both exogenous and endogenous attention can affect second-order contrast sensitivity, but that endogenous attention is more flexible, benefitting performance under different conditions.

Temporal Frequency Modulates Reaction Time Responses to First-Order and Second-Order Motion

Science.gov (United States)

Hutchinson, Claire V.; Ledgeway, Tim

2010-01-01

This study investigated the effect of temporal frequency and modulation depth on reaction times for discriminating the direction of first-order (luminance-defined) and second-order (contrast-defined) motion, equated for visibility using equal multiples of direction-discrimination threshold. Results showed that reaction times were heavily…

Exact calculation of three-body contact interaction to second order

International Nuclear Information System (INIS)

Kaiser, N.

2012-01-01

For a system of fermions with a three-body contact interaction the second-order contributions to the energy per particle anti E(k f ) are calculated exactly. The three-particle scattering amplitude in the medium is derived in closed analytical form from the corresponding two-loop rescattering diagram. We compare the (genuine) second-order three-body contribution to anti E(k f )∝k f 10 with the second-order term due to the density-dependent effective two-body interaction, and find that the latter term dominates. The results of the present study are of interest for nuclear many-body calculations where chiral three-nucleon forces are treated beyond leading order via a density-dependent effective two-body interaction. (orig.)

Generalized Cross-Gramian for Linear Systems

DEFF Research Database (Denmark)

Shaker, Hamid Reza

2012-01-01

The cross-gramian is a well-known matrix with embedded controllability and observability information. The cross-gramian is related to the Hankel operator and the Hankel singular values of a linear square system and it has several interesting properties. These properties make the cross...... square symmetric systems, the ordinary cross-gramian does not exist. To cope with this problem, a new generalized cross-gramian is introduced in this paper. In contrast to the ordinary cross-gramian, the generalized cross-gramian can be easily obtained for general linear systems and therefore can be used...

Mixed hyperbolic-second-order-parabolic formulations of general relativity

International Nuclear Information System (INIS)

Paschalidis, Vasileios

2008-01-01

Two new formulations of general relativity are introduced. The first one is a parabolization of the Arnowitt-Deser-Misner formulation and is derived by the addition of combinations of the constraints and their derivatives to the right-hand side of the Arnowitt-Deser-Misner evolution equations. The desirable property of this modification is that it turns the surface of constraints into a local attractor because the constraint propagation equations become second-order parabolic independently of the gauge conditions employed. This system may be classified as mixed hyperbolic--second-order parabolic. The second formulation is a parabolization of the Kidder-Scheel-Teukolsky formulation and is a manifestly mixed strongly hyperbolic--second-order-parabolic set of equations, bearing thus resemblance to the compressible Navier-Stokes equations. As a first test, a stability analysis of flat space is carried out and it is shown that the first modification exponentially damps and smoothes all constraint-violating modes. These systems provide a new basis for constructing schemes for long-term and stable numerical integration of the Einstein field equations.

Time-dependent quantum transport through an interacting quantum dot beyond sequential tunneling: second-order quantum rate equations

International Nuclear Information System (INIS)

Dong, B; Ding, G H; Lei, X L

2015-01-01

A general theoretical formulation for the effect of a strong on-site Coulomb interaction on the time-dependent electron transport through a quantum dot under the influence of arbitrary time-varying bias voltages and/or external fields is presented, based on slave bosons and the Keldysh nonequilibrium Green's function (GF) techniques. To avoid the difficulties of computing double-time GFs, we generalize the propagation scheme recently developed by Croy and Saalmann to combine the auxiliary-mode expansion with the celebrated Lacroix's decoupling approximation in dealing with the second-order correlated GFs and then establish a closed set of coupled equations of motion, called second-order quantum rate equations (SOQREs), for an exact description of transient dynamics of electron correlated tunneling. We verify that the stationary solution of our SOQREs is able to correctly describe the Kondo effect on a qualitative level. Moreover, a comparison with other methods, such as the second-order von Neumann approach and Hubbard-I approximation, is performed. As illustrations, we investigate the transient current behaviors in response to a step voltage pulse and a harmonic driving voltage, and linear admittance as well, in the cotunneling regime. (paper)

Aspects of second-order analysis of structured inhomogeneous spatio-temporal processes

DEFF Research Database (Denmark)

Møller, Jesper; Ghorbani, Mohammad

2012-01-01

Statistical methodology for spatio-temporal point processes is in its infancy. We consider second-order analysis based on pair correlation functions and K-functions for general inhomogeneous spatio-temporal point processes and for inhomogeneous spatio-temporal Cox processes. Assuming spatio......-temporal separability of the intensity function, we clarify different meanings of second-order spatio-temporal separability. One is second-order spatio-temporal independence and relates to log-Gaussian Cox processes with an additive covariance structure of the underlying spatio-temporal Gaussian process. Another...... concerns shot-noise Cox processes with a separable spatio-temporal covariance density. We propose diagnostic procedures for checking hypotheses of second-order spatio-temporal separability, which we apply on simulated and real data....

A porous flow model of flank eruptions on Mt. Etna: second-order perturbation theory

Directory of Open Access Journals (Sweden)

N. Cenni

1997-06-01

Full Text Available A porous flow model for magma migration from a deep source within a volcanic edifice is developed. The model is based on the assumption that an isotropic and homogeneous system of fractures allows magma migration from one localized feeding dyke up to the surface of the volcano. The maximum level that magma can reach within the volcano (i.e., the «free surface» of magma, where fluid pressure equals the atmospheric pressure is reproduced through a second-order perturbation approach to the non-linear equations governing the migration of incompressible fluids through a porous medium. The perturbation parameter is found to depend on the ratio of the volumic discharge rate at the source (m3/s divided by the product of the hydraulic conductivity of the medium (m1/s times the square of the source depth. The second-order corrections for the free surface of Mt. Etna are found to be small but not negligible; from the comparison between first-order and second-order free surfaces it appears that the former is higher near the summit, slightly lower at intermediate altitudes and slightly higher far away from the axis of the volcano. Flank eruptions in the southern sector are found to be located in regions where the topography is actually lower than the theoretical free surface of magma. In this sector, modulations in the eruption site density correlate well with even minor differences between free surface and topography. In the northern and western sectors similar good fits are found, while the NE rift and the eastern sector seem to require mechanisms or structures respectively favouring and inhibiting magma migration.

Variability and Variation in Second Language Acquisition Orders : A Dynamic Reevaluation

NARCIS (Netherlands)

Lowie, Wander; Verspoor, Marjolijn

2015-01-01

The traditional morpheme order studies in second language acquisition have tried to demonstrate the existence of a fixed order of acquisition of English morphemes, regardless of the second language learner's background. Such orders have been taken as evidence of the preprogrammed nature of language

Pyrolytic Graphite as a Tunable Second order Neutron Filter

International Nuclear Information System (INIS)

Adib, M.

2009-01-01

A study has been carried out on the neutron transmission through pyrolytic graphite (PG) crystals in order to check its applicability as an efficient tunable second order neutron filter. The neutron transmission have been calculated as a function of neutron wavelengths in the range from 0.01 nm up to 0.7 nm at various PG mosaic spread, thickness and orientation of its c-axis with respect to the beam direction The Computer package Graphite has been used to provide the required calculation. It was shown that highly aligned (10 FWHM on mosaic spread) PG crystal ∼2 cm thick, may be tuned for optimum scattering of 2 second order neutrons within some favorable wavelength intervals in the range between 0.112 and 0.425 nm by adjusting the crystal in an appropriate orientation. .However, a less quality and thinner PG was found to almost eliminate 2 second order neutrons at only tuned values of wavelength corresponding to the poison of the triple intersection points of the curves (hkl) ± and (00l)

Second order guiding-center Vlasov–Maxwell equations

DEFF Research Database (Denmark)

Madsen, Jens

2010-01-01

Second order gyrogauge invariant guiding-center coordinates with strong E×B-flow are derived using the Lie transformation method. The corresponding Poisson bracket structure and equations of motion are obtained. From a variational principle the explicit Vlasov–Maxwell equations are derived...

Second order bounce back boundary condition for the lattice Boltzmann fluid simulation

International Nuclear Information System (INIS)

Kim, In Chan

2000-01-01

A new bounce back boundary method of the second order in error is proposed for the lattice Boltzmann fluid simulation. This new method can be used for the arbitrarily irregular lattice geometry of a non-slip boundary. The traditional bounce back boundary condition for the lattice Boltzmann simulation is of the first order in error. Since the lattice Boltzmann method is the second order scheme by itself, a boundary technique of the second order has been desired to replace the first order bounce back method. This study shows that, contrary to the common belief that the bounce back boundary condition is unilaterally of the first order, the second order bounce back boundary condition can be realized. This study also shows that there exists a generalized bounce back technique that can be characterized by a single interpolation parameter. The second order bounce back method can be obtained by proper selection of this parameter in accordance with the detailed lattice geometry of the boundary. For an illustrative purpose, the transient Couette and the plane Poiseuille flows are solved by the lattice Boltzmann simulation with various boundary conditions. The results show that the generalized bounce back method yields the second order behavior in the error of the solution, provided that the interpolation parameter is properly selected. Coupled with its intuitive nature and the ease of implementation, the bounce back method can be as good as any second order boundary method

The contribution of second-order processes to (3He, n) calculations

International Nuclear Information System (INIS)

Brissaud, I.

1978-01-01

The reactions 90 Zr, 116 Cd( 3 He, n) have been analysed by adding two second-order processes to the usual one-step DWBA calculations: sequential stripping or inelastic transition followed by double stripping. These second-order contributions increase the cross sections, especially for 90 Zr, and improve the shape of the angular distribution for 116 Cd. It shows that such second-order processes cannot be omitted in the analysis of two-particle stripping reactions. (author)

Construction of a Smooth Lyapunov Function for the Robust and Exact Second-Order Differentiator

Directory of Open Access Journals (Sweden)

Tonametl Sanchez

2016-01-01

Full Text Available Differentiators play an important role in (continuous feedback control systems. In particular, the robust and exact second-order differentiator has shown some very interesting properties and it has been used successfully in sliding mode control, in spite of the lack of a Lyapunov based procedure to design its gains. As contribution of this paper, we provide a constructive method to determine a differentiable Lyapunov function for such a differentiator. Moreover, the Lyapunov function is used to provide a procedure to design the differentiator’s parameters. Also, some sets of such parameters are provided. The determination of the positive definiteness of the Lyapunov function and negative definiteness of its derivative is converted to the problem of solving a system of inequalities linear in the parameters of the Lyapunov function candidate and also linear in the gains of the differentiator, but bilinear in both.

Sound dispersion in a spin-1 Ising system near the second-order phase transition point

International Nuclear Information System (INIS)

Erdem, Ryza; Keskin, Mustafa

2003-01-01

Sound dispersion relation is derived for a spin-1 Ising system and its behaviour near the second-order phase transition point or the critical point is analyzed. The method used is a combination of molecular field approximation and Onsager theory of irreversible thermodynamics. If we assume a linear coupling of sound wave with the order parameter fluctuations in the system, we find that the dispersion which is the relative sound velocity change with frequency behaves as ω 0 ε 0 , where ω is the sound frequency and ε the temperature distance from the critical point. In the ordered region, one also observes a frequency-dependent velocity or dispersion minimum which is shifted from the corresponding attenuation maxima. These phenomena are in good agreement with the calculations of sound velocity in other magnetic systems such as magnetic metals, magnetic insulators, and magnetic semiconductors

Second-Rate Coverage of Second-Order Elections: Czech and Slovak Elections to the EP in the Media

Directory of Open Access Journals (Sweden)

Jan Kovář

2010-12-01

Full Text Available Elections to the European Parliament (EP are considered second-order national elections (SOE. The SOE model suggests that there is a qualitative difference between different types of elections depending on the perception of what is at stake. Compared to first order elections, in second order elections there is less at stake because they do not determine the composition of government. Given that voters behave differently in second-order elections, the question arises: do the media also consider second-order elections less interesting and therefore devote to them less coverage? The media play a crucial role in informing citizens about such events as elections; they function as intermediaries between the electorate and the political arena. However, little is known about how EU issues are covered in the media, particularly in the new EU member states. Conducting a content analysis and applying the second-order election model, this paper analyses TV news coverage of the 2004 and 2009 European elections in the Czech Republic and Slovakia in a comparative fashion. The findings are discussed in the light of existing research literature on the EU’s legitimacy as well as its alleged democratic and communication deficit, not least because the EU relies on the media in strengthening (albeit indirectly its legitimacy by increasing citizen awareness of its activities.

Variability and Variation in Second Language Acquisition Orders: A Dynamic Reevaluation

Science.gov (United States)

Lowie, Wander; Verspoor, Marjolijn

2015-01-01

The traditional morpheme order studies in second language acquisition have tried to demonstrate the existence of a fixed order of acquisition of English morphemes, regardless of the second language learner's background. Such orders have been taken as evidence of the preprogrammed nature of language acquisition. This article argues for a…

A cell-centred finite volume method for the Poisson problem on non-graded quadtrees with second order accurate gradients

Science.gov (United States)

Batty, Christopher

2017-02-01

This paper introduces a two-dimensional cell-centred finite volume discretization of the Poisson problem on adaptive Cartesian quadtree grids which exhibits second order accuracy in both the solution and its gradients, and requires no grading condition between adjacent cells. At T-junction configurations, which occur wherever resolution differs between neighboring cells, use of the standard centred difference gradient stencil requires that ghost values be constructed by interpolation. To properly recover second order accuracy in the resulting numerical gradients, prior work addressing block-structured grids and graded trees has shown that quadratic, rather than linear, interpolation is required; the gradients otherwise exhibit only first order convergence, which limits potential applications such as fluid flow. However, previous schemes fail or lose accuracy in the presence of the more complex T-junction geometries arising in the case of general non-graded quadtrees, which place no restrictions on the resolution of neighboring cells. We therefore propose novel quadratic interpolant constructions for this case that enable second order convergence by relying on stencils oriented diagonally and applied recursively as needed. The method handles complex tree topologies and large resolution jumps between neighboring cells, even along the domain boundary, and both Dirichlet and Neumann boundary conditions are supported. Numerical experiments confirm the overall second order accuracy of the method in the L∞ norm.

An exactly solvable model for first- and second-order transitions

International Nuclear Information System (INIS)

Klushin, L I; Skvortsov, A M; Gorbunov, A A

1998-01-01

The possibility of an exact analytical description of first-order and second-order transitions is demonstrated using a specific microscopic model. Predictions using the exactly calculated partition function are compared with those based on the Landau and Yang-Lee approaches. The model employed is an adsorbed polymer chain with an arbitrary number of links and an external force applied to its end, for which the variation of the partition function with the adsorption interaction parameter and the magnitude of the applied force is calculated. In the thermodynamic limit, the system has one isotropic and two anisotropic, ordered phases, each of which is characterized by two order parameters and between which first-order and second-order transitions occur and a bicritical point exists. The Landau free energy is found exactly as a function of each order parameter separately and, near the bicritical point, as a function of both of them simultaneously. An exact analytical formula is found for the distribution of the complex zeros of the partition function in first-order and second-order phase transitions. Hypotheses concerning the way in which the free energy and the positions of the complex zeros scale with the number of particles N in the system are verified. (reviews of topical problems)

Three-dimensional linear peeling-ballooning theory in magnetic fusion devices

Energy Technology Data Exchange (ETDEWEB)

Weyens, T., E-mail: tweyens@fis.uc3m.es; Sánchez, R.; García, L. [Departamento de Física, Universidad Carlos III de Madrid, Madrid 28911 (Spain); Loarte, A.; Huijsmans, G. [ITER Organization, Route de Vinon sur Verdon, 13067 Saint Paul Lez Durance (France)

2014-04-15

Ideal magnetohydrodynamics theory is extended to fully 3D magnetic configurations to investigate the linear stability of intermediate to high n peeling-ballooning modes, with n the toroidal mode number. These are thought to be important for the behavior of edge localized modes and for the limit of the size of the pedestal that governs the high confinement H-mode. The end point of the derivation is a set of coupled second order ordinary differential equations with appropriate boundary conditions that minimize the perturbed energy and that can be solved to find the growth rate of the perturbations. This theory allows of the evaluation of 3D effects on edge plasma stability in tokamaks such as those associated with the toroidal ripple due to the finite number of toroidal field coils, the application of external 3D fields for elm control, local modification of the magnetic field in the vicinity of ferromagnetic components such as the test blanket modules in ITER, etc.

A second order penalized direct forcing for hybrid Cartesian/immersed boundary flow simulations

International Nuclear Information System (INIS)

Introini, C.; Belliard, M.; Fournier, C.

2014-01-01

In this paper, we propose a second order penalized direct forcing method to deal with fluid-structure interaction problems involving complex static or time-varying geometries. As this work constitutes a first step toward more complicated problems, our developments are restricted to Dirichlet boundary condition in purely hydraulic context. The proposed method belongs to the class of immersed boundary techniques and consists in immersing the physical domain in a Cartesian fictitious one of simpler geometry on fixed grids. A penalized forcing term is added to the momentum equation to take the boundary conditions around/inside the obstacles into account. This approach avoids the tedious task of re-meshing and allows us to use fast and accurate numerical schemes. In contrary, as the immersed boundary is described by a set of Lagrangian points that does not generally coincide with those of the Eulerian grid, numerical procedures are required to reconstruct the velocity field near the immersed boundary. Here, we develop a second order linear interpolation scheme and we compare it to a simpler model of order one. As far as the governing equations are concerned, we use a particular fractional-step method in which the penalized forcing term is distributed both in prediction and correction equations. The accuracy of the proposed method is assessed through 2-D numerical experiments involving static and rotating solids. We show in particular that the numerical rate of convergence of our method is quasi-quadratic. (authors)

Non-critically phase-matched second harmonic generation and third order nonlinearity in organic crystal glucuronic acid γ-lactone

Science.gov (United States)

Saripalli, Ravi Kiran; Katturi, Naga Krishnakanth; Soma, Venugopal Rao; Bhat, H. L.; Elizabeth, Suja

2017-12-01

The linear, second order, and third order nonlinear optical properties of glucuronic acid γ-lactone single crystals were investigated. The optic axes and principal dielectric axes were identified through optical conoscopy and the principal refractive indices were obtained using the Brewster's angle method. Conic sections were observed which is perceived to be due to spontaneous non-collinear phase matching. The direction of collinear phase matching was determined and the deff evaluated in this direction was 0.71 pm/V. Open and closed aperture Z-scan measurements with femtosecond pulses revealed high third order nonlinearity in the form of self-defocusing, two-photon absorption, as well as saturable absorption.

On nonlinear differential equation with exact solutions having various pole orders

International Nuclear Information System (INIS)

Kudryashov, N.A.

2015-01-01

We consider a nonlinear ordinary differential equation having solutions with various movable pole order on the complex plane. We show that the pole order of exact solution is determined by values of parameters of the equation. Exact solutions in the form of the solitary waves for the second order nonlinear differential equation are found taking into account the method of the logistic function. Exact solutions of differential equations are discussed and analyzed

Generalized second-order Coulomb phase shift functions

International Nuclear Information System (INIS)

Rosendorff, S.

1982-01-01

Some specific properties and the evaluation of the generalized second-order Coulomb phase shift functions (two-dimensional integrals of four spherical cylinder functions) are discussed. The dependence on the three momenta k 1 ,k-bar,k 2 , corresponding to the final, intermediate, and initial states is illustrated

Conservation properties of numerical integration methods for systems of ordinary differential equations

Science.gov (United States)

Rosenbaum, J. S.

1976-01-01

If a system of ordinary differential equations represents a property conserving system that can be expressed linearly (e.g., conservation of mass), it is then desirable that the numerical integration method used conserve the same quantity. It is shown that both linear multistep methods and Runge-Kutta methods are 'conservative' and that Newton-type methods used to solve the implicit equations preserve the inherent conservation of the numerical method. It is further shown that a method used by several authors is not conservative.

Synchronization from Second Order Network Connectivity Statistics

Science.gov (United States)

Zhao, Liqiong; Beverlin, Bryce; Netoff, Theoden; Nykamp, Duane Q.

2011-01-01

We investigate how network structure can influence the tendency for a neuronal network to synchronize, or its synchronizability, independent of the dynamical model for each neuron. The synchrony analysis takes advantage of the framework of second order networks, which defines four second order connectivity statistics based on the relative frequency of two-connection network motifs. The analysis identifies two of these statistics, convergent connections, and chain connections, as highly influencing the synchrony. Simulations verify that synchrony decreases with the frequency of convergent connections and increases with the frequency of chain connections. These trends persist with simulations of multiple models for the neuron dynamics and for different types of networks. Surprisingly, divergent connections, which determine the fraction of shared inputs, do not strongly influence the synchrony. The critical role of chains, rather than divergent connections, in influencing synchrony can be explained by their increasing the effective coupling strength. The decrease of synchrony with convergent connections is primarily due to the resulting heterogeneity in firing rates. PMID:21779239

Reduced-cost second-order algebraic-diagrammatic construction method for excitation energies and transition moments

Science.gov (United States)

Mester, Dávid; Nagy, Péter R.; Kállay, Mihály

2018-03-01

A reduced-cost implementation of the second-order algebraic-diagrammatic construction [ADC(2)] method is presented. We introduce approximations by restricting virtual natural orbitals and natural auxiliary functions, which results, on average, in more than an order of magnitude speedup compared to conventional, density-fitting ADC(2) algorithms. The present scheme is the successor of our previous approach [D. Mester, P. R. Nagy, and M. Kállay, J. Chem. Phys. 146, 194102 (2017)], which has been successfully applied to obtain singlet excitation energies with the linear-response second-order coupled-cluster singles and doubles model. Here we report further methodological improvements and the extension of the method to compute singlet and triplet ADC(2) excitation energies and transition moments. The various approximations are carefully benchmarked, and conservative truncation thresholds are selected which guarantee errors much smaller than the intrinsic error of the ADC(2) method. Using the canonical values as reference, we find that the mean absolute error for both singlet and triplet ADC(2) excitation energies is 0.02 eV, while that for oscillator strengths is 0.001 a.u. The rigorous cutoff parameters together with the significantly reduced operation count and storage requirements allow us to obtain accurate ADC(2) excitation energies and transition properties using triple-ζ basis sets for systems of up to one hundred atoms.

Investigation of the Impact of Different Terms in the Second Order Hamiltonian on Excitation Energies of Valence and Rydberg States.

Science.gov (United States)

Tajti, Attila; Szalay, Péter G

2016-11-08

Describing electronically excited states of molecules accurately poses a challenging problem for theoretical methods. Popular second order techniques like Linear Response CC2 (CC2-LR), Partitioned Equation-of-Motion MBPT(2) (P-EOM-MBPT(2)), or Equation-of-Motion CCSD(2) (EOM-CCSD(2)) often produce results that are controversial and are ill-balanced with their accuracy on valence and Rydberg type states. In this study, we connect the theory of these methods and, to investigate the origin of their different behavior, establish a series of intermediate variants. The accuracy of these on excitation energies of singlet valence and Rydberg electronic states is benchmarked on a large sample against high-accuracy Linear Response CC3 references. The results reveal the role of individual terms of the second order similarity transformed Hamiltonian, and the reason for the bad performance of CC2-LR in the description of Rydberg states. We also clarify the importance of the T̂ 1 transformation employed in the CC2 procedure, which is found to be very small for vertical excitation energies.

Gravitational waves from global second order phase transitions

Energy Technology Data Exchange (ETDEWEB)

Jr, John T. Giblin [Department of Physics, Kenyon College, 201 North College Rd, Gambier, OH 43022 (United States); Price, Larry R.; Siemens, Xavier; Vlcek, Brian, E-mail: giblinj@kenyon.edu, E-mail: larryp@caltech.edu, E-mail: siemens@gravity.phys.uwm.edu, E-mail: bvlcek@uwm.edu [Center for Gravitation and Cosmology, Department of Physics, University of Wisconsin-Milwaukee, P.O. Box 413, Milwaukee, WI 53201 (United States)

2012-11-01

Global second-order phase transitions are expected to produce scale-invariant gravitational wave spectra. In this manuscript we explore the dynamics of a symmetry-breaking phase transition using lattice simulations. We explicitly calculate the stochastic gravitational wave background produced during the transition and subsequent self-ordering phase. We comment on this signal as it compares to the scale-invariant spectrum produced during inflation.

Linear and quasi-linear equations of parabolic type

CERN Document Server

Ladyženskaja, O A; Ural′ceva, N N; Uralceva, N N

1968-01-01

Equations of parabolic type are encountered in many areas of mathematics and mathematical physics, and those encountered most frequently are linear and quasi-linear parabolic equations of the second order. In this volume, boundary value problems for such equations are studied from two points of view: solvability, unique or otherwise, and the effect of smoothness properties of the functions entering the initial and boundary conditions on the smoothness of the solutions.

Contrast gain control in first- and second-order motion perception.

Science.gov (United States)

Lu, Z L; Sperling, G

1996-12-01

A novel pedestal-plus-test paradigm is used to determine the nonlinear gain-control properties of the first-order (luminance) and the second-order (texture-contrast) motion systems, that is, how these systems' responses to motion stimuli are reduced by pedestals and other masking stimuli. Motion-direction thresholds were measured for test stimuli consisting of drifting luminance and texture-contrast-modulation stimuli superimposed on pedestals of various amplitudes. (A pedestal is a static sine-wave grating of the same type and same spatial frequency as the moving test grating.) It was found that first-order motion-direction thresholds are unaffected by small pedestals, but at pedestal contrasts above 1-2% (5-10 x pedestal threshold), motion thresholds increase proportionally to pedestal amplitude (a Weber law). For first-order stimuli, pedestal masking is specific to the spatial frequency of the test. On the other hand, motion-direction thresholds for texture-contrast stimuli are independent of pedestal amplitude (no gain control whatever) throughout the accessible pedestal amplitude range (from 0 to 40%). However, when baseline carrier contrast increases (with constant pedestal modulation amplitude), motion thresholds increase, showing that gain control in second-order motion is determined not by the modulator (as in first-order motion) but by the carrier. Note that baseline contrast of the carrier is inherently independent of spatial frequency of the modulator. The drastically different gain-control properties of the two motion systems and prior observations of motion masking and motion saturation are all encompassed in a functional theory. The stimulus inputs to both first- and second-order motion process are normalized by feedforward, shunting gain control. The different properties arise because the modulator is used to control the first-order gain and the carrier is used to control the second-order gain.

High-order upwind schemes for the wave equation on overlapping grids: Maxwell's equations in second-order form

Science.gov (United States)

Angel, Jordan B.; Banks, Jeffrey W.; Henshaw, William D.

2018-01-01

High-order accurate upwind approximations for the wave equation in second-order form on overlapping grids are developed. Although upwind schemes are well established for first-order hyperbolic systems, it was only recently shown by Banks and Henshaw [1] how upwinding could be incorporated into the second-order form of the wave equation. This new upwind approach is extended here to solve the time-domain Maxwell's equations in second-order form; schemes of arbitrary order of accuracy are formulated for general curvilinear grids. Taylor time-stepping is used to develop single-step space-time schemes, and the upwind dissipation is incorporated by embedding the exact solution of a local Riemann problem into the discretization. Second-order and fourth-order accurate schemes are implemented for problems in two and three space dimensions, and overlapping grids are used to treat complex geometry and problems with multiple materials. Stability analysis of the upwind-scheme on overlapping grids is performed using normal mode theory. The stability analysis and computations confirm that the upwind scheme remains stable on overlapping grids, including the difficult case of thin boundary grids when the traditional non-dissipative scheme becomes unstable. The accuracy properties of the scheme are carefully evaluated on a series of classical scattering problems for both perfect conductors and dielectric materials in two and three space dimensions. The upwind scheme is shown to be robust and provide high-order accuracy.

On differential operators generating iterative systems of linear ODEs of maximal symmetry algebra

Science.gov (United States)

Ndogmo, J. C.

2017-06-01

Although every iterative scalar linear ordinary differential equation is of maximal symmetry algebra, the situation is different and far more complex for systems of linear ordinary differential equations, and an iterative system of linear equations need not be of maximal symmetry algebra. We illustrate these facts by examples and derive families of vector differential operators whose iterations are all linear systems of equations of maximal symmetry algebra. Some consequences of these results are also discussed.

Massless second-order tetradic spin-3 and higher-helicity bosons

Energy Technology Data Exchange (ETDEWEB)

Aragone, C; La Roche, H [Universidad Simon Bolivar, Caracas (Venezuela) Dept. de Fisica

1982-11-21

The unique, uniform, second-order formulation of massless bosons of helicity >=3 is presented here in terms of tetradic fields. The actions we find are shown to coincide both with the first-order (tetradic) formulation of Vasiliev and with the symmetric second-order description of Fronsdal. We carefully analyse the gravitational coupling of the spin-3 field and find that tetradic spin-3 matter presents the same translational consistency problem as symmetric matter does. Furthermore, in the curved tetradic case the generalized Lorentz invariance can be restored by the addition of nominimal terms.

High-order quantum algorithm for solving linear differential equations

International Nuclear Information System (INIS)

Berry, Dominic W

2014-01-01

Linear differential equations are ubiquitous in science and engineering. Quantum computers can simulate quantum systems, which are described by a restricted type of linear differential equations. Here we extend quantum simulation algorithms to general inhomogeneous sparse linear differential equations, which describe many classical physical systems. We examine the use of high-order methods (where the error over a time step is a high power of the size of the time step) to improve the efficiency. These provide scaling close to Δt 2 in the evolution time Δt. As with other algorithms of this type, the solution is encoded in amplitudes of the quantum state, and it is possible to extract global features of the solution. (paper)

A Second-Order Conditionally Linear Mixed Effects Model with Observed and Latent Variable Covariates

Science.gov (United States)

Harring, Jeffrey R.; Kohli, Nidhi; Silverman, Rebecca D.; Speece, Deborah L.

2012-01-01

A conditionally linear mixed effects model is an appropriate framework for investigating nonlinear change in a continuous latent variable that is repeatedly measured over time. The efficacy of the model is that it allows parameters that enter the specified nonlinear time-response function to be stochastic, whereas those parameters that enter in a…

High-Order Sparse Linear Predictors for Audio Processing

DEFF Research Database (Denmark)

Giacobello, Daniele; van Waterschoot, Toon; Christensen, Mads Græsbøll

2010-01-01

Linear prediction has generally failed to make a breakthrough in audio processing, as it has done in speech processing. This is mostly due to its poor modeling performance, since an audio signal is usually an ensemble of different sources. Nevertheless, linear prediction comes with a whole set...... of interesting features that make the idea of using it in audio processing not far fetched, e.g., the strong ability of modeling the spectral peaks that play a dominant role in perception. In this paper, we provide some preliminary conjectures and experiments on the use of high-order sparse linear predictors...... in audio processing. These predictors, successfully implemented in modeling the short-term and long-term redundancies present in speech signals, will be used to model tonal audio signals, both monophonic and polyphonic. We will show how the sparse predictors are able to model efficiently the different...

Expanded porphyrins as third order non-linear optical materials ...

Indian Academy of Sciences (India)

WINTEC

function correlations ... An understanding of the structure–function corre- lations of these expanded porphyrins is an important first step for ... where χ (2) and χ (3) are the quadratic χ (2) (first- order) and χ (3) cubic (second-order) susceptibilities.

Stability problems for linear hyperbolic systems

International Nuclear Information System (INIS)

Eckhoff, K.S.

1975-05-01

The stability properties for the trivial solution of a general linear hyperbolic system of partial differential equations of the first order are studied. It is shown that results may be obtained by studying the stability properties of certain systems of ordinary differential equations which can be constructed from the hyperbolic system (the so-called transport equations). In some cases the associated stability problem for the transport equations can in fact be shown to be equivalent to the stability problem for the hyperbolic system, but in general the transport equations will only give the necessary conditions for stability. (Auth.)

Consensus of second-order multi-agent dynamic systems with quantized data

Energy Technology Data Exchange (ETDEWEB)

Guan, Zhi-Hong, E-mail: zhguan@mail.hust.edu.cn [Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, 430074 (China); Meng, Cheng [Department of Control Science and Engineering, Huazhong University of Science and Technology, Wuhan, 430074 (China); Liao, Rui-Quan [Petroleum Engineering College,Yangtze University, Jingzhou, 420400 (China); Zhang, Ding-Xue, E-mail: zdx7773@163.com [Petroleum Engineering College,Yangtze University, Jingzhou, 420400 (China)

2012-01-09

The consensus problem of second-order multi-agent systems with quantized link is investigated in this Letter. Some conditions are derived for the quantized consensus of the second-order multi-agent systems by the stability theory. Moreover, a result characterizing the relationship between the eigenvalues of the Laplacians matrix and the quantized consensus is obtained. Examples are given to illustrate the theoretical analysis. -- Highlights: ► A second-order multi-agent model with quantized data is proposed. ► Two sufficient and necessary conditions are obtained. ► The relationship between the eigenvalues of the Laplacians matrix and the quantized consensus is discovered.

Application of ordinary kriging for interpolation of micro-structured technical surfaces

International Nuclear Information System (INIS)

Raid, Indek; Kusnezowa, Tatjana; Seewig, Jörg

2013-01-01

Kriging is an interpolation technique used in geostatistics. In this paper we present kriging applied in the field of three-dimensional optical surface metrology. Technical surfaces are not always optically cooperative, meaning that measurements of technical surfaces contain invalid data points because of different effects. These data points need to be interpolated to obtain a complete area in order to fulfil further processing. We present an elementary type of kriging, known as ordinary kriging, and apply it to interpolate measurements of different technical surfaces containing different kinds of realistic defects. The result of the interpolation with kriging is compared to six common interpolation techniques: nearest neighbour, natural neighbour, inverse distance to a power, triangulation with linear interpolation, modified Shepard's method and radial basis function. In order to quantify the results of different interpolations, the topographies are compared to defect-free reference topographies. Kriging is derived from a stochastic model that suggests providing an unbiased, linear estimation with a minimized error variance. The estimation with kriging is based on a preceding statistical analysis of the spatial structure of the surface. This comprises the choice and adaptation of specific models of spatial continuity. In contrast to common methods, kriging furthermore considers specific anisotropy in the data and adopts the interpolation accordingly. The gained benefit requires some additional effort in preparation and makes the overall estimation more time-consuming than common methods. However, the adaptation to the data makes this method very flexible and accurate. (paper)

Improved system blind identification based on second-order ...

Indian Academy of Sciences (India)

An improved system blind identification method based on second- order cyclostationary statistics and the properties of group delay, has been ... In the last decade, there has been considerable research on achieving blind identification.

Concordance measures and second order stochastic dominance-portfolio efficiency analysis

Czech Academy of Sciences Publication Activity Database

Kopa, Miloš; Tichý, T.

2012-01-01

Roč. 15, č. 4 (2012), s. 110-120 ISSN 1212-3609 R&D Projects: GA ČR(CZ) GBP402/12/G097 Institutional support: RVO:67985556 Keywords : dependency * concordance * portfolio selection * second order stochastic dominance Subject RIV: BB - Applied Statistics, Operational Research Impact factor: 0.633, year: 2012 http://library.utia.cas.cz/separaty/2013/E/kopa-concordance measures and second order stochastic dominance- portfolio efficiency analysis.pdf

On a class of fourth order linear recurrence equations

Directory of Open Access Journals (Sweden)

Sui-Sun Cheng

1984-01-01

Full Text Available This paper is concerned with sequences that satisfy a class of fourth order linear recurrence equations. Basic properties of such sequences are derived. In addition, we discuss the oscillatory and nonoscillatory behavior of such sequences.

An integral for second-order multiple scattering perturbation theory

International Nuclear Information System (INIS)

Hoffman, G.G.

1997-01-01

This paper presents the closed form evaluation of a six-dimensional integral. The integral arises in the application to many-electron systems of a multiple scattering perturbation expansion at second order when formulated in fourier space. The resulting function can be used for the calculation of both the electron density and the effective one-electron potential in an SCF calculations. The closed form expression derived here greatly facilitates these calculations. In addition, the evaluated integral can be used for the computation of second-order corrections to the open-quotes optimized Thomas-Fermi theory.close quotes 10 refs., 2 figs

Stability and bifurcation of numerical discretization of a second-order delay differential equation with negative feedback

International Nuclear Information System (INIS)

Ding Xiaohua; Su Huan; Liu Mingzhu

2008-01-01

The paper analyzes a discrete second-order, nonlinear delay differential equation with negative feedback. The characteristic equation of linear stability is solved, as a function of two parameters describing the strength of the feedback and the damping in the autonomous system. The existence of local Hopf bifurcations is investigated, and the direction and stability of periodic solutions bifurcating from the Hopf bifurcation of the discrete model are determined by the Hopf bifurcation theory of discrete system. Finally, some numerical simulations are performed to illustrate the analytical results found

Influence of temperature rise distribution in second harmonic generation crystal on intensity distributions of output second harmonic wave

International Nuclear Information System (INIS)

Li Wei; Feng Guoying; Li Gang; Huang Yu; Zhang Qiuhui

2009-01-01

Second-harmonic generation (SHG) of high-intensity laser with an SHG crystal for type I angle phase matching has been studied by the use of a split-step algorithm based on the fast Fourier transform and a fourth-order Runge-Kutta (R-K) integrator. The transverse walk-off effect, diffraction, the second-order and the third-order nonlinear effects have been taken into consideration. Influences of a temperature rise distribution of the SHG crystal on the refractive indices of ordinary wave and extraordinary wave have been discussed. The rules of phase mismatching quantity, intensity distribution of output beam and frequency conversion efficiency varying with the temperature rise distribution of the SHG crystal have been analyzed quantitatively. The calculated results indicate that in a high power frequency conversion system, the temperature rise distribution of SHG crystal would result in the phase mismatching of fundamental and harmonic waves, leading to the variation of intensity distribution of the output beam and the decrease of the conversion efficiency. (authors)

Second Order Sliding Mode Controller Design for Pneumatic Artificial Muscle

Directory of Open Access Journals (Sweden)

Ammar Al-Jodah

2018-01-01

Full Text Available In this paper, first and second order sliding mode controllers are designed for a single link robotic arm actuated by two Pneumatic Artificial Muscles (PAMs. A new mathematical model for the arm has been developed based on the model of large scale pneumatic muscle actuator model. Uncertainty in parameters has been presented and tested for the two controllers. The simulation results of the second-order sliding mode controller proves to have a low tracking error and chattering effect as compared to the first order one. The verification has been done by using MATLAB and Simulink software.

Solution of second order supersymmetrical intertwining relations in Minkowski plane

Energy Technology Data Exchange (ETDEWEB)

Ioffe, M. V., E-mail: m.ioffe@spbu.ru; Kolevatova, E. V., E-mail: e.v.kolev@yandex.ru [Saint Petersburg State University, 7/9 Universitetskaya nab., St. Petersburg 199034 (Russian Federation); Nishnianidze, D. N., E-mail: cutaisi@yahoo.com [Saint Petersburg State University, 7/9 Universitetskaya nab., St. Petersburg 199034 (Russian Federation); Akaki Tsereteli State University, 4600 Kutaisi, Georgia (United States)

2016-08-15

Supersymmetrical (SUSY) intertwining relations are generalized to the case of quantum Hamiltonians in Minkowski space. For intertwining operators (supercharges) of second order in derivatives, the intertwined Hamiltonians correspond to completely integrable systems with the symmetry operators of fourth order in momenta. In terms of components, the intertwining relations correspond to the system of nonlinear differential equations which are solvable with the simplest—constant—ansatzes for the “metric” matrix in second order part of the supercharges. The corresponding potentials are built explicitly both for diagonalizable and nondiagonalizable form of “metric” matrices, and their properties are discussed.

The numerical solution of linear multi-term fractional differential equations: systems of equations

Science.gov (United States)

Edwards, John T.; Ford, Neville J.; Simpson, A. Charles

2002-11-01

In this paper, we show how the numerical approximation of the solution of a linear multi-term fractional differential equation can be calculated by reduction of the problem to a system of ordinary and fractional differential equations each of order at most unity. We begin by showing how our method applies to a simple class of problems and we give a convergence result. We solve the Bagley Torvik equation as an example. We show how the method can be applied to a general linear multi-term equation and give two further examples.

Accelerating distributed average consensus by exploring the information of second-order neighbors

Energy Technology Data Exchange (ETDEWEB)

Yuan Deming [School of Automation, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu (China); Xu Shengyuan, E-mail: syxu02@yahoo.com.c [School of Automation, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu (China); Zhao Huanyu [School of Automation, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu (China); Chu Yuming [Department of Mathematics, Huzhou Teacher' s College, Huzhou 313000, Zhejiang (China)

2010-05-17

The problem of accelerating distributed average consensus by using the information of second-order neighbors in both the discrete- and continuous-time cases is addressed in this Letter. In both two cases, when the information of second-order neighbors is used in each iteration, the network will converge with a speed faster than the algorithm only using the information of first-order neighbors. Moreover, the problem of using partial information of second-order neighbors is considered, and the edges are not chosen randomly from second-order neighbors. In the continuous-time case, the edges are chosen by solving a convex optimization problem which is formed by using the convex relaxation method. In the discrete-time case, for small network the edges are chosen optimally via the brute force method. Finally, simulation examples are provided to demonstrate the effectiveness of the proposed algorithm.

Transport coefficients in second-order non-conformal viscous hydrodynamics

International Nuclear Information System (INIS)

Ryblewski, Radoslaw

2015-01-01

Based on the exact solution of Boltzmann kinetic equation in the relaxation-time approximation, the precision of the two most recent formulations of relativistic second-order non-conformal viscous hydrodynamics (14-moment approximation and causal Chapman-Enskog method), standard Israel-Stewart theory, and anisotropic hydrodynamics framework, in the simple case of one-dimensional Bjorken expansion, is tested. It is demonstrated that the failure of Israel-Stewart theory in reproducing exact solutions of the Boltzmann kinetic equation occurs due to neglecting and/or choosing wrong forms of some of the second-order transport coefficients. In particular, the importance of shear-bulk couplings in the evolution equations for dissipative quantities is shown. One finds that, in the case of the bulk viscous pressure correction, such coupling terms are as important as the corresponding first-order Navier-Stokes term and must be included in order to obtain, at least qualitative, overall agreement with the kinetic theory. (paper)

Numerical integration of asymptotic solutions of ordinary differential equations

Science.gov (United States)

Thurston, Gaylen A.

1989-01-01

Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.

Communication: Surface-facilitated softening of ordinary and vapor-deposited glasses

Science.gov (United States)

Cubeta, Ulyana; Bhattacharya, Deepanjan; Sadtchenko, Vlad

2017-08-01

A common distinction between the ordinary glasses formed by melt cooling and the stable amorphous films formed by vapor deposition is the apparent mechanism of their devitrification. Using quasi-adiabatic, fast scanning calorimetry that is capable of heating rates in excess of 105 K s-1, we have investigated the softening kinetics of micrometer-scale, ordinary glass films of methylbenzene and 2-propanol. At the limit of high heating rates, the transformation mechanism of ordinary glasses is identical to that of their stable vapor-deposited counterparts. In both cases, softening is likely to begin at the sample surface and progress into its bulk via a transformation front. Furthermore, such a surface-facilitated mechanism complies with zero-order, Arrhenius rate law. The activation energy barriers for the softening transformation imply that the kinetics must be defined, at least in part, by the initial thermodynamic and structural state of the samples.

Cascading second-order nonlinear processes in a lithium niobate-on-insulator microdisk.

Science.gov (United States)

Liu, Shijie; Zheng, Yuanlin; Chen, Xianfeng

2017-09-15

Whispering-gallery-mode (WGM) microcavities are very important in both fundamental science and practical applications, among which on-chip second-order nonlinear microresonators play an important role in integrated photonic functionalities. Here we demonstrate resonant second-harmonic generation (SHG) and cascaded third-harmonic generation (THG) in a lithium niobate-on-insulator (LNOI) microdisk resonator. Efficient SHG in the visible range was obtained with only several mW input powers at telecom wavelengths. THG was also observed through a cascading process, which reveals simultaneous phase matching and strong mode coupling in the resonator. Cascading of second-order nonlinear processes gives rise to an effectively large third-order nonlinearity, which makes on-chip second-order nonlinear microresonators a promising frequency converter for integrated nonlinear photonics.

Second order gauge invariant measure of a tidally deformed black hole

Energy Technology Data Exchange (ETDEWEB)

Ahmadi, Nahid, E-mail: nahmadi@ut.ac.ir [Department of Physics, University of Tehran, Kargar Avenue North, Tehran 14395-547 (Iran, Islamic Republic of)

2012-08-01

In this paper, a Lagrangian perturbation theory for the second order treatment of small disturbances of the event horizon in Schwarzchild black holes is introduced. The issue of gauge invariance in the context of general relativistic theory is also discussed. The developments of this paper is a logical continuation of the calculations presented in [1], in which the first order coordinate dependance of the intrinsic and exterinsic geometry of the horizon is examined and the first order gauge invariance of the intrinsic geometry of the horizon is shown. In context of second order perturbation theory, It is shown that the rate of the expansion of the congruence of the horizon generators is invariant under a second order reparametrization; so it can be considered as a measure of tidal perturbation. A generally non-vanishing expression for this observable, which accomodates tidal perturbations and implies nonlinear response of the horizon, is also presented.

A Numerical Scheme for Ordinary Differential Equations Having Time Varying and Nonlinear Coefficients Based on the State Transition Matrix

Science.gov (United States)

Bartels, Robert E.

2002-01-01

A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.

Constrained generalized mechanics. The second-order case

International Nuclear Information System (INIS)

Tapia, V.

1985-01-01

The Dirac formalism for constrained systems is developed for systems described by a Lagrangian depending on up to a second-order time derivatives of the generalized co-ordinates (accelerations). It turns out that for a Lagrangian of this kind differing by a total time derivative from a Lagrangian depending on only up to first-order time-derivatives of the generalized co-ordinates (velocities), both classical mechanics at the Lagrangian level are the same; at the Hamiltonian level the two classical mechanics differ conceptually even when the solutions to both sets of Hamiltonian equations of motion are the same

Large optical second-order nonlinearity of poled WO3-TeO2 glass.

Science.gov (United States)

Tanaka, K; Narazaki, A; Hirao, K

2000-02-15

Second-harmonic generation, one of the second-order nonlinear optical properties of thermally and electrically poled WO>(3)-TeO>(2) glasses, has been examined. We poled glass samples with two thicknesses (0.60 and 0.86 mm) at various temperatures to explore the effects of external electric field strength and poling temperature on second-order nonlinearity. The dependence of second-harmonic intensity on the poling temperature is maximum at a specific poling temperature. A second-order nonlinear susceptibility of 2.1 pm/V was attained for the 0.60-mm-thick glass poled at 250 degrees C. This value is fairly large compared with those for poled silica and tellurite glasses reported thus far. We speculate that the large third-order nonlinear susceptibility of WO>(3)- TeO>(2) glasses gives rise to the large second-order nonlinearity by means of a X((2)) = 3X((3)) E(dc) process.

First- and Second-Order Full-Differential in Edge Analysis of Images

Directory of Open Access Journals (Sweden)

Dong-Mei Pu

2014-01-01

mathematics. We propose and reformulate them with a uniform definition framework. Based on our observation and analysis with the difference, we propose an algorithm to detect the edge from image. Experiments on Corel5K and PASCAL VOC 2007 are done to show the difference between the first order and the second order. After comparison with Canny operator and the proposed first-order differential, the main result is that the second-order differential has the better performance in analysis of changes of the context of images with good selection of control parameter.

The Poisson equation at second order in relativistic cosmology

International Nuclear Information System (INIS)

Hidalgo, J.C.; Christopherson, Adam J.; Malik, Karim A.

2013-01-01

We calculate the relativistic constraint equation which relates the curvature perturbation to the matter density contrast at second order in cosmological perturbation theory. This relativistic ''second order Poisson equation'' is presented in a gauge where the hydrodynamical inhomogeneities coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the framework of General Relativity. We then derive expressions that can be used as the initial conditions of N-body codes for structure formation which probe the observable signature of primordial non-Gaussianity in the statistics of the evolved matter density field

Combined First and Second Order Total Variation Inpainting using Split Bregman

KAUST Repository

Papafitsoros, Konstantinos

2013-07-12

In this article we discuss the implementation of the combined first and second order total variation inpainting that was introduced by Papafitsoros and Schdönlieb. We describe the algorithm we use (split Bregman) in detail, and we give some examples that indicate the difference between pure first and pure second order total variation inpainting.

Combined First and Second Order Total Variation Inpainting using Split Bregman

KAUST Repository

Papafitsoros, Konstantinos; Schoenlieb, Carola Bibiane; Sengul, Bati

2013-01-01

In this article we discuss the implementation of the combined first and second order total variation inpainting that was introduced by Papafitsoros and Schdönlieb. We describe the algorithm we use (split Bregman) in detail, and we give some examples that indicate the difference between pure first and pure second order total variation inpainting.

Numerical analysis of non-linear vibrations of a fractionally damped cylindrical shell under the conditions of combinational internal resonance

Directory of Open Access Journals (Sweden)

Rossikhin Yury A.

2018-01-01

Full Text Available Non-linear damped vibrations of a cylindrical shell embedded into a fractional derivative medium are investigated for the case of the combinational internal resonance, resulting in modal interaction, using two different numerical methods with further comparison of the results obtained. The damping properties of the surrounding medium are described by the fractional derivative Kelvin-Voigt model utilizing the Riemann-Liouville fractional derivatives. Within the first method, the generalized displacements of a coupled set of nonlinear ordinary differential equations of the second order are estimated using numerical solution of nonlinear multi-term fractional differential equations by the procedure based on the reduction of the problem to a system of fractional differential equations. According to the second method, the amplitudes and phases of nonlinear vibrations are estimated from the governing nonlinear differential equations describing amplitude-and-phase modulations for the case of the combinational internal resonance. A good agreement in results is declared.

A connection between the asymptotic iteration method and the continued fractions formalism

International Nuclear Information System (INIS)

Matamala, A.R.; Gutierrez, F.A.; Diaz-Valdes, J.

2007-01-01

In this work, we show that there is a connection between the asymptotic iteration method (a method to solve second order linear ordinary differential equations) and the older method of continued fractions to solve differential equations

Second Order Sliding Mode Controller Design for Pneumatic Artificial Muscle

OpenAIRE

Ammar Al-Jodah; Laith Khames

2018-01-01

In this paper, first and second order sliding mode controllers are designed for a single link robotic arm actuated by two Pneumatic Artificial Muscles (PAMs). A new mathematical model for the arm has been developed based on the model of large scale pneumatic muscle actuator model. Uncertainty in parameters has been presented and tested for the two controllers. The simulation results of the second-order sliding mode controller proves to have a low tracking error and chattering effect as compar...

Growth of meromorphic solutions of higher-order linear differential equations

Directory of Open Access Journals (Sweden)

Wenjuan Chen

2009-01-01

Full Text Available In this paper, we investigate the higher-order linear differential equations with meromorphic coefficients. We improve and extend a result of M.S. Liu and C.L. Yuan, by using the estimates for the logarithmic derivative of a transcendental meromorphic function due to Gundersen, and the extended Winman-Valiron theory which proved by J. Wang and H.X. Yi. In addition, we also consider the nonhomogeneous linear differential equations.

Second-order Born effect in coplanar doubly symmetric (e,2e) collisions for sodium

Energy Technology Data Exchange (ETDEWEB)

Wang, Yang; Jiao, Liguang [Center for Theoretical Atomic and Molecular Physics, Academy of Fundamental and Interdisciplinary Sciences, Harbin Institute of Technology, Harbin 150080 (China); Zhou, Yajun, E-mail: yajunzhou2003@yahoo.com.cn [Center for Theoretical Atomic and Molecular Physics, Academy of Fundamental and Interdisciplinary Sciences, Harbin Institute of Technology, Harbin 150080 (China)

2012-06-18

The second-order distorted wave Born approximation (DWBA) method is employed to investigate the triple differential cross sections (TDCS) of coplanar doubly symmetric (e,2e) collisions for alkali target sodium at excess energies of 6–60 eV. Comparing with the first-order DWBA calculations, the inclusion of second-order Born term in the scattering amplitude improves the degree of agreement with experiments, especially for backward scattering region of TDCS. This indicates the present second-order Born term is capable to give a reasonable correction to DWBA model in studying coplanar symmetric (e,2e) problems in low and intermediate energy range. -- Highlights: ► We consider second-order Born effect in (e,2e) collisions for sodium. ► Our second-order term gives a correct description on the multi scattering process. ► Our second-order DWBA model improves the agreement between theory and experiment.

Effect of Second-Order Hydrodynamics on a Floating Offshore Wind Turbine

Energy Technology Data Exchange (ETDEWEB)

Roald, L.; Jonkman, J.; Robertson, A.

2014-05-01

The design of offshore floating wind turbines uses design codes that can simulate the entire coupled system behavior. At the present, most codes include only first-order hydrodynamics, which induce forces and motions varying with the same frequency as the incident waves. Effects due to second- and higher-order hydrodynamics are often ignored in the offshore industry, because the forces induced typically are smaller than the first-order forces. In this report, first- and second-order hydrodynamic analysis used in the offshore oil and gas industry is applied to two different wind turbine concepts--a spar and a tension leg platform.

Optimization of an intracavity Q-switched solid-state second order Raman laser

Science.gov (United States)

Chen, Zhiqiong; Fu, Xihong; Peng, Hangyu; Zhang, Jun; Qin, Li; Ning, Yongqiang

2017-01-01

In this paper, the model of an intracavity Q-switched second order Raman laser is established, the characteristics of the output 2nd Stokes are simulated. The dynamic balance mechanism among intracavity conversion rates of stimulated emission, first order Raman and second order Raman is obtained. Finally, optimization solutions for increasing output 2nd Stokes pulse energy are proposed.

PID control of second-order systems with hysteresis

NARCIS (Netherlands)

Jayawardhana, Bayu; Logemann, Hartmut; Ryan, Eugene P.

2008-01-01

The efficacy of proportional, integral and derivative (PID) control for set point regulation and disturbance rejection is investigated in a context of second-order systems with hysteretic components. Two basic structures are studied: in the first, the hysteretic component resides (internally) in the

Effect of second-order and fully nonlinear wave kinematics on a tension-leg-platform wind turbine in extreme wave conditions

DEFF Research Database (Denmark)

Pegalajar Jurado, Antonio Manuel; Borg, Michael; Robertson, Amy

2017-01-01

In this study, we assess the impact of different wave kinematics models on the dynamic response of a tension-leg-platform wind turbine. Aero-hydro-elastic simulations of the floating wind turbine are carried out employing linear, second-order, and fully nonlinear kinematics using the Morison equa...... damping coefficients in the model by a more detailed, customizable definition of the user-defined numerical damping....

Numerical solution of ordinary differential equations. For classical, relativistic and nano systems

International Nuclear Information System (INIS)

Greenspan, D.

2006-01-01

An up-to-date survey on numerical solutions with theory, intuition and applications. Ordinary differential equations (ODE) play a significant role in mathematics, physics and engineering sciences, and thus are part of relevant college and university courses. Many problems, however, both traditional and modern, do not possess exact solutions, and must be treated numerically. Usually this is done with software packages, but for this to be efficient requires a sound understanding of the mathematics involved. This work meets the need for an affordable textbook that helps in understanding numerical solutions of ODE. Carefully structured by an experienced textbook author, it provides a survey of ODE for various applications, both classical and modern, including such special applications as relativistic and nano systems. The examples are carefully explained and compiled into an algorithm, each of which is presented generically, independent of a specific programming language, while each chapter is rounded off with exercises. The text meets the demands of MA200 courses and of the newly created Numerical Solution of Differential Equations courses, making it ideal for both students and lecturers in physics, mathematics, mechanical engineering, electrical engineering, as well as for physicists, mathematicians, engineers, and electrical engineers. From the Contents - Euler's Method - Runge-Kutta Methods - The Method of Taylor Expansions - Large Second Order Systems with Application to Nano Systems - Completely Conservative, Covariant Numerical Methodology - Instability - Numerical Solution of Tridiagonal Linear Algebraic Systems and Related Nonlinear Systems - Approximate Solution of Boundary Value Problems - Special Relativistic Motion - Special Topics - Appendix: Basic Matrix Operations - Bibliography. (orig.) (orig.)

[Application of ordinary Kriging method in entomologic ecology].

Science.gov (United States)

Zhang, Runjie; Zhou, Qiang; Chen, Cuixian; Wang, Shousong

2003-01-01

Geostatistics is a statistic method based on regional variables and using the tool of variogram to analyze the spatial structure and the patterns of organism. In simulating the variogram within a great range, though optimal simulation cannot be obtained, the simulation method of a dialogue between human and computer can be used to optimize the parameters of the spherical models. In this paper, the method mentioned above and the weighted polynomial regression were utilized to simulate the one-step spherical model, the two-step spherical model and linear function model, and the available nearby samples were used to draw on the ordinary Kriging procedure, which provided a best linear unbiased estimate of the constraint of the unbiased estimation. The sum of square deviation between the estimating and measuring values of varying theory models were figured out, and the relative graphs were shown. It was showed that the simulation based on the two-step spherical model was the best simulation, and the one-step spherical model was better than the linear function model.

A unified approach to fixed-order controller design via linear matrix inequalities

Directory of Open Access Journals (Sweden)

T. Iwasaki

1995-01-01

Full Text Available We consider the design of fixed-order (or low-order linear controllers which meet certain performance and/or robustness specifications. The following three problems are considered; covariance control as a nominal performance problem,-stabilization as a robust stabilization problem, and robust L∞ control problem as a robust performance problem. All three control problems are converted to a single linear algebra problem of solving a linear matrix inequality (LMI of the type BGC+(BGCT+Q<0 for the unknown matrix G. Thus this paper addresses the fixed-order controller design problem in a unified way. Necessary and sufficient conditions for the existence of a fixed-order controller which satisfies the design specifications for each problem are derived, and an explicit controller formula is given. In any case, the resulting problem is shown to be a search for a (structured positive definite matrix X such that X∈1 and X−1∈2 where 1 and 2 are convex sets defined by LMIs. Computational aspects of the nonconvex LMI problem are discussed.

Successive substitution one-leg hybrid P-stable LMM for initial value ...

African Journals Online (AJOL)

This paper derives P-stable successive substitution one-leg hybrid linear multistep methods for the numerical solution of second order initial value problems in ordinary differential equations without explicit first order derivative. The methods are demonstrated by a numerical example also considered by Fatunla, et al (1997) ...

Global weighted estimates for second-order nondivergence elliptic ...

Indian Academy of Sciences (India)

Fengping Yao

2018-03-21

Mar 21, 2018 ... One of the key a priori estimates in the theory of second-order elliptic .... It is well known that the maximal functions satisfy strong p–p .... Here we prove the following auxiliary result, which will be a crucial ingredient in the proof.

Validity of second order analysis of superdense matter

International Nuclear Information System (INIS)

Bowers, R.L.; Gleeson, A.M.; Pedigo, R.D.

1975-01-01

The limitations of relativistic calculations of the properties of superdense matter obtained from strictly second order terms is discussed. Extension of the model to overcome these limitations leads to serious complications which can only be overcome by a fully self-consistent treatment. (U.S.)

Existence of solutions for second-order evolution inclusions

Directory of Open Access Journals (Sweden)

Nikolaos S. Papageorgiou

1994-01-01

Full Text Available In this paper we examine second-order nonlinear evolution inclusions and prove two existence theorems; one with a convex-valued orientor field and the other with a nonconvex-valued field. An example of a hyperbolic partial differential inclusion is also presented.

Oscillation of solutions of some higher order linear differential equations

Directory of Open Access Journals (Sweden)

Hong-Yan Xu

2009-11-01

Full Text Available In this paper, we deal with the order of growth and the hyper order of solutions of higher order linear differential equations $$f^{(k}+B_{k-1}f^{(k-1}+\\cdots+B_1f'+B_0f=F$$ where $B_j(z (j=0,1,\\ldots,k-1$ and $F$ are entire functions or polynomials. Some results are obtained which improve and extend previous results given by Z.-X. Chen, J. Wang, T.-B. Cao and C.-H. Li.

Enhancing second-order conditioning with lesions of the basolateral amygdala.

Science.gov (United States)

Holland, Peter C

2016-04-01

Because the occurrence of primary reinforcers in natural environments is relatively rare, conditioned reinforcement plays an important role in many accounts of behavior, including pathological behaviors such as the abuse of alcohol or drugs. As a result of pairing with natural or drug reinforcers, initially neutral cues acquire the ability to serve as reinforcers for subsequent learning. Accepting a major role for conditioned reinforcement in everyday learning is complicated by the often-evanescent nature of this phenomenon in the laboratory, especially when primary reinforcers are entirely absent from the test situation. Here, I found that under certain conditions, the impact of conditioned reinforcement could be extended by lesions of the basolateral amygdala (BLA). Rats received first-order Pavlovian conditioning pairings of 1 visual conditioned stimulus (CS) with food prior to receiving excitotoxic or sham lesions of the BLA, and first-order pairings of another visual CS with food after that surgery. Finally, each rat received second-order pairings of a different auditory cue with each visual first-order CS. As in prior studies, relative to sham-lesioned control rats, lesioned rats were impaired in their acquisition of second-order conditioning to the auditory cue paired with the first-order CS that was trained after surgery. However, lesioned rats showed enhanced and prolonged second-order conditioning to the auditory cue paired with the first-order CS that was trained before amygdala damage was made. Implications for an enhanced role for conditioned reinforcement by drug-related cues after drug-induced alterations in neural plasticity are discussed. (PsycINFO Database Record (c) 2016 APA, all rights reserved).

On Application of Second Order Sliding Mode Control to Electro-Hydraulic Systems

DEFF Research Database (Denmark)

Schmidt, Lasse; Andersen, Torben Ole; Pedersen, Henrik C.

2014-01-01

This paper discusses the application of second order mode controls to hydraulic valve-cylinder drives with a special focus on the limitations resulting from nonlinear dynamic effects in flow control valves. Second order sliding mode algorithms appear highly attractive in the successive...

Second-order Hydrodynamics in QCD at NLO arXiv

CERN Document Server

Ghiglieri, Jacopo; Teaney, Derek

We compute the hydrodynamic relaxation times $\\tau_\\pi$ and $\\tau_j$ for hot QCD at next-to-leading order in the coupling by using kinetic theory. We show that certain dimensionless ratios of second-order to first-order transport coefficients obey bounds which apply whenever a kinetic theory description is possible; the computed values lie somewhat above these bounds. Strongly coupled theories with holographic duals strongly violate these bounds, highlighting their distance from a quasiparticle description.

About a definition of metric over an abelian linearly ordered group

Directory of Open Access Journals (Sweden)

Bice Cavallo

2012-06-01

Full Text Available A G-metric over an abelian linearly ordered group G = (G,⊙,≤ is a binary operation, d G , verifying suitable properties. We consider a particular G metric derived by the group operation ⊙ and the total weak order ≤, and show that it provides a base for the order topology associated to G.

Exploiting Small Leakages in Masks to Turn a Second-Order Attack into a First-Order Attack and Improved Rotating Substitution Box Masking with Linear Code Cosets

Science.gov (United States)

DeTrano, Alexander; Karimi, Naghmeh; Karri, Ramesh; Guo, Xiaofei; Carlet, Claude; Guilley, Sylvain

2015-01-01

Masking countermeasures, used to thwart side-channel attacks, have been shown to be vulnerable to mask-extraction attacks. State-of-the-art mask-extraction attacks on the Advanced Encryption Standard (AES) algorithm target S-Box recomputation schemes but have not been applied to scenarios where S-Boxes are precomputed offline. We propose an attack targeting precomputed S-Boxes stored in nonvolatile memory. Our attack targets AES implemented in software protected by a low entropy masking scheme and recovers the masks with 91% success rate. Recovering the secret key requires fewer power traces (in fact, by at least two orders of magnitude) compared to a classical second-order attack. Moreover, we show that this attack remains viable in a noisy environment or with a reduced number of leakage points. Eventually, we specify a method to enhance the countermeasure by selecting a suitable coset of the masks set. PMID:26491717

A Fifth Order Hybrid Linear Multistep method For the Direct Solution ...

African Journals Online (AJOL)

A linear multistep hybrid method (LMHM)with continuous coefficients isconsidered and directly applied to solve third order initial and boundary value problems (IBVPs). The continuous method is used to obtain Multiple Finite Difference Methods (MFDMs) (each of order 5) which are combined as simultaneous numerical ...

Second-Order Controllability of Multi-Agent Systems with Multiple Leaders

International Nuclear Information System (INIS)

Liu Bo; Han Xiao; Shi Yun-Tao; Su Hou-Sheng

2016-01-01

This paper proposes a new second-order continuous-time multi-agent model and analyzes the controllability of second-order multi-agent system with multiple leaders based on the asymmetric topology. This paper considers the more general case: velocity coupling topology is different from location coupling topology. Some sufficient and necessary conditions are presented for the controllability of the system with multiple leaders. In addition, the paper studies the controllability of the system with velocity damping gain. Simulation results are given to illustrate the correctness of theoretical results. (paper)

Practical considerations for a second-order directional hearing aid microphone system

Science.gov (United States)

Thompson, Stephen C.

2003-04-01

First-order directional microphone systems for hearing aids have been available for several years. Such a system uses two microphones and has a theoretical maximum free-field directivity index (DI) of 6.0 dB. A second-order microphone system using three microphones could provide a theoretical increase in free-field DI to 9.5 dB. These theoretical maximum DI values assume that the microphones have exactly matched sensitivities at all frequencies of interest. In practice, the individual microphones in the hearing aid always have slightly different sensitivities. For the small microphone separation necessary to fit in a hearing aid, these sensitivity matching errors degrade the directivity from the theoretical values, especially at low frequencies. This paper shows that, for first-order systems the directivity degradation due to sensitivity errors is relatively small. However, for second-order systems with practical microphone sensitivity matching specifications, the directivity degradation below 1 kHz is not tolerable. A hybrid order directive system is proposed that uses first-order processing at low frequencies and second-order directive processing at higher frequencies. This hybrid system is suggested as an alternative that could provide improved directivity index in the frequency regions that are important to speech intelligibility.

VODE, Variable Coefficient Ordinary Differential Equations (ODE) Solver

International Nuclear Information System (INIS)

Brown, P.N.; Hindmarsh, A.C.; Byrne, G.D.

2002-01-01

1 - Description of program or function: VODE is a package of subroutines for the numerical solution of the initial-value problem for systems of first-order ordinary differential equations. The package can be used for either stiff or non-stiff systems. In the stiff case, the Jacobian matrix is treated as full or banded. An algorithm is included for saving and reusing the Jacobian matrix under certain conditions. If storage is limited, this option may be suppressed. 2 - Method of solution - VODE uses the variable-order, variable- coefficient Adams-Moulton method for non-stiff systems and the variable-order, fixed-leading-coefficient Backward Differentiation Formula (BDF) method for stiff systems

Discrete second order trajectory generator with nonlinear constraints

NARCIS (Netherlands)

Morselli, R.; Zanasi, R.; Stramigioli, Stefano

2005-01-01

A discrete second order trajectory generator for motion control systems is presented. The considered generator is a nonlinear system which receives as input a raw reference signal and provides as output a smooth reference signal satisfying nonlinear constraints on the output derivatives as UM-(x) ≤

The qualitative theory of ordinary differential equations an introduction

CERN Document Server

Brauer, Fred

1989-01-01

""This is a very good book ... with many well-chosen examples and illustrations."" - American Mathematical MonthlyThis highly regarded text presents a self-contained introduction to some important aspects of modern qualitative theory for ordinary differential equations. It is accessible to any student of physical sciences, mathematics or engineering who has a good knowledge of calculus and of the elements of linear algebra. In addition, algebraic results are stated as needed; the less familiar ones are proved either in the text or in appendixes.The topics covered in the first three chapters a

The second-order interference of two independent single-mode He-Ne lasers

Science.gov (United States)

Liu, Jianbin; Le, Mingnan; Bai, Bin; Wang, Wentao; Chen, Hui; Zhou, Yu; Li, Fu-li; Xu, Zhuo

2015-09-01

The second-order spatial and temporal interference patterns with two independent single-mode continuous-wave He-Ne lasers are observed when these two lasers are incident to two adjacent input ports of a 1:1 non-polarizing beam splitter, respectively. Two-photon interference based on the superposition principle in Feynman's path integral theory is employed to interpret the experimental results. The conditions to observe the second-order interference pattern with two independent single-mode continuous-wave lasers are discussed. It is concluded that frequency stability is important to observe the second-order interference pattern with two independent light beams.

Zeroth-order design report for the next linear collider. Volume 2

International Nuclear Information System (INIS)

Raubenheimer, T.O.

1996-05-01

This Zeroth-Order Design Report (ZDR) for the Next Linear Collider (NLC) has been completed as a feasibility study for a TeV-scale linear collider that incorporates a room-temperature accelerator powered by rf microwaves at 11.424 GHz--similar to that presently used in the SLC, but at four times the rf frequency. The purpose of this study is to examine the complete systems of such a collider, to understand how the parts fit together, and to make certain that every required piece has been included. The ''design'' presented here is not fully engineered in any sense, but to be assured that the NLC can be built, attention has been given to a number of critical components and issues that present special challenges. More engineering and development of a number of mechanical and electrical systems remain to be done, but the conclusion of this study is that indeed the NLC is technically feasible and can be expected to reach the performance levels required to perform research at the TeV energy scale. Volume II covers the following: collimation systems; IP switch and big bend; final focus; the interaction region; multiple bunch issues; control systems; instrumentation; machine protection systems; NLC reliability considerations; NLC conventional facilities. Also included are four appendices on the following topics: An RF power source upgrade to the NLC; a second interaction region for gamma-gamma, gamma-electron; ground motion: theory and measurement; and beam-based feedback: theory and implementation

Zeroth-order design report for the next linear collider. Volume 2

Energy Technology Data Exchange (ETDEWEB)

Raubenheimer, T.O. [ed.

1996-05-01

This Zeroth-Order Design Report (ZDR) for the Next Linear Collider (NLC) has been completed as a feasibility study for a TeV-scale linear collider that incorporates a room-temperature accelerator powered by rf microwaves at 11.424 GHz--similar to that presently used in the SLC, but at four times the rf frequency. The purpose of this study is to examine the complete systems of such a collider, to understand how the parts fit together, and to make certain that every required piece has been included. The ``design`` presented here is not fully engineered in any sense, but to be assured that the NLC can be built, attention has been given to a number of critical components and issues that present special challenges. More engineering and development of a number of mechanical and electrical systems remain to be done, but the conclusion of this study is that indeed the NLC is technically feasible and can be expected to reach the performance levels required to perform research at the TeV energy scale. Volume II covers the following: collimation systems; IP switch and big bend; final focus; the interaction region; multiple bunch issues; control systems; instrumentation; machine protection systems; NLC reliability considerations; NLC conventional facilities. Also included are four appendices on the following topics: An RF power source upgrade to the NLC; a second interaction region for gamma-gamma, gamma-electron; ground motion: theory and measurement; and beam-based feedback: theory and implementation.

Second order finite-difference ghost-point multigrid methods for elliptic problems with discontinuous coefficients on an arbitrary interface

Science.gov (United States)

Coco, Armando; Russo, Giovanni

2018-05-01

In this paper we propose a second-order accurate numerical method to solve elliptic problems with discontinuous coefficients (with general non-homogeneous jumps in the solution and its gradient) in 2D and 3D. The method consists of a finite-difference method on a Cartesian grid in which complex geometries (boundaries and interfaces) are embedded, and is second order accurate in the solution and the gradient itself. In order to avoid the drop in accuracy caused by the discontinuity of the coefficients across the interface, two numerical values are assigned on grid points that are close to the interface: a real value, that represents the numerical solution on that grid point, and a ghost value, that represents the numerical solution extrapolated from the other side of the interface, obtained by enforcing the assigned non-homogeneous jump conditions on the solution and its flux. The method is also extended to the case of matrix coefficient. The linear system arising from the discretization is solved by an efficient multigrid approach. Unlike the 1D case, grid points are not necessarily aligned with the normal derivative and therefore suitable stencils must be chosen to discretize interface conditions in order to achieve second order accuracy in the solution and its gradient. A proper treatment of the interface conditions will allow the multigrid to attain the optimal convergence factor, comparable with the one obtained by Local Fourier Analysis for rectangular domains. The method is robust enough to handle large jump in the coefficients: order of accuracy, monotonicity of the errors and good convergence factor are maintained by the scheme.

Linear models of coregionalization for multivariate lattice data: Order-dependent and order-free cMCARs.

Science.gov (United States)

MacNab, Ying C

2016-08-01

This paper concerns with multivariate conditional autoregressive models defined by linear combination of independent or correlated underlying spatial processes. Known as linear models of coregionalization, the method offers a systematic and unified approach for formulating multivariate extensions to a broad range of univariate conditional autoregressive models. The resulting multivariate spatial models represent classes of coregionalized multivariate conditional autoregressive models that enable flexible modelling of multivariate spatial interactions, yielding coregionalization models with symmetric or asymmetric cross-covariances of different spatial variation and smoothness. In the context of multivariate disease mapping, for example, they facilitate borrowing strength both over space and cross variables, allowing for more flexible multivariate spatial smoothing. Specifically, we present a broadened coregionalization framework to include order-dependent, order-free, and order-robust multivariate models; a new class of order-free coregionalized multivariate conditional autoregressives is introduced. We tackle computational challenges and present solutions that are integral for Bayesian analysis of these models. We also discuss two ways of computing deviance information criterion for comparison among competing hierarchical models with or without unidentifiable prior parameters. The models and related methodology are developed in the broad context of modelling multivariate data on spatial lattice and illustrated in the context of multivariate disease mapping. The coregionalization framework and related methods also present a general approach for building spatially structured cross-covariance functions for multivariate geostatistics. © The Author(s) 2016.

A new numerical approximation of the fractal ordinary differential equation

Science.gov (United States)

Atangana, Abdon; Jain, Sonal

2018-02-01

The concept of fractal medium is present in several real-world problems, for instance, in the geological formation that constitutes the well-known subsurface water called aquifers. However, attention has not been quite devoted to modeling for instance, the flow of a fluid within these media. We deem it important to remind the reader that the concept of fractal derivative is not to represent the fractal sharps but to describe the movement of the fluid within these media. Since this class of ordinary differential equations is highly complex to solve analytically, we present a novel numerical scheme that allows to solve fractal ordinary differential equations. Error analysis of the method is also presented. Application of the method and numerical approximation are presented for fractal order differential equation. The stability and the convergence of the numerical schemes are investigated in detail. Also some exact solutions of fractal order differential equations are presented and finally some numerical simulations are presented.

Second-home electricity consumption

International Nuclear Information System (INIS)

Andersen, Frits M.; Christensen, Morten S.; Jensen, Ole Michael; Kofoed, Niels-Ulrik; Morthorst, Poul Erik

2008-01-01

In Denmark, electricity consumption in first and second homes has developed quite differently. Since 1990, electricity consumption in ordinary residences has grown moderately, while consumption in weekend and second homes has increased considerably. In turn, this development has been blamed on a growing number of luxury cottages, new legislation permitting senior citizens to have their permanent address in their second home and a growing number of electric appliances. In order to examine the growing electricity consumption in second homes and to estimate future demand, a multidisciplinary study combining top-down and bottom-up analyses was conducted, i.e., combining models using aggregated economic parameters and feasibility studies using technical parameters, respectively. The top-down estimation showed that changes in electricity consumption in second homes correlate to changes in income. The bottom-up estimation showed that consumption was mainly affected by the frequency with which second homes were used in the winter time. This indicates that additional second homes, increased full-time use and intensified use of electric appliances are the main reasons for the observed increases in electricity consumption. Luxury tourism use and senior citizens' that use a few per cent of the second homes as their home contribute to a minor degree to the overall increase of electricity consumption. Scenarios show that this development may accelerate with increased leisure time, increased use and more permanent settlement in second homes

Comparison of some biased estimation methods (including ordinary subset regression) in the linear model

Science.gov (United States)

Sidik, S. M.

1975-01-01

Ridge, Marquardt's generalized inverse, shrunken, and principal components estimators are discussed in terms of the objectives of point estimation of parameters, estimation of the predictive regression function, and hypothesis testing. It is found that as the normal equations approach singularity, more consideration must be given to estimable functions of the parameters as opposed to estimation of the full parameter vector; that biased estimators all introduce constraints on the parameter space; that adoption of mean squared error as a criterion of goodness should be independent of the degree of singularity; and that ordinary least-squares subset regression is the best overall method.

Synchronization from second order network connectivity statistics

Directory of Open Access Journals (Sweden)

Liqiong eZhao

2011-07-01

Full Text Available We investigate how network structure can influence the tendency for a neuronal network to synchronize, or its synchronizability, independent of the dynamical model for each neuron. The synchrony analysis takes advantage of the framework of second order networks (SONETs, which defines four second order connectivity statistics based on the relative frequency of two-connection network motifs. The analysis identifies two of these statistics, convergent connections and chain connections, as highly influencing the synchrony. Simulations verify that synchrony decreases with the frequency of convergent connections and increases with the frequency of chain connections. These trends persist with simulations of multiple models for the neuron dynamics and for different types of networks. Surprisingly, divergent connections, which determine the fraction of shared inputs, do not strongly influence the synchrony. The critical role of chains, rather than divergent connections, in influencing synchrony can be explained by a pool and redistribute mechanism. The pooling of many inputs averages out independent fluctuations, amplifying weak correlations in the inputs. With increased chain connections, neurons with many inputs tend to have many outputs. Hence, chains ensure that the amplified correlations in the neurons with many inputs are redistributed throughout the network, enhancing the development of synchrony across the network.

A unified approach to fixed-order controller design via linear matrix inequalities

Directory of Open Access Journals (Sweden)

Iwasaki T.

1995-01-01

Full Text Available We consider the design of fixed-order (or low-order linear controllers which meet certain performance and/or robustness specifications. The following three problems are considered; covariance control as a nominal performance problem, 𝒬 -stabilization as a robust stabilization problem, and robust L ∞ control problem as a robust performance problem. All three control problems are converted to a single linear algebra problem of solving a linear matrix inequality (LMI of the type B G C + ( B G C T + Q < 0 for the unknown matrix G . Thus this paper addresses the fixed-order controller design problem in a unified way. Necessary and sufficient conditions for the existence of a fixed-order controller which satisfies the design specifications for each problem are derived, and an explicit controller formula is given. In any case, the resulting problem is shown to be a search for a (structured positive definite matrix X such that X ∈ 𝒞 1 and X − 1 ∈ 𝒞 2 where 𝒞 1 and 𝒞 2 are convex sets defined by LMIs. Computational aspects of the nonconvex LMI problem are discussed.

Parameters and Fractional Differentiation Orders Estimation for Linear Continuous-Time Non-Commensurate Fractional Order Systems

KAUST Repository

Belkhatir, Zehor; Laleg-Kirati, Taous-Meriem

2017-01-01

This paper proposes a two-stage estimation algorithm to solve the problem of joint estimation of the parameters and the fractional differentiation orders of a linear continuous-time fractional system with non-commensurate orders. The proposed algorithm combines the modulating functions and the first-order Newton methods. Sufficient conditions ensuring the convergence of the method are provided. An error analysis in the discrete case is performed. Moreover, the method is extended to the joint estimation of smooth unknown input and fractional differentiation orders. The performance of the proposed approach is illustrated with different numerical examples. Furthermore, a potential application of the algorithm is proposed which consists in the estimation of the differentiation orders of a fractional neurovascular model along with the neural activity considered as input for this model.

Parameters and Fractional Differentiation Orders Estimation for Linear Continuous-Time Non-Commensurate Fractional Order Systems

KAUST Repository

Belkhatir, Zehor

2017-05-31

This paper proposes a two-stage estimation algorithm to solve the problem of joint estimation of the parameters and the fractional differentiation orders of a linear continuous-time fractional system with non-commensurate orders. The proposed algorithm combines the modulating functions and the first-order Newton methods. Sufficient conditions ensuring the convergence of the method are provided. An error analysis in the discrete case is performed. Moreover, the method is extended to the joint estimation of smooth unknown input and fractional differentiation orders. The performance of the proposed approach is illustrated with different numerical examples. Furthermore, a potential application of the algorithm is proposed which consists in the estimation of the differentiation orders of a fractional neurovascular model along with the neural activity considered as input for this model.

First-order system least-squares for second-order elliptic problems with discontinuous coefficients: Further results

Energy Technology Data Exchange (ETDEWEB)

Bloechle, B.; Manteuffel, T.; McCormick, S.; Starke, G.

1996-12-31

Many physical phenomena are modeled as scalar second-order elliptic boundary value problems with discontinuous coefficients. The first-order system least-squares (FOSLS) methodology is an alternative to standard mixed finite element methods for such problems. The occurrence of singularities at interface corners and cross-points requires that care be taken when implementing the least-squares finite element method in the FOSLS context. We introduce two methods of handling the challenges resulting from singularities. The first method is based on a weighted least-squares functional and results in non-conforming finite elements. The second method is based on the use of singular basis functions and results in conforming finite elements. We also share numerical results comparing the two approaches.

Second-order phase transition in gφ42 theory

International Nuclear Information System (INIS)

Ganbold, G.; Efimov, G.V.

1993-08-01

We have suggested a regular scheme for calculating systematically the leading term and next corrections to it up to the fourth order for the effective potential in the scalar φ 4 2 theory. The obtained results give evidence in favour of a second-order phase transition at (g/2πm 2 ) crit ≅ 0.9 in the theory under consideration. (author). 18 refs, 1 fig

1/N-expansion of the non-linear σ-model: The first three orders

International Nuclear Information System (INIS)

Flyvbjerg, H.; Varsted, S.

1990-01-01

The two-point function of the O(N)-symmetric non-linear σ-model is expanded in 1/N, keeping terms of three leading orders. The mass gap and the magnetic susceptibility are obtained from the two-point function. They are evaluated on square lattices for N=3 and N=4. The systematic errors of the 1/N-series truncated after the first, second, or third term are found by using recent high precision Monte Carlo results as bench marks. For all three truncations, we find systematic errors which are smaller than the expected magnitude of neglected terms, both for the mass gap and for the susceptibility. This result is uniform in the inverse coupling β, and valid for N as small as 3. We conclude that the 1/N-series approach the exact results as rapidly as one could ever hope for. (orig.)

Second Order Optimality in Markov Decision Chains

Czech Academy of Sciences Publication Activity Database

Sladký, Karel

2017-01-01

Roč. 53, č. 6 (2017), s. 1086-1099 ISSN 0023-5954 R&D Projects: GA ČR GA15-10331S Institutional support: RVO:67985556 Keywords : Markov decision chains * second order optimality * optimalilty conditions for transient, discounted and average models * policy and value iterations Subject RIV: BB - Applied Statistics, Operational Research OBOR OECD: Statistics and probability Impact factor: 0.379, year: 2016 http://library.utia.cas.cz/separaty/2017/E/sladky-0485146.pdf

Physical origin of third order non-linear optical response of porphyrin nanorods

International Nuclear Information System (INIS)

Mongwaketsi, N.; Khamlich, S.; Pranaitis, M.; Sahraoui, B.; Khammar, F.; Garab, G.; Sparrow, R.; Maaza, M.

2012-01-01

The non-linear optical properties of porphyrin nanorods were studied using Z-scan, Second and Third harmonic generation techniques. We investigated in details the heteroaggregate behaviour formation of [H 4 TPPS 4 ] 2- and [SnTPyP] 2+ mixture by means of the UV-VIS spectroscopy and aggregates structure and morphology by transmission electron microscopy. The porphyrin nanorods under investigation were synthesized by self assembly and molecular recognition method. They have been optimized in view of future application in the construction of the light harvesting system. The focus of this study was geared towards understanding the influence of the type of solvent used on these porphyrins nanorods using spectroscopic and microscopic techniques. Highlights: ► We synthesized porphyrin nanorods by self assembly and molecular recognition method. ► TEM images confirmed solid cylindrical shapes. ► UV-VIS spectroscopy showed the decrease in the absorbance peaks of the precursors. ► The enhanced third-order nonlinearities were observed.

Hybrid approximations via second order combined dynamic derivatives on time scales

Directory of Open Access Journals (Sweden)

Qin Sheng

2007-09-01

Full Text Available This article focuses on the approximation of conventional second order derivative via the combined (diamond-$\\alpha$ dynamic derivative on time scales with necessary smoothness conditions embedded. We will show the constraints under which the second order dynamic derivative provides a consistent approximation to the conventional second derivative; the cases where the dynamic derivative approximates the derivative only via a proper modification of the existing formula; and the situations in which the dynamic derivative can never approximate consistently even with the help of available structure correction methods. Constructive error analysis will be given via asymptotic expansions for practical hybrid modeling and computational applications.

CMB anisotropies at all orders: the non-linear Sachs-Wolfe formula

International Nuclear Information System (INIS)

Roldan, Omar

2017-01-01

We obtain the non-linear generalization of the Sachs-Wolfe + integrated Sachs-Wolfe (ISW) formula describing the CMB temperature anisotropies. Our formula is valid at all orders in perturbation theory, is also valid in all gauges and includes scalar, vector and tensor modes. A direct consequence of our results is that the maps of the logarithmic temperature anisotropies are much cleaner than the usual CMB maps, because they automatically remove many secondary anisotropies. This can for instance, facilitate the search for primordial non-Gaussianity in future works. It also disentangles the non-linear ISW from other effects. Finally, we provide a method which can iteratively be used to obtain the lensing solution at the desired order.

CMB anisotropies at all orders: the non-linear Sachs-Wolfe formula

Energy Technology Data Exchange (ETDEWEB)

Roldan, Omar, E-mail: oaroldan@if.ufrj.br [Instituto de Física, Universidade Federal do Rio de Janeiro, 21941-972, Rio de Janeiro, RJ (Brazil)

2017-08-01

We obtain the non-linear generalization of the Sachs-Wolfe + integrated Sachs-Wolfe (ISW) formula describing the CMB temperature anisotropies. Our formula is valid at all orders in perturbation theory, is also valid in all gauges and includes scalar, vector and tensor modes. A direct consequence of our results is that the maps of the logarithmic temperature anisotropies are much cleaner than the usual CMB maps, because they automatically remove many secondary anisotropies. This can for instance, facilitate the search for primordial non-Gaussianity in future works. It also disentangles the non-linear ISW from other effects. Finally, we provide a method which can iteratively be used to obtain the lensing solution at the desired order.

7 CFR 28.407 - Good Ordinary Color.

Science.gov (United States)

2010-01-01

... 7 Agriculture 2 2010-01-01 2010-01-01 false Good Ordinary Color. 28.407 Section 28.407 Agriculture..., TESTING, AND STANDARDS Standards Official Cotton Standards of the United States for the Color Grade of American Upland Cotton § 28.407 Good Ordinary Color. Good Ordinary Color is color which is within the range...

Estimating epidemic arrival times using linear spreading theory

Science.gov (United States)

Chen, Lawrence M.; Holzer, Matt; Shapiro, Anne

2018-01-01

We study the dynamics of a spatially structured model of worldwide epidemics and formulate predictions for arrival times of the disease at any city in the network. The model is composed of a system of ordinary differential equations describing a meta-population susceptible-infected-recovered compartmental model defined on a network where each node represents a city and the edges represent the flight paths connecting cities. Making use of the linear determinacy of the system, we consider spreading speeds and arrival times in the system linearized about the unstable disease free state and compare these to arrival times in the nonlinear system. Two predictions are presented. The first is based upon expansion of the heat kernel for the linearized system. The second assumes that the dominant transmission pathway between any two cities can be approximated by a one dimensional lattice or a homogeneous tree and gives a uniform prediction for arrival times independent of the specific network features. We test these predictions on a real network describing worldwide airline traffic.

Second-order domain derivative of normal-dependent boundary integrals

KAUST Repository

Balzer, Jonathan

2010-01-01

Numerous reconstruction tasks in (optical) surface metrology allow for a variational formulation. The occurring boundary integrals may be interpreted as shape functions. The paper is concerned with the second-order analysis of such functions. Shape

An Alternating Direction Method for Convex Quadratic Second-Order Cone Programming with Bounded Constraints

Directory of Open Access Journals (Sweden)

Xuewen Mu

2015-01-01

quadratic programming over second-order cones and a bounded set. At each iteration, we only need to compute the metric projection onto the second-order cones and the projection onto the bound set. The result of convergence is given. Numerical results demonstrate that our method is efficient for the convex quadratic second-order cone programming problems with bounded constraints.

The Semianalytical Solutions for Stiff Systems of Ordinary Differential Equations by Using Variational Iteration Method and Modified Variational Iteration Method with Comparison to Exact Solutions

Directory of Open Access Journals (Sweden)

Mehmet Tarik Atay

2013-01-01

Full Text Available The Variational Iteration Method (VIM and Modified Variational Iteration Method (MVIM are used to find solutions of systems of stiff ordinary differential equations for both linear and nonlinear problems. Some examples are given to illustrate the accuracy and effectiveness of these methods. We compare our results with exact results. In some studies related to stiff ordinary differential equations, problems were solved by Adomian Decomposition Method and VIM and Homotopy Perturbation Method. Comparisons with exact solutions reveal that the Variational Iteration Method (VIM and the Modified Variational Iteration Method (MVIM are easier to implement. In fact, these methods are promising methods for various systems of linear and nonlinear stiff ordinary differential equations. Furthermore, VIM, or in some cases MVIM, is giving exact solutions in linear cases and very satisfactory solutions when compared to exact solutions for nonlinear cases depending on the stiffness ratio of the stiff system to be solved.

Projective Synchronization of N-Dimensional Chaotic Fractional-Order Systems via Linear State Error Feedback Control

Directory of Open Access Journals (Sweden)

Baogui Xin

2012-01-01

Full Text Available Based on linear feedback control technique, a projective synchronization scheme of N-dimensional chaotic fractional-order systems is proposed, which consists of master and slave fractional-order financial systems coupled by linear state error variables. It is shown that the slave system can be projectively synchronized with the master system constructed by state transformation. Based on the stability theory of linear fractional order systems, a suitable controller for achieving synchronization is designed. The given scheme is applied to achieve projective synchronization of chaotic fractional-order financial systems. Numerical simulations are given to verify the effectiveness of the proposed projective synchronization scheme.

Description and use of LSODE, the Livermore Solver for Ordinary Differential Equations

Science.gov (United States)

Radhakrishnan, Krishnan; Hindmarsh, Alan C.

1993-01-01

LSODE, the Livermore Solver for Ordinary Differential Equations, is a package of FORTRAN subroutines designed for the numerical solution of the initial value problem for a system of ordinary differential equations. It is particularly well suited for 'stiff' differential systems, for which the backward differentiation formula method of orders 1 to 5 is provided. The code includes the Adams-Moulton method of orders 1 to 12, so it can be used for nonstiff problems as well. In addition, the user can easily switch methods to increase computational efficiency for problems that change character. For both methods a variety of corrector iteration techniques is included in the code. Also, to minimize computational work, both the step size and method order are varied dynamically. This report presents complete descriptions of the code and integration methods, including their implementation. It also provides a detailed guide to the use of the code, as well as an illustrative example problem.

Accelerating transient simulation of linear reduced order models.

Energy Technology Data Exchange (ETDEWEB)

Thornquist, Heidi K.; Mei, Ting; Keiter, Eric Richard; Bond, Brad

2011-10-01

Model order reduction (MOR) techniques have been used to facilitate the analysis of dynamical systems for many years. Although existing model reduction techniques are capable of providing huge speedups in the frequency domain analysis (i.e. AC response) of linear systems, such speedups are often not obtained when performing transient analysis on the systems, particularly when coupled with other circuit components. Reduced system size, which is the ostensible goal of MOR methods, is often insufficient to improve transient simulation speed on realistic circuit problems. It can be shown that making the correct reduced order model (ROM) implementation choices is crucial to the practical application of MOR methods. In this report we investigate methods for accelerating the simulation of circuits containing ROM blocks using the circuit simulator Xyce.

Second-harmonic generation circular dichroism spectroscopy from tripod-like chiral molecular films

International Nuclear Information System (INIS)

Wang Xiao-Ou; Chen Li-An; Chen Li-Xue; Sun Xiu-Dong; Li Jun-Qing; Li Chun-Fei

2010-01-01

The second-harmonic generation (SHG) circular dichroism in the light of reflection from chiral films of tripod-like chiral molecules is investigated. The expressions of the second-harmonic generation circular dichroism are derived from our presented three-coupled-oscillator model for the tripod-like chiral molecules. Spectral dependence of the circular dichroism of SHG from film surface composed of tripod-like chiral molecules is simulated numerically and analysed. Influence of chiral parameters on the second-harmonic generation circular dichroism spectrum in chiral films is studied. The result shows that the second-harmonic generation circular dichroism is a sensitive method of detecting chirality compared with the ordinary circular dichroism in linear optics. All of our work indicates that the classical molecular models are very effective to explain the second-harmonic generation circular dichroism of chiral molecular system. The classical molecular model theory can give us a clear physical picture and brings us very instructive information about the link between the molecular configuration and the nonlinear processes

Comparison of second-order impact line shifts

International Nuclear Information System (INIS)

Griem, H.R.; Iglesias, C.A.; Boercker, D.B.

1991-01-01

The second-order impact shifts in hydrogen obtained from the Baranger formalism are compared with those from a kinetic theory approach. The resulting Δn=0 contributions to the shift from the two theories are shown to be identical, except for the neglect of electron-electron correlations in the Baranger formalism. It is also shown that some care is required in taking the classical limit for the perturbing electrons, or else the shift from Δn=0 interactions vanishes

Advanced Mathematics Online: Assessing Particularities in the Online Delivery of a Second Linear Algebra Course

Science.gov (United States)

Montiel, Mariana; Bhatti, Uzma

2010-01-01

This article presents an overview of some issues that were confronted when delivering an online second Linear Algebra course (assuming a previous Introductory Linear Algebra course) to graduate students enrolled in a Secondary Mathematics Education program. The focus is on performance in one particular aspect of the course: "change of basis" and…

Asymptotic integration of a linear fourth order differential equation of Poincaré type

Directory of Open Access Journals (Sweden)

Anibal Coronel

2015-11-01

Full Text Available This article deals with the asymptotic behavior of nonoscillatory solutions of fourth order linear differential equation where the coefficients are perturbations of constants. We define a change of variable and deduce that the new variable satisfies a third order nonlinear differential equation. We assume three hypotheses. The first hypothesis is related to the constant coefficients and set up that the characteristic polynomial associated with the fourth order linear equation has simple and real roots. The other two hypotheses are related to the behavior of the perturbation functions and establish asymptotic integral smallness conditions of the perturbations. Under these general hypotheses, we obtain four main results. The first two results are related to the application of a fixed point argument to prove that the nonlinear third order equation has a unique solution. The next result concerns with the asymptotic behavior of the solutions of the nonlinear third order equation. The fourth main theorem is introduced to establish the existence of a fundamental system of solutions and to precise the formulas for the asymptotic behavior of the linear fourth order differential equation. In addition, we present an example to show that the results introduced in this paper can be applied in situations where the assumptions of some classical theorems are not satisfied.

Second Order Sliding Mode Control of the Coupled Tanks System

Directory of Open Access Journals (Sweden)

Fayiz Abu Khadra

2015-01-01

Full Text Available Four classes of second order sliding mode controllers (2-SMC have been successfully applied to regulate the liquid level in the second tank of a coupled tanks system. The robustness of these classes of 2-SMC is investigated and their performances are compared with a first order controller to show the merits of these controllers. The effectiveness of these controllers is verified through computer simulations. Comparison between the controllers is based on the time domain performance measures such as rise time, settling time, and the integral absolute error. Results showed that controllers are able to regulate the liquid level with small differences in their performance.

SECOND-ORDER CYBERNETICS, SEMIOTICS AND THE ART

Directory of Open Access Journals (Sweden)

Niculae V. Mihaita

2011-04-01

Full Text Available We take into consideration the concept of second order cybernetics and Pierce‘s approach of semiotics fundamentals. I am also an observer, experimenter and mental interpreter of metasigns given to the audience by Eugene Ionesco‘s absurd theatre. The interpreting of signs meaning is determinate by the context. From Semiotics ‗point of view, the objects I‘m studying (The Love Poem Lucifer or Evening Star, the short play Foursome and the most known, The Chairs gives me a lot of information about differences or NOT between actors, positive and negative interactions and become knowledge when I see them as signs. Second order cybernetics brings to the semiotics the idea of closure of structural coupling, interpretation and language [Soren, Cybersemiotics, 2008]. Them, the objects chosen are, for EXPERIMENTER, the YOYO in figure 1, and signifies the OBJECT of recursion. Boje [Boje, David, 2005] redefines antenarrative communication more holistically as an enactive phenomenon, and makes connections between varieties of disciplines in order to find out how antenarratives help us understand communication in the world. Instead of the finite event of producing an artifact, betting is a process and an end in itself, through which the practitioners might gain self-awareness. By synthesizing enactive-thinking in virtual space and the practice of communicating we appeal for valuable insights into the creative mind, challenging scholars and practitioners alike. Drawing contributions as above ideograms are useful for practicing cyberneticians, statisticians, researchers and academics, Informational Statistics applications [Mihaita, 2010] explores the ways in which liberal arts writers seek to involve, create and engage with new and diverse audiences from beginners encountering and participating in the work unexpectedly, to professionals from other disciplines and members of particular communities. Taking into consideration the Second-order Cybernetics�

A Solution to the Fundamental Linear Fractional Order Differential Equation

Science.gov (United States)

Hartley, Tom T.; Lorenzo, Carl F.

1998-01-01

This paper provides a solution to the fundamental linear fractional order differential equation, namely, (sub c)d(sup q, sub t) + ax(t) = bu(t). The impulse response solution is shown to be a series, named the F-function, which generalizes the normal exponential function. The F-function provides the basis for a qth order "fractional pole". Complex plane behavior is elucidated and a simple example, the inductor terminated semi- infinite lossy line, is used to demonstrate the theory.

From Ordinary Differential Equations to Structural Causal Models: the deterministic case

NARCIS (Netherlands)

Mooij, J.M.; Janzing, D.; Schölkopf, B.; Nicholson, A.; Smyth, P.

2013-01-01

We show how, and under which conditions, the equilibrium states of a first-order Ordinary Differential Equation (ODE) system can be described with a deterministic Structural Causal Model (SCM). Our exposition sheds more light on the concept of causality as expressed within the framework of

A Separation Algorithm for Sources with Temporal Structure Only Using Second-order Statistics

Directory of Open Access Journals (Sweden)

J.G. Wang

2013-09-01

Full Text Available Unlike conventional blind source separation (BSS deals with independent identically distributed (i.i.d. sources, this paper addresses the separation from mixtures of sources with temporal structure, such as linear autocorrelations. Many sequential extraction algorithms have been reported, resulting in inevitable cumulated errors introduced by the deflation scheme. We propose a robust separation algorithm to recover original sources simultaneously, through a joint diagonalizer of several average delayed covariance matrices at positions of the optimal time delay and its integers. The proposed algorithm is computationally simple and efficient, since it is based on the second-order statistics only. Extensive simulation results confirm the validity and high performance of the algorithm. Compared with related extraction algorithms, its separation signal-to-noise rate for a desired source can reach 20dB higher, and it seems rather insensitive to the estimation error of the time delay.

New Ideas: Ordinary is Extraordinary

Directory of Open Access Journals (Sweden)

Antonio Jose

2004-06-01

Full Text Available Abstract With the initial issue of this journal, a new challenge has been offered tothe world of sports nutrition: initiate "team oriented" research and clinical trials in order to make dynamic progress in terms of understandingand applying nutrition principals to the field of competitive sports. It is our further challenge that these teams think "outside the box" in terms of their approach to elucidating new concepts through which nutritional interventions might play a role in the regulation of muscle growth and repair, athletic performance and endurance, and mental acuity. What was once thought of as extraordinary might now be approached as ordinary, if the correct composition of "teams" were formed.

Stanley Cavell, Classical Hollywood and the Constitution of the Ordinary (With Notes on Billy Wilder

Directory of Open Access Journals (Sweden)

Tatjana Jukić

2016-04-01

Full Text Available When in his Tanner lectures Stanley Cavell sets out to define Ordinary Language Philosophy or – rather – to explain how it demarcates philosophy as such, he takes up psychoanalytic literary criticism in order to articulate the terms of this task. Yet the constitution of the ordinary, in Cavell, is never quite accessed from within psychoanalysis-cum-literature alone; instead, it takes another relation, that of psychoanalysis and literature to classical Hollywood, for Cavell to address the ordinary in terms of its constitution. I propose to discuss this complex using two films by Billy Wilder as a passageway to Cavell’s analytic procedure.

Synthesis of models for order-sorted first-order theories using linear algebra and constraint solving

Directory of Open Access Journals (Sweden)

Salvador Lucas

2015-12-01

Full Text Available Recent developments in termination analysis for declarative programs emphasize the use of appropriate models for the logical theory representing the program at stake as a generic approach to prove termination of declarative programs. In this setting, Order-Sorted First-Order Logic provides a powerful framework to represent declarative programs. It also provides a target logic to obtain models for other logics via transformations. We investigate the automatic generation of numerical models for order-sorted first-order logics and its use in program analysis, in particular in termination analysis of declarative programs. We use convex domains to give domains to the different sorts of an order-sorted signature; we interpret the ranked symbols of sorted signatures by means of appropriately adapted convex matrix interpretations. Such numerical interpretations permit the use of existing algorithms and tools from linear algebra and arithmetic constraint solving to synthesize the models.

Pap-smear Classification Using Efficient Second Order Neural Network Training Algorithms

DEFF Research Database (Denmark)

Ampazis, Nikolaos; Dounias, George; Jantzen, Jan

2004-01-01

In this paper we make use of two highly efficient second order neural network training algorithms, namely the LMAM (Levenberg-Marquardt with Adaptive Momentum) and OLMAM (Optimized Levenberg-Marquardt with Adaptive Momentum), for the construction of an efficient pap-smear test classifier. The alg......In this paper we make use of two highly efficient second order neural network training algorithms, namely the LMAM (Levenberg-Marquardt with Adaptive Momentum) and OLMAM (Optimized Levenberg-Marquardt with Adaptive Momentum), for the construction of an efficient pap-smear test classifier...

Existence of positive solutions for systems of second order multi ...

Indian Academy of Sciences (India)

Multi-point boundary value problems (BVPs) for ordinary differential or difference ... of a curved beam having a constant or varying cross section, three layer ..... Applying Theorem 2.5 to (3.1) and (3.2), we obtain that T has a fixed point (u, v).

3D Laser Scanning Assisted by Ordinary Plane Mirror for Non-direct Viewing Area

Directory of Open Access Journals (Sweden)

ZHANG Fan

2017-12-01

Full Text Available Terrestrial 3D laser scanning is one of principal methods to get the geometric information of object surface,and the integrity of the scanned object is a basic requirement in data acquisition. In order to solve the missing point cloud problem due to the scanning dead angle caused by confined working space,this paper proposes a method using ordinary plane mirror to obtain laser scanning data for non-direct viewing area according to the plane mirror reflection principle,analyzes the influence mechanism of the ordinary plane mirror on the propagation path and distance of laser beam,deduces the coordinate equation of the object point corresponding to the image point reflected by ordinary plane mirror in laser scanning. Given the laser scanning characteristic,this paper introduces a mirror reflection system included target balls and ordinary plane mirror,and expounds the system construction,system calibration and constructing method of system coordinate system. The feasibility and precision of the method are verified by experiments.

Learning to fear a second-order stimulus following vicarious learning.

Science.gov (United States)

Reynolds, Gemma; Field, Andy P; Askew, Chris

2017-04-01

Vicarious fear learning refers to the acquisition of fear via observation of the fearful responses of others. The present study aims to extend current knowledge by exploring whether second-order vicarious fear learning can be demonstrated in children. That is, whether vicariously learnt fear responses for one stimulus can be elicited in a second stimulus associated with that initial stimulus. Results demonstrated that children's (5-11 years) fear responses for marsupials and caterpillars increased when they were seen with fearful faces compared to no faces. Additionally, the results indicated a second-order effect in which fear-related learning occurred for other animals seen together with the fear-paired animal, even though the animals were never observed with fearful faces themselves. Overall, the findings indicate that for children in this age group vicariously learnt fear-related responses for one stimulus can subsequently be observed for a second stimulus without it being experienced in a fear-related vicarious learning event. These findings may help to explain why some individuals do not recall involvement of a traumatic learning episode in the development of their fear of a specific stimulus.

Second Order Impulsive Retarded Differential Inclusions with Nonlocal Conditions

Directory of Open Access Journals (Sweden)

Hernán R. Henríquez

2014-01-01

Full Text Available In this work we establish some existence results for abstract second order Cauchy problems modeled by a retarded differential inclusion involving nonlocal and impulsive conditions. Our results are obtained by using fixed point theory for the measure of noncompactness.

Second-home electricity consumption

Energy Technology Data Exchange (ETDEWEB)

Andersen, Frits M.; Morthorst, Poul Erik [Risoe, Systems Analysis Department, Technical University of Denmark, P.O. Box 49, DK-4000 Roskilde (Denmark); Christensen, Morten S.; Kofoed, Niels-Ulrik [Esbensen Consulting Engineers, Carl Jacobsens Vej 25D, DK-2500 Valby (Denmark); Jensen, Ole Michael [Danish Building Research Institute, Department of Energy and Environment, University of Aalborg, P.O. Box 119, DK-2970 Horsholm (Denmark)

2008-01-15

In Denmark, electricity consumption in first and second homes has developed quite differently. Since 1990, electricity consumption in ordinary residences has grown moderately, while consumption in weekend and second homes has increased considerably. In turn, this development has been blamed on a growing number of luxury cottages, new legislation permitting senior citizens to have their permanent address in their second home and a growing number of electric appliances. In order to examine the growing electricity consumption in second homes and to estimate future demand, a multidisciplinary study combining top-down and bottom-up analyses was conducted, i.e., combining models using aggregated economic parameters and feasibility studies using technical parameters, respectively. The top-down estimation showed that changes in electricity consumption in second homes correlate to changes in income. The bottom-up estimation showed that consumption was mainly affected by the frequency with which second homes were used in the winter time. This indicates that additional second homes, increased full-time use and intensified use of electric appliances are the main reasons for the observed increases in electricity consumption. Luxury tourism use and senior citizens' that use a few per cent of the second homes as their home contribute to a minor degree to the overall increase of electricity consumption. Scenarios show that this development may accelerate with increased leisure time, increased use and more permanent settlement in second homes. (author)

Optimal Control Strategies in a Two Dimensional Differential Game Using Linear Equation under a Perturbed System

Directory of Open Access Journals (Sweden)

Musa Danjuma SHEHU

2008-06-01

Full Text Available This paper lays emphasis on formulation of two dimensional differential games via optimal control theory and consideration of control systems whose dynamics is described by a system of Ordinary Differential equation in the form of linear equation under the influence of two controls U(. and V(.. Base on this, strategies were constructed. Hence we determine the optimal strategy for a control say U(. under a perturbation generated by the second control V(. within a given manifold M.

Second-order contributions to relativistic time delay in the parametrized post-Newtonian formalism

International Nuclear Information System (INIS)

Richter, G.W.; Matzner, R.A.

1983-01-01

Using a parametrized expansion of the solar metric to second order in the Newtonian potential, we calculate the relativistic delay in the round-trip travel time of a radar signal reflected from a nearby planet. We find that one second-order contribution to the delay is on the order of ten nanoseconds, which is comparable to the uncertainties in present-day experiments involving the Viking spacecraft

Adaptive Second-Order Total Variation: An Approach Aware of Slope Discontinuities

KAUST Repository

Lenzen, Frank; Becker, Florian; Lellmann, Jan

2013-01-01

Total variation (TV) regularization, originally introduced by Rudin, Osher and Fatemi in the context of image denoising, has become widely used in the field of inverse problems. Two major directions of modifications of the original approach were proposed later on. The first concerns adaptive variants of TV regularization, the second focuses on higher-order TV models. In the present paper, we combine the ideas of both directions by proposing adaptive second-order TV models, including one anisotropic model. Experiments demonstrate that introducing adaptivity results in an improvement of the reconstruction error. © 2013 Springer-Verlag.

Efficient method for time-domain simulation of the linear feedback systems containing fractional order controllers.

Science.gov (United States)

Merrikh-Bayat, Farshad

2011-04-01

One main approach for time-domain simulation of the linear output-feedback systems containing fractional-order controllers is to approximate the transfer function of the controller with an integer-order transfer function and then perform the simulation. In general, this approach suffers from two main disadvantages: first, the internal stability of the resulting feedback system is not guaranteed, and second, the amount of error caused by this approximation is not exactly known. The aim of this paper is to propose an efficient method for time-domain simulation of such systems without facing the above mentioned drawbacks. For this purpose, the fractional-order controller is approximated with an integer-order transfer function (possibly in combination with the delay term) such that the internal stability of the closed-loop system is guaranteed, and then the simulation is performed. It is also shown that the resulting approximate controller can effectively be realized by using the proposed method. Some formulas for estimating and correcting the simulation error, when the feedback system under consideration is subjected to the unit step command or the unit step disturbance, are also presented. Finally, three numerical examples are studied and the results are compared with the Oustaloup continuous approximation method. Copyright © 2011 ISA. Published by Elsevier Ltd. All rights reserved.

A Solution Method for Linear and Geometrically Nonlinear MDOF Systems with Random Properties subject to Random Excitation

DEFF Research Database (Denmark)

Micaletti, R. C.; Cakmak, A. S.; Nielsen, Søren R. K.

structural properties. The resulting state-space formulation is a system of ordinary stochastic differential equations with random coefficient and deterministic initial conditions which are subsequently transformed into ordinary stochastic differential equations with deterministic coefficients and random......A method for computing the lower-order moments of randomly-excited multi-degree-of-freedom (MDOF) systems with random structural properties is proposed. The method is grounded in the techniques of stochastic calculus, utilizing a Markov diffusion process to model the structural system with random...... initial conditions. This transformation facilitates the derivation of differential equations which govern the evolution of the unconditional statistical moments of response. Primary consideration is given to linear systems and systems with odd polynomial nonlinearities, for in these cases...

Nonlinear second order evolution inclusions with noncoercive viscosity term

Science.gov (United States)

Papageorgiou, Nikolaos S.; Rădulescu, Vicenţiu D.; Repovš, Dušan D.

2018-04-01

In this paper we deal with a second order nonlinear evolution inclusion, with a nonmonotone, noncoercive viscosity term. Using a parabolic regularization (approximation) of the problem and a priori bounds that permit passing to the limit, we prove that the problem has a solution.

An optimal PID controller via LQR for standard second order plus time delay systems.

Science.gov (United States)

Srivastava, Saurabh; Misra, Anuraag; Thakur, S K; Pandit, V S

2016-01-01

An improved tuning methodology of PID controller for standard second order plus time delay systems (SOPTD) is developed using the approach of Linear Quadratic Regulator (LQR) and pole placement technique to obtain the desired performance measures. The pole placement method together with LQR is ingeniously used for SOPTD systems where the time delay part is handled in the controller output equation instead of characteristic equation. The effectiveness of the proposed methodology has been demonstrated via simulation of stable open loop oscillatory, over damped, critical damped and unstable open loop systems. Results show improved closed loop time response over the existing LQR based PI/PID tuning methods with less control effort. The effect of non-dominant pole on the stability and robustness of the controller has also been discussed. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Real-effectiveness medicine--pursuing the best effectiveness in the ordinary care of patients.

Science.gov (United States)

Malmivaara, Antti

2013-03-01

Clinical know-how and skills as well as up-to-date scientific evidence are cornerstones for providing effective treatment for patients. However, in order to improve the effectiveness of treatment in ordinary practice, also appropriate documentation of care at the health care units and benchmarking based on this documentation are needed. This article presents the new concept of real-effectiveness medicine (REM) which pursues the best effectiveness of patient care in the real-world setting. In order to reach the goal, four layers of information are utilized: 1) good medical know-how and skills combined with the patient view, 2) up-to-date scientific evidence, 3) continuous documentation of performance in ordinary settings, and 4) benchmarking between providers. The new framework is suggested for clinicians, organizations, policy-makers, and researchers.

The second workshop on phase separation with ordering

International Nuclear Information System (INIS)

Osamura, K.; Furusaka, M.

1993-04-01

The second workshop on phase separation with ordering was held at the seminar room of Booster, National Laboratory for High Energy Physics, KEK, Tsukuba, in March 16-17, 1992. 31 participants attended this meeting. The structure and its dynamical change were discussed mainly in the experimental viewpoint, and the theories have been developed and the results of simulation were reported. (J.P.N.) 115 refs

Comparison of Second-Order Loads on a Tension-Leg Platform for Wind Turbines: Preprint

Energy Technology Data Exchange (ETDEWEB)

Gueydon, S.; Wuillaume, P.; Jonkman, J.; Robertson, A.; Platt, A.

2015-03-01

The first objective of this work is to compare the two floating offshore wind turbine simulation packages {DIFFRAC+aNySIM} and {WAMIT+FAST}. The focus is on second-order wave loads, and so first- and second-order wave loads are applied to a structure sequentially for a detailed comparison and a more precise analysis of the effects of the second-order loads. aNySIM does not have the capability to model flexible bodies, and so the simulations performed in this tool are done assuming a rigid body. FAST also assumes that the platform is rigid, but can account for the flexibility of the tower. The second objective is to study the effects of the second-order loads on the response of a TLP floating wind turbine. The flexibility of the tower must be considered for this investigation, and therefore only FAST is used.

Dynamic Uncertainty for Compensated Second-Order Systems

Directory of Open Access Journals (Sweden)

Clemens Elster

2010-08-01

Full Text Available The compensation of LTI systems and the evaluation of the according uncertainty is of growing interest in metrology. Uncertainty evaluation in metrology ought to follow specific guidelines, and recently two corresponding uncertainty evaluation schemes have been proposed for FIR and IIR filtering. We employ these schemes to compare an FIR and an IIR approach for compensating a second-order LTI system which has relevance in metrology. Our results suggest that the FIR approach is superior in the sense that it yields significantly smaller uncertainties when real-time evaluation of uncertainties is desired.

A probabilistic approach to second order variational inequalities with ...

Indian Academy of Sciences (India)

R. Narasimhan (Krishtel eMaging) 1461 1996 Oct 15 13:05:22

. TX 78712, USA. ‡CMI, Université de Provence, 39, Rue F. J. Curie, 13 453 Marseille, France. Email: mkg@math.iisc.ernet.in; mrinal@ece.utexas.edu. MS received 5 April 2002; revised 8 May 2003. Abstract. We study a class of second order ...

Effect of Second-Order and Fully Nonlinear Wave Kinematics on a Tension-Leg-Platform Wind Turbine in Extreme Wave Conditions

Energy Technology Data Exchange (ETDEWEB)

Robertson, Amy N [National Renewable Energy Laboratory (NREL), Golden, CO (United States); Jonkman, Jason [National Renewable Energy Laboratory (NREL), Golden, CO (United States); Pegalajar-Jurado, Antonio [Technical University of Denmark; Borg, Michael [Technical University of Denmark; Bredmose, Henrik [Technical University of Denmark

2017-06-03

In this study, we assess the impact of different wave kinematics models on the dynamic response of a tension-leg-platform wind turbine. Aero-hydro-elastic simulations of the floating wind turbine are carried out employing linear, second-order, and fully nonlinear kinematics using the Morison equation for the hydrodynamic forcing. The wave kinematics are computed from either theoretical or measured signals of free-surface elevation. The numerical results from each model are compared to results from wave basin tests on a scaled prototype. The comparison shows that sub and superharmonic responses can be introduced by second-order and fully nonlinear wave kinematics. The response at the wave frequency range is better reproduced when kinematics are generated from the measured surface elevation. In the future, the numerical response may be further improved by replacing the global, constant damping coefficients in the model by a more detailed, customizable definition of the user-defined numerical damping.

Multireference second order perturbation theory with a simplified treatment of dynamical correlation.

Science.gov (United States)

Xu, Enhua; Zhao, Dongbo; Li, Shuhua

2015-10-13

A multireference second order perturbation theory based on a complete active space configuration interaction (CASCI) function or density matrix renormalized group (DMRG) function has been proposed. This method may be considered as an approximation to the CAS/A approach with the same reference, in which the dynamical correlation is simplified with blocked correlated second order perturbation theory based on the generalized valence bond (GVB) reference (GVB-BCPT2). This method, denoted as CASCI-BCPT2/GVB or DMRG-BCPT2/GVB, is size consistent and has a similar computational cost as the conventional second order perturbation theory (MP2). We have applied it to investigate a number of problems of chemical interest. These problems include bond-breaking potential energy surfaces in four molecules, the spectroscopic constants of six diatomic molecules, the reaction barrier for the automerization of cyclobutadiene, and the energy difference between the monocyclic and bicyclic forms of 2,6-pyridyne. Our test applications demonstrate that CASCI-BCPT2/GVB can provide comparable results with CASPT2 (second order perturbation theory based on the complete active space self-consistent-field wave function) for systems under study. Furthermore, the DMRG-BCPT2/GVB method is applicable to treat strongly correlated systems with large active spaces, which are beyond the capability of CASPT2.

Second-order temporal interference of two independent light beams at an asymmetrical beam splitter

International Nuclear Information System (INIS)