On the nucleon-nucleon potential obtained from non-linear coupling

International Nuclear Information System (INIS)

El Ghabaty, S.S.

1975-07-01

The static limit of a pseudoscalar symmetric meson theory of nuclear forces is examined. The Born-Oppenheimer potential is determined for the case of two very heavy nucleons exchanging pseudoscalar isovector pions with non-linear coupling. It is found that the non-linear terms induced by the γ 5 coupling are cancelled by the additional pion-nucleon coupling of the non-linear sigma model. The nucleon-nucleon potential thus obtained is the same as the Yukava potential except for strength at different separations between the two nucleons

Equations for the non linear evolution of the resistive tearing modes in toroidal plasmas

International Nuclear Information System (INIS)

Edery, D.; Pellat, R.; Soule, J.L.

1979-09-01

Following the tokamak ordering, we simplify the resistive MHD equations in toroidal geometry. We obtain a closed system of non linear equations for two scalar potentials of the magnetic and velocity fields and for plasma density and temperature. If we expand these equations in the inverse of aspect ratio they are exact to the two first orders. Our formalism should correctly describe the mode coupling by curvature effects /1/ and the toroidal displacement of magnetic surfaces /2/. It provides a natural extension of the well known cylindrical model /3/ and is now being solved on computer

Infinite sets of conservation laws for linear and non-linear field equations

International Nuclear Information System (INIS)

Niederle, J.

1984-01-01

The work was motivated by a desire to understand group theoretically the existence of an infinite set of conservation laws for non-interacting fields and to carry over these conservation laws to the case of interacting fields. The relation between an infinite set of conservation laws of a linear field equation and the enveloping algebra of its space-time symmetry group was established. It is shown that in the case of the Korteweg-de Vries (KdV) equation to each symmetry of the corresponding linear equation delta sub(o)uxxx=u sub() determined by an element of the enveloping algebra of the space translation algebra, there corresponds a symmetry of the full KdV equation

Dynamical symmetries of semi-linear Schrodinger and diffusion equations

International Nuclear Information System (INIS)

Stoimenov, Stoimen; Henkel, Malte

2005-01-01

Conditional and Lie symmetries of semi-linear 1D Schrodinger and diffusion equations are studied if the mass (or the diffusion constant) is considered as an additional variable. In this way, dynamical symmetries of semi-linear Schrodinger equations become related to the parabolic and almost-parabolic subalgebras of a three-dimensional conformal Lie algebra (conf 3 ) C . We consider non-hermitian representations and also include a dimensionful coupling constant of the non-linearity. The corresponding representations of the parabolic and almost-parabolic subalgebras of (conf 3 ) C are classified and the complete list of conditionally invariant semi-linear Schrodinger equations is obtained. Possible applications to the dynamical scaling behaviour of phase-ordering kinetics are discussed

Effects of toroidal coupling on the non-linear evolution of tearing modes and on the stochastisation of the magnetic field topology in plasmas

International Nuclear Information System (INIS)

Edery, D.; Pellat, R.; Soule, J.L.

1981-01-01

The resistive MHD equations have been handled in toroidal geometry following the tokamak ordering, in order to obtain a simplified set of non-linear equations. This system of equations is compact, closed, consistent and exact to the first two orders in the expansion in the inverse aspect ratio. Studies of the non-linear evolution of tearing modes in the real geometry of tokamak discharges are now in progress, and quite significant results have been obtained from the numerical code REVE of Fontenay based on our above model. From the analytical results, strong linear coupling between neighbouring modes is expected as is demonstrated by the numerical results in the linear, and non-linear regimes. Moreover, coupling exhibits a stochastic structure of the magnetic field lines, the threshold of which is seen to be easily computed by a simple analytical criterion. (orig.)

Efficient solution of the non-linear Reynolds equation for compressible fluid using the finite element method

DEFF Research Database (Denmark)

Larsen, Jon Steffen; Santos, Ilmar

2015-01-01

An efficient finite element scheme for solving the non-linear Reynolds equation for compressible fluid coupled to compliant structures is presented. The method is general and fast and can be used in the analysis of airfoil bearings with simplified or complex foil structure models. To illustrate...

Variational iteration method for solving coupled-KdV equations

International Nuclear Information System (INIS)

Assas, Laila M.B.

2008-01-01

In this paper, the He's variational iteration method is applied to solve the non-linear coupled-KdV equations. This method is based on the use of Lagrange multipliers for identification of optimal value of a parameter in a functional. This technique provides a sequence of functions which converge to the exact solution of the coupled-KdV equations. This procedure is a powerful tool for solving coupled-KdV equations

Quantum osp-invariant non-linear Schroedinger equation

International Nuclear Information System (INIS)

Kulish, P.P.

1985-04-01

The generalizations of the non-linear Schroedinger equation (NS) associated with the orthosymplectic superalgebras are formulated. The simplest osp(1/2)-NS model is solved by the quantum inverse scattering method on a finite interval under periodic boundary conditions as well as on the wholeline in the case of a finite number of excitations. (author)

Linearized holographic isotropization at finite coupling

Energy Technology Data Exchange (ETDEWEB)

Atashi, Mahdi; Fadafan, Kazem Bitaghsir [Shahrood University of Technology, Physics Department (Iran, Islamic Republic of); Jafari, Ghadir [Institute for Research in Fundamental Sciences (IPM), School of Physics, Tehran (Iran, Islamic Republic of)

2017-06-15

We study holographic isotropization of an anisotropic homogeneous non-Abelian strongly coupled plasma in the presence of Gauss-Bonnet corrections. It was verified before that one can linearize Einstein's equations around the final black hole background and simplify the complicated setup. Using this approach, we study the expectation value of the boundary stress tensor. Although we consider small values of the Gauss-Bonnet coupling constant, it is found that finite coupling leads to significant increasing of the thermalization time. By including higher order corrections in linearization, we extend the results to study the effect of the Gauss-Bonnet coupling on the entropy production on the event horizon. (orig.)

Oscillation and non-oscillation criterion for Riemann–Weber type half-linear differential equations

Directory of Open Access Journals (Sweden)

Petr Hasil

2016-08-01

Full Text Available By the combination of the modified half-linear Prüfer method and the Riccati technique, we study oscillatory properties of half-linear differential equations. Taking into account the transformation theory of half-linear equations and using some known results, we show that the analysed equations in the Riemann–Weber form with perturbations in both terms are conditionally oscillatory. Within the process, we identify the critical oscillation values of their coefficients and, consequently, we decide when the considered equations are oscillatory and when they are non-oscillatory. As a direct corollary of our main result, we solve the so-called critical case for a certain type of half-linear non-perturbed equations.

Dissipative behavior of some fully non-linear KdV-type equations

Science.gov (United States)

Brenier, Yann; Levy, Doron

2000-03-01

The KdV equation can be considered as a special case of the general equation u t+f(u) x-δg(u xx) x=0, δ>0, where f is non-linear and g is linear, namely f( u)= u2/2 and g( v)= v. As the parameter δ tends to 0, the dispersive behavior of the KdV equation has been throughly investigated (see, e.g., [P.G. Drazin, Solitons, London Math. Soc. Lect. Note Ser. 85, Cambridge University Press, Cambridge, 1983; P.D. Lax, C.D. Levermore, The small dispersion limit of the Korteweg-de Vries equation, III, Commun. Pure Appl. Math. 36 (1983) 809-829; G.B. Whitham, Linear and Nonlinear Waves, Wiley/Interscience, New York, 1974] and the references therein). We show through numerical evidence that a completely different, dissipative behavior occurs when g is non-linear, namely when g is an even concave function such as g( v)=-∣ v∣ or g( v)=- v2. In particular, our numerical results hint that as δ→0 the solutions strongly converge to the unique entropy solution of the formal limit equation, in total contrast with the solutions of the KdV equation.

Regression of non-linear coupling of noise in LIGO detectors

Science.gov (United States)

Da Silva Costa, C. F.; Billman, C.; Effler, A.; Klimenko, S.; Cheng, H.-P.

2018-03-01

In 2015, after their upgrade, the advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) detectors started acquiring data. The effort to improve their sensitivity has never stopped since then. The goal to achieve design sensitivity is challenging. Environmental and instrumental noise couple to the detector output with different, linear and non-linear, coupling mechanisms. The noise regression method we use is based on the Wiener–Kolmogorov filter, which uses witness channels to make noise predictions. We present here how this method helped to determine complex non-linear noise couplings in the output mode cleaner and in the mirror suspension system of the LIGO detector.

Fourier imaging of non-linear structure formation

Energy Technology Data Exchange (ETDEWEB)

Brandbyge, Jacob; Hannestad, Steen, E-mail: jacobb@phys.au.dk, E-mail: sth@phys.au.dk [Department of Physics and Astronomy, University of Aarhus, Ny Munkegade 120, DK-8000 Aarhus C (Denmark)

2017-04-01

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in ΛCDM models. From N -body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important, and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.

Fourier imaging of non-linear structure formation

International Nuclear Information System (INIS)

Brandbyge, Jacob; Hannestad, Steen

2017-01-01

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in ΛCDM models. From N -body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important, and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.

From the hypergeometric differential equation to a non-linear Schrödinger one

International Nuclear Information System (INIS)

Plastino, A.; Rocca, M.C.

2015-01-01

We show that the q-exponential function is a hypergeometric function. Accordingly, it obeys the hypergeometric differential equation. We demonstrate that this differential equation can be transformed into a non-linear Schrödinger equation (NLSE). This NLSE exhibits both similarities and differences vis-a-vis the Nobre–Rego-Monteiro–Tsallis one. - Highlights: • We show that the q-exponential is a hypergeometric function. • It thus obeys the hypergeometric differential equation (HDE). • We show that the HDE can be cast as a non-linear Schrödinger equation. • This is different from the Nobre, Rego-Monteiro, Tsallis one.

On non-linear dynamics of coupled 1+1DOF versus 1+1/2DOF Electro-Mechanical System

DEFF Research Database (Denmark)

Darula, Radoslav; Sorokin, Sergey

2014-01-01

The electro-mechanical systems (EMS) are used from nano-/micro-scale (NEMS/MEMS) up to macro-scale applications. From mathematical view point, they are modelled with the second order differential equation (or a set of equations) for mechanical system, which is nonlinearly coupled with the second...... or the first order differential equation (or a set of equations) for electrical system, depending on properties of the electrical circuit. For the sake of brevity, we assume a 1DOF mechanical system, coupled to 1 or 1/2DOF electrical system (depending whether the capacitance is, or is not considered......). In the paper, authors perform a parametric study to identify operation regimes, where the capacitance term contributes to the non-linear behaviour of the coupled system. To accomplish this task, the classical method of multiple scales is used. The parametric study allows us to assess for which applications...

Non self-similar collapses described by the non-linear Schroedinger equation

International Nuclear Information System (INIS)

Berge, L.; Pesme, D.

1992-01-01

We develop a rapid method in order to find the contraction rates of the radially symmetric collapsing solutions of the nonlinear Schroedinger equation defined for space dimensions exceeding a threshold value. We explicitly determine the asymptotic behaviour of these latter solutions by solving the non stationary linear problem relative to the nonlinear Schroedinger equation. We show that the self-similar states associated with the collapsing solutions are characterized by a spatial extent which is bounded from the top by a cut-off radius

An inhomogeneous wave equation and non-linear Diophantine approximation

DEFF Research Database (Denmark)

Beresnevich, V.; Dodson, M. M.; Kristensen, S.

2008-01-01

A non-linear Diophantine condition involving perfect squares and arising from an inhomogeneous wave equation on the torus guarantees the existence of a smooth solution. The exceptional set associated with the failure of the Diophantine condition and hence of the existence of a smooth solution...

Non-linear mixed-effects pharmacokinetic/pharmacodynamic modelling in NLME using differential equations

DEFF Research Database (Denmark)

Tornøe, Christoffer Wenzel; Agersø, Henrik; Madsen, Henrik

2004-01-01

The standard software for non-linear mixed-effect analysis of pharmacokinetic/phar-macodynamic (PK/PD) data is NONMEM while the non-linear mixed-effects package NLME is an alternative as tong as the models are fairly simple. We present the nlmeODE package which combines the ordinary differential...... equation (ODE) solver package odesolve and the non-Linear mixed effects package NLME thereby enabling the analysis of complicated systems of ODEs by non-linear mixed-effects modelling. The pharmacokinetics of the anti-asthmatic drug theophylline is used to illustrate the applicability of the nlme...

GDTM-Padé technique for the non-linear differential-difference equation

Directory of Open Access Journals (Sweden)

Lu Jun-Feng

2013-01-01

Full Text Available This paper focuses on applying the GDTM-Padé technique to solve the non-linear differential-difference equation. The bell-shaped solitary wave solution of Belov-Chaltikian lattice equation is considered. Comparison between the approximate solutions and the exact ones shows that this technique is an efficient and attractive method for solving the differential-difference equations.

A discrete homotopy perturbation method for non-linear Schrodinger equation

Directory of Open Access Journals (Sweden)

H. A. Wahab

2015-12-01

Full Text Available A general analysis is made by homotopy perturbation method while taking the advantages of the initial guess, appearance of the embedding parameter, different choices of the linear operator to the approximated solution to the non-linear Schrodinger equation. We are not dependent upon the Adomian polynomials and find the linear forms of the components without these calculations. The discretised forms of the nonlinear Schrodinger equation allow us whether to apply any numerical technique on the discritisation forms or proceed for perturbation solution of the problem. The discretised forms obtained by constructed homotopy provide the linear parts of the components of the solution series and hence a new discretised form is obtained. The general discretised form for the NLSE allows us to choose any initial guess and the solution in the closed form.

Could solitons be adiabatic invariants attached to certain non linear equations

International Nuclear Information System (INIS)

Lochak, P.

1984-01-01

Arguments are given to support the claim that solitons should be the adiabatic invariants associated to certain non linear partial differential equations; a precise mathematical form of this conjecture is then stated. As a particular case of the conjecture, the Korteweg-de Vries equation is studied. (Auth.)

Study of the linear and non-linear coupling of the LH wave to the tokamak plasmas

International Nuclear Information System (INIS)

Preynas, M.

2012-10-01

In order to achieve long pulse operation with a tokamak, additional heating and current drive systems are necessary. High frequency antennas, which deliver several megawatts of power to the plasma, are currently used in several tokamaks. Moreover, a good control of the coupling of the wave launched by the antenna to the edge plasma is required to optimize the efficiency of heating and current drive LH systems. However, non-linear effects which depend on the level of injected power in the plasma strongly damage the coupling of the LH wave at particular edge parameters (density and temperature profiles). Results presented in the manuscript deal with the study of the linear and non-linear coupling of the LH wave to the plasma. In the framework of the commissioning of the Passive Active Multijunction antenna in 2009 on the Tore Supra tokamak aiming at validating the LH system suggested for ITER, the characterisation of its coupling properties was realized from low power experiments. The experimental results, which are compared with the linear coupling code ALOHA, have validated the theoretical predictions of good coupling at edge plasma density around the cut-off density. Besides, the ponderomotive effect is clearly identified as responsible for the deterioration in the coupling of the wave, which is measured under particular edge plasma conditions. A theoretical model combining the coupling of the LH wave with the ponderomotive force is suggested to explain the experimental observations. Thus, a new full wave code (named PICCOLO-2D) was developed and results from simulations validate the working hypothesis of the contribution of the ponderomotive effect in the non-linear observations of LHCD coupling on Tore Supra. (author)

Approximate Forward Difference Equations for the Lower Order Non-Stationary Statistics of Geometrically Non-Linear Systems subject to Random Excitation

DEFF Research Database (Denmark)

Köylüoglu, H. U.; Nielsen, Søren R. K.; Cakmak, A. S.

Geometrically non-linear multi-degree-of-freedom (MDOF) systems subject to random excitation are considered. New semi-analytical approximate forward difference equations for the lower order non-stationary statistical moments of the response are derived from the stochastic differential equations...... of motion, and, the accuracy of these equations is numerically investigated. For stationary excitations, the proposed method computes the stationary statistical moments of the response from the solution of non-linear algebraic equations....

Unique solvability of a non-linear non-local boundary-value problem for systems of non-linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Dilna, N.; Rontó, András

2010-01-01

Roč. 60, č. 3 (2010), s. 327-338 ISSN 0139-9918 R&D Projects: GA ČR(CZ) GA201/06/0254 Institutional research plan: CEZ:AV0Z10190503 Keywords : non-linear boundary value-problem * functional differential equation * non-local condition * unique solvability * differential inequality Subject RIV: BA - General Mathematics Impact factor: 0.316, year: 2010 http://link.springer.com/article/10.2478%2Fs12175-010-0015-9

Classical electromagnetic non-minimal coupling for spin 3/2 fields

International Nuclear Information System (INIS)

Villanueva, V. M.; Obregon, O.; Nieto, J. A.

1996-01-01

We obtain a non-minimal electromagnetic coupling for spin 3/2 particles from linearized N=2 Supergravity. This coupling coincides with the one found by Ferrara et al. by demanding g=2 at the tree level. Linearized Einstein field equations plus interaction terms are obtained by squaring the Rarita-Schwinger with this non-minimal coupling by using generalized Poisson brackets

Anti-symmetrically fused model and non-linear integral equations in the three-state Uimin-Sutherland model

International Nuclear Information System (INIS)

Fujii, Akira; Kluemper, Andreas

1999-01-01

We derive the non-linear integral equations determining the free energy of the three-state pure bosonic Uimin-Sutherland model. In order to find a complete set of auxiliary functions, the anti-symmetric fusion procedure is utilized. We solve the non-linear integral equations numerically and see that the low-temperature behavior coincides with that predicted by conformal field theory. The magnetization and magnetic susceptibility are also calculated by means of the non-linear integral equation

Non-linear analysis of wave progagation using transform methods and plates and shells using integral equations

Science.gov (United States)

Pipkins, Daniel Scott

Two diverse topics of relevance in modern computational mechanics are treated. The first involves the modeling of linear and non-linear wave propagation in flexible, lattice structures. The technique used combines the Laplace Transform with the Finite Element Method (FEM). The procedure is to transform the governing differential equations and boundary conditions into the transform domain where the FEM formulation is carried out. For linear problems, the transformed differential equations can be solved exactly, hence the method is exact. As a result, each member of the lattice structure is modeled using only one element. In the non-linear problem, the method is no longer exact. The approximation introduced is a spatial discretization of the transformed non-linear terms. The non-linear terms are represented in the transform domain by making use of the complex convolution theorem. A weak formulation of the resulting transformed non-linear equations yields a set of element level matrix equations. The trial and test functions used in the weak formulation correspond to the exact solution of the linear part of the transformed governing differential equation. Numerical results are presented for both linear and non-linear systems. The linear systems modeled are longitudinal and torsional rods and Bernoulli-Euler and Timoshenko beams. For non-linear systems, a viscoelastic rod and Von Karman type beam are modeled. The second topic is the analysis of plates and shallow shells under-going finite deflections by the Field/Boundary Element Method. Numerical results are presented for two plate problems. The first is the bifurcation problem associated with a square plate having free boundaries which is loaded by four, self equilibrating corner forces. The results are compared to two existing numerical solutions of the problem which differ substantially. non-linear model are compared to those

Fitting and forecasting coupled dark energy in the non-linear regime

Energy Technology Data Exchange (ETDEWEB)

Casas, Santiago; Amendola, Luca; Pettorino, Valeria; Vollmer, Adrian [Institut für Theoretische Physik, Ruprecht-Karls-Universität Heidelberg, Philosophenweg 16, Heidelberg, 69120 Germany (Germany); Baldi, Marco, E-mail: casas@thphys.uni-heidelberg.de, E-mail: l.amendola@thphys.uni-heidelberg.de, E-mail: mail@marcobaldi.it, E-mail: v.pettorino@thphys.uni-heidelberg.de, E-mail: vollmer@thphys.uni-heidelberg.de [Dipartimento di Fisica e Astronomia, Alma Mater Studiorum Università di Bologna, viale Berti Pichat, 6/2, Bologna, I-40127 Italy (Italy)

2016-01-01

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range 0z=–1.6 and wave modes below 0k=1 h/Mpc. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, β{sup 2}, with precision smaller than 4% and all other cosmological parameters better than 1%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications.

Fitting and forecasting coupled dark energy in the non-linear regime

International Nuclear Information System (INIS)

Casas, Santiago; Amendola, Luca; Pettorino, Valeria; Vollmer, Adrian; Baldi, Marco

2016-01-01

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range 0z=–1.6 and wave modes below 0k=1 h/Mpc. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, β 2 , with precision smaller than 4% and all other cosmological parameters better than 1%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications

Linear supermultiplets and non-holomorphic gauge coupling functions

International Nuclear Information System (INIS)

Binetruy, P.; Grimm, R.; Girardi, G.

1991-04-01

The general couplings of linear multiplets, including Chern-Simons forms, to chiral matter as well as to the standard supergravity-matter system are constructed. Insisting on a canonically normalised Einstein term in particular the appearance of non-holomorphic gauge couplings are discussed and duality transformations in full generality are performed. The implications of these structures for the effective description of sigma model anomalies are presented with and without coupling to supergravity, following recent proposals of Derendinger, Ferrara, Kounnas and Zwirner and of Cardoso and Ovrut. (author) 14 refs

Equations of motion for a (non-linear) scalar field model as derived from the field equations

International Nuclear Information System (INIS)

Kaniel, S.; Itin, Y.

2006-01-01

The problem of derivation of the equations of motion from the field equations is considered. Einstein's field equations have a specific analytical form: They are linear in the second order derivatives and quadratic in the first order derivatives of the field variables. We utilize this particular form and propose a novel algorithm for the derivation of the equations of motion from the field equations. It is based on the condition of the balance between the singular terms of the field equation. We apply the algorithm to a non-linear Lorentz invariant scalar field model. We show that it results in the Newton law of attraction between the singularities of the field moved on approximately geodesic curves. The algorithm is applicable to the N-body problem of the Lorentz invariant field equations. (Abstract Copyright [2006], Wiley Periodicals, Inc.)

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.

A New Theory of Non-Linear Thermo-Elastic Constitutive Equation of Isotropic Hyperelastic Materials

Science.gov (United States)

Li, Chen; Liao, Yufei

2018-03-01

Considering the influence of temperature and strain variables on materials. According to the relationship of conjugate stress-strain, a complete and irreducible non-linear constitutive equation of isotropic hyperelastic materials is derived and the constitutive equations of 16 types of isotropic hyperelastic materials are given we study the transformation methods and routes of 16 kinds of constitutive equations and the study proves that transformation of two forms of constitutive equation. As an example of application, the non-linear thermo-elastic constitutive equation of isotropic hyperelastic materials is combined with the natural vulcanized rubber experimental data in the existing literature base on MATLAB, The results show that the fitting accuracy is satisfactory.

Separation-induced boundary layer transition: Modeling with a non-linear eddy-viscosity model coupled with the laminar kinetic energy equation

International Nuclear Information System (INIS)

Vlahostergios, Z.; Yakinthos, K.; Goulas, A.

2009-01-01

We present an effort to model the separation-induced transition on a flat plate with a semi-circular leading edge, using a cubic non-linear eddy-viscosity model combined with the laminar kinetic energy. A non-linear model, compared to a linear one, has the advantage to resolve the anisotropic behavior of the Reynolds-stresses in the near-wall region and it provides a more accurate expression for the generation of turbulence in the transport equation of the turbulence kinetic energy. Although in its original formulation the model is not able to accurately predict the separation-induced transition, the inclusion of the laminar kinetic energy increases its accuracy. The adoption of the laminar kinetic energy by the non-linear model is presented in detail, together with some additional modifications required for the adaption of the laminar kinetic energy into the basic concepts of the non-linear eddy-viscosity model. The computational results using the proposed combined model are shown together with the ones obtained using an isotropic linear eddy-viscosity model, which adopts also the laminar kinetic energy concept and in comparison with the existing experimental data.

New exact solutions of the Tzitzéica-type equations in non-linear optics using the expa function method

Science.gov (United States)

Hosseini, K.; Ayati, Z.; Ansari, R.

2018-04-01

One specific class of non-linear evolution equations, known as the Tzitzéica-type equations, has received great attention from a group of researchers involved in non-linear science. In this article, new exact solutions of the Tzitzéica-type equations arising in non-linear optics, including the Tzitzéica, Dodd-Bullough-Mikhailov and Tzitzéica-Dodd-Bullough equations, are obtained using the expa function method. The integration technique actually suggests a useful and reliable method to extract new exact solutions of a wide range of non-linear evolution equations.

Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series

Science.gov (United States)

Gnoffo, Peter A.

2015-01-01

Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.

Calculations of stationary solutions for the non linear viscous resistive MHD equations in slab geometry

International Nuclear Information System (INIS)

Edery, D.

1983-11-01

The reduced system of the non linear resistive MHD equations is used in the 2-D one helicity approximation in the numerical computations of stationary tearing modes. The critical magnetic Raynolds number S (S=tausub(r)/tausub(H) where tausub(R) and tausub(H) are respectively the characteristic resistive and hydro magnetic times) and the corresponding linear solution are computed as a starting approximation for the full non linear equations. These equations are then treated numerically by an iterative procedure which is shown to be rapidly convergent. A numerical application is given in the last part of this paper

Neutron stars in non-linear coupling models

International Nuclear Information System (INIS)

Taurines, Andre R.; Vasconcellos, Cesar A.Z.; Malheiro, Manuel; Chiapparini, Marcelo

2001-01-01

We present a class of relativistic models for nuclear matter and neutron stars which exhibits a parameterization, through mathematical constants, of the non-linear meson-baryon couplings. For appropriate choices of the parameters, it recovers current QHD models found in the literature: Walecka, ZM and ZM3 models. We have found that the ZM3 model predicts a very small maximum neutron star mass, ∼ 0.72M s un. A strong similarity between the results of ZM-like models and those with exponential couplings is noted. Finally, we discuss the very intense scalar condensates found in the interior of neutron stars which may lead to negative effective masses. (author)

Neutron stars in non-linear coupling models

Energy Technology Data Exchange (ETDEWEB)

Taurines, Andre R.; Vasconcellos, Cesar A.Z. [Rio Grande do Sul Univ., Porto Alegre, RS (Brazil); Malheiro, Manuel [Universidade Federal Fluminense, Niteroi, RJ (Brazil); Chiapparini, Marcelo [Universidade do Estado, Rio de Janeiro, RJ (Brazil)

2001-07-01

We present a class of relativistic models for nuclear matter and neutron stars which exhibits a parameterization, through mathematical constants, of the non-linear meson-baryon couplings. For appropriate choices of the parameters, it recovers current QHD models found in the literature: Walecka, ZM and ZM3 models. We have found that the ZM3 model predicts a very small maximum neutron star mass, {approx} 0.72M{sub s}un. A strong similarity between the results of ZM-like models and those with exponential couplings is noted. Finally, we discuss the very intense scalar condensates found in the interior of neutron stars which may lead to negative effective masses. (author)

Linear differential equations to solve nonlinear mechanical problems: A novel approach

OpenAIRE

Nair, C. Radhakrishnan

2004-01-01

Often a non-linear mechanical problem is formulated as a non-linear differential equation. A new method is introduced to find out new solutions of non-linear differential equations if one of the solutions of a given non-linear differential equation is known. Using the known solution of the non-linear differential equation, linear differential equations are set up. The solutions of these linear differential equations are found using standard techniques. Then the solutions of the linear differe...

On the equivalence between particular types of Navier-Stokes and non-linear Schroedinger equations

International Nuclear Information System (INIS)

Dietrich, K.; Vautherin, D.

1985-01-01

We derive a Schroedinger equation equivalent to the Navier-Stokes equation in the special case of constant kinematic viscosities. This equation contains a non-linear term similar to that proposed by Kostin for a quantum description of friction [fr

Linear determining equations for differential constraints

International Nuclear Information System (INIS)

Kaptsov, O V

1998-01-01

A construction of differential constraints compatible with partial differential equations is considered. Certain linear determining equations with parameters are used to find such differential constraints. They generalize the classical determining equations used in the search for admissible Lie operators. As applications of this approach equations of an ideal incompressible fluid and non-linear heat equations are discussed

Improved harmonic balance approach to periodic solutions of non-linear jerk equations

International Nuclear Information System (INIS)

Wu, B.S.; Lim, C.W.; Sun, W.P.

2006-01-01

An analytical approximate approach for determining periodic solutions of non-linear jerk equations involving third-order time-derivative is presented. This approach incorporates salient features of both Newton's method and the method of harmonic balance. By appropriately imposing the method of harmonic balance to the linearized equation, the approach requires only one or two iterations to predict very accurate analytical approximate solutions for a large range of initial velocity amplitude. One typical example is used to verify and illustrate the usefulness and effectiveness of the proposed approach

A three operator split-step method covering a larger set of non-linear partial differential equations

Science.gov (United States)

Zia, Haider

2017-06-01

This paper describes an updated exponential Fourier based split-step method that can be applied to a greater class of partial differential equations than previous methods would allow. These equations arise in physics and engineering, a notable example being the generalized derivative non-linear Schrödinger equation that arises in non-linear optics with self-steepening terms. These differential equations feature terms that were previously inaccessible to model accurately with low computational resources. The new method maintains a 3rd order error even with these additional terms and models the equation in all three spatial dimensions and time. The class of non-linear differential equations that this method applies to is shown. The method is fully derived and implementation of the method in the split-step architecture is shown. This paper lays the mathematical ground work for an upcoming paper employing this method in white-light generation simulations in bulk material.

Quantum linear Boltzmann equation

International Nuclear Information System (INIS)

Vacchini, Bassano; Hornberger, Klaus

2009-01-01

We review the quantum version of the linear Boltzmann equation, which describes in a non-perturbative fashion, by means of scattering theory, how the quantum motion of a single test particle is affected by collisions with an ideal background gas. A heuristic derivation of this Lindblad master equation is presented, based on the requirement of translation-covariance and on the relation to the classical linear Boltzmann equation. After analyzing its general symmetry properties and the associated relaxation dynamics, we discuss a quantum Monte Carlo method for its numerical solution. We then review important limiting forms of the quantum linear Boltzmann equation, such as the case of quantum Brownian motion and pure collisional decoherence, as well as the application to matter wave optics. Finally, we point to the incorporation of quantum degeneracies and self-interactions in the gas by relating the equation to the dynamic structure factor of the ambient medium, and we provide an extension of the equation to include internal degrees of freedom.

An analytical study of non-linear behaviour of coupled 2+2x0.5 DOF electro-magneto-mechanical system by a method of multiple scales

DEFF Research Database (Denmark)

Darula, Radoslav; Sorokin, Sergey

2013-01-01

An electro-magneto-mechanical system combines three physical domains - a mechanical structure, a magnetic field and an electric circuit. The interaction between these domains is analysed for a structure with two degrees of freedom (translational and rotational) and two electrical circuits. Each...... electrical circuit is described by a differential equation of the 1st order, which is considered to contribute to the coupled system by 0.5 DOF. The electrical and mechanical systems are coupled via a magnetic circuit, which is inherently non-linear, due to a non-linear nature of the electro-magnetic force...

Nodal methods with non linear feedback for the three dimensional resolution of the diffusion's multigroup equations

International Nuclear Information System (INIS)

Ferri, A.A.

1986-01-01

Nodal methods applied in order to calculate the power distribution in a nuclear reactor core are presented. These methods have received special attention, because they yield accurate results in short computing times. Present nodal schemes contain several unknowns per node and per group. In the methods presented here, non linear feedback of the coupling coefficients has been applied to reduce this number to only one unknown per node and per group. The resulting algorithm is a 7- points formula, and the iterative process has proved stable in the response matrix scheme. The intranodal flux shape is determined by partial integration of the diffusion equations over two of the coordinates, leading to a set of three coupled one-dimensional equations. These can be solved by using a polynomial approximation or by integration (analytic solution). The tranverse net leakage is responsible for the coupling between the spatial directions, and two alternative methods are presented to evaluate its shape: direct parabolic approximation and local model expansion. Numerical results, which include the IAEA two-dimensional benchmark problem illustrate the efficiency of the developed methods. (M.E.L.) [es

Quantum Non-Markovian Langevin Equations and Transport Coefficients

International Nuclear Information System (INIS)

Sargsyan, V.V.; Antonenko, N.V.; Kanokov, Z.; Adamian, G.G.

2005-01-01

Quantum diffusion equations featuring explicitly time-dependent transport coefficients are derived from generalized non-Markovian Langevin equations. Generalized fluctuation-dissipation relations and analytic expressions for calculating the friction and diffusion coefficients in nuclear processes are obtained. The asymptotic behavior of the transport coefficients and correlation functions for a damped harmonic oscillator that is linearly coupled in momentum to a heat bath is studied. The coupling to a heat bath in momentum is responsible for the appearance of the diffusion coefficient in coordinate. The problem of regression of correlations in quantum dissipative systems is analyzed

Optimal linear-quadratic control of coupled parabolic-hyperbolic PDEs

Science.gov (United States)

Aksikas, I.; Moghadam, A. Alizadeh; Forbes, J. F.

2017-10-01

This paper focuses on the optimal control design for a system of coupled parabolic-hypebolic partial differential equations by using the infinite-dimensional state-space description and the corresponding operator Riccati equation. Some dynamical properties of the coupled system of interest are analysed to guarantee the existence and uniqueness of the solution of the linear-quadratic (LQ)-optimal control problem. A state LQ-feedback operator is computed by solving the operator Riccati equation, which is converted into a set of algebraic and differential Riccati equations, thanks to the eigenvalues and the eigenvectors of the parabolic operator. The results are applied to a non-isothermal packed-bed catalytic reactor. The LQ-optimal controller designed in the early portion of the paper is implemented for the original nonlinear model. Numerical simulations are performed to show the controller performances.

Coupling and reduction of the HAWC equations

DEFF Research Database (Denmark)

Nim, E.

2001-01-01

This report contains a description of a general method for coupling and reduction of the so-called HAWC equations, which constitute the basis equations of motion of the aeroelastic model HAWC used widely by research institutes and industrial companies formore than the ten years. The principal aim....... In addition, the method enables the reduction of the number of degrees of freedom of the structure in order to increase the calculation efficiency and improve thecondition of the system.......This report contains a description of a general method for coupling and reduction of the so-called HAWC equations, which constitute the basis equations of motion of the aeroelastic model HAWC used widely by research institutes and industrial companies formore than the ten years. The principal aim...... of the work has been to enable the modelling wind turbines with large displacements of the blades in order to predict phenomena caused by geometric non-linear effects. However, the method can also be applied tomodel the nacelle/shaft structure of a turbine more detailed than the present HAWC model...

Solvable linear potentials in the Dirac equation

International Nuclear Information System (INIS)

Dominguez-Adame, F.; Gonzalez, M.A.

1990-01-01

The Dirac equation for some linear potentials leading to Schroedinger-like oscillator equations for the upper and lower components of the Dirac spinor have been solved. Energy levels for the bound states appear in pairs, so that both particles and antiparticles may be bound with the same energy. For weak coupling, the spacing between levels is proportional to the coupling constant while in the strong limit those levels are depressed compared to the nonrelativistic ones

Linear representation of algebras with non-associative operations which are satisfy in the balanced functional equations

International Nuclear Information System (INIS)

Ehsani, Amir

2015-01-01

Algebras with a pair of non-associative binary operations (f, g) which are satisfy in the balanced quadratic functional equations with four object variables considered. First, we obtain a linear representation for the operations, of this kind of binary algebras (A,f,g), over an abelian group (A, +) and then we generalize the linear representation of operations, to an algebra (A,F) with non-associative binary operations which are satisfy in the balanced quadratic functional equations with four object variables. (paper)

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

Value distribution of meromorphic solutions of homogeneous and non-homogeneous complex linear differential-difference equations

Directory of Open Access Journals (Sweden)

Luo Li-Qin

2016-01-01

Full Text Available In this paper, we investigate the value distribution of meromorphic solutions of homogeneous and non-homogeneous complex linear differential-difference equations, and obtain the results on the relations between the order of the solutions and the convergence exponents of the zeros, poles, a-points and small function value points of the solutions, which show the relations in the case of non-homogeneous equations are sharper than the ones in the case of homogeneous equations.

A non linear half space problem for radiative transfer equations. Application to the Rosseland approximation

International Nuclear Information System (INIS)

Sentis, R.

1984-07-01

The radiative transfer equations may be approximated by a non linear diffusion equation (called Rosseland equation) when the mean free paths of the photons are small with respect to the size of the medium. Some technical assomptions are made, namely about the initial conditions, to avoid any problem of initial layer terms

A linear evolution for non-linear dynamics and correlations in realistic nuclei

International Nuclear Information System (INIS)

Levin, E.; Lublinsky, M.

2004-01-01

A new approach to high energy evolution based on a linear equation for QCD generating functional is developed. This approach opens a possibility for systematic study of correlations inside targets, and, in particular, inside realistic nuclei. Our results are presented as three new equations. The first one is a linear equation for QCD generating functional (and for scattering amplitude) that sums the 'fan' diagrams. For the amplitude this equation is equivalent to the non-linear Balitsky-Kovchegov equation. The second equation is a generalization of the Balitsky-Kovchegov non-linear equation to interactions with realistic nuclei. It includes a new correlation parameter which incorporates, in a model-dependent way, correlations inside the nuclei. The third equation is a non-linear equation for QCD generating functional (and for scattering amplitude) that in addition to the 'fan' diagrams sums the Glauber-Mueller multiple rescatterings

On the Cauchy problem for a Sobolev-type equation with quadratic non-linearity

International Nuclear Information System (INIS)

Aristov, Anatoly I

2011-01-01

We investigate the asymptotic behaviour as t→∞ of the solution of the Cauchy problem for a Sobolev-type equation with quadratic non-linearity and develop ideas used by I. A. Shishmarev and other authors in the study of classical and Sobolev-type equations. Conditions are found under which it is possible to consider the case of an arbitrary dimension of the spatial variable.

Computing with linear equations and matrices

International Nuclear Information System (INIS)

Churchhouse, R.F.

1983-01-01

Systems of linear equations and matrices arise in many disciplines. The equations may accurately represent conditions satisfied by a system or, more likely, provide an approximation to a more complex system of non-linear or differential equations. The system may involve a few or many thousand unknowns and each individual equation may involve few or many of them. Over the past 50 years a vast literature on methods for solving systems of linear equations and the associated problems of finding the inverse or eigenvalues of a matrix has been produced. These lectures cover those methods which have been found to be most useful for dealing with such types of problem. References are given where appropriate and attention is drawn to the possibility of improved methods for use on vector and parallel processors. (orig.)

Painlevйe analysis and integrability of two-coupled non-linear ...

Indian Academy of Sciences (India)

the Painlevйe property. In this case the system is expected to be integrable. In recent years more attention is paid to the study of coupled non-linear oscilla- ... Painlevйe analysis. To be self-contained, in §2 we briefly outline the salient features.

Analytical solutions of coupled-mode equations for microring ...

Indian Academy of Sciences (India)

equivalent to waveguide and single microring coupled system. The 3 × 3 coupled system is equivalent to waveguide and double microring coupled system. In this paper, we adopt a novel approach for obtaining coupled-mode equations for linearly distributed and circularly distributed multiwaveguide systems with different ...

Fourth order Douglas implicit scheme for solving three dimension reaction diffusion equation with non-linear source term

Science.gov (United States)

Hasnain, Shahid; Saqib, Muhammad; Mashat, Daoud Suleiman

2017-07-01

This research paper represents a numerical approximation to non-linear three dimension reaction diffusion equation with non-linear source term from population genetics. Since various initial and boundary value problems exist in three dimension reaction diffusion phenomena, which are studied numerically by different numerical methods, here we use finite difference schemes (Alternating Direction Implicit and Fourth Order Douglas Implicit) to approximate the solution. Accuracy is studied in term of L2, L∞ and relative error norms by random selected grids along time levels for comparison with analytical results. The test example demonstrates the accuracy, efficiency and versatility of the proposed schemes. Numerical results showed that Fourth Order Douglas Implicit scheme is very efficient and reliable for solving 3-D non-linear reaction diffusion equation.

Fourth order Douglas implicit scheme for solving three dimension reaction diffusion equation with non-linear source term

Directory of Open Access Journals (Sweden)

Shahid Hasnain

2017-07-01

Full Text Available This research paper represents a numerical approximation to non-linear three dimension reaction diffusion equation with non-linear source term from population genetics. Since various initial and boundary value problems exist in three dimension reaction diffusion phenomena, which are studied numerically by different numerical methods, here we use finite difference schemes (Alternating Direction Implicit and Fourth Order Douglas Implicit to approximate the solution. Accuracy is studied in term of L2, L∞ and relative error norms by random selected grids along time levels for comparison with analytical results. The test example demonstrates the accuracy, efficiency and versatility of the proposed schemes. Numerical results showed that Fourth Order Douglas Implicit scheme is very efficient and reliable for solving 3-D non-linear reaction diffusion equation.

Perfect observables for the hierarchical non-linear O(N)-invariant σ-model

International Nuclear Information System (INIS)

Wieczerkowski, C.; Xylander, Y.

1995-05-01

We compute moving eigenvalues and the eigenvectors of the linear renormalization group transformation for observables along the renormalized trajectory of the hierarchical non-linear O(N)-invariant σ-model by means of perturbation theory in the running coupling constant. Moving eigenvectors are defined as solutions to a Callan-Symanzik type equation. (orig.)

On the prolongation structure and Backlund transformation for new non-linear Klein-Gordon equations

International Nuclear Information System (INIS)

Roy Chowdhury, A.; Mukherjee, J.

1986-07-01

We have considered the complete integrability of two nonlinear equations which are some kind of extensions of usual Sine-Gordon and Sinh-Gordon equations. The first one is of non-autonomous version of Sinh-Gordon system and the second is closely related to the usual Sine-Gordon theory. The first problem indicates how (x,t) dependent non-linear equations can be treated in the prolongation theory and how a Backlund map can be constructed. The second one is a variation of the usual Sine-Gordon equation and suggests that there may be other equations (similar to Sine-Gordon) which are completely integrable. In both cases we have been able to construct the Lax pair. We then construct an auto-Backlund map by following the idea of Konno and Wadati, for the generation of multisolution states. (author)

Non-linear corrections to the time-covariance function derived from a multi-state chemical master equation.

Science.gov (United States)

Scott, M

2012-08-01

The time-covariance function captures the dynamics of biochemical fluctuations and contains important information about the underlying kinetic rate parameters. Intrinsic fluctuations in biochemical reaction networks are typically modelled using a master equation formalism. In general, the equation cannot be solved exactly and approximation methods are required. For small fluctuations close to equilibrium, a linearisation of the dynamics provides a very good description of the relaxation of the time-covariance function. As the number of molecules in the system decrease, deviations from the linear theory appear. Carrying out a systematic perturbation expansion of the master equation to capture these effects results in formidable algebra; however, symbolic mathematics packages considerably expedite the computation. The authors demonstrate that non-linear effects can reveal features of the underlying dynamics, such as reaction stoichiometry, not available in linearised theory. Furthermore, in models that exhibit noise-induced oscillations, non-linear corrections result in a shift in the base frequency along with the appearance of a secondary harmonic.

Numerical simulation of electro-osmotic consolidation coupling non-linear variation of soil parameters

Science.gov (United States)

Wu, Hui; Hu, Liming; Wen, Qingbo

2017-06-01

Electro-osmotic consolidation is an effective method for soft ground improvement. A main limitation of previous numerical models on this technique is the ignorance of the non-linear variation of soil parameters. In the present study, a multi-field numerical model is developed with the consideration of the non-linear variation of soil parameters during electro-osmotic consolidation process. The numerical simulations on an axisymmetric model indicated that the non-linear variation of soil parameters showed remarkable impact on the development of the excess pore water pressure and degree of consolidation. A field experiment with complex geometry, boundary conditions, electrode configuration and voltage application was further simulated with the developed numerical model. The comparison between field and numerical data indicated that the numerical model coupling of the non-linear variation of soil parameters gave more reasonable results. The developed numerical model is capable to analyze engineering cases with complex operating conditions.

Stimulated Raman scattering and ion dynamics: the role of Langmuir wave non-linearities

International Nuclear Information System (INIS)

Bonnaud, G.; Pesme, D.

1988-02-01

The non-linear evolution of stimulated Raman scattering by coupling of the SRS-driven Langmuir waves to ion acoustic waves is studied numerically, in a homogeneous density laser-irradiated plasma. The coupled wave amplitude behaviour is represented either by envelope equations or by complete wave-like equations. The various physical phenomena which are involved are described. This preliminary work has been presented at the 17th Anomalous Absorption Conference, held in last May, in Lake Tahoe City (USA) [fr

Non-Linear Multi-Physics Analysis and Multi-Objective Optimization in Electroheating Applications

Czech Academy of Sciences Publication Activity Database

di Barba, P.; Doležel, Ivo; Mognaschi, M. E.; Savini, A.; Karban, P.

2014-01-01

Roč. 50, č. 2 (2014), s. 7016604-7016604 ISSN 0018-9464 Institutional support: RVO:61388998 Keywords : coupled multi-physics problems * finite element method * non-linear equations Subject RIV: JA - Electronics ; Optoelectronics, Electrical Engineering Impact factor: 1.386, year: 2014

Non-linear neutron star oscillations viewed as deviations from an equilibrium state

International Nuclear Information System (INIS)

Sperhake, U

2002-01-01

A numerical technique is presented which facilitates the evolution of non-linear neutron star oscillations with a high accuracy essentially independent of the oscillation amplitude. We apply this technique to radial neutron star oscillations in a Lagrangian formulation and demonstrate the superior performance of the new scheme compared with 'conventional' techniques. The key feature of our approach is to describe the evolution in terms of deviations from an equilibrium configuration. In contrast to standard perturbation analysis we keep all higher order terms in the evolution equations and thus obtain a fully non-linear description. The advantage of our scheme lies in the elimination of background terms from the equations and the associated numerical errors. The improvements thus achieved will be particularly significant in the study of mildly non-linear effects where the amplitude of the dynamic signal is small compared with the equilibrium values but large enough to warrant non-linear effects. We apply the new technique to the study of non-linear coupling of Eigenmodes and non-linear effects in the oscillations of marginally stable neutron stars. We find non-linear effects in low amplitude oscillations to be particularly pronounced in the range of modes with vanishing frequency which typically mark the onset of instability. (author)

Synchronization of Two Non-Identical Coupled Exciters in a Non-Resonant Vibrating System of Linear Motion. Part I: Theoretical Analysis

Directory of Open Access Journals (Sweden)

Chunyu Zhao

2009-01-01

Full Text Available In this paper an analytical approach is proposed to study the feature of frequency capture of two non-identical coupled exciters in a non-resonant vibrating system. The electromagnetic torque of an induction motor in the quasi-steady-state operation is derived. With the introduction of two perturbation small parameters to average angular velocity of two exciters and their phase difference, we deduce the Equation of Frequency Capture by averaging two motion equations of two exciters over their average period. It converts the synchronization problem of two exciters into that of existence and stability of zero solution for the Equation of Frequency Capture. The conditions of implementing frequency capture and that of stabilizing synchronous operation of two motors have been derived. The concept of torque of frequency capture is proposed to physically explain the peculiarity of self-synchronization of the two exciters. An interesting conclusion is reached that the moments of inertia of the two exciters in the Equation of Frequency Capture reduce and there is a coupling moment of inertia between the two exciters. The reduction of moments of inertia and the coupling moment of inertia have an effect on the stability of synchronous operation.

Weak ωNN coupling in the non-linear chiral model

International Nuclear Information System (INIS)

Shmatikov, M.

1988-01-01

In the non-linear chiral model with the soliton solution stabilized by the ω-meson field the weak ωNN coupling constants are calculated. Applying the vector dominance model for the isoscalar current the constant of the isoscalar P-odd ωNN interaction h ω (0) =0 is obtained while the constant of the isovector (of the Lagrangian of the ωNN interaction proves to be h ω (1) ≅ 1.0x10 -7

TBA equations for the mass gap in the O(2r) non-linear σ-models

International Nuclear Information System (INIS)

Balog, Janos; Hegedues, Arpad

2005-01-01

We propose TBA integral equations for 1-particle states in the O(n) non-linear σ-model for even n. The equations are conjectured on the basis of the analytic properties of the large volume asymptotics of the problem, which is explicitly constructed starting from Luscher's asymptotic formula. For small volumes the mass gap values computed numerically from the TBA equations agree very well with results of three-loop perturbation theory calculations, providing support for the validity of the proposed TBA system

Fast and local non-linear evolution of steep wave-groups on deep water: A comparison of approximate models to fully non-linear simulations

International Nuclear Information System (INIS)

Adcock, T. A. A.; Taylor, P. H.

2016-01-01

The non-linear Schrödinger equation and its higher order extensions are routinely used for analysis of extreme ocean waves. This paper compares the evolution of individual wave-packets modelled using non-linear Schrödinger type equations with packets modelled using fully non-linear potential flow models. The modified non-linear Schrödinger Equation accurately models the relatively large scale non-linear changes to the shape of wave-groups, with a dramatic contraction of the group along the mean propagation direction and a corresponding extension of the width of the wave-crests. In addition, as extreme wave form, there is a local non-linear contraction of the wave-group around the crest which leads to a localised broadening of the wave spectrum which the bandwidth limited non-linear Schrödinger Equations struggle to capture. This limitation occurs for waves of moderate steepness and a narrow underlying spectrum

(2,0)-Super-Yang-Mills coupled to non-linear {sigma}-model

Energy Technology Data Exchange (ETDEWEB)

Goes-Negrao, M.S.; Penna-Firme, A.B. [Centro Brasileiro de Pesquisas Fisicas (CBPF), Rio de Janeiro, RJ (Brazil); Negrao, M.R. [Universidade Federal, Rio de Janeiro, RJ (Brazil). Inst. de Fisica

1999-07-01

Considering a class of (2,0)-super yang-Mills multiplets that accommodate a pair of independent gauge potentials in connection with a single symmetry group, we present here their coupling to ordinary matter to non-linear {sigma}-models in (2,0)-superspace. The dynamics and the coupling of the gauge potentials are discussed and the interesting feature that comes out is a sort of chirality for one of the gauge potentials are discussed and the interesting feature that comes out is a sort of chirality for one of the gauge potentials once light-cone coordinates are chosen. (author)

(2,0)-Super-Yang-Mills coupled to non-linear σ-model

International Nuclear Information System (INIS)

Goes-Negrao, M.S.; Penna-Firme, A.B.; Negrao, M.R.

1999-07-01

Considering a class of (2,0)-super yang-Mills multiplets that accommodate a pair of independent gauge potentials in connection with a single symmetry group, we present here their coupling to ordinary matter to non-linear σ-models in (2,0)-superspace. The dynamics and the coupling of the gauge potentials are discussed and the interesting feature that comes out is a sort of chirality for one of the gauge potentials are discussed and the interesting feature that comes out is a sort of chirality for one of the gauge potentials once light-cone coordinates are chosen. (author)

Linear dynamic coupling in geared rotor systems

Science.gov (United States)

David, J. W.; Mitchell, L. D.

1986-01-01

The effects of high frequency oscillations caused by the gear mesh, on components of a geared system that can be modeled as rigid discs are analyzed using linear dynamic coupling terms. The coupled, nonlinear equations of motion for a disc attached to a rotating shaft are presented. The results of a trial problem analysis show that the inclusion of the linear dynamic coupling terms can produce significant changes in the predicted response of geared rotor systems, and that the produced sideband responses are greater than the unbalanced response. The method is useful in designing gear drives for heavy-lift helicopters, industrial speed reducers, naval propulsion systems, and heavy off-road equipment.

Optimal control of coupled parabolic-hyperbolic non-autonomous PDEs: infinite-dimensional state-space approach

Science.gov (United States)

Aksikas, I.; Moghadam, A. Alizadeh; Forbes, J. F.

2018-04-01

This paper deals with the design of an optimal state-feedback linear-quadratic (LQ) controller for a system of coupled parabolic-hypebolic non-autonomous partial differential equations (PDEs). The infinite-dimensional state space representation and the corresponding operator Riccati differential equation are used to solve the control problem. Dynamical properties of the coupled system of interest are analysed to guarantee the existence and uniqueness of the solution of the LQ-optimal control problem and also to guarantee the exponential stability of the closed-loop system. Thanks to the eigenvalues and eigenfunctions of the parabolic operator and also the fact that the hyperbolic-associated operator Riccati differential equation can be converted to a scalar Riccati PDE, an algorithm to solve the LQ control problem has been presented. The results are applied to a non-isothermal packed-bed catalytic reactor. The LQ optimal controller designed in the early portion of the paper is implemented for the original non-linear model. Numerical simulations are performed to show the controller performances.

Linear perturbations of a self-similar solution of hydrodynamics with non-linear heat conduction

International Nuclear Information System (INIS)

Dubois-Boudesocque, Carine

2000-01-01

The stability of an ablative flow, where a shock wave is located upstream a thermal front, is of importance in inertial confinement fusion. The present model considers an exact self-similar solution to the hydrodynamic equations with non-linear heat conduction for a semi-infinite slab. For lack of an analytical solution, a high resolution numerical procedure is devised, which couples a finite difference method with a relaxation algorithm using a two-domain pseudo-spectral method. Stability of this solution is studied by introducing linear perturbation method within a Lagrangian-Eulerian framework. The initial and boundary value problem is solved by a splitting of the equations between a hyperbolic system and a parabolic equation. The boundary conditions of the hyperbolic system are treated, in the case of spectral methods, according to Thompson's approach. The parabolic equation is solved by an influence matrix method. These numerical procedures have been tested versus exact solutions. Considering a boundary heat flux perturbation, the space-time evolution of density, velocity and temperature are shown. (author) [fr

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.

Non-linear coupling of the lower hybrid grill in ASDEX

International Nuclear Information System (INIS)

Petrzilka, V.A.

1991-01-01

Computations of the reflection coefficient based on a non-linear lower hybrid (LH) coupling theory are presented and compared with the measurements of the reflection coefficient of the ASDEX tokamak LH grill, where powers up to 4 kW/cm 2 have been launched. This high LH power density modifies the electron density in front of the grill because of ponderomotive forces. Thus, the coupling and the power reflection coefficient change. To explain the observed saturation of the growth of the reflection coefficient with power, it is necessary to take into account some heating of the plasma in front of the grill by the transmitted LH power, which also leads to a poloidally inhomogeneous edge electron density. (author). Letter-to-the-editor. 14 refs, 13 figs

Non-linear coupling of the lower hybrid grill in ASDEX

Energy Technology Data Exchange (ETDEWEB)

Petrzilka, V A [Ceskoslovenska Akademie Ved, Prague (Czechoslovakia). Ustav Fyziky Plazmatu; Leuterer, F; Soeldner, F X; Giannone, L.; Schubert, R [Association Euratom-Max-Planck-Institut fuer Plasmaphysik, Garching (Germany, F.R.)

1991-09-01

Computations of the reflection coefficient based on a non-linear lower hybrid (LH) coupling theory are presented and compared with the measurements of the reflection coefficient of the ASDEX tokamak LH grill, where powers up to 4 kW/cm{sup 2} have been launched. This high LH power density modifies the electron density in front of the grill because of ponderomotive forces. Thus, the coupling and the power reflection coefficient change. To explain the observed saturation of the growth of the reflection coefficient with power, it is necessary to take into account some heating of the plasma in front of the grill by the transmitted LH power, which also leads to a poloidally inhomogeneous edge electron density. (author). Letter-to-the-editor. 14 refs, 13 figs.

Exponential synchronization of two nonlinearly non-delayed and delayed coupled complex dynamical networks

International Nuclear Information System (INIS)

Zheng Song

2012-01-01

In this paper, the exponential synchronization between two nonlinearly coupled complex networks with non-delayed and delayed coupling is investigated with Lyapunov-Krasovskii-type functionals. Based on the stability analysis of the impulsive differential equation and the linear matrix inequality, sufficient delay-dependent conditions for exponential synchronization are derived, and a linear impulsive controller and simple updated laws are also designed. Furthermore, the coupling matrices need not be symmetric or irreducible. Numerical examples are presented to verify the effectiveness and correctness of the synchronization criteria obtained.

KAM for the non-linear Schroedinger equation

CERN Document Server

Eliasson, L H

2006-01-01

We consider the $d$-dimensional nonlinear Schr\\"o\\-dinger equation under periodic boundary conditions:-i\\dot u=\\Delta u+V(x)*u+\\ep|u|^2u;\\quad u=u(t,x),\\;x\\in\\T^dwhere $V(x)=\\sum \\hat V(a)e^{i\\sc{a,x}}$ is an analytic function with $\\hat V$ real. (This equation is a popular model for the `real' NLS equation, where instead of the convolution term $V*u$ we have the potential term $Vu$.) For $\\ep=0$ the equation is linear and has time--quasi-periodic solutions $u$,u(t,x)=\\sum_{s\\in \\AA}\\hat u_0(a)e^{i(|a|^2+\\hat V(a))t}e^{i\\sc{a,x}}, \\quad 0<|\\hat u_0(a)|\\le1,where $\\AA$ is any finite subset of $\\Z^d$. We shall treat $\\omega_a=|a|^2+\\hat V(a)$, $a\\in\\AA$, as free parameters in some domain $U\\subset\\R^{\\AA}$. This is a Hamiltonian system in infinite degrees of freedom, degenerate but with external parameters, and we shall describe a KAM-theory which, in particular, will have the following consequence: \\smallskip {\\it If $|\\ep|$ is sufficiently small, then there is a large subset $U'$ of $U$ such that for all $...

Infinite sets of conservation laws for linear and nonlinear field equations

International Nuclear Information System (INIS)

Mickelsson, J.

1984-01-01

The relation between an infinite set of conservation laws of a linear field equation and the enveloping algebra of the space-time symmetry group is established. It is shown that each symmetric element of the enveloping algebra of the space-time symmetry group of a linear field equation generates a one-parameter group of symmetries of the field equation. The cases of the Maxwell and Dirac equations are studied in detail. Then it is shown that (at least in the sense of a power series in the 'coupling constant') the conservation laws of the linear case can be deformed to conservation laws of a nonlinear field equation which is obtained from the linear one by adding a nonlinear term invariant under the group of space-time symmetries. As an example, our method is applied to the Korteweg-de Vries equation and to the massless Thirring model. (orig.)

Stationary solution of the Fokker-Planck equation for linearly coupled motion in an electron storage ring

International Nuclear Information System (INIS)

Chao, A.W.; Lee, M.J.

1975-09-01

Effects upon longitudinal bunch shape in a storage ring due to linear and nonlinear potential can be calculated by finding the stationary solution to the Fokker-Planck equation for the particle distribution. Effects upon transverse bunch shape of a stored electron beam due to photon emissions and damping can be calculated by this method. It has been found that this method can also be used for a case in which the transverse modes of oscillation are coupled to the energy deviation δ. Examples of lattice elements which produce linear coupling between these oscillations are skew quadrupole magnets and solenoid magnets. For the linearly coupled case the stationary solution has been found to be given by exp (ΣΣA/sub ij/ x/sub i/x/sub j/) with x/sub i/ the canonical variables (x,p/sub x/, y, p/sub y/, δ, p/sub δ/) and A /sub ij/ some constants. The solution for the values of A /sub ij/'s will be described in this report. It will be shown that this solution can be expressed in a compact form. For simple cases, this form of solution leads directly to analytic expressions for the values of A /sub ij/'s and the bunch shape can be calculated by integrating the distribution function over some of the coordinates; for the more complex cases, it can be conveniently adapted as an algorithm for numerical evaluation. 16 refs

Linear and non-linear optics of condensed matter

International Nuclear Information System (INIS)

McLean, T.P.

1977-01-01

Part I - Linear optics: 1. General introduction. 2. Frequency dependence of epsilon(ω, k vector). 3. Wave-vector dependence of epsilon(ω, k vector). 4. Tensor character of epsilon(ω, k vector). Part II - Non-linear optics: 5. Introduction. 6. A classical theory of non-linear response in one dimension. 7. The generalization to three dimensions. 8. General properties of the polarizability tensors. 9. The phase-matching condition. 10. Propagation in a non-linear dielectric. 11. Second harmonic generation. 12. Coupling of three waves. 13. Materials and their non-linearities. 14. Processes involving energy exchange with the medium. 15. Two-photon absorption. 16. Stimulated Raman effect. 17. Electro-optic effects. 18. Limitations of the approach presented here. (author)

Convergence of hybrid methods for solving non-linear partial ...

African Journals Online (AJOL)

This paper is concerned with the numerical solution and convergence analysis of non-linear partial differential equations using a hybrid method. The solution technique involves discretizing the non-linear system of PDE to obtain a corresponding non-linear system of algebraic difference equations to be solved at each time ...

Mathematical problems in non-linear Physics: some results

International Nuclear Information System (INIS)

1979-01-01

The basic results presented in this report are the following: 1) Characterization of the range and Kernel of the variational derivative. 2) Determination of general conservation laws in linear evolution equations, as well as bounds for the number of polynomial conserved densities in non-linear evolution equations in two independent variables of even order. 3) Construction of the most general evolution equation which has a given family of conserved densities. 4) Regularity conditions for the validity of the Lie invariance method. 5) A simple class of perturbations in non-linear wave equations. 6) Soliton solutions in generalized KdV equations. (author)

Rate equation analysis and non-Hermiticity in coupled semiconductor laser arrays

Science.gov (United States)

Gao, Zihe; Johnson, Matthew T.; Choquette, Kent D.

2018-05-01

Optically coupled semiconductor laser arrays are described by coupled rate equations. The coupled mode equations and carrier densities are included in the analysis, which inherently incorporate the carrier-induced nonlinearities including gain saturation and amplitude-phase coupling. We solve the steady-state coupled rate equations and consider the cavity frequency detuning and the individual laser pump rates as the experimentally controlled variables. We show that the carrier-induced nonlinearities play a critical role in the mode control, and we identify gain contrast induced by cavity frequency detuning as a unique mechanism for mode control. Photon-mediated energy transfer between cavities is also discussed. Parity-time symmetry and exceptional points in this system are studied. Unbroken parity-time symmetry can be achieved by judiciously combining cavity detuning and unequal pump rates, while broken symmetry lies on the boundary of the optical locking region. Exceptional points are identified at the intersection between broken symmetry and unbroken parity-time symmetry.

Some results on the neutron transport and the coupling of equations

International Nuclear Information System (INIS)

Bal, G.

1997-01-01

Neutron transport in nuclear reactors is well modeled by the linear Boltzmann transport equation. Its resolution is relatively easy but very expensive. To achieve whole core calculations, one has to consider simpler models, such as diffusion or homogeneous transport equations. However, the solutions may become inaccurate in particular situations (as accidents for instance). That is the reason why we wish to solve the equations on small area accurately and more coarsely on the remaining part of the core. It is than necessary to introduce some links between different discretizations or modelizations. In this note, we give some results on the coupling of different discretizations of all degrees of freedom of the integral-differential neutron transport equation (two degrees for the angular variable, on for the energy component, and two or three degrees for spatial position respectively in 2D (cylindrical symmetry) and 3D). Two chapters are devoted to the coupling of discrete ordinates methods (for angular discretization). The first one is theoretical and shows the well posing of the coupled problem, whereas the second one deals with numerical applications of practical interest (the results have been obtained from the neutron transport code developed at the R and D, which has been modified for introducing the coupling). Next, we present the nodal scheme RTN0, used for the spatial discretization. We show well posing results for the non-coupled and the coupled problems. At the end, we deal with the coupling of energy discretizations for the multigroup equations obtained by homogenization. Some theoretical results of the discretization of the velocity variable (well-posing of problems), which do not deal directly with the purposes of coupling, are presented in the annexes. (author)

Non linear photons: a non singular cosmological solution

International Nuclear Information System (INIS)

Alves, G.A.

1986-01-01

The validity of equivalence principle as principle of minimum coupling between field interactions, is discussed. The non minimum coupling between vector field and gravitational field, and some consequences of this coupling are analysed. Starting from spherical symmetry metric, the coupled field equations, obtaining exact solutions and interpreting these solutions, are solved. (M.C.K.) [pt

A novel mixed-synchronization phenomenon in coupled Chua's circuits via non-fragile linear control

International Nuclear Information System (INIS)

Wang Jun-Wei; Ma Qing-Hua; Zeng Li

2011-01-01

Dynamical variables of coupled nonlinear oscillators can exhibit different synchronization patterns depending on the designed coupling scheme. In this paper, a non-fragile linear feedback control strategy with multiplicative controller gain uncertainties is proposed for realizing the mixed-synchronization of Chua's circuits connected in a drive-response configuration. In particular, in the mixed-synchronization regime, different state variables of the response system can evolve into complete synchronization, anti-synchronization and even amplitude death simultaneously with the drive variables for an appropriate choice of scaling matrix. Using Lyapunov stability theory, we derive some sufficient criteria for achieving global mixed-synchronization. It is shown that the desired non-fragile state feedback controller can be constructed by solving a set of linear matrix inequalities (LMIs). Numerical simulations are also provided to demonstrate the effectiveness of the proposed control approach. (general)

Non-linear mixing in coupled photonic crystal nanobeam cavities due to cross-coupling opto-mechanical mechanisms

Energy Technology Data Exchange (ETDEWEB)

Ramos, Daniel, E-mail: daniel.ramos@csic.es; Frank, Ian W.; Deotare, Parag B.; Bulu, Irfan; Lončar, Marko [School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138 (United States)

2014-11-03

We investigate the coupling between mechanical and optical modes supported by coupled, freestanding, photonic crystal nanobeam cavities. We show that localized cavity modes for a given gap between the nanobeams provide weak optomechanical coupling with out-of-plane mechanical modes. However, we show that the coupling can be significantly increased, more than an order of magnitude for the symmetric mechanical mode, due to optical resonances that arise from the interaction of the localized cavity modes with standing waves formed by the reflection from thesubstrate. Finally, amplification of motion for the symmetric mode has been observed and attributed to the strong optomechanical interaction of our hybrid system. The amplitude of these self-sustained oscillations is large enough to put the system into a non-linear oscillation regime where a mixing between the mechanical modes is experimentally observed and theoretically explained.

Non-linear seismic analysis of structures coupled with fluid

International Nuclear Information System (INIS)

Descleve, P.; Derom, P.; Dubois, J.

1983-01-01

This paper presents a method to calculate non-linear structure behaviour under horizontal and vertical seismic excitation, making possible the full non-linear seismic analysis of a reactor vessel. A pseudo forces method is used to introduce non linear effects and the problem is solved by superposition. Two steps are used in the method: - Linear calculation of the complete model. - Non linear analysis of thin shell elements and calculation of seismic induced pressure originating from linear and non linear effects, including permanent loads and thermal stresses. Basic aspects of the mathematical formulation are developed. It has been applied to axi-symmetric shell element using a Fourier series solution. For the fluid interaction effect, a comparison is made with a dynamic test. In an example of application, the displacement and pressure time history are given. (orig./GL)

Space-time coupled spectral/hp least-squares finite element formulation for the incompressible Navier-Stokes equations

International Nuclear Information System (INIS)

Pontaza, J.P.; Reddy, J.N.

2004-01-01

We consider least-squares finite element models for the numerical solution of the non-stationary Navier-Stokes equations governing viscous incompressible fluid flows. The paper presents a formulation where the effects of space and time are coupled, resulting in a true space-time least-squares minimization procedure, as opposed to a space-time decoupled formulation where a least-squares minimization procedure is performed in space at each time step. The formulation is first presented for the linear advection-diffusion equation and then extended to the Navier-Stokes equations. The formulation has no time step stability restrictions and is spectrally accurate in both space and time. To allow the use of practical C 0 element expansions in the resulting finite element model, the Navier-Stokes equations are expressed as an equivalent set of first-order equations by introducing vorticity as an additional independent variable and the least-squares method is used to develop the finite element model of the governing equations. High-order element expansions are used to construct the discrete model. The discrete model thus obtained is linearized by Newton's method, resulting in a linear system of equations with a symmetric positive definite coefficient matrix that is solved in a fully coupled manner by a preconditioned conjugate gradient method in matrix-free form. Spectral convergence of the L 2 least-squares functional and L 2 error norms in space-time is verified using a smooth solution to the two-dimensional non-stationary incompressible Navier-Stokes equations. Numerical results are presented for impulsively started lid-driven cavity flow, oscillatory lid-driven cavity flow, transient flow over a backward-facing step, and flow around a circular cylinder; the results demonstrate the predictive capability and robustness of the proposed formulation. Even though the space-time coupled formulation is emphasized, we also present the formulation and numerical results for least

Experimental characterization and modelling of non-linear coupling of the lower hybrid current drive power on Tore Supra

Science.gov (United States)

Preynas, M.; Goniche, M.; Hillairet, J.; Litaudon, X.; Ekedahl, A.; Colas, L.

2013-01-01

To achieve steady-state operation on future fusion devices, in particular on ITER, the coupling of the lower hybrid wave must be optimized on a wide range of edge conditions. However, under some specific conditions, deleterious effects on the lower hybrid current drive (LHCD) coupling are sometimes observed on Tore Supra. In this way, dedicated LHCD experiments have been performed using the LHCD system of Tore Supra, composed of two different conceptual designs of launcher: the fully active multi-junction (FAM) and the new passive active multi-junction (PAM) antennas. A non-linear interaction between the electron density and the electric field has been characterized in a thin plasma layer in front of the two LHCD antennas. The resulting dependence of the power reflection coefficient (RC) with the LHCD power is not predicted by the standard linear theory of the LH wave coupling. A theoretical model is suggested to describe the non-linear wave-plasma interaction induced by the ponderomotive effect and implemented in a new full wave LHCD code, PICCOLO-2D (ponderomotive effect in a coupling code of lower hybrid wave-2D). The code self-consistently treats the wave propagation in the antenna vicinity and its interaction with the local edge plasma density. The simulation reproduces very well the occurrence of a non-linear behaviour in the coupling observed in the LHCD experiments. The important differences and trends between the FAM and the PAM antennas, especially a larger increase in RC for the FAM, are also reproduced by the PICCOLO-2D simulation. The working hypothesis of the contribution of the ponderomotive effect in the non-linear observations of LHCD coupling is therefore validated through this comprehensive modelling for the first time on the FAM and PAM antennas on Tore Supra.

Experimental characterization and modelling of non-linear coupling of the lower hybrid current drive power on Tore Supra

International Nuclear Information System (INIS)

Preynas, M.; Goniche, M.; Hillairet, J.; Litaudon, X.; Ekedahl, A.; Colas, L.

2013-01-01

To achieve steady-state operation on future fusion devices, in particular on ITER, the coupling of the lower hybrid wave must be optimized on a wide range of edge conditions. However, under some specific conditions, deleterious effects on the lower hybrid current drive (LHCD) coupling are sometimes observed on Tore Supra. In this way, dedicated LHCD experiments have been performed using the LHCD system of Tore Supra, composed of two different conceptual designs of launcher: the fully active multi-junction (FAM) and the new passive active multi-junction (PAM) antennas. A non-linear interaction between the electron density and the electric field has been characterized in a thin plasma layer in front of the two LHCD antennas. The resulting dependence of the power reflection coefficient (RC) with the LHCD power is not predicted by the standard linear theory of the LH wave coupling. A theoretical model is suggested to describe the non-linear wave–plasma interaction induced by the ponderomotive effect and implemented in a new full wave LHCD code, PICCOLO-2D (ponderomotive effect in a coupling code of lower hybrid wave-2D). The code self-consistently treats the wave propagation in the antenna vicinity and its interaction with the local edge plasma density. The simulation reproduces very well the occurrence of a non-linear behaviour in the coupling observed in the LHCD experiments. The important differences and trends between the FAM and the PAM antennas, especially a larger increase in RC for the FAM, are also reproduced by the PICCOLO-2D simulation. The working hypothesis of the contribution of the ponderomotive effect in the non-linear observations of LHCD coupling is therefore validated through this comprehensive modelling for the first time on the FAM and PAM antennas on Tore Supra. (paper)

Engineering equations for characterizing non-linear laser intensity propagation in air with loss.

Science.gov (United States)

Karr, Thomas; Stotts, Larry B; Tellez, Jason A; Schmidt, Jason D; Mansell, Justin D

2018-02-19

The propagation of high peak-power laser beams in real atmospheres will be affected at long range by both linear and nonlinear effects contained therein. Arguably, J. H. Marburger is associated with the mathematical characterization of this phenomenon. This paper provides a validated set of engineering equations for characterizing the self-focusing distance from a laser beam propagating through non-turbulent air with, and without, loss as well as three source configurations: (1) no lens, (2) converging lens and (3) diverging lens. The validation was done against wave-optics simulation results. Some validated equations follow Marburger completely, but others do not, requiring modification of the original theory. Our results can provide a guide for numerical simulations and field experiments.

Coupled-channel equations and off-shell transformations in many-body scattering

International Nuclear Information System (INIS)

Cattapan, G.; Vanzani, V.

1977-01-01

The general structure and the basic features of several many-body coupled-channel integral equations, obtained by means of the channel coupling array device, are studied in a systematic way. Particular attention is paid to the employment of symmetric transition operators. The connection between different formulations has been clarified and the role played by some off-shell transformations for many-body transition operators has been discussed. Specific choices of the coupling scheme are considered and the corresponding coupled equations are compared with similar equations previously derived. Several sets of linear relations between transition operators have also been presented and used in a three-body context to derive uncoupled integral equations with connected kernel

Impact of quadratic non-linearity on the dynamics of periodic solutions of a wave equation

International Nuclear Information System (INIS)

Kolesov, Andrei Yu; Rozov, Nikolai Kh

2002-01-01

For the non-linear telegraph equation with homogeneous Dirichlet or Neumann conditions at the end-points of a finite interval the question of the existence and the stability of time-periodic solutions bifurcating from the zero equilibrium state is considered. The dynamics of these solutions under a change of the diffusion coefficient (that is, the coefficient of the second derivative with respect to the space variable) is investigated. For the Dirichlet boundary conditions it is shown that this dynamics substantially depends on the presence - or the absence - of quadratic terms in the non-linearity. More precisely, it is shown that a quadratic non-linearity results in the occurrence, under an unbounded decrease of diffusion, of an infinite sequence of bifurcations of each periodic solution. En route, the related issue of the limits of applicability of Yu.S. Kolesov's method of quasinormal forms to the construction of self-oscillations in singularly perturbed hyperbolic boundary value problems is studied

The linear parameters and the decoupling matrix for linearly coupled motion in 6 dimensional phase space

International Nuclear Information System (INIS)

Parzen, G.

1997-01-01

It will be shown that starting from a coordinate system where the 6 phase space coordinates are linearly coupled, one can go to a new coordinate system, where the motion is uncoupled, by means of a linear transformation. The original coupled coordinates and the new uncoupled coordinates are related by a 6 x 6 matrix, R. It will be shown that of the 36 elements of the 6 x 6 decoupling matrix R, only 12 elements are independent. A set of equations is given from which the 12 elements of R can be computed form the one period transfer matrix. This set of equations also allows the linear parameters, the β i , α i , i = 1, 3, for the uncoupled coordinates, to be computed from the one period transfer matrix

Hartree Fock-type equations in relativistic quantum electrodynamics with non-linear gauge fixing

International Nuclear Information System (INIS)

Dietz, K.; Hess, B.A.

1990-08-01

Relativistic mean-field equations are obtained by minimizing the effective energy obtained from the gauge-invariant energy density by eliminating electro-magnetic degrees of freedom in certain characteristic non-linear gauges. It is shown that by an appropriate choice of gauge many-body correlations, e.g. screening, three-body 'forces' etc. can be included already at the mean-field level. The many-body perturbation theory built on the latter is then expected to show improved 'convergence'. (orig.)

Linear vs non-linear QCD evolution: from HERA data to LHC phenomenology

CERN Document Server

Albacete, J L; Quiroga-Arias, P; Rojo, J

2012-01-01

The very precise combined HERA data provides a testing ground in which the relevance of novel QCD regimes, other than the successful linear DGLAP evolution, in small-x inclusive DIS data can be ascertained. We present a study of the dependence of the AAMQS fits, based on the running coupling BK non-linear evolution equations (rcBK), on the fitted dataset. This allows for the identification of the kinematical region where rcBK accurately describes the data, and thus for the determination of its applicability boundary. We compare the rcBK results with NNLO DGLAP fits, obtained with the NNPDF methodology with analogous kinematical cuts. Further, we explore the impact on LHC phenomenology of applying stringent kinematical cuts to the low-x HERA data in a DGLAP fit.

Coupled radiative transfer equation and diffusion approximation model for photon migration in turbid medium with low-scattering and non-scattering regions

International Nuclear Information System (INIS)

Tarvainen, Tanja; Vauhkonen, Marko; Kolehmainen, Ville; Arridge, Simon R; Kaipio, Jari P

2005-01-01

In this paper, a coupled radiative transfer equation and diffusion approximation model is extended for light propagation in turbid medium with low-scattering and non-scattering regions. The light propagation is modelled with the radiative transfer equation in sub-domains in which the assumptions of the diffusion approximation are not valid. The diffusion approximation is used elsewhere in the domain. The two equations are coupled through their boundary conditions and they are solved simultaneously using the finite element method. The streamline diffusion modification is used to avoid the ray-effect problem in the finite element solution of the radiative transfer equation. The proposed method is tested with simulations. The results of the coupled model are compared with the finite element solutions of the radiative transfer equation and the diffusion approximation and with results of Monte Carlo simulation. The results show that the coupled model can be used to describe photon migration in turbid medium with low-scattering and non-scattering regions more accurately than the conventional diffusion model

Generalized non-linear Schroedinger hierarchy

International Nuclear Information System (INIS)

Aratyn, H.; Gomes, J.F.; Zimerman, A.H.

1994-01-01

The importance in studying the completely integrable models have became evident in the last years due to the fact that those models present an algebraic structure extremely rich, providing the natural scenery for solitons description. Those models can be described through non-linear differential equations, pseudo-linear operators (Lax formulation), or a matrix formulation. The integrability implies in the existence of a conservation law associated to each of degree of freedom. Each conserved charge Q i can be associated to a Hamiltonian, defining a time evolution related to to a time t i through the Hamilton equation ∂A/∂t i =[A,Q i ]. Particularly, for a two-dimensions field theory, infinite degree of freedom exist, and consequently infinite conservation laws describing the time evolution in space of infinite times. The Hamilton equation defines a hierarchy of models which present a infinite set of conservation laws. This paper studies the generalized non-linear Schroedinger hierarchy

An implicit meshless scheme for the solution of transient non-linear Poisson-type equations

KAUST Repository

Bourantas, Georgios

2013-07-01

A meshfree point collocation method is used for the numerical simulation of both transient and steady state non-linear Poisson-type partial differential equations. Particular emphasis is placed on the application of the linearization method with special attention to the lagging of coefficients method and the Newton linearization method. The localized form of the Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are performed for regular nodal distributions, stressing the positivity conditions that make the resulting system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through representative and well-established benchmark problems. © 2013 Elsevier Ltd.

An implicit meshless scheme for the solution of transient non-linear Poisson-type equations

KAUST Repository

Bourantas, Georgios; Burganos, Vasilis N.

2013-01-01

A meshfree point collocation method is used for the numerical simulation of both transient and steady state non-linear Poisson-type partial differential equations. Particular emphasis is placed on the application of the linearization method with special attention to the lagging of coefficients method and the Newton linearization method. The localized form of the Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are performed for regular nodal distributions, stressing the positivity conditions that make the resulting system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through representative and well-established benchmark problems. © 2013 Elsevier Ltd.

Two-dimensional differential transform method for solving linear and non-linear Schroedinger equations

International Nuclear Information System (INIS)

Ravi Kanth, A.S.V.; Aruna, K.

2009-01-01

In this paper, we propose a reliable algorithm to develop exact and approximate solutions for the linear and nonlinear Schroedinger equations. The approach rest mainly on two-dimensional differential transform method which is one of the approximate methods. The method can easily be applied to many linear and nonlinear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. Several illustrative examples are given to demonstrate the effectiveness of the present method.

Linear integral equations and soliton systems

International Nuclear Information System (INIS)

Quispel, G.R.W.

1983-01-01

A study is presented of classical integrable dynamical systems in one temporal and one spatial dimension. The direct linearizations are given of several nonlinear partial differential equations, for example the Korteweg-de Vries equation, the modified Korteweg-de Vries equation, the sine-Gordon equation, the nonlinear Schroedinger equation, and the equation of motion for the isotropic Heisenberg spin chain; the author also discusses several relations between these equations. The Baecklund transformations of these partial differential equations are treated on the basis of a singular transformation of the measure (or equivalently of the plane-wave factor) occurring in the corresponding linear integral equations, and the Baecklund transformations are used to derive the direct linearization of a chain of so-called modified partial differential equations. Finally it is shown that the singular linear integral equations lead in a natural way to the direct linearizations of various nonlinear difference-difference equations. (Auth.)

Algebraic inversion of the Dirac equation for the vector potential in the non-Abelian case

International Nuclear Information System (INIS)

Inglis, S M; Jarvis, P D

2012-01-01

We study the Dirac equation for spinor wavefunctions minimally coupled to an external field, from the perspective of an algebraic system of linear equations for the vector potential. By analogy with the method in electromagnetism, which has been well-studied, and leads to classical solutions of the Maxwell–Dirac equations, we set up the formalism for non-Abelian gauge symmetry, with the SU(2) group and the case of four-spinor doublets. An extended isospin-charge conjugation operator is defined, enabling the hermiticity constraint on the gauge potential to be imposed in a covariant fashion, and rendering the algebraic system tractable. The outcome is an invertible linear equation for the non-Abelian vector potential in terms of bispinor current densities. We show that, via application of suitable extended Fierz identities, the solution of this system for the non-Abelian vector potential is a rational expression involving only Pauli scalar and Pauli triplet, Lorentz scalar, vector and axial vector current densities, albeit in the non-closed form of a Neumann series. (paper)

Stability of non-linear constitutive formulations for viscoelastic fluids

CERN Document Server

Siginer, Dennis A

2014-01-01

Stability of Non-linear Constitutive Formulations for Viscoelastic Fluids provides a complete and up-to-date view of the field of constitutive equations for flowing viscoelastic fluids, in particular on their non-linear behavior, the stability of these constitutive equations that is their predictive power, and the impact of these constitutive equations on the dynamics of viscoelastic fluid flow in tubes. This book gives an overall view of the theories and attendant methodologies developed independently of thermodynamic considerations as well as those set within a thermodynamic framework to derive non-linear rheological constitutive equations for viscoelastic fluids. Developments in formulating Maxwell-like constitutive differential equations as well as single integral constitutive formulations are discussed in the light of Hadamard and dissipative type of instabilities.

Resummation of the 1/N-expansion of the non-linear σ-model by Dyson-Schwinger equations

International Nuclear Information System (INIS)

Drouffe, J.M.; Flyvbjerg, H.

1988-02-01

Dyson-Schwinger equations for the O(N)-symmetric non-linear σ-model are derived and expanded in 1/N. A closed set of equations is obtained by keeping only the leading term and the first correction term in this expansion. These equations are solved numerically in 2 dimensions on square lattices of sizes 50x50 and 100x100. Results for the magnetic susceptibility and the mass gap are compared with predictions of the ordinary 1/N-expansion and with Monte Carlo results. The results obtained with the Dyson-Schwinger equations show the same scaling behavior as found in the Monte Carlo results. This is not the behavior predicted by the perturbative renormalization group. (orig.)

Novel method for solution of coupled radial Schrödinger equations

International Nuclear Information System (INIS)

Ershov, S. N.; Vaagen, J. S.; Zhukov, M. V.

2011-01-01

One of the major problems in numerical solution of coupled differential equations is the maintenance of linear independence for different sets of solution vectors. A novel method for solution of radial Schrödinger equations is suggested. It consists of rearrangement of coupled equations in a way that is appropriate to avoid usual numerical instabilities associated with components of the wave function in their classically forbidden regions. Applications of the new method for nuclear structure calculations within the hyperspherical harmonics approach are given.

Numerical solution of two-dimensional non-linear partial differential ...

African Journals Online (AJOL)

linear partial differential equations using a hybrid method. The solution technique involves discritizing the non-linear system of partial differential equations (PDEs) to obtain a corresponding nonlinear system of algebraic difference equations to be ...

Linear versus non-linear supersymmetry, in general

Energy Technology Data Exchange (ETDEWEB)

Ferrara, Sergio [Theoretical Physics Department, CERN,CH-1211 Geneva 23 (Switzerland); INFN - Laboratori Nazionali di Frascati,Via Enrico Fermi 40, I-00044 Frascati (Italy); Department of Physics and Astronomy, UniversityC.L.A.,Los Angeles, CA 90095-1547 (United States); Kallosh, Renata [SITP and Department of Physics, Stanford University,Stanford, California 94305 (United States); Proeyen, Antoine Van [Institute for Theoretical Physics, Katholieke Universiteit Leuven,Celestijnenlaan 200D, B-3001 Leuven (Belgium); Wrase, Timm [Institute for Theoretical Physics, Technische Universität Wien,Wiedner Hauptstr. 8-10, A-1040 Vienna (Austria)

2016-04-12

We study superconformal and supergravity models with constrained superfields. The underlying version of such models with all unconstrained superfields and linearly realized supersymmetry is presented here, in addition to the physical multiplets there are Lagrange multiplier (LM) superfields. Once the equations of motion for the LM superfields are solved, some of the physical superfields become constrained. The linear supersymmetry of the original models becomes non-linearly realized, its exact form can be deduced from the original linear supersymmetry. Known examples of constrained superfields are shown to require the following LM’s: chiral superfields, linear superfields, general complex superfields, some of them are multiplets with a spin.

Linear versus non-linear supersymmetry, in general

International Nuclear Information System (INIS)

Ferrara, Sergio; Kallosh, Renata; Proeyen, Antoine Van; Wrase, Timm

2016-01-01

We study superconformal and supergravity models with constrained superfields. The underlying version of such models with all unconstrained superfields and linearly realized supersymmetry is presented here, in addition to the physical multiplets there are Lagrange multiplier (LM) superfields. Once the equations of motion for the LM superfields are solved, some of the physical superfields become constrained. The linear supersymmetry of the original models becomes non-linearly realized, its exact form can be deduced from the original linear supersymmetry. Known examples of constrained superfields are shown to require the following LM’s: chiral superfields, linear superfields, general complex superfields, some of them are multiplets with a spin.

Relativistic mean-field theory for unstable nuclei with non-linear σ and ω terms

International Nuclear Information System (INIS)

Sugahara, Y.; Toki, H.

1994-01-01

We search for a new parameter set for the description of stable as well as unstable nuclei in the wide mass range within the relativistic mean-field theory. We include a non-linear ω self-coupling term in addition to the non-linear σ self-coupling terms, the necessity of which is suggested by the relativistic Brueckner-Hartree-Fock (RBHF) theory of nuclear matter. We find two parameter sets, one of which is for nuclei above Z=20 and the other for nuclei below that. The calculated results agree very well with the existing data for finite nuclei. The parameter set for the heavy nuclei provides the equation of state of nuclear matter similar to the one of the RBHF theory. ((orig.))

Identification of a time-varying point source in a system of two coupled linear diffusion-advection- reaction equations: application to surface water pollution

International Nuclear Information System (INIS)

Hamdi, Adel

2009-01-01

This paper deals with the identification of a point source (localization of its position and recovering the history of its time-varying intensity function) that constitutes the right-hand side of the first equation in a system of two coupled 1D linear transport equations. Assuming that the source intensity function vanishes before reaching the final control time, we prove the identifiability of the sought point source from recording the state relative to the second coupled transport equation at two observation points framing the source region. Note that at least one of the two observation points should be strategic. We establish an identification method that uses these records to identify the source position as the root of a continuous and strictly monotonic function. Whereas the source intensity function is recovered using a recursive formula without any need of an iterative process. Some numerical experiments on a variant of the surface water pollution BOD–OD coupled model are presented

Ion-acoustic cnoidal wave and associated non-linear ion flux in dusty plasma

Energy Technology Data Exchange (ETDEWEB)

Jain, S. L. [Poornima Group of Institution, Sitapura, Jaipur 302022 (India); Tiwari, R. S. [Regional College for Education, Research and Technology, Jaipur 302022 (India); Mishra, M. K. [Department of Physics, University of Rajasthan, Jaipur 302004 (India)

2012-10-15

Using reductive perturbation method with appropriate boundary conditions, coupled evolution equations for first and second order potentials are derived for ion-acoustic waves in a collisionless, un-magnetized plasma consisting of hot isothermal electrons, cold ions, and massive mobile charged dust grains. The boundary conditions give rise to renormalization term, which enable us to eliminate secular contribution in higher order terms. Determining the non secular solution of these coupled equations, expressions for wave phase velocity and averaged non-linear ion flux associated with ion-acoustic cnoidal wave are obtained. Variation of the wave phase velocity and averaged non-linear ion flux as a function of modulus (k{sup 2}) dependent wave amplitude are numerically examined for different values of dust concentration, charge on dust grains, and mass ratio of dust grains with plasma ions. It is found that for a given amplitude, the presence of positively (negatively) charged dust grains in plasma decreases (increases) the wave phase velocity. This behavior is more pronounced with increase in dust concentrations or increase in charge on dust grains or decrease in mass ratio of dust grains. The averaged non-linear ion flux associated with wave is positive (negative) for negatively (positively) charged dust grains in the plasma and increases (decreases) with modulus (k{sup 2}) dependent wave amplitude. For given amplitude, it increases (decreases) as dust concentration or charge of negatively (positively) charged dust grains increases in the plasma.

Analysis of fractional non-linear diffusion behaviors based on Adomian polynomials

Directory of Open Access Journals (Sweden)

Wu Guo-Cheng

2017-01-01

Full Text Available A time-fractional non-linear diffusion equation of two orders is considered to investigate strong non-linearity through porous media. An equivalent integral equation is established and Adomian polynomials are adopted to linearize non-linear terms. With the Taylor expansion of fractional order, recurrence formulae are proposed and novel numerical solutions are obtained to depict the diffusion behaviors more accurately. The result shows that the method is suitable for numerical simulation of the fractional diffusion equations of multi-orders.

Transport equations, Level Set and Eulerian mechanics. Application to fluid-structure coupling

International Nuclear Information System (INIS)

Maitre, E.

2008-11-01

My works were devoted to numerical analysis of non-linear elliptic-parabolic equations, to neutron transport equation and to the simulation of fabrics draping. More recently I developed an Eulerian method based on a level set formulation of the immersed boundary method to deal with fluid-structure coupling problems arising in bio-mechanics. Some of the more efficient algorithms to solve the neutron transport equation make use of the splitting of the transport operator taking into account its characteristics. In the present work we introduced a new algorithm based on this splitting and an adaptation of minimal residual methods to infinite dimensional case. We present the case where the velocity space is of dimension 1 (slab geometry) and 2 (plane geometry) because the splitting is simpler in the former

On the removal of boundary errors caused by Runge-Kutta integration of non-linear partial differential equations

Science.gov (United States)

Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.

1994-01-01

It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.

Non-linear partial differential equations an algebraic view of generalized solutions

CERN Document Server

Rosinger, Elemer E

1990-01-01

A massive transition of interest from solving linear partial differential equations to solving nonlinear ones has taken place during the last two or three decades. The availability of better computers has often made numerical experimentations progress faster than the theoretical understanding of nonlinear partial differential equations. The three most important nonlinear phenomena observed so far both experimentally and numerically, and studied theoretically in connection with such equations have been the solitons, shock waves and turbulence or chaotical processes. In many ways, these phenomen

E11 and the non-linear dual graviton

Science.gov (United States)

Tumanov, Alexander G.; West, Peter

2018-04-01

The non-linear dual graviton equation of motion as well as the duality relation between the gravity and dual gravity fields are found in E theory by carrying out E11 variations of previously found equations of motion. As a result the equations of motion in E theory have now been found at the full non-linear level up to, and including, level three, which contains the dual graviton field. When truncated to contain fields at levels three and less, and the spacetime is restricted to be the familiar eleven dimensional space time, the equations are equivalent to those of eleven dimensional supergravity.

Linear superposition solutions to nonlinear wave equations

International Nuclear Information System (INIS)

Liu Yu

2012-01-01

The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K(2,2,1), the Oliver water wave equation, and the k(n, n) equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed

Analytical solutions of time-fractional models for homogeneous Gardner equation and non-homogeneous differential equations

Directory of Open Access Journals (Sweden)

Olaniyi Samuel Iyiola

2014-09-01

Full Text Available In this paper, we obtain analytical solutions of homogeneous time-fractional Gardner equation and non-homogeneous time-fractional models (including Buck-master equation using q-Homotopy Analysis Method (q-HAM. Our work displays the elegant nature of the application of q-HAM not only to solve homogeneous non-linear fractional differential equations but also to solve the non-homogeneous fractional differential equations. The presence of the auxiliary parameter h helps in an effective way to obtain better approximation comparable to exact solutions. The fraction-factor in this method gives it an edge over other existing analytical methods for non-linear differential equations. Comparisons are made upon the existence of exact solutions to these models. The analysis shows that our analytical solutions converge very rapidly to the exact solutions.

Dynamics of atom-field probability amplitudes in a coupled cavity system with Kerr non-linearity

Energy Technology Data Exchange (ETDEWEB)

Priyesh, K. V.; Thayyullathil, Ramesh Babu [Department of Physics, Cochin University of Science and Technology, Cochin (India)

2014-01-28

We have investigated the dynamics of two cavities coupled together via photon hopping, filled with Kerr non-linear medium and each containing a two level atom in it. The evolution of various atom (field) state probabilities of the coupled cavity system in two excitation sub space are obtained numerically. Detailed analysis has been done by taking different initial conditions of the system, with various coupling strengths and by varying the susceptibility of the medium. The role of susceptibility factor, on the dynamics atom field probability has been examined. In a coupled cavity system with strong photon hopping it is found that the susceptibility factor modifies the behaviour of probability amplitudes.

Adaptive discontinuous Galerkin methods for non-linear reactive flows

CERN Document Server

Uzunca, Murat

2016-01-01

The focus of this monograph is the development of space-time adaptive methods to solve the convection/reaction dominated non-stationary semi-linear advection diffusion reaction (ADR) equations with internal/boundary layers in an accurate and efficient way. After introducing the ADR equations and discontinuous Galerkin discretization, robust residual-based a posteriori error estimators in space and time are derived. The elliptic reconstruction technique is then utilized to derive the a posteriori error bounds for the fully discrete system and to obtain optimal orders of convergence. As coupled surface and subsurface flow over large space and time scales is described by (ADR) equation the methods described in this book are of high importance in many areas of Geosciences including oil and gas recovery, groundwater contamination and sustainable use of groundwater resources, storing greenhouse gases or radioactive waste in the subsurface.

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.

Inhomogeneous linear equation in Rota-Baxter algebra

OpenAIRE

Pietrzkowski, Gabriel

2014-01-01

We consider a complete filtered Rota-Baxter algebra of weight $\\lambda$ over a commutative ring. Finding the unique solution of a non-homogeneous linear algebraic equation in this algebra, we generalize Spitzer's identity in both commutative and non-commutative cases. As an application, considering the Rota-Baxter algebra of power series in one variable with q-integral as the Rota-Baxter operator, we show certain Eulerian identities.

Finite element procedures for coupled linear analysis of heat transfer, fluid and solid mechanics

Science.gov (United States)

Sutjahjo, Edhi; Chamis, Christos C.

1993-01-01

Coupled finite element formulations for fluid mechanics, heat transfer, and solid mechanics are derived from the conservation laws for energy, mass, and momentum. To model the physics of interactions among the participating disciplines, the linearized equations are coupled by combining domain and boundary coupling procedures. Iterative numerical solution strategy is presented to solve the equations, with the partitioning of temporal discretization implemented.

New approach to solve fully fuzzy system of linear equations using ...

Indian Academy of Sciences (India)

This paper proposes two new methods to solve fully fuzzy system of linear equations. The fuzzy system has been converted to a crisp system of linear equations by using single and double parametric form of fuzzy numbers to obtain the non-negative solution. Double parametric form of fuzzy numbers is defined and applied ...

Smooth, cusped, and discontinuous traveling waves in the periodic fluid resonance equation

Science.gov (United States)

Kruse, Matthew Thomas

The principal motivation for this dissertation is to extend the study of small amplitude high frequency wave propagation in solutions for hyperbolic conservation laws begun by A. Majda and R. Rosales in 1984. It was then that Majda and Rosales obtained equations governing the leading order wave amplitudes of resonantly interacting weakly nonlinear high frequency wave trains in the compressible Euler equations. The equations were obtained through systematic application of multiple scales and result in a pair of nonlinear acoustic wave equations coupled through a convolution operator. The extended solutions satisfy a pair of inviscid Burgers' equations coupled via a spatial convolution operator. Since then, many mathematicians have used this technique to extend the time validity of solutions to systems of equations other than the Euler equations and have arrived at similar nonlinear non-local systems. This work attempts to look at some of the basic features of the linear and nonlinear coupled and decoupled non- local equations, offering some analytic solutions and numerical insight into the phenomena associated with these equations. We do so by examining a single non-local linear equation, and then a single equation coupling a Burgers' nonlinearity with a linear convolution operator. The linear case is completely solvable. Analytic solutions are provided along with numerical results showing the fundamental properties of the linear non- local equations. In the nonlinear case some analytic solutions, including steady state profiles and traveling wave solutions, are provided along with a battery of numerical simulations. Evidence indicates the existence of attractors for solutions of the single equation with a single mode kernel. Provided resonant interaction takes place, the profile of the attractor is uniquely dependent on the kernel alone. Hamiltonian equations are obtained for both the linear and nonlinear equations with the condition that the resonant kernel must

Linear and non-linear calculations of the hose instability in the ion-focused regime

International Nuclear Information System (INIS)

Buchanan, H.L.

1982-01-01

A simple model is adopted to study the hose instability of an intense relativistic electron beam in a partially neutralized, low density ion channel (ion focused regime). Equations of motion for the beam and the channel are derived and linearized to obtain an approximate dispersion relation. The non-linear equations of motion are then solved numerically and the results compared to linearized data

Some results on the neutron transport and the coupling of equations; Quelques resultats sur le transport neutronique et le couplage d`equations

Energy Technology Data Exchange (ETDEWEB)

Bal, G. [Electricite de France (EDF), Direction des Etudes et Recherches, 92 - Clamart (France)

1997-12-31

Neutron transport in nuclear reactors is well modeled by the linear Boltzmann transport equation. Its resolution is relatively easy but very expensive. To achieve whole core calculations, one has to consider simpler models, such as diffusion or homogeneous transport equations. However, the solutions may become inaccurate in particular situations (as accidents for instance). That is the reason why we wish to solve the equations on small area accurately and more coarsely on the remaining part of the core. It is than necessary to introduce some links between different discretizations or modelizations. In this note, we give some results on the coupling of different discretizations of all degrees of freedom of the integral-differential neutron transport equation (two degrees for the angular variable, on for the energy component, and two or three degrees for spatial position respectively in 2D (cylindrical symmetry) and 3D). Two chapters are devoted to the coupling of discrete ordinates methods (for angular discretization). The first one is theoretical and shows the well posing of the coupled problem, whereas the second one deals with numerical applications of practical interest (the results have been obtained from the neutron transport code developed at the R and D, which has been modified for introducing the coupling). Next, we present the nodal scheme RTN0, used for the spatial discretization. We show well posing results for the non-coupled and the coupled problems. At the end, we deal with the coupling of energy discretizations for the multigroup equations obtained by homogenization. Some theoretical results of the discretization of the velocity variable (well-posing of problems), which do not deal directly with the purposes of coupling, are presented in the annexes. (author). 34 refs.

Applicability of refined Born approximation to non-linear equations

International Nuclear Information System (INIS)

Rayski, J.

1990-01-01

A computational method called ''Refined Born Approximation'', formerly applied exclusively to linear problems, is shown to be successfully applicable also to non-linear problems enabling me to compute bifurcations and other irregular solutions which cannot be obtained by the standard perturbation procedures. (author)

Non-linear auto-regressive models for cross-frequency coupling in neural time series

Science.gov (United States)

Tallot, Lucille; Grabot, Laetitia; Doyère, Valérie; Grenier, Yves; Gramfort, Alexandre

2017-01-01

We address the issue of reliably detecting and quantifying cross-frequency coupling (CFC) in neural time series. Based on non-linear auto-regressive models, the proposed method provides a generative and parametric model of the time-varying spectral content of the signals. As this method models the entire spectrum simultaneously, it avoids the pitfalls related to incorrect filtering or the use of the Hilbert transform on wide-band signals. As the model is probabilistic, it also provides a score of the model “goodness of fit” via the likelihood, enabling easy and legitimate model selection and parameter comparison; this data-driven feature is unique to our model-based approach. Using three datasets obtained with invasive neurophysiological recordings in humans and rodents, we demonstrate that these models are able to replicate previous results obtained with other metrics, but also reveal new insights such as the influence of the amplitude of the slow oscillation. Using simulations, we demonstrate that our parametric method can reveal neural couplings with shorter signals than non-parametric methods. We also show how the likelihood can be used to find optimal filtering parameters, suggesting new properties on the spectrum of the driving signal, but also to estimate the optimal delay between the coupled signals, enabling a directionality estimation in the coupling. PMID:29227989

Stable explicit coupling of the Yee scheme with a linear current model in fluctuating magnetized plasmas

International Nuclear Information System (INIS)

Silva, Filipe da; Pinto, Martin Campos; Després, Bruno; Heuraux, Stéphane

2015-01-01

This work analyzes the stability of the Yee scheme for non-stationary Maxwell's equations coupled with a linear current model with density fluctuations. We show that the usual procedure may yield unstable scheme for physical situations that correspond to strongly magnetized plasmas in X-mode (TE) polarization. We propose to use first order clustered discretization of the vectorial product that gives back a stable coupling. We validate the schemes on some test cases representative of direct numerical simulations of X-mode in a magnetic fusion plasma including turbulence

A less-constrained (2,0) super-Yang-Mills model: the coupling to non-linear σ-models

International Nuclear Information System (INIS)

Almeida, C.A.S.; Doria, R.M.

1990-01-01

Considering a class of (2,0) super-Yang-Mills multiplets characterised by the appearance of a pair of independent gauge potentials, we present here their coupling to non-linear σ-models in (2,0)-superspace. Contrary to the case of the coupling to (2,0) matter superfields, the extra gauge potential present in the Yang-Mills sector does not decouple from the theory in the case one gauges isometry groups of σ-models. (author)

Solution of Large Systems of Linear Equations with Quadratic or Non-Quadratic Matrices and Deconvoiution of Spectra

Energy Technology Data Exchange (ETDEWEB)

Nygaard, K

1967-12-15

The numerical deconvolution of spectra is equivalent to the solution of a (large) system of linear equations with a matrix which is not necessarily a square matrix. The demand that the square sum of the residual errors shall be minimum is not in general sufficient to ensure a unique or 'sound' solution. Therefore other demands which may include the demand for minimum square errors are introduced which lead to 'sound' and 'non-oscillatory' solutions irrespective of the shape of the original matrix and of the determinant of the matrix of the normal equations.

A Lie-admissible method of integration of Fokker-Planck equations with non-linear coefficients (exact and numerical solutions)

International Nuclear Information System (INIS)

Fronteau, J.; Combis, P.

1984-08-01

A Lagrangian method is introduced for the integration of non-linear Fokker-Planck equations. Examples of exact solutions obtained in this way are given, and also the explicit scheme used for the computation of numerical solutions. The method is, in addition, shown to be of a Lie-admissible type

Positive solution of non-square fully Fuzzy linear system of equation in general form using least square method

Directory of Open Access Journals (Sweden)

Reza Ezzati

2014-08-01

Full Text Available In this paper, we propose the least square method for computing the positive solution of a non-square fully fuzzy linear system. To this end, we use Kaffman' arithmetic operations on fuzzy numbers \\cite{17}. Here, considered existence of exact solution using pseudoinverse, if they are not satisfy in positive solution condition, we will compute fuzzy vector core and then we will obtain right and left spreads of positive fuzzy vector by introducing constrained least squares problem. Using our proposed method, non-square fully fuzzy linear system of equations always has a solution. Finally, we illustrate the efficiency of proposed method by solving some numerical examples.

A 2+1 non-isospectral discrete integrable system and its discrete integrable coupling system

International Nuclear Information System (INIS)

Yu Fajun; Zhang Hongqing

2006-01-01

In this Letter by considering a (2+1)-dimensional discrete non-isospectral linear problem, a new (2+1)-dimensional integrable lattice hierarchy is constructed. It shows that generalization of the Blaszak-Marciniak lattice hierarchy can be obtained as a reduction. Then an extended algebraic system X-bar of X is presented, from which the integrable coupling system of the (2+1)-dimensional discrete non-isospectral Blaszak-Marciniak lattice equations are obtained

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.

Nonoscillation of half-linear dynamic equations

Czech Academy of Sciences Publication Activity Database

Matucci, S.; Řehák, Pavel

2010-01-01

Roč. 60, č. 5 (2010), s. 1421-1429 ISSN 0898-1221 R&D Projects: GA AV ČR KJB100190701 Grant - others:GA ČR(CZ) GA201/07/0145 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear dynamic equation * time scale * (non)oscillation * Riccati technique Subject RIV: BA - General Mathematics Impact factor: 1.472, year: 2010 http://www.sciencedirect.com/science/article/pii/S0898122110004384

New conditions on synchronization of networks of linearly coupled dynamical systems with non-Lipschitz right-hand sides.

Science.gov (United States)

Liu, Bo; Lu, Wenlian; Chen, Tianping

2012-01-01

In this paper, we study synchronization of networks of linearly coupled dynamical systems. The node dynamics of the network can be very general, which may not satisfy the QUAD condition. We derive sufficient conditions for synchronization, which can be regarded as extensions of previous results. These results can be employed to networks of coupled systems, of which, in particular, the node dynamics have non-Lipschitz or even discontinuous right-hand sides. We also give several corollaries where the synchronization of some specific non-QUAD systems can be deduced. As an application, we propose a scheme to realize synchronization of coupled switching systems via coupling the signals which drive the switchings. Examples with numerical simulations are also provided to illustrate the theoretical results. Copyright © 2011 Elsevier Ltd. All rights reserved.

Control of Coherent Instabilities by Linear Coupling

CERN Document Server

Cappi, R; Möhl, D

2001-01-01

One of the main challenges in the design of high-energy colliders is the very high luminosity necessary to provide significant event rates. This imposes strong constraints to achieve and preserve beams of high brightness, i.e. intensity to emittance ratio, all along the injector chain. Amongst the phenomena that can blow up and even destroy the beam are transverse coherent instabilities. Two methods are widely used to damp these instabilities. The first one is Landau damping by non-linearities. The second consists in using an electronic feedback system. However, non-linearities are harmful to single-particle motion due to resonance phenomena, and powerful wideband feedback systems are expensive. It is shown in this paper that linear coupling is a further method that can be used to damp transverse coherent instabilities. The theory of collective motion is outlined, including the coupling of instability rise and damping rates, chromaticity and Landau damping. Experimental results obtained at the CERN PS are rep...

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.

Isomorphism of Intransitive Linear Lie Equations

Directory of Open Access Journals (Sweden)

Jose Miguel Martins Veloso

2009-11-01

Full Text Available We show that formal isomorphism of intransitive linear Lie equations along transversal to the orbits can be extended to neighborhoods of these transversal. In analytic cases, the word formal is dropped from theorems. Also, we associate an intransitive Lie algebra with each intransitive linear Lie equation, and from the intransitive Lie algebra we recover the linear Lie equation, unless of formal isomorphism. The intransitive Lie algebra gives the structure functions introduced by É. Cartan.

Nonlinear coupled equations for electrochemical cells as developed by the general equation for nonequilibrium reversible-irreversible coupling.

Science.gov (United States)

Bedeaux, Dick; Kjelstrup, Signe; Öttinger, Hans Christian

2014-09-28

We show how the Butler-Volmer and Nernst equations, as well as Peltier effects, are contained in the general equation for nonequilibrium reversible and irreversible coupling, GENERIC, with a unique definition of the overpotential. Linear flux-force relations are used to describe the transport in the homogeneous parts of the electrochemical system. For the electrode interface, we choose nonlinear flux-force relationships. We give the general thermodynamic basis for an example cell with oxygen electrodes and electrolyte from the solid oxide fuel cell. In the example cell, there are two activated chemical steps coupled also to thermal driving forces at the surface. The equilibrium exchange current density obtains contributions from both rate-limiting steps. The measured overpotential is identified at constant temperature and stationary states, in terms of the difference in electrochemical potential of products and reactants. Away from these conditions, new terms appear. The accompanying energy flux out of the surface, as well as the heat generation at the surface are formulated, adding to the general thermodynamic basis.

Nonlinear coupled equations for electrochemical cells as developed by the general equation for nonequilibrium reversible-irreversible coupling

Science.gov (United States)

Bedeaux, Dick; Kjelstrup, Signe; Öttinger, Hans Christian

2014-09-01

We show how the Butler-Volmer and Nernst equations, as well as Peltier effects, are contained in the general equation for nonequilibrium reversible and irreversible coupling, GENERIC, with a unique definition of the overpotential. Linear flux-force relations are used to describe the transport in the homogeneous parts of the electrochemical system. For the electrode interface, we choose nonlinear flux-force relationships. We give the general thermodynamic basis for an example cell with oxygen electrodes and electrolyte from the solid oxide fuel cell. In the example cell, there are two activated chemical steps coupled also to thermal driving forces at the surface. The equilibrium exchange current density obtains contributions from both rate-limiting steps. The measured overpotential is identified at constant temperature and stationary states, in terms of the difference in electrochemical potential of products and reactants. Away from these conditions, new terms appear. The accompanying energy flux out of the surface, as well as the heat generation at the surface are formulated, adding to the general thermodynamic basis.

Non-linear soil-structure interaction

International Nuclear Information System (INIS)

Wolf, J.P.

1984-01-01

The basic equation of motion to analyse the interaction of a non-linear structure and an irregular soil with the linear unbounded soil is formulated in the time domain. The contribution of the unbounded soil involves convolution integrals of the dynamic-stiffness coefficients in the time domain and the corresponding motions. As another possibility, a flexibility formulation fot the contribution of the unbounded soil using the dynamic-flexibility coefficients in the time domain, together with the direct-stiffness method for the structure and the irregular soil can be applied. As an example of a non-linear soil-structure-interaction analysis, the partial uplift of the basemat of a structure is examined. (Author) [pt

Performance prediction of gas turbines by solving a system of non-linear equations

Energy Technology Data Exchange (ETDEWEB)

Kaikko, J

1998-09-01

This study presents a novel method for implementing the performance prediction of gas turbines from the component models. It is based on solving the non-linear set of equations that corresponds to the process equations, and the mass and energy balances for the engine. General models have been presented for determining the steady state operation of single components. Single and multiple shad arrangements have been examined with consideration also being given to heat regeneration and intercooling. Emphasis has been placed upon axial gas turbines of an industrial scale. Applying the models requires no information of the structural dimensions of the gas turbines. On comparison with the commonly applied component matching procedures, this method incorporates several advantages. The application of the models for providing results is facilitated as less attention needs to be paid to calculation sequences and routines. Solving the set of equations is based on zeroing co-ordinate functions that are directly derived from the modelling equations. Therefore, controlling the accuracy of the results is easy. This method gives more freedom for the selection of the modelling parameters since, unlike for the matching procedures, exchanging these criteria does not itself affect the algorithms. Implicit relationships between the variables are of no significance, thus increasing the freedom for the modelling equations as well. The mathematical models developed in this thesis will provide facilities to optimise the operation of any major gas turbine configuration with respect to the desired process parameters. The computational methods used in this study may also be adapted to any other modelling problems arising in industry. (orig.) 36 refs.

EDITORIAL: Non-linear and non-Gaussian cosmological perturbations Non-linear and non-Gaussian cosmological perturbations

Science.gov (United States)

Sasaki, Misao; Wands, David

2010-06-01

In recent years there has been a resurgence of interest in the study of non-linear perturbations of cosmological models. This has been the result of both theoretical developments and observational advances. New theoretical challenges arise at second and higher order due to mode coupling and the need to develop new gauge-invariant variables beyond first order. In particular, non-linear interactions lead to deviations from a Gaussian distribution of primordial perturbations even if initial vacuum fluctuations are exactly Gaussian. These non-Gaussianities provide an important probe of models for the origin of structure in the very early universe. We now have a detailed picture of the primordial distribution of matter from surveys of the cosmic microwave background, notably NASA's WMAP satellite. The situation will continue to improve with future data from the ESA Planck satellite launched in 2009. To fully exploit these data cosmologists need to extend non-linear cosmological perturbation theory beyond the linear theory that has previously been sufficient on cosmological scales. Another recent development has been the realization that large-scale structure, revealed in high-redshift galaxy surveys, could also be sensitive to non-linearities in the primordial curvature perturbation. This focus section brings together a collection of invited papers which explore several topical issues in this subject. We hope it will be of interest to theoretical physicists and astrophysicists alike interested in understanding and interpreting recent developments in cosmological perturbation theory and models of the early universe. Of course it is only an incomplete snapshot of a rapidly developing field and we hope the reader will be inspired to read further work on the subject and, perhaps, fill in some of the missing pieces. This focus section is dedicated to the memory of Lev Kofman (1957-2009), an enthusiastic pioneer of inflationary cosmology and non-Gaussian perturbations.

Existence of entire solutions of some non-linear differential-difference equations.

Science.gov (United States)

Chen, Minfeng; Gao, Zongsheng; Du, Yunfei

2017-01-01

In this paper, we investigate the admissible entire solutions of finite order of the differential-difference equations [Formula: see text] and [Formula: see text], where [Formula: see text], [Formula: see text] are two non-zero polynomials, [Formula: see text] is a polynomial and [Formula: see text]. In addition, we investigate the non-existence of entire solutions of finite order of the differential-difference equation [Formula: see text], where [Formula: see text], [Formula: see text] are two non-constant polynomials, [Formula: see text], m , n are positive integers and satisfy [Formula: see text] except for [Formula: see text], [Formula: see text].

Non-linear finite element analysis in structural mechanics

CERN Document Server

Rust, Wilhelm

2015-01-01

This monograph describes the numerical analysis of non-linearities in structural mechanics, i.e. large rotations, large strain (geometric non-linearities), non-linear material behaviour, in particular elasto-plasticity as well as time-dependent behaviour, and contact. Based on that, the book treats stability problems and limit-load analyses, as well as non-linear equations of a large number of variables. Moreover, the author presents a wide range of problem sets and their solutions. The target audience primarily comprises advanced undergraduate and graduate students of mechanical and civil engineering, but the book may also be beneficial for practising engineers in industry.

Diagonalizing quadratic bosonic operators by non-autonomous flow equations

CERN Document Server

Bach, Volker

2016-01-01

The authors study a non-autonomous, non-linear evolution equation on the space of operators on a complex Hilbert space. They specify assumptions that ensure the global existence of its solutions and allow them to derive its asymptotics at temporal infinity. They demonstrate that these assumptions are optimal in a suitable sense and more general than those used before. The evolution equation derives from the Brocketâe"Wegner flow that was proposed to diagonalize matrices and operators by a strongly continuous unitary flow. In fact, the solution of the non-linear flow equation leads to a diagonalization of Hamiltonian operators in boson quantum field theory which are quadratic in the field.

Linear and non-linear stability analysis for finite difference discretizations of high-order Boussinesq equations

DEFF Research Database (Denmark)

Fuhrman, David R.; Bingham, Harry B.; Madsen, Per A.

2004-01-01

of rotational and irrotational formulations in two horizontal dimensions provides evidence that the irrotational formulation has significantly better stability properties when the deep-water non-linearity is high, particularly on refined grids. Computation of matrix pseudospectra shows that the system is only...... insight into the numerical behaviour of this rather complicated system of non-linear PDEs....

Localized chaoticity in two linearly coupled inverted double-well ...

African Journals Online (AJOL)

Two linearly coupled inverted double-well oscillators for a fixed energy and varying coupling strength were studied. The dynamics yielded a chaotic system in which the Poincare surface was characterised by two non-mixing regions, one of regular motion and the other region that became chaotic as the coupling increased.

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

Chaos and bifurcation of a flexible rotor supported by porous squeeze couple stress fluid film journal bearings with non-linear suspension

International Nuclear Information System (INIS)

Chang-Jian, C.-W.; Chen, C.-K.

2008-01-01

This study presents a dynamic analysis of a flexible rotor supported by two porous squeeze couple stress fluid film journal bearings with non-linear suspension. The dynamics of the rotor center and bearing center are studied. The analysis of the rotor-bearing system is investigated under the assumptions of non-Newtonian fluid and a short bearing approximation. The spatial displacements in the horizontal and vertical directions are considered for various non-dimensional speed ratios. The dynamic equations are solved using the Runge-Kutta method. The analysis methods employed in this study is inclusive of the dynamic trajectories of the rotor center and bearing center, power spectra, Poincare maps and bifurcation diagrams. The maximum Lyapunov exponent analysis is also used to identify the onset of chaotic motion. The numerical results show that the stability of the system varies with the non-dimensional speed ratios, the non-dimensional parameter l* and the permeability. The modeling results thus obtained by using the method proposed in this paper can be employed to predict the stability of the rotor-bearing system and the undesirable behavior of the rotor and bearing center can be avoided

Students’ difficulties in solving linear equation problems

Science.gov (United States)

Wati, S.; Fitriana, L.; Mardiyana

2018-03-01

A linear equation is an algebra material that exists in junior high school to university. It is a very important material for students in order to learn more advanced mathematics topics. Therefore, linear equation material is essential to be mastered. However, the result of 2016 national examination in Indonesia showed that students’ achievement in solving linear equation problem was low. This fact became a background to investigate students’ difficulties in solving linear equation problems. This study used qualitative descriptive method. An individual written test on linear equation tasks was administered, followed by interviews. Twenty-one sample students of grade VIII of SMPIT Insan Kamil Karanganyar did the written test, and 6 of them were interviewed afterward. The result showed that students with high mathematics achievement donot have difficulties, students with medium mathematics achievement have factual difficulties, and students with low mathematics achievement have factual, conceptual, operational, and principle difficulties. Based on the result there is a need of meaningfulness teaching strategy to help students to overcome difficulties in solving linear equation problems.

Non-linear Matter Spectra in Coupled Quintessence

CERN Document Server

Saracco, F; Tetradis, N; Pettorino, V; Robbers, G

2010-01-01

We consider cosmologies in which a dark-energy scalar field interacts with cold dark matter. The growth of perturbations is followed beyond the linear level by means of the time-renormalization-group method, which is extended to describe a multi-component matter sector. Even in the absence of the extra interaction, a scale-dependent bias is generated as a consequence of the different initial conditions for baryons and dark matter after decoupling. The effect is greatly enhanced by the extra coupling and can be at the percent level in the range of scales of baryonic acoustic oscillations. We compare our results with N-body simulations, finding very good agreement.

Selective Linear or Quadratic Optomechanical Coupling via Measurement

Directory of Open Access Journals (Sweden)

Michael R. Vanner

2011-11-01

Full Text Available The ability to engineer both linear and nonlinear coupling with a mechanical resonator is an important goal for the preparation and investigation of macroscopic mechanical quantum behavior. In this work, a measurement based scheme is presented where linear or square mechanical-displacement coupling can be achieved using the optomechanical interaction that is linearly proportional to the mechanical position. The resulting square-displacement measurement strength is compared to that attainable in the dispersive case that has a direct interaction with the mechanical-displacement squared. An experimental protocol and parameter set are discussed for the generation and observation of non-Gaussian states of motion of the mechanical element.

Effective electrodiffusion equation for non-uniform nanochannels.

Science.gov (United States)

Marini Bettolo Marconi, Umberto; Melchionna, Simone; Pagonabarraga, Ignacio

2013-06-28

We derive a one-dimensional formulation of the Planck-Nernst-Poisson equation to describe the dynamics of a symmetric binary electrolyte in channels whose section is nanometric and varies along the axial direction. The approach is in the spirit of the Fick-Jacobs diffusion equation and leads to a system of coupled equations for the partial densities which depends on the charge sitting at the walls in a non-trivial fashion. We consider two kinds of non-uniformities, those due to the spatial variation of charge distribution and those due to the shape variation of the pore and report one- and three-dimensional solutions of the electrokinetic equations.

Numerical simulation of a nonlinear coupled fluid-structure problem. Application to the design of naval nuclear propulsion structures; Modelisation et simulation numerique d'un probleme couple fluide/structure non lineaire: application au dimensionnement de structures nucleaires de propulsion navale

Energy Technology Data Exchange (ETDEWEB)

Sigrist, J.F

2004-11-15

The present work deals with the numerical simulation of a coupled fluid/structure problem with fluid free surface. A generic coupled fluid/structure system is defined, on which a linear problem (modal analysis) and a non-linear problem (temporal analysis) are stated. In the linear case, a strong coupled method is used. It is based on a finite element approach of the structure problem and a finite or a boundary element approach of the fluid problem. The coupled problem is formulated in terms of pressure and displacement, leading to a non-symmetric problem which is solved with an appropriate algorithm. In the non-linear case, the structure problem is described with non-linear equations of motion, whereas the fluid problem is modeled with the Stokes equations. The numerical resolution of the coupled problem is based on a weak coupling procedure. The fluid problem is solved with a finite volume technique, using a moving mesh technique to adjust the structure motion, a VOF method for the description of the free surface and the PISO algorithm for the time integration. The structure problem is solved with a finite element technique, using an explicit/implicit time integration algorithm. A procedure is developed in order to handle the coupling in space (fluid forces and structure displacement exchanges between fluid and structure mesh, fluid re-meshing) and in time (staggered explicit algorithm, dynamic filtering of numerical oscillations). The non linear coupled problem is solved using a CFD code, whose use for FSI problem is validated with a benchmark presented in this work. A comparison is proposed between numerical results and analytical solution for two elementary fluid problems. The validation process can be applied for any CFD numerical code. A numerical study is then proposed on the generic coupled case in order to describe the fluid/structure interaction phenomenon (added mass, displaced mass, mode coupling, influence of structural non-linearity). An industrial

CPDES2: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in two dimensions

Science.gov (United States)

Anderson, D. V.; Koniges, A. E.; Shumaker, D. E.

1988-11-01

Many physical problems require the solution of coupled partial differential equations on two-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils and allows for general couplings between all of the component PDE's and it automatically generates the matrix structures needed to perform the algorithm. The resulting sparse matrix equation is solved by either the preconditioned conjugate gradient (CG) method or by the preconditioned biconjugate gradient (BCG) algorithm. An arbitrary number of component equations are permitted only limited by available memory. In the sub-band representation used, we generate an algorithm that is written compactly in terms of indirect indices which is vectorizable on some of the newer scientific computers.

CPDES3: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions

Science.gov (United States)

Anderson, D. V.; Koniges, A. E.; Shumaker, D. E.

1988-11-01

Many physical problems require the solution of coupled partial differential equations on three-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils and allows for general couplings between all of the component PDE's and it automatically generates the matrix structures needed to perform the algorithm. The resulting sparse matrix equation is solved by either the preconditioned conjugate gradient (CG) method or by the preconditioned biconjugate gradient (BCG) algorithm. An arbitrary number of component equations are permitted only limited by available memory. In the sub-band representation used, we generate an algorithm that is written compactly in terms of indirect induces which is vectorizable on some of the newer scientific computers.

Correct Linearization of Einstein's Equations

Directory of Open Access Journals (Sweden)

Rabounski D.

2006-06-01

Full Text Available Regularly Einstein's equations can be reduced to a wave form (linearly dependent from the second derivatives of the space metric in the absence of gravitation, the space rotation and Christoffel's symbols. As shown here, the origin of the problem is that one uses the general covariant theory of measurement. Here the wave form of Einstein's equations is obtained in the terms of Zelmanov's chronometric invariants (physically observable projections on the observer's time line and spatial section. The obtained equations depend on solely the second derivatives even if gravitation, the space rotation and Christoffel's symbols. The correct linearization proves: the Einstein equations are completely compatible with weak waves of the metric.

Diffusive limits for linear transport equations

International Nuclear Information System (INIS)

Pomraning, G.C.

1992-01-01

The authors show that the Hibert and Chapman-Enskog asymptotic treatments that reduce the nonlinear Boltzmann equation to the Euler and Navier-Stokes fluid equations have analogs in linear transport theory. In this linear setting, these fluid limits are described by diffusion equations, involving familiar and less familiar diffusion coefficients. Because of the linearity extant, one can carry out explicitly the initial and boundary layer analyses required to obtain asymptotically consistent initial and boundary conditions for the diffusion equations. In particular, the effects of boundary curvature and boundary condition variation along the surface can be included in the boundary layer analysis. A brief review of heuristic (nonasymptotic) diffusion description derivations is also included in our discussion

Spectral theories for linear differential equations

International Nuclear Information System (INIS)

Sell, G.R.

1976-01-01

The use of spectral analysis in the study of linear differential equations with constant coefficients is not only a fundamental technique but also leads to far-reaching consequences in describing the qualitative behaviour of the solutions. The spectral analysis, via the Jordan canonical form, will not only lead to a representation theorem for a basis of solutions, but will also give a rather precise statement of the (exponential) growth rates of various solutions. Various attempts have been made to extend this analysis to linear differential equations with time-varying coefficients. The most complete such extensions is the Floquet theory for equations with periodic coefficients. For time-varying linear differential equations with aperiodic coefficients several authors have attempted to ''extend'' the Foquet theory. The precise meaning of such an extension is itself a problem, and we present here several attempts in this direction that are related to the general problem of extending the spectral analysis of equations with constant coefficients. The main purpose of this paper is to introduce some problems of current research. The primary problem we shall examine occurs in the context of linear differential equations with almost periodic coefficients. We call it ''the Floquet problem''. (author)

Comparison of linear and non-linear models for the adsorption of fluoride onto geo-material: limonite.

Science.gov (United States)

Sahin, Rubina; Tapadia, Kavita

2015-01-01

The three widely used isotherms Langmuir, Freundlich and Temkin were examined in an experiment using fluoride (F⁻) ion adsorption on a geo-material (limonite) at four different temperatures by linear and non-linear models. Comparison of linear and non-linear regression models were given in selecting the optimum isotherm for the experimental results. The coefficient of determination, r², was used to select the best theoretical isotherm. The four Langmuir linear equations (1, 2, 3, and 4) are discussed. Langmuir isotherm parameters obtained from the four Langmuir linear equations using the linear model differed but they were the same when using the nonlinear model. Langmuir-2 isotherm is one of the linear forms, and it had the highest coefficient of determination (r² = 0.99) compared to the other Langmuir linear equations (1, 3 and 4) in linear form, whereas, for non-linear, Langmuir-4 fitted best among all the isotherms because it had the highest coefficient of determination (r² = 0.99). The results showed that the non-linear model may be a better way to obtain the parameters. In the present work, the thermodynamic parameters show that the absorption of fluoride onto limonite is both spontaneous (ΔG 0). Scanning electron microscope and X-ray diffraction images also confirm the adsorption of F⁻ ion onto limonite. The isotherm and kinetic study reveals that limonite can be used as an adsorbent for fluoride removal. In future we can develop new technology for fluoride removal in large scale by using limonite which is cost-effective, eco-friendly and is easily available in the study area.

Canonical structure of evolution equations with non-linear ...

Indian Academy of Sciences (India)

The dispersion produced is compensated by non-linear effects resulting in the formation of exponentially localized .... determining the values of Lagrange's multipliers αis. We postulate that a slightly .... c3 «w2x -v. (36). To include the effect of the secondary constraint c3 in the total Hamiltonian H we modify. (33) as. 104.

A non-linear optimal Discontinuous Petrov-Galerkin method for stabilising the solution of the transport equation

International Nuclear Information System (INIS)

Merton, S. R.; Smedley-Stevenson, R. P.; Pain, C. C.; Buchan, A. G.; Eaton, M. D.

2009-01-01

This paper describes a new Non-Linear Discontinuous Petrov-Galerkin (NDPG) method and application to the one-speed Boltzmann Transport Equation (BTE) for space-time problems. The purpose of the method is to remove unwanted oscillations in the transport solution which occur in the vicinity of sharp flux gradients, while improving computational efficiency and numerical accuracy. This is achieved by applying artificial dissipation in the solution gradient direction, internal to an element using a novel finite element (FE) Riemann approach. The amount of dissipation added acts internal to each element. This is done using a gradient-informed scaling of the advection velocities in the stabilisation term. This makes the method in its most general form non-linear. The method is designed to be independent of angular expansion framework. This is demonstrated for the both discrete ordinates (S N ) and spherical harmonics (P N ) descriptions of the angular variable. Results show the scheme performs consistently well in demanding time dependent and multi-dimensional radiation transport problems. (authors)

approximate controllability of a non-autonomous differential equation

Indian Academy of Sciences (India)

53

for a non-autonomous functional differential equation using the theory of linear ... approximate controllability of various functional differential equations in abstract ...... the operator A(t) and into the requirement that x(t) ∈ D(A) for all t ≥ 0.

CPDS3, Coupled 3-D Partial Differential Equation Solution

International Nuclear Information System (INIS)

Anderson, D.V.; Koniges, A.E.; Shumaker, D.E.

1992-01-01

1 - Description of program or function: CPDES3 solves the linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils, permits general couplings between all of the component PDE's, and automatically generates the matrix structures needed to perform the algorithm. 2 - Method of solution: The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned bi-conjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used. 3 - Restrictions on the complexity of the problem: The discretization of the coupled three-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective

CPDES2, Coupled 2-D Partial Differential Equation Solution

International Nuclear Information System (INIS)

Anderson, D.V.; Koniges, A.E.; Shumaker, D.E.

1992-01-01

1 - Description of program or function: CPDES2 solves the linear asymmetric equations arising from coupled partial differential equations in two dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils, permits general coupling between all of the component PDE's, and automatically generates the matrix structures needed to perform the algorithm. 2 - Method of solution: The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned bi-conjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used. 3 - Restrictions on the complexity of the problem: The discretization of the coupled two-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective

KAM for the non-linear Schroedinger equation a short presentation

CERN Document Server

Eliasson, H L

2006-01-01

We consider the $d$-dimensional nonlinear Schr\\"o\\-dinger equation under periodic boundary conditions:-i\\dot u=\\Delta u+V(x)*u+\\ep \\frac{\\p F}{\\p \\bar u}(x,u,\\bar u) ;\\quad u=u(t,x),\\;x\\in\\T^dwhere $V(x)=\\sum \\hat V(a)e^{i\\sc{a,x}}$ is an analytic function with $\\hat V$ real and $F$ is a real analytic function in $\\Re u$, $\\Im u$ and $x$. (This equation is a popular model for the `real' NLS equation, where instead of the convolution term $V*u$ we have the potential term $Vu$.) For $\\ep=0$ the equation is linear and has time--quasi-periodic solutions $u$,u(t,x)=\\sum_{s\\in \\AA}\\hat u_0(a)e^{i(|a|^2+\\hat V(a))t}e^{i\\sc{a,x}}, \\quad 0<|\\hat u_0(a)|\\le1,where $\\AA$ is any finite subset of $\\Z^d$. We shall treat $\\omega_a=|a|^2+\\hat V(a)$, $a\\in\\AA$, as free parameters in some domain $U\\subset\\R^{\\AA}$. This is a Hamiltonian system in infinite degrees of freedom, degenerate but with external parameters, and we shall describe a KAM-theory which, in particular, will have the following consequence: \\smallskip {\\it ...

An Etude in non-linear Dyson-Schwinger Equations

International Nuclear Information System (INIS)

Kreimer, Dirk; Yeats, Karen

2006-01-01

We show how to use the Hopf algebra structure of quantum field theory to derive nonperturbative results for the short-distance singular sector of a renormalizable quantum field theory in a simple but generic example. We discuss renormalized Green functions G R (α,L) in such circumstances which depend on a single scale L=lnq 2 /μ 2 and start from an expansion in the scale G R (α,L)=1+-bar k γ k (α)L k . We derive recursion relations between the γ k which make full use of the renormalization group. We then show how to determine the Green function by the use of a Mellin transform on suitable integral kernels. We exhibit our approach in an example for which we find a functional equation relating weak and strong coupling expansions

Solving the Coupled System Improves Computational Efficiency of the Bidomain Equations

KAUST Repository

Southern, J.A.

2009-10-01

The bidomain equations are frequently used to model the propagation of cardiac action potentials across cardiac tissue. At the whole organ level, the size of the computational mesh required makes their solution a significant computational challenge. As the accuracy of the numerical solution cannot be compromised, efficiency of the solution technique is important to ensure that the results of the simulation can be obtained in a reasonable time while still encapsulating the complexities of the system. In an attempt to increase efficiency of the solver, the bidomain equations are often decoupled into one parabolic equation that is computationally very cheap to solve and an elliptic equation that is much more expensive to solve. In this study, the performance of this uncoupled solution method is compared with an alternative strategy in which the bidomain equations are solved as a coupled system. This seems counterintuitive as the alternative method requires the solution of a much larger linear system at each time step. However, in tests on two 3-D rabbit ventricle benchmarks, it is shown that the coupled method is up to 80% faster than the conventional uncoupled method-and that parallel performance is better for the larger coupled problem.

Solving the Coupled System Improves Computational Efficiency of the Bidomain Equations

KAUST Repository

Southern, J.A.; Plank, G.; Vigmond, E.J.; Whiteley, J.P.

2009-01-01

The bidomain equations are frequently used to model the propagation of cardiac action potentials across cardiac tissue. At the whole organ level, the size of the computational mesh required makes their solution a significant computational challenge. As the accuracy of the numerical solution cannot be compromised, efficiency of the solution technique is important to ensure that the results of the simulation can be obtained in a reasonable time while still encapsulating the complexities of the system. In an attempt to increase efficiency of the solver, the bidomain equations are often decoupled into one parabolic equation that is computationally very cheap to solve and an elliptic equation that is much more expensive to solve. In this study, the performance of this uncoupled solution method is compared with an alternative strategy in which the bidomain equations are solved as a coupled system. This seems counterintuitive as the alternative method requires the solution of a much larger linear system at each time step. However, in tests on two 3-D rabbit ventricle benchmarks, it is shown that the coupled method is up to 80% faster than the conventional uncoupled method-and that parallel performance is better for the larger coupled problem.

Lie algebras and linear differential equations.

Science.gov (United States)

Brockett, R. W.; Rahimi, A.

1972-01-01

Certain symmetry properties possessed by the solutions of linear differential equations are examined. For this purpose, some basic ideas from the theory of finite dimensional linear systems are used together with the work of Wei and Norman on the use of Lie algebraic methods in differential equation theory.

Non-linear Structures in the Non-critical NSR String

International Nuclear Information System (INIS)

Hamada, K.; Ishikawa, H.

1996-01-01

We investigate the Ward identities of the W ∞ symmetry in the super-Liouville theory coupled to the super-conformal matter of central charge c M =1-2(p-q) 2 /pq. The theory is classified into two chiralities. For the positive chirality, all gravitationally dressed scaling operators are generated from the q-1 repeatedly. After fixing the normalizations of the dressed scaling operators, we find that the Ward identities are expressed in the form of the usual W q algebra constraints as in the bosonic case: W n (k+1) τ=0, (k=1,..,q-1; nεZ≥1- k), where the equations for even and odd n come from the currents in the NS- and the R-sector respectively. The non-linear terms come from the anomalous contributions at the boundaries of moduli space. The negative chirality is defined by interchanging the roles of p and q. Then we get the W p algebra constraints. (orig.)

Heterotic sigma models and non-linear strings

International Nuclear Information System (INIS)

Hull, C.M.

1986-01-01

The two-dimensional supersymmetric non-linear sigma models are examined with respect to the heterotic string. The paper was presented at the workshop on :Supersymmetry and its applications', Cambridge, United Kingdom, 1985. The non-linear sigma model with Wess-Zumino-type term, the coupling of the fermionic superfields to the sigma model, super-conformal invariance, and the supersymmetric string, are all discussed. (U.K.)

Linear q-nonuniform difference equations

International Nuclear Information System (INIS)

Bangerezako, Gaspard

2010-01-01

We introduce basic concepts of q-nonuniform differentiation and integration and study linear q-nonuniform difference equations and systems, as well as their application in q-nonuniform difference linear control systems. (author)

Dark matter as a non-linear effect of gravitation

International Nuclear Information System (INIS)

Maia, M.D.; Capistrano, A.J.S.

2006-01-01

The rotation curves of stars in disk galaxies are calculated with the Newtonian law of motion applied to a scalar potential derived from the geodesic equation, only, under the slow motion condition, the so-called Nearly Newtonian Gravity (NNG). A nearly Newtonian gravitational potential, Φ NN = -1/2 c 2 (1+g 44 ), is obtained, characterized by an exact solution of Einsteins equations, with the non-linear effects present in the component g 44 . This gravitational field lies somewhere between General Relativity and Newtonian Gravity. Therefore, Einsteins equations and the equivalence principle are preserved, but the general covariance is broken. The resulting curves are remarkably close to the observed rotation curves in spiral galaxies, suggesting that a substantial component of dark matter may be explained by the non-linearity of Einsteins equations. (author)

Non-Linear Fibres for Widely Tunable Femtosecond Fibre Lasers

DEFF Research Database (Denmark)

Pedersen, Martin Erland Vestergaard

and numerically. For the intermodal four-wave mixing experiment an alternative version of the Generalised Non-Linear Schrödinger Equation is derived, which includes the correct dispersion of the transverse field. It is observed that the alternative version of the Generalised Non-Linear Schrödinger Equation......, as opposed to the commonly used version, is able to reproduce the intermodal four-wave mixing experiment. The relation between the intramodal self-phase modulation and the intramodal Raman effect is determined from experimental measurements on a number of step-index fibres. The Raman fraction is found...

Exact non-Markovian master equations for multiple qubit systems: Quantum-trajectory approach

Science.gov (United States)

Chen, Yusui; You, J. Q.; Yu, Ting

2014-11-01

A wide class of exact master equations for a multiple qubit system can be explicitly constructed by using the corresponding exact non-Markovian quantum-state diffusion equations. These exact master equations arise naturally from the quantum decoherence dynamics of qubit system as a quantum memory coupled to a collective colored noisy source. The exact master equations are also important in optimal quantum control, quantum dissipation, and quantum thermodynamics. In this paper, we show that the exact non-Markovian master equation for a dissipative N -qubit system can be derived explicitly from the statistical average of the corresponding non-Markovian quantum trajectories. We illustrated our general formulation by an explicit construction of a three-qubit system coupled to a non-Markovian bosonic environment. This multiple qubit master equation offers an accurate time evolution of quantum systems in various domains, and paves the way to investigate the memory effect of an open system in a non-Markovian regime without any approximation.

Increase of nonlinear signal distortions due to linear mode coupling in space division multiplexed systems

DEFF Research Database (Denmark)

Kutluyarov, Ruslan V.; Bagmanov, Valeriy Kh; Antonov, Vyacheslav V.

2017-01-01

This paper is focused on the analysis of linear and nonlinear mode coupling in space division multiplexed (SDM) optical communications over step-index fiber in few-mode regime. Linear mode coupling is caused by the fiber imperfections, while the nonlinear coupling is caused by the Kerr......-nonlinearities. Therefore, we use the system of generalized coupled nonlinear Schrödinger equations (GCNLSE) to describe the signal propagation. We analytically show that the presence of linear mode coupling may cause increasing of the nonlinear signal distortions. For the detailed study we solve GCNLSE numerically...... for the standard step index fiber at the wavelength of 850 nm in the basis of spatial modes with helical phase front (vortex modes) and for a special kind of few-mode fiber with enlarged core, providing propagation of five spatial modes at 1550 nm. Simulation results confirm that the linear mode coupling may lead...

Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.

Energy Technology Data Exchange (ETDEWEB)

Preston, Leiph [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-08-01

Explosions within the earth nonlinearly deform the local media, but at typical seismological observation distances, the seismic waves can be considered linear. Although nonlinear algorithms can simulate explosions in the very near field well, these codes are computationally expensive and inaccurate at propagating these signals to great distances. A linearized wave propagation code, coupled to a nonlinear code, provides an efficient mechanism to both accurately simulate the explosion itself and to propagate these signals to distant receivers. To this end we have coupled Sandia's nonlinear simulation algorithm CTH to a linearized elastic wave propagation code for 2-D axisymmetric media (axiElasti) by passing information from the nonlinear to the linear code via time-varying boundary conditions. In this report, we first develop the 2-D axisymmetric elastic wave equations in cylindrical coordinates. Next we show how we design the time-varying boundary conditions passing information from CTH to axiElasti, and finally we demonstrate the coupling code via a simple study of the elastic radius.

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.

A fast linearized conservative finite element method for the strongly coupled nonlinear fractional Schrödinger equations

Science.gov (United States)

Li, Meng; Gu, Xian-Ming; Huang, Chengming; Fei, Mingfa; Zhang, Guoyu

2018-04-01

In this paper, a fast linearized conservative finite element method is studied for solving the strongly coupled nonlinear fractional Schrödinger equations. We prove that the scheme preserves both the mass and energy, which are defined by virtue of some recursion relationships. Using the Sobolev inequalities and then employing the mathematical induction, the discrete scheme is proved to be unconditionally convergent in the sense of L2-norm and H α / 2-norm, which means that there are no any constraints on the grid ratios. Then, the prior bound of the discrete solution in L2-norm and L∞-norm are also obtained. Moreover, we propose an iterative algorithm, by which the coefficient matrix is independent of the time level, and thus it leads to Toeplitz-like linear systems that can be efficiently solved by Krylov subspace solvers with circulant preconditioners. This method can reduce the memory requirement of the proposed linearized finite element scheme from O (M2) to O (M) and the computational complexity from O (M3) to O (Mlog ⁡ M) in each iterative step, where M is the number of grid nodes. Finally, numerical results are carried out to verify the correction of the theoretical analysis, simulate the collision of two solitary waves, and show the utility of the fast numerical solution techniques.

Equation of state of strongly coupled plasma mixtures

International Nuclear Information System (INIS)

DeWitt, H.E.

1984-01-01

Thermodynamic properties of strongly coupled (high density) plasmas of mixtures of light elements have been obtained by Monte Carlo simulations. For an assumed uniform charge background the equation of state of ionic mixtures is a simple extension of the one-component plasma EOS. More realistic electron screening effects are treated in linear response theory and with an appropriate electron dielectric function. Results have been obtained for the ionic pair distribution functions, and for the electric microfield distribution

The non-linear ion trap. Part 5. Nature of non-linear resonances and resonant ion ejection

Science.gov (United States)

Franzen, J.

1994-01-01

The superposition of higher order multipole fields on the basic quadrupole field in ion traps generates a non-harmonic oscillator system for the ions. Fourier analyses of simulated secular oscillations in non-linear ion traps, therefore, not only reveal the sideband frequencies, well-known from the Mathieu theory, but additionally a commonwealth of multipole-specific overtones (or higher harmonics), and corresponding sidebands of overtones. Non-linear resonances occur when the overtone frequencies match sideband frequencies. It can be shown that in each of the resonance conditions, not just one overtone matches one sideband, instead, groups of overtones match groups of sidebands. The generation of overtones is studied by Fourier analysis of computed ion oscillations in the direction of thez axis. Even multipoles (octopole, dodecapole, etc.) generate only odd orders of higher harmonics (3, 5, etc.) of the secular frequency, explainable by the symmetry with regard to the planez = 0. In contrast, odd multipoles (hexapole, decapole, etc.) generate all orders of higher harmonics. For all multipoles, the lowest higher harmonics are found to be strongest. With multipoles of higher orders, the strength of the overtones decreases weaker with the order of the harmonics. Forz direction resonances in stationary trapping fields, the function governing the amplitude growth is investigated by computer simulations. The ejection in thez direction, as a function of timet, follows, at least in good approximation, the equation wheren is the order of multipole, andC is a constant. This equation is strictly valid for the electrically applied dipole field (n = 1), matching the secular frequency or one of its sidebands, resulting in a linear increase of the amplitude. It is valid also for the basic quadrupole field (n = 2) outside the stability area, giving an exponential increase. It is at least approximately valid for the non-linear resonances by weak superpositions of all higher odd

Basic linear partial differential equations

CERN Document Server

Treves, Francois

1975-01-01

Focusing on the archetypes of linear partial differential equations, this text for upper-level undergraduates and graduate students features most of the basic classical results. The methods, however, are decidedly nontraditional: in practically every instance, they tend toward a high level of abstraction. This approach recalls classical material to contemporary analysts in a language they can understand, as well as exploiting the field's wealth of examples as an introduction to modern theories.The four-part treatment covers the basic examples of linear partial differential equations and their

Magnetohydrodynamic flow of Carreau fluid over a convectively heated surface in the presence of non-linear radiation

Energy Technology Data Exchange (ETDEWEB)

Khan, Masood [Department of Mathematics, Quaid-i-Azam University, Islamabad 44000 (Pakistan); Hashim, E-mail: hashim_alik@yahoo.com [Department of Mathematics, Quaid-i-Azam University, Islamabad 44000 (Pakistan); Hussain, M. [Department of Sciences and Humanities, National University of Computer and Emerging Sciences, Islamabad 44000 (Pakistan); Azam, M. [Department of Mathematics, Quaid-i-Azam University, Islamabad 44000 (Pakistan)

2016-08-15

This paper presents a study of the magnetohydrodynamic (MHD) boundary layer flow of a non-Newtonian Carreau fluid over a convectively heated surface. The analysis of heat transfer is further performed in the presence of non-linear thermal radiation. The appropriate transformations are employed to bring the governing equations into dimensionless form. The numerical solutions of the partially coupled non-linear ordinary differential equations are obtained by using the Runge-Kutta Fehlberg integration scheme. The influence of non-dimensional governing parameters on the velocity, temperature, local skin friction coefficient and local Nusselt number is studied and discussed with the help of graphs and tables. Results proved that there is significant decrease in the velocity and the corresponding momentum boundary layer thickness with the growth in the magnetic parameter. However, a quite the opposite is true for the temperature and the corresponding thermal boundary layer thickness. - Highlights: • We investigated the Magnetohydrodynamic flow of Carreau constitutive fluid model. • Impact of non-linear thermal radiation is further taken into account. • Runge-Kutta Fehlberg method is employed to obtain the numerical solutions. • Fluid velocity is higher in case of hydromagnetic flow in comparison with hydrodynamic flow. • The local Nusselt number is a decreasing function of the thermal radiation parameter.

Non-linear instability analysis of the two-dimensional Navier-Stokes equation: The Taylor-Green vortex problem

Science.gov (United States)

Sengupta, Tapan K.; Sharma, Nidhi; Sengupta, Aditi

2018-05-01

An enstrophy-based non-linear instability analysis of the Navier-Stokes equation for two-dimensional (2D) flows is presented here, using the Taylor-Green vortex (TGV) problem as an example. This problem admits a time-dependent analytical solution as the base flow, whose instability is traced here. The numerical study of the evolution of the Taylor-Green vortices shows that the flow becomes turbulent, but an explanation for this transition has not been advanced so far. The deviation of the numerical solution from the analytical solution is studied here using a high accuracy compact scheme on a non-uniform grid (NUC6), with the fourth-order Runge-Kutta method. The stream function-vorticity (ψ, ω) formulation of the governing equations is solved here in a periodic square domain with four vortices at t = 0. Simulations performed at different Reynolds numbers reveal that numerical errors in computations induce a breakdown of symmetry and simultaneous fragmentation of vortices. It is shown that the actual physical instability is triggered by the growth of disturbances and is explained by the evolution of disturbance mechanical energy and enstrophy. The disturbance evolution equations have been traced by looking at (a) disturbance mechanical energy of the Navier-Stokes equation, as described in the work of Sengupta et al., "Vortex-induced instability of an incompressible wall-bounded shear layer," J. Fluid Mech. 493, 277-286 (2003), and (b) the creation of rotationality via the enstrophy transport equation in the work of Sengupta et al., "Diffusion in inhomogeneous flows: Unique equilibrium state in an internal flow," Comput. Fluids 88, 440-451 (2013).

A symbiotic approach to fluid equations and non-linear flux-driven simulations of plasma dynamics

Science.gov (United States)

Halpern, Federico

2017-10-01

The fluid framework is ubiquitous in studies of plasma transport and stability. Typical forms of the fluid equations are motivated by analytical work dating several decades ago, before computer simulations were indispensable, and can be, therefore, not optimal for numerical computation. We demonstrate a new first-principles approach to obtaining manifestly consistent, skew-symmetric fluid models, ensuring internal consistency and conservation properties even in discrete form. Mass, kinetic, and internal energy become quadratic (and always positive) invariants of the system. The model lends itself to a robust, straightforward discretization scheme with inherent non-linear stability. A simpler, drift-ordered form of the equations is obtained, and first results of their numerical implementation as a binary framework for bulk-fluid global plasma simulations are demonstrated. This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, Theory Program, under Award No. DE-FG02-95ER54309.

AAMQS: a non-linear QCD description of new HERA data at small-x

CERN Document Server

Quiroga-Arias, Paloma; Armesto, Nestor; Milhano, Jose Guilherme; Salgado, Carlos A

2011-01-01

We present a global analysis of available data on inclusive structure functions measured in electron-proton scattering at small values of Bjorken-x, including the latest data from the combined HERA analysis on reduced cross sections. Our approach relies on the dipole formulation of DIS together with the use of the non-linear running coupling BK equation for the description of the small-x dynamics. With the resulting parametrization we are able to describe the latest FL data measured by the H1 collaboration. Further, we discuss the kinematical domain where significant deviations from NLO-DGLAP should be expected and the ability of non-linnear physics to account for such deviations.

The linear-non-linear frontier for the Goldstone Higgs

International Nuclear Information System (INIS)

Gavela, M.B.; Saa, S.; Kanshin, K.; Machado, P.A.N.

2016-01-01

The minimal SO(5)/SO(4) σ-model is used as a template for the ultraviolet completion of scenarios in which the Higgs particle is a low-energy remnant of some high-energy dynamics, enjoying a (pseudo) Nambu-Goldstone-boson ancestry. Varying the σ mass allows one to sweep from the perturbative regime to the customary non-linear implementations. The low-energy benchmark effective non-linear Lagrangian for bosons and fermions is obtained, determining as well the operator coefficients including linear corrections. At first order in the latter, three effective bosonic operators emerge which are independent of the explicit soft breaking assumed. The Higgs couplings to vector bosons and fermions turn out to be quite universal: the linear corrections are proportional to the explicit symmetry-breaking parameters. Furthermore, we define an effective Yukawa operator which allows a simple parametrization and comparison of different heavy-fermion ultraviolet completions. In addition, one particular fermionic completion is explored in detail, obtaining the corresponding leading low-energy fermionic operators. (orig.)

The linear-non-linear frontier for the Goldstone Higgs

Energy Technology Data Exchange (ETDEWEB)

Gavela, M.B.; Saa, S. [IFT-UAM/CSIC, Universidad Autonoma de Madrid, Departamento de Fisica Teorica y Instituto de Fisica Teorica, Madrid (Spain); Kanshin, K. [Universita di Padova, Dipartimento di Fisica e Astronomia ' G. Galilei' , Padua (Italy); INFN, Padova (Italy); Machado, P.A.N. [IFT-UAM/CSIC, Universidad Autonoma de Madrid, Departamento de Fisica Teorica y Instituto de Fisica Teorica, Madrid (Spain); Fermi National Accelerator Laboratory, Theoretical Physics Department, Batavia, IL (United States)

2016-12-15

The minimal SO(5)/SO(4) σ-model is used as a template for the ultraviolet completion of scenarios in which the Higgs particle is a low-energy remnant of some high-energy dynamics, enjoying a (pseudo) Nambu-Goldstone-boson ancestry. Varying the σ mass allows one to sweep from the perturbative regime to the customary non-linear implementations. The low-energy benchmark effective non-linear Lagrangian for bosons and fermions is obtained, determining as well the operator coefficients including linear corrections. At first order in the latter, three effective bosonic operators emerge which are independent of the explicit soft breaking assumed. The Higgs couplings to vector bosons and fermions turn out to be quite universal: the linear corrections are proportional to the explicit symmetry-breaking parameters. Furthermore, we define an effective Yukawa operator which allows a simple parametrization and comparison of different heavy-fermion ultraviolet completions. In addition, one particular fermionic completion is explored in detail, obtaining the corresponding leading low-energy fermionic operators. (orig.)

Dark energy cosmology with generalized linear equation of state

International Nuclear Information System (INIS)

Babichev, E; Dokuchaev, V; Eroshenko, Yu

2005-01-01

Dark energy with the usually used equation of state p = wρ, where w const 0 ), where the constants α and ρ 0 are free parameters. This non-homogeneous linear equation of state provides the description of both hydrodynamically stable (α > 0) and unstable (α < 0) fluids. In particular, the considered cosmological model describes the hydrodynamically stable dark (and phantom) energy. The possible types of cosmological scenarios in this model are determined and classified in terms of attractors and unstable points by using phase trajectories analysis. For the dark energy case, some distinctive types of cosmological scenarios are possible: (i) the universe with the de Sitter attractor at late times, (ii) the bouncing universe, (iii) the universe with the big rip and with the anti-big rip. In the framework of a linear equation of state the universe filled with a phantom energy, w < -1, may have either the de Sitter attractor or the big rip

Vibration suppression in ultrasonic machining described by non-linear differential equations

International Nuclear Information System (INIS)

Kamel, M. M.; El-Ganaini, W. A. A.; Hamed, Y. S.

2009-01-01

Vibrations are usually undesired phenomena as they may cause damage or destruction of the system. However, sometimes they are desirable, as in ultrasonic machining (USM). In such case, the problem is a complicated one, as it is required to reduce the vibration of the machine head and have reasonable amplitude for the tool. In the present work, the coupling of two non-linear oscillators of the tool holder and tool representing ultrasonic cutting process is investigated. This leads to a two-degree-of-freedom system subjected to multi-external excitation force. The aim of this work is to control the tool holder behavior at simultaneous primary and internal resonance condition and have high amplitude for the tool. Multiple scale perturbation method is applied to obtain a solution up to the second order approximations. Other different resonance cases are reported and studied numerically. The stability of the system is investigated applying both phase-plane and frequency response techniques. The effects of the different parameters of the tool on the system behavior are studied numerically. Comparison with the available published work is reported

Comparison of Linear and Non-linear Regression Analysis to Determine Pulmonary Pressure in Hyperthyroidism.

Science.gov (United States)

Scarneciu, Camelia C; Sangeorzan, Livia; Rus, Horatiu; Scarneciu, Vlad D; Varciu, Mihai S; Andreescu, Oana; Scarneciu, Ioan

2017-01-01

This study aimed at assessing the incidence of pulmonary hypertension (PH) at newly diagnosed hyperthyroid patients and at finding a simple model showing the complex functional relation between pulmonary hypertension in hyperthyroidism and the factors causing it. The 53 hyperthyroid patients (H-group) were evaluated mainly by using an echocardiographical method and compared with 35 euthyroid (E-group) and 25 healthy people (C-group). In order to identify the factors causing pulmonary hypertension the statistical method of comparing the values of arithmetical means is used. The functional relation between the two random variables (PAPs and each of the factors determining it within our research study) can be expressed by linear or non-linear function. By applying the linear regression method described by a first-degree equation the line of regression (linear model) has been determined; by applying the non-linear regression method described by a second degree equation, a parabola-type curve of regression (non-linear or polynomial model) has been determined. We made the comparison and the validation of these two models by calculating the determination coefficient (criterion 1), the comparison of residuals (criterion 2), application of AIC criterion (criterion 3) and use of F-test (criterion 4). From the H-group, 47% have pulmonary hypertension completely reversible when obtaining euthyroidism. The factors causing pulmonary hypertension were identified: previously known- level of free thyroxin, pulmonary vascular resistance, cardiac output; new factors identified in this study- pretreatment period, age, systolic blood pressure. According to the four criteria and to the clinical judgment, we consider that the polynomial model (graphically parabola- type) is better than the linear one. The better model showing the functional relation between the pulmonary hypertension in hyperthyroidism and the factors identified in this study is given by a polynomial equation of second

Modulational Instability in Linearly Coupled Asymmetric Dual-Core Fibers

Directory of Open Access Journals (Sweden)

Arjunan Govindarajan

2017-06-01

Full Text Available We investigate modulational instability (MI in asymmetric dual-core nonlinear directional couplers incorporating the effects of the differences in effective mode areas and group velocity dispersions, as well as phase- and group-velocity mismatches. Using coupled-mode equations for this system, we identify MI conditions from the linearization with respect to small perturbations. First, we compare the MI spectra of the asymmetric system and its symmetric counterpart in the case of the anomalous group-velocity dispersion (GVD. In particular, it is demonstrated that the increase of the inter-core linear-coupling coefficient leads to a reduction of the MI gain spectrum in the asymmetric coupler. The analysis is extended for the asymmetric system in the normal-GVD regime, where the coupling induces and controls the MI, as well as for the system with opposite GVD signs in the two cores. Following the analytical consideration of the MI, numerical simulations are carried out to explore nonlinear development of the MI, revealing the generation of periodic chains of localized peaks with growing amplitudes, which may transform into arrays of solitons.

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...

S2SA preconditioning for the Sn equations with strictly non negative spatial discretization

International Nuclear Information System (INIS)

Bruss, D. E.; Morel, J. E.; Ragusa, J. C.

2013-01-01

Preconditioners based upon sweeps and diffusion-synthetic acceleration have been constructed and applied to the zeroth and first spatial moments of the 1-D S n transport equation using a strictly non negative nonlinear spatial closure. Linear and nonlinear preconditioners have been analyzed. The effectiveness of various combinations of these preconditioners are compared. In one dimension, nonlinear sweep preconditioning is shown to be superior to linear sweep preconditioning, and DSA preconditioning using nonlinear sweeps in conjunction with a linear diffusion equation is found to be essentially equivalent to nonlinear sweeps in conjunction with a nonlinear diffusion equation. The ability to use a linear diffusion equation has important implications for preconditioning the S n equations with a strictly non negative spatial discretization in multiple dimensions. (authors)

A spline-based non-linear diffeomorphism for multimodal prostate registration.

Science.gov (United States)

Mitra, Jhimli; Kato, Zoltan; Martí, Robert; Oliver, Arnau; Lladó, Xavier; Sidibé, Désiré; Ghose, Soumya; Vilanova, Joan C; Comet, Josep; Meriaudeau, Fabrice

2012-08-01

This paper presents a novel method for non-rigid registration of transrectal ultrasound and magnetic resonance prostate images based on a non-linear regularized framework of point correspondences obtained from a statistical measure of shape-contexts. The segmented prostate shapes are represented by shape-contexts and the Bhattacharyya distance between the shape representations is used to find the point correspondences between the 2D fixed and moving images. The registration method involves parametric estimation of the non-linear diffeomorphism between the multimodal images and has its basis in solving a set of non-linear equations of thin-plate splines. The solution is obtained as the least-squares solution of an over-determined system of non-linear equations constructed by integrating a set of non-linear functions over the fixed and moving images. However, this may not result in clinically acceptable transformations of the anatomical targets. Therefore, the regularized bending energy of the thin-plate splines along with the localization error of established correspondences should be included in the system of equations. The registration accuracies of the proposed method are evaluated in 20 pairs of prostate mid-gland ultrasound and magnetic resonance images. The results obtained in terms of Dice similarity coefficient show an average of 0.980±0.004, average 95% Hausdorff distance of 1.63±0.48 mm and mean target registration and target localization errors of 1.60±1.17 mm and 0.15±0.12 mm respectively. Copyright © 2012 Elsevier B.V. All rights reserved.

The Importance of Non-Linearity on Turbulent Fluxes

DEFF Research Database (Denmark)

Rokni, Masoud

2007-01-01

Two new non-linear models for the turbulent heat fluxes are derived and developed from the transport equation of the scalar passive flux. These models are called as non-linear eddy diffusivity and non-linear scalar flux. The structure of these models is compared with the exact solution which...... is derived from the Cayley-Hamilton theorem and contains a three term-basis plus a non-linear term due to scalar fluxes. In order to study the performance of the model itself, all other turbulent quantities are taken from a DNS channel flow data-base and thus the error source has been minimized. The results...... are compared with the DNS channel flow and good agreement is achieved. It has been shown that the non-linearity parts of the models are important to capture the true path of the streamwise scalar fluxes. It has also been shown that one of model constant should have negative sign rather than positive, which had...

General coupled mode theory in non-Hermitian waveguides.

Science.gov (United States)

Xu, Jing; Chen, Yuntian

2015-08-24

In the presence of loss and gain, the coupled mode equation on describing the mode hybridization of various waveguides or cavities, or cavities coupled to waveguides becomes intrinsically non-Hermitian. In such non-Hermitian waveguides, the standard coupled mode theory fails. We generalize the coupled mode theory with a properly defined inner product based on reaction conservation. We apply our theory to the non-Hermitian parity-time symmetric waveguides, and obtain excellent agreement with results obtained by finite element fullwave simulations. The theory presented here is typically formulated in space to study coupling between waveguides, which can be transformed into time domain by proper reformulation to study coupling between non-Hermitian resonators. Our theory has the strength of studying non-Hermitian optical systems with inclusion of the full vector fields, thus is useful to study and design non-Hermitian devices that support asymmetric and even nonreciprocal light propagations.

Krylov Subspace Methods for Complex Non-Hermitian Linear Systems. Thesis

Science.gov (United States)

Freund, Roland W.

1991-01-01

We consider Krylov subspace methods for the solution of large sparse linear systems Ax = b with complex non-Hermitian coefficient matrices. Such linear systems arise in important applications, such as inverse scattering, numerical solution of time-dependent Schrodinger equations, underwater acoustics, eddy current computations, numerical computations in quantum chromodynamics, and numerical conformal mapping. Typically, the resulting coefficient matrices A exhibit special structures, such as complex symmetry, or they are shifted Hermitian matrices. In this paper, we first describe a Krylov subspace approach with iterates defined by a quasi-minimal residual property, the QMR method, for solving general complex non-Hermitian linear systems. Then, we study special Krylov subspace methods designed for the two families of complex symmetric respectively shifted Hermitian linear systems. We also include some results concerning the obvious approach to general complex linear systems by solving equivalent real linear systems for the real and imaginary parts of x. Finally, numerical experiments for linear systems arising from the complex Helmholtz equation are reported.

A Non-linear Stochastic Model for an Office Building with Air Infiltration

DEFF Research Database (Denmark)

Thavlov, Anders; Madsen, Henrik

2015-01-01

This paper presents a non-linear heat dynamic model for a multi-room office building with air infiltration. Several linear and non-linear models, with and without air infiltration, are investigated and compared. The models are formulated using stochastic differential equations and the model...

Linear measure functional differential equations with infinite delay

OpenAIRE

Monteiro, G. (Giselle Antunes); Slavík, A.

2014-01-01

We use the theory of generalized linear ordinary differential equations in Banach spaces to study linear measure functional differential equations with infinite delay. We obtain new results concerning the existence, uniqueness, and continuous dependence of solutions. Even for equations with a finite delay, our results are stronger than the existing ones. Finally, we present an application to functional differential equations with impulses.

Localization of the eigenvalues of linear integral equations with applications to linear ordinary differential equations.

Science.gov (United States)

Sloss, J. M.; Kranzler, S. K.

1972-01-01

The equivalence of a considered integral equation form with an infinite system of linear equations is proved, and the localization of the eigenvalues of the infinite system is expressed. Error estimates are derived, and the problems of finding upper bounds and lower bounds for the eigenvalues are solved simultaneously.

Bayesian analysis of non-linear differential equation models with application to a gut microbial ecosystem.

Science.gov (United States)

Lawson, Daniel J; Holtrop, Grietje; Flint, Harry

2011-07-01

Process models specified by non-linear dynamic differential equations contain many parameters, which often must be inferred from a limited amount of data. We discuss a hierarchical Bayesian approach combining data from multiple related experiments in a meaningful way, which permits more powerful inference than treating each experiment as independent. The approach is illustrated with a simulation study and example data from experiments replicating the aspects of the human gut microbial ecosystem. A predictive model is obtained that contains prediction uncertainty caused by uncertainty in the parameters, and we extend the model to capture situations of interest that cannot easily be studied experimentally. Copyright © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Structure Learning in Stochastic Non-linear Dynamical Systems

Science.gov (United States)

Morris, R. D.; Smelyanskiy, V. N.; Luchinsky, D. G.

2005-12-01

A great many systems can be modeled in the non-linear dynamical systems framework, as x˙ = f(x) + ξ(t), where f(x) is the potential function for the system, and ξ(t) is the driving noise. Modeling the potential using a set of basis functions, we derive the posterior for the basis coefficients. A more challenging problem is to determine the set of basis functions that are required to model a particular system. We show that using the Bayesian Information Criteria (BIC) to rank models, and the beam search technique, that we can accurately determine the structure of simple non-linear dynamical system models, and the structure of the coupling between non-linear dynamical systems where the individual systems are known. This last case has important ecological applications, for example in predator-prey systems, where the very structure of the coupling between predator-prey pairs can have great ecological significance.

Saturation and linear transport equation

International Nuclear Information System (INIS)

Kutak, K.

2009-03-01

We show that the GBW saturation model provides an exact solution to the one dimensional linear transport equation. We also show that it is motivated by the BK equation considered in the saturated regime when the diffusion and the splitting term in the diffusive approximation are balanced by the nonlinear term. (orig.)

On Coupled System of Navier-Stokes Equations and Temperature

African Journals Online (AJOL)

Dr. Anthony Peter

ABSTRACT. This paper deals with the coupled system of Navier-Stokes equations and temperature (Thermohydraulics) in a strip in the class of spatially non-decaying (infinite-energy) solutions belonging to the properly chosen uniformly local Sobolev spaces. The global well-posedness and dissipativity of the Navier- ...

Coupling Integrable Couplings of an Equation Hierarchy

International Nuclear Information System (INIS)

Wang Hui; Xia Tie-Cheng

2013-01-01

Based on a kind of Lie algebra G proposed by Zhang, one isospectral problem is designed. Under the framework of zero curvature equation, a new kind of integrable coupling of an equation hierarchy is generated using the methods proposed by Ma and Gao. With the help of variational identity, we get the Hamiltonian structure of the hierarchy. (general)

Quasi-linear equation for magnetoplasma oscillations in the weakly relativistic approximation

International Nuclear Information System (INIS)

Rizzato, F.B.

1985-01-01

Some limitations which are present in the dynamical equations for collisionless plasmas are discussed. Some elementary corrections to the linear theories are obtained in a heuristic form, which directly lead to the so-called quasi-linear theories in its non-relativistic and relativistic forms. The effect of the relativistic variation of the gyrofrequency on the diffusion coefficient is examined in a typically perturbative approximation. (author)

Non linear evolution of plasma waves excited to mode conversion at the vicinity of plasma resonance. Application to experiments of ionosphere modification

International Nuclear Information System (INIS)

Cros, Brigitte

1989-01-01

This research thesis reports the study of the non linear evolution of plasma waves excited by mode conversion in a non homogeneous, non collisional, and free-of-external-magnetic-field plasma. Experiments performed in the microwave domain in a plasma created by means of a multi-polar device show that the evolution of plasma waves displays a transition between a non linear quasi-steady regime and a stochastic regime when the power of incident electromagnetic waves or plasma gradient length is increased. These regimes are characterized through a numerical resolution of Zakharov equations which describe the coupled evolution of plasma wave envelope and low frequency density perturbations [fr

Linearized pseudo-Einstein equations on the Heisenberg group

Science.gov (United States)

Barletta, Elisabetta; Dragomir, Sorin; Jacobowitz, Howard

2017-02-01

We study the pseudo-Einstein equation R11bar = 0 on the Heisenberg group H1 = C × R. We consider first order perturbations θɛ =θ0 + ɛ θ and linearize the pseudo-Einstein equation about θ0 (the canonical Tanaka-Webster flat contact form on H1 thought of as a strictly pseudoconvex CR manifold). If θ =e2uθ0 the linearized pseudo-Einstein equation is Δb u - 4 | Lu|2 = 0 where Δb is the sublaplacian of (H1 ,θ0) and L bar is the Lewy operator. We solve the linearized pseudo-Einstein equation on a bounded domain Ω ⊂H1 by applying subelliptic theory i.e. existence and regularity results for weak subelliptic harmonic maps. We determine a solution u to the linearized pseudo-Einstein equation, possessing Heisenberg spherical symmetry, and such that u(x) → - ∞ as | x | → + ∞.

STABLE STATIONARY STATES OF NON-LOCAL INTERACTION EQUATIONS

KAUST Repository

FELLNER, KLEMENS

2010-12-01

In this paper, we are interested in the large-time behaviour of a solution to a non-local interaction equation, where a density of particles/individuals evolves subject to an interaction potential and an external potential. It is known that for regular interaction potentials, stable stationary states of these equations are generically finite sums of Dirac masses. For a finite sum of Dirac masses, we give (i) a condition to be a stationary state, (ii) two necessary conditions of linear stability w.r.t. shifts and reallocations of individual Dirac masses, and (iii) show that these linear stability conditions imply local non-linear stability. Finally, we show that for regular repulsive interaction potential Wε converging to a singular repulsive interaction potential W, the Dirac-type stationary states ρ̄ ε approximate weakly a unique stationary state ρ̄ ∈ L∞. We illustrate our results with numerical examples. © 2010 World Scientific Publishing Company.

Solving Nonlinear Coupled Differential Equations

Science.gov (United States)

Mitchell, L.; David, J.

1986-01-01

Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.

Electron non-linearities in Langmuir waves with application to beat-wave experiments

International Nuclear Information System (INIS)

Bell, A.R.; Gibbon, P.

1988-01-01

Non-linear Langmuir waves are examined in the context of the beat-wave accelerator. With a background of immobile ions the waves in one dimension are subject to the relativistic non-linearity of Rosenbluth, M.N. and Liu, C.S., Phys. Rev. Lett., 1972, 29, 701. In two or three dimensions, other electron non-linearities occur which involve electric and magnetic fields. The quasi-linear equations for these non-linearities are developed and solved numerically in a geometry representative of laser-driven beat waves. (author)

Non-linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents

International Nuclear Information System (INIS)

Hoelzl, M; Merkel, P; Lackner, K; Strumberger, E; Huijsmans, G T A; Aleynikova, K; Liu, F; Atanasiu, C; Nardon, E; Fil, A; McAdams, R; Chapman, I

2014-01-01

The dynamics of large scale plasma instabilities can be strongly influenced by the mutual interaction with currents flowing in conducting vessel structures. Especially eddy currents caused by time-varying magnetic perturbations and halo currents flowing directly from the plasma into the walls are important. The relevance of a resistive wall model is directly evident for Resistive Wall Modes (RWMs) or Vertical Displacement Events (VDEs). However, also the linear and non-linear properties of most other large-scale instabilities may be influenced significantly by the interaction with currents in conducting structures near the plasma. The understanding of halo currents arising during disruptions and VDEs, which are a serious concern for ITER as they may lead to strong asymmetric forces on vessel structures, could also benefit strongly from these non-linear modeling capabilities. Modeling the plasma dynamics and its interaction with wall currents requires solving the magneto-hydrodynamic (MHD) equations in realistic toroidal X-point geometry consistently coupled with a model for the vacuum region and the resistive conducting structures. With this in mind, the non-linear finite element MHD code JOREK [1, 2] has been coupled [3] with the resistive wall code STARWALL [4], which allows us to include the effects of eddy currents in 3D conducting structures in non-linear MHD simulations. This article summarizes the capabilities of the coupled JOREK-STARWALL system and presents benchmark results as well as first applications to non-linear simulations of RWMs, VDEs, disruptions triggered by massive gas injection, and Quiescent H-Mode. As an outlook, the perspectives for extending the model to halo currents are described

A Homotopy-Perturbation analysis of the non-linear contaminant ...

African Journals Online (AJOL)

In this research work, a Homotopy-perturbation analysis of a non –linear contaminant flow equation with an initial continuous point source is provided. The equation is characterized by advection, diffusion and adsorption. We assume that the adsorption term is modeled by Freudlich Isotherm. We provide an approximation of ...

Linear causal modeling with structural equations

CERN Document Server

Mulaik, Stanley A

2009-01-01

Emphasizing causation as a functional relationship between variables that describe objects, Linear Causal Modeling with Structural Equations integrates a general philosophical theory of causation with structural equation modeling (SEM) that concerns the special case of linear causal relations. In addition to describing how the functional relation concept may be generalized to treat probabilistic causation, the book reviews historical treatments of causation and explores recent developments in experimental psychology on studies of the perception of causation. It looks at how to perceive causal

Quantization of non-Hamiltonian physical systems

International Nuclear Information System (INIS)

Bolivar, A.O.

1998-09-01

We propose a general method of quantization of non-Hamiltonian physical systems. Applying it, for example, to a dissipative system coupled to a thermal reservoir described by the Fokker-Planck equation, we are able to obtain the Caldeira-Leggett master equation, the non-linear Schroedinger-Langevin equation and Caldirola-Kanai equation (with an additional term), as particular cases. (author)

Non-linear aeroelastic prediction for aircraft applications

Science.gov (United States)

de C. Henshaw, M. J.; Badcock, K. J.; Vio, G. A.; Allen, C. B.; Chamberlain, J.; Kaynes, I.; Dimitriadis, G.; Cooper, J. E.; Woodgate, M. A.; Rampurawala, A. M.; Jones, D.; Fenwick, C.; Gaitonde, A. L.; Taylor, N. V.; Amor, D. S.; Eccles, T. A.; Denley, C. J.

2007-05-01

in this domain. This is set within the context of a generic industrial process and the requirements of UK and US aeroelastic qualification. A range of test cases, from simple small DOF cases to full aircraft, have been used to evaluate and validate the non-linear methods developed and to make comparison with the linear methods in everyday use. These have focused mainly on aerodynamic non-linearity, although some results for structural non-linearity are also presented. The challenges associated with time domain (coupled computational fluid dynamics-computational structural model (CFD-CSM)) methods have been addressed through the development of grid movement, fluid-structure coupling, and control surface movement technologies. Conclusions regarding the accuracy and computational cost of these are presented. The computational cost of time-domain methods, despite substantial improvements in efficiency, remains high. However, significant advances have been made in reduced order methods, that allow non-linear behaviour to be modelled, but at a cost comparable with that of the regular linear methods. Of particular note is a method based on Hopf bifurcation that has reached an appropriate maturity for deployment on real aircraft configurations, though only limited results are presented herein. Results are also presented for dynamically linearised CFD approaches that hold out the possibility of non-linear results at a fraction of the cost of time coupled CFD-CSM methods. Local linearisation approaches (higher order harmonic balance and continuation method) are also presented; these have the advantage that no prior assumption of the nature of the aeroelastic instability is required, but currently these methods are limited to low DOF problems and it is thought that these will not reach a level of maturity appropriate to real aircraft problems for some years to come. Nevertheless, guidance on the most likely approaches has been derived and this forms the basis for ongoing

A coupling method for a cardiovascular simulation model which includes the Kalman filter.

Science.gov (United States)

Hasegawa, Yuki; Shimayoshi, Takao; Amano, Akira; Matsuda, Tetsuya

2012-01-01

Multi-scale models of the cardiovascular system provide new insight that was unavailable with in vivo and in vitro experiments. For the cardiovascular system, multi-scale simulations provide a valuable perspective in analyzing the interaction of three phenomenons occurring at different spatial scales: circulatory hemodynamics, ventricular structural dynamics, and myocardial excitation-contraction. In order to simulate these interactions, multiscale cardiovascular simulation systems couple models that simulate different phenomena. However, coupling methods require a significant amount of calculation, since a system of non-linear equations must be solved for each timestep. Therefore, we proposed a coupling method which decreases the amount of calculation by using the Kalman filter. In our method, the Kalman filter calculates approximations for the solution to the system of non-linear equations at each timestep. The approximations are then used as initial values for solving the system of non-linear equations. The proposed method decreases the number of iterations required by 94.0% compared to the conventional strong coupling method. When compared with a smoothing spline predictor, the proposed method required 49.4% fewer iterations.

Modelling non-linear effects of dark energy

Science.gov (United States)

Bose, Benjamin; Baldi, Marco; Pourtsidou, Alkistis

2018-04-01

We investigate the capabilities of perturbation theory in capturing non-linear effects of dark energy. We test constant and evolving w models, as well as models involving momentum exchange between dark energy and dark matter. Specifically, we compare perturbative predictions at 1-loop level against N-body results for four non-standard equations of state as well as varying degrees of momentum exchange between dark energy and dark matter. The interaction is modelled phenomenologically using a time dependent drag term in the Euler equation. We make comparisons at the level of the matter power spectrum and the redshift space monopole and quadrupole. The multipoles are modelled using the Taruya, Nishimichi and Saito (TNS) redshift space spectrum. We find perturbation theory does very well in capturing non-linear effects coming from dark sector interaction. We isolate and quantify the 1-loop contribution coming from the interaction and from the non-standard equation of state. We find the interaction parameter ξ amplifies scale dependent signatures in the range of scales considered. Non-standard equations of state also give scale dependent signatures within this same regime. In redshift space the match with N-body is improved at smaller scales by the addition of the TNS free parameter σv. To quantify the importance of modelling the interaction, we create mock data sets for varying values of ξ using perturbation theory. This data is given errors typical of Stage IV surveys. We then perform a likelihood analysis using the first two multipoles on these sets and a ξ=0 modelling, ignoring the interaction. We find the fiducial growth parameter f is generally recovered even for very large values of ξ both at z=0.5 and z=1. The ξ=0 modelling is most biased in its estimation of f for the phantom w=‑1.1 case.

Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation

Directory of Open Access Journals (Sweden)

Hamidreza Rezazadeh

2014-05-01

Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.

Scattering of massless lumps and non-local charges in the two-dimensional classical non-linear sigma-model

International Nuclear Information System (INIS)

Luescher, M.; Pohlmeyer, K.

1977-09-01

Finite energy solutions of the field equations of the non-linear sigma-model are shown to decay asymptotically into massless lumps. By means of a linear eigenvalue problem connected with the field equations we then find an infinite set of dynamical conserved charges. They, however, do not provide sufficient information to decode the complicated scattering of lumps. (orig.) [de

Dynamical and geometric aspects of Hamilton-Jacobi and linearized Monge-Ampère equations VIASM 2016

CERN Document Server

Tran, Hung

2017-01-01

Consisting of two parts, the first part of this volume is an essentially self-contained exposition of the geometric aspects of local and global regularity theory for the Monge–Ampère and linearized Monge–Ampère equations. As an application, we solve the second boundary value problem of the prescribed affine mean curvature equation, which can be viewed as a coupling of the latter two equations. Of interest in its own right, the linearized Monge–Ampère equation also has deep connections and applications in analysis, fluid mechanics and geometry, including the semi-geostrophic equations in atmospheric flows, the affine maximal surface equation in affine geometry and the problem of finding Kahler metrics of constant scalar curvature in complex geometry. Among other topics, the second part provides a thorough exposition of the large time behavior and discounted approximation of Hamilton–Jacobi equations, which have received much attention in the last two decades, and a new approach to the subject, the n...

Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods

Science.gov (United States)

Faraji, Shirin; Matsika, Spiridoula; Krylov, Anna I.

2018-01-01

We report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states' character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. The numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches.

Systems of Inhomogeneous Linear Equations

Science.gov (United States)

Scherer, Philipp O. J.

Many problems in physics and especially computational physics involve systems of linear equations which arise e.g. from linearization of a general nonlinear problem or from discretization of differential equations. If the dimension of the system is not too large standard methods like Gaussian elimination or QR decomposition are sufficient. Systems with a tridiagonal matrix are important for cubic spline interpolation and numerical second derivatives. They can be solved very efficiently with a specialized Gaussian elimination method. Practical applications often involve very large dimensions and require iterative methods. Convergence of Jacobi and Gauss-Seidel methods is slow and can be improved by relaxation or over-relaxation. An alternative for large systems is the method of conjugate gradients.

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.

Diffusion phenomenon for linear dissipative wave equations

KAUST Repository

Said-Houari, Belkacem

2012-01-01

In this paper we prove the diffusion phenomenon for the linear wave equation. To derive the diffusion phenomenon, a new method is used. In fact, for initial data in some weighted spaces, we prove that for {equation presented} decays with the rate {equation presented} [0,1] faster than that of either u or v, where u is the solution of the linear wave equation with initial data {equation presented} [0,1], and v is the solution of the related heat equation with initial data v 0 = u 0 + u 1. This result improves the result in H. Yang and A. Milani [Bull. Sci. Math. 124 (2000), 415-433] in the sense that, under the above restriction on the initial data, the decay rate given in that paper can be improved by t -γ/2. © European Mathematical Society.

Numerical solution of non-linear dual-phase-lag bioheat transfer equation within skin tissues.

Science.gov (United States)

Kumar, Dinesh; Kumar, P; Rai, K N

2017-11-01

This paper deals with numerical modeling and simulation of heat transfer in skin tissues using non-linear dual-phase-lag (DPL) bioheat transfer model under periodic heat flux boundary condition. The blood perfusion is assumed temperature-dependent which results in non-linear DPL bioheat transfer model in order to predict more accurate results. A numerical method of line which is based on finite difference and Runge-Kutta (4,5) schemes, is used to solve the present non-linear problem. Under specific case, the exact solution has been obtained and compared with the present numerical scheme, and we found that those are in good agreement. A comparison based on model selection criterion (AIC) has been made among non-linear DPL models when the variation of blood perfusion rate with temperature is of constant, linear and exponential type with the experimental data and it has been found that non-linear DPL model with exponential variation of blood perfusion rate is closest to the experimental data. In addition, it is found that due to absence of phase-lag phenomena in Pennes bioheat transfer model, it achieves steady state more quickly and always predict higher temperature than thermal and DPL non-linear models. The effect of coefficient of blood perfusion rate, dimensionless heating frequency and Kirchoff number on dimensionless temperature distribution has also been analyzed. The whole analysis is presented in dimensionless form. Copyright © 2017 Elsevier Inc. All rights reserved.

A new linearized equation for servo valve in hydraulic control systems

International Nuclear Information System (INIS)

Kim, Tae Hyung; Lee, Ill Yeong

2002-01-01

In the procedure of the hydraulic control system analysis, a linearized approximate equation described by the first order term of Taylor's series has been widely used. Such a linearized equation is effective just near the operating point. And, as of now, there are no general standards on how to determine the operating point of a servo valve in the process of applying the linearized equation. So, in this study, a new linearized equation for valve characteristics is proposed as a modified form of the existing linearized equation. And, a method for selecting an optimal operating point is proposed for the new linearized equation. The effectiveness of the new linearized equation is confirmed through numerical simulations and experiments for a model hydraulic control system

Spinor-electron wave guided modes in coupled quantum wells structures by solving the Dirac equation

International Nuclear Information System (INIS)

Linares, Jesus; Nistal, Maria C.

2009-01-01

A quantum analysis based on the Dirac equation of the propagation of spinor-electron waves in coupled quantum wells, or equivalently coupled electron waveguides, is presented. The complete optical wave equations for Spin-Up (SU) and Spin-Down (SD) spinor-electron waves in these electron guides couplers are derived from the Dirac equation. The relativistic amplitudes and dispersion equations of the spinor-electron wave-guided modes in a planar quantum coupler formed by two coupled quantum wells, or equivalently by two coupled slab electron waveguides, are exactly derived. The main outcomes related to the spinor modal structure, such as the breaking of the non-relativistic degenerate spin states, the appearance of phase shifts associated with the spin polarization and so on, are shown.

Hypocoercivity for linear kinetic equations conserving mass

KAUST Repository

Dolbeault, Jean; Mouhot, Clé ment; Schmeiser, Christian

2015-01-01

We develop a new method for proving hypocoercivity for a large class of linear kinetic equations with only one conservation law. Local mass conservation is assumed at the level of the collision kernel, while transport involves a confining potential, so that the solution relaxes towards a unique equilibrium state. Our goal is to evaluate in an appropriately weighted $ L^2$ norm the exponential rate of convergence to the equilibrium. The method covers various models, ranging from diffusive kinetic equations like Vlasov-Fokker-Planck equations, to scattering models or models with time relaxation collision kernels corresponding to polytropic Gibbs equilibria, including the case of the linear Boltzmann model. In this last case and in the case of Vlasov-Fokker-Planck equations, any linear or superlinear growth of the potential is allowed. - See more at: http://www.ams.org/journals/tran/2015-367-06/S0002-9947-2015-06012-7/#sthash.ChjyK6rc.dpuf

Hypocoercivity for linear kinetic equations conserving mass

KAUST Repository

Dolbeault, Jean

2015-02-03

We develop a new method for proving hypocoercivity for a large class of linear kinetic equations with only one conservation law. Local mass conservation is assumed at the level of the collision kernel, while transport involves a confining potential, so that the solution relaxes towards a unique equilibrium state. Our goal is to evaluate in an appropriately weighted $ L^2$ norm the exponential rate of convergence to the equilibrium. The method covers various models, ranging from diffusive kinetic equations like Vlasov-Fokker-Planck equations, to scattering models or models with time relaxation collision kernels corresponding to polytropic Gibbs equilibria, including the case of the linear Boltzmann model. In this last case and in the case of Vlasov-Fokker-Planck equations, any linear or superlinear growth of the potential is allowed. - See more at: http://www.ams.org/journals/tran/2015-367-06/S0002-9947-2015-06012-7/#sthash.ChjyK6rc.dpuf

Non-linear thermal fluctuations in a diode

NARCIS (Netherlands)

Kampen, N.G. van

As an example of non-linear noise the fluctuations in a circuit consisting of a diode and a condenser C are studied. From the master equation for this system the following results are derived. 1. (i) The equilibrium distribution of the voltage is rigorously Gaussian, the average voltage being

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.

Plasma heating by non-linear wave-Plasma interaction | Echi ...

African Journals Online (AJOL)

We simulate the non-linear interaction of waves with magnetized tritium plasma with the aim of determining the parameter values that characterize the response of the plasma. The wave-plasma interaction has a non-conservative Hamiltonian description. The resulting system of Hamilton's equations is integrated numerically ...

Non-linear coupling of drift modes in a quadrupole

International Nuclear Information System (INIS)

Elliott, J.A.; Sandeman, J.C.; Tessema, G.Y.

1990-01-01

We report continuing experimental studies of non-linear interactions of drift waves, with direct evidence of a growth saturation mechanism by transfer of energy to lower frequency modes. Wave launching experiments show that the decay rate of drift waves can be strongly amplitude dependent. (author) 9 refs., 5 figs

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.

Non-linear analysis of skew thin plate by finite difference method

International Nuclear Information System (INIS)

Kim, Chi Kyung; Hwang, Myung Hwan

2012-01-01

This paper deals with a discrete analysis capability for predicting the geometrically nonlinear behavior of skew thin plate subjected to uniform pressure. The differential equations are discretized by means of the finite difference method which are used to determine the deflections and the in-plane stress functions of plates and reduced to several sets of linear algebraic simultaneous equations. For the geometrically non-linear, large deflection behavior of the plate, the non-linear plate theory is used for the analysis. An iterative scheme is employed to solve these quasi-linear algebraic equations. Several problems are solved which illustrate the potential of the method for predicting the finite deflection and stress. For increasing lateral pressures, the maximum principal tensile stress occurs at the center of the plate and migrates toward the corners as the load increases. It was deemed important to describe the locations of the maximum principal tensile stress as it occurs. The load-deflection relations and the maximum bending and membrane stresses for each case are presented and discussed

Problem on primary radiation filtration effect on form of coupling equations during X-ray fluorescence analysis

International Nuclear Information System (INIS)

Lavrent'ev, Yu.G.; Kuznetsova, A.I.

1976-01-01

A simplified method for x-ray fluorescence analysis is given. It is shown that the system of coupling equations with constant coefficients and with the number of equations equal to the number of unknown elements allows to obtain the same accuracy of the analysis as with the considerably more complex equations with variable coefficients which take into account the filtration of the primary radiation in a direct way. The system can even be more simplified by using linear equations with constant coefficients. In order to test these systems and to compare them with known coupling equations experimental data for the determination of zirconium and niobium from 16 artificial preparations with fillers of variable composition are presented. The calculation of the absorption of the secondary as well as the primary radiation by means of the proposed equations with constant coefficients is sufficiently good

Remarks on the derivation of the governing equations for the dynamics of a nonlinear beam to a non ideal shaft coupling

Energy Technology Data Exchange (ETDEWEB)

Fenili, André; Lopes Rebello da Fonseca Brasil, Reyolando Manoel [Universidade Federal do ABC (UFABC), Centro de Engenharia, Modelagem e Ciências Sociais Aplicadas (CECS) / Aerospace Engineering Santo André, São Paulo (Brazil); Balthazar, José M., E-mail: jmbaltha@gmail.com [Universidade Federal do ABC (UFABC), Centro de Engenharia, Modelagem e Ciências Sociais Aplicadas (CECS) / Aerospace Engineering Santo André, São Paulo, Brazil and Universidade Estadual Paulista, Faculdade de Engenharia Mec and #x00E (Brazil); Francisco, Cayo Prado Fernandes [Universidade Federal do ABC (UFABC), Centro de Engenharia, Modelagem e Ciências Sociais Aplicadas (CECS) / Aerospace Engineering Santo André, São Paulo, Brazil and Instituto de Aeronáutica e Espaço, Departamento de (Brazil)

2014-12-10

We derive nonlinear governing equations without assuming that the beam is inextensible. The derivation couples the equations that govern a weak electric motor, which is used to rotate the base of the beam, to those that govern the motion of the beam. The system is considered non-ideal in the sense that the response of the motor to an applied voltage and the motion of the beam must be obtained interactively. The moment that the motor exerts on the base of the beam cannot be determined without solving for the motion of the beam.

Remarks on the derivation of the governing equations for the dynamics of a nonlinear beam to a non ideal shaft coupling

International Nuclear Information System (INIS)

Fenili, André; Lopes Rebello da Fonseca Brasil, Reyolando Manoel; Balthazar, José M.; Francisco, Cayo Prado Fernandes

2014-01-01

We derive nonlinear governing equations without assuming that the beam is inextensible. The derivation couples the equations that govern a weak electric motor, which is used to rotate the base of the beam, to those that govern the motion of the beam. The system is considered non-ideal in the sense that the response of the motor to an applied voltage and the motion of the beam must be obtained interactively. The moment that the motor exerts on the base of the beam cannot be determined without solving for the motion of the beam

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)

Invariant imbedding equations for linear scattering problems

International Nuclear Information System (INIS)

Apresyan, L.

1988-01-01

A general form of the invariant imbedding equations is investigated for the linear problem of scattering by a bounded scattering volume. The conditions for the derivability of such equations are described. It is noted that the possibility of the explicit representation of these equations for a sphere and for a layer involves the separation of variables in the unperturbed wave equation

A non-differentiable solution for the local fractional telegraph equation

Directory of Open Access Journals (Sweden)

Li Jie

2017-01-01

Full Text Available In this paper, we consider the linear telegraph equations with local fractional derivative. The local fractional Laplace series expansion method is used to handle the local fractional telegraph equation. The analytical solution with the non-differentiable graphs is discussed in detail. The proposed method is efficient and accurate.

Non-Linear Dynamics of Saturn's Rings

Science.gov (United States)

Esposito, L. W.

2016-12-01

Non-linear processes can explain why Saturn's rings are so active and dynamic. Ring systems differ from simple linear systems in two significant ways: 1. They are systems of granular material: where particle-to-particle collisions dominate; thus a kinetic, not a fluid description needed. Stresses are strikingly inhomogeneous and fluctuations are large compared to equilibrium. 2. They are strongly forced by resonances: which drive a non-linear response, that push the system across thresholds that lead to persistent states. Some of this non-linearity is captured in a simple Predator-Prey Model: Periodic forcing from the moon causes streamline crowding; This damps the relative velocity. About a quarter phase later, the aggregates stir the system to higher relative velocity and the limit cycle repeats each orbit, with relative velocity ranging from nearly zero to a multiple of the orbit average. Summary of Halo Results: A predator-prey model for ring dynamics produces transient structures like `straw' that can explain the halo morphology and spectroscopy: Cyclic velocity changes cause perturbed regions to reach higher collision speeds at some orbital phases, which preferentially removes small regolith particles; surrounding particles diffuse back too slowly to erase the effect: this gives the halo morphology; this requires energetic collisions (v ≈ 10m/sec, with throw distances about 200km, implying objects of scale R ≈ 20km).Transform to Duffing Eqn : With the coordinate transformation, z = M2/3, the Predator-Prey equations can be combined to form a single second-order differential equation with harmonic resonance forcing.Ring dynamics and history implications: Moon-triggered clumping explains both small and large particles at resonances. We calculate the stationary size distribution using a cell-to-cell mapping procedure that converts the phase-plane trajectories to a Markov chain. Approximating it as an asymmetric random walk with reflecting boundaries

Integrable coupling system of fractional soliton equation hierarchy

Energy Technology Data Exchange (ETDEWEB)

Yu Fajun, E-mail: yfajun@163.co [College of Maths and Systematic Science, Shenyang Normal University, Shenyang 110034 (China)

2009-10-05

In this Letter, we consider the derivatives and integrals of fractional order and present a class of the integrable coupling system of the fractional order soliton equations. The fractional order coupled Boussinesq and KdV equations are the special cases of this class. Furthermore, the fractional AKNS soliton equation hierarchy is obtained.

Two hierarchies of multi-component Kaup-Newell equations and theirs integrable couplings

International Nuclear Information System (INIS)

Zhu Fubo; Ji Jie; Zhang Jianbin

2008-01-01

Two hierarchies of multi-component Kaup-Newell equations are derived from an arbitrary order matrix spectral problem, including positive non-isospectral Kaup-Newell hierarchy and negative non-isospectral Kaup-Newell hierarchy. Moreover, new integrable couplings of the resulting Kaup-Newell soliton hierarchies are constructed by enlarging the associated matrix spectral problem

Linear Einstein equations and Kerr-Schild maps

International Nuclear Information System (INIS)

Gergely, Laszlo A

2002-01-01

We prove that given a solution of the Einstein equations g ab for the matter field T ab , an autoparallel null vector field l a and a solution (l a l c , T ac ) of the linearized Einstein equation on the given background, the Kerr-Schild metric g ac + λl a l c (λ arbitrary constant) is an exact solution of the Einstein equation for the energy-momentum tensor T ac + λT ac + λ 2 l (a T c)b l b . The mixed form of the Einstein equation for Kerr-Schild metrics with autoparallel null congruence is also linear. Some more technical conditions hold when the null congruence is not autoparallel. These results generalize previous theorems for vacuum due to Xanthopoulos and for flat seed spacetime due to Guerses and Guersey

Coupling of discrete ordinates methods by transmission of boundary conditions in solving the neutron transport equation in slab geometry; Couplage de discretisations aux ordonnees discretes d`equations de transport 1D par passage de conditions frontieres

Energy Technology Data Exchange (ETDEWEB)

Bal, G. [Departement MMN, Service IMA, Direction des Etudes et Recherches, Electricite de France (EDF), 92 - Clamart (France)

1995-10-01

Neutron transport in nuclear reactors is quite well modelled by the linear Boltzmann transport equation. Its solution is relatively easy, but unfortunately too expensive to achieve whole core computations. Thus, we have to simplify it, for example by homogenizing some physical characteristics. However, the solution may then be inaccurate. Moreover, in strongly homogeneous areas, the error may be too big. Then we would like to deal with such an inconvenient by solving the equation accurately on this area, but more coarsely away from it, so that the computation is not too expensive. This problem is the subject of a thesis. We present here some results obtained for slab geometry. The couplings between the fine and coarse discretization regions could be conceived in a number of approaches. Here, we only deal with the coupling at crossing the interface between two sub-domains. In the first section, we present the coupling of discrete ordinate methods for solving the homogeneous, isotropic and mono-kinetic equation. Coupling operators are defined and shown to be optimal. The second and the third sections are devoted to an extension of the previous results when the equation is non-homogeneous, anisotropic and multigroup (under some restrictive assumptions). Some numerical results are given in the case of isotropic and mono-kinetic equations. (author) 15 refs.

Uniqueness of non-linear ground states for fractional Laplacians in R

DEFF Research Database (Denmark)

Frank, Rupert L.; Lenzmann, Enno

2013-01-01

We prove uniqueness of ground state solutions Q = Q(|x|) ≥ 0 of the non-linear equation (−Δ)sQ+Q−Qα+1=0inR,where 0 fractional Laplacian in one dimension. In particular, we answer affirmatively an open question...... recently raised by Kenig–Martel–Robbiano and we generalize (by completely different techniques) the specific uniqueness result obtained by Amick and Toland for s=12 and α = 1 in [5] for the Benjamin–Ono equation. As a technical key result in this paper, we show that the associated linearized operator L...... + = (−Δ) s +1−(α+1)Q α is non-degenerate; i.e., its kernel satisfies ker L + = span{Q′}. This result about L + proves a spectral assumption, which plays a central role for the stability of solitary waves and blowup analysis for non-linear dispersive PDEs with fractional Laplacians, such as the generalized...

Solving large sets of coupled equations iteratively by vector processing on the CYBER 205 computer

International Nuclear Information System (INIS)

Tolsma, L.D.

1985-01-01

The set of coupled linear second-order differential equations which has to be solved for the quantum-mechanical description of inelastic scattering of atomic and nuclear particles can be rewritten as an equivalent set of coupled integral equations. When some type of functions is used as piecewise analytic reference solutions, the integrals that arise in this set can be evaluated analytically. The set of integral equations can be solved iteratively. For the results mentioned an inward-outward iteration scheme has been applied. A concept of vectorization of coupled-channel Fortran programs, based on this integral method, is presented for the use on the Cyber 205 computer. It turns out that, for two heavy ion nuclear scattering test cases, this vector algorithm gives an overall speed-up of about a factor of 2 to 3 compared to a highly optimized scalar algorithm for a one vector pipeline computer

Equations of motion for massive spin 2 field coupled to gravity

International Nuclear Information System (INIS)

Buchbinder, I.L.; Gitman, D.M.; Krykhtin, V.A.; Pershin, V.D.

2000-01-01

We investigate the problems of consistency and causality for the equations of motion describing massive spin two field in external gravitational and massless scalar dilaton fields in arbitrary spacetime dimension. From the field theoretical point of view we consider a general classical action with non-minimal couplings and find gravitational and dilaton background on which this action describes a theory consistent with the flat space limit. In the case of pure gravitational background all field components propagate causally. We show also that the massive spin two field can be consistently described in arbitrary background by means of the lagrangian representing an infinite series in the inverse mass. Within string theory we obtain equations of motion for the massive spin two field coupled to gravity from the requirement of quantum Weyl invariance of the corresponding two-dimensional sigma-model. In the lowest order in α' we demonstrate that these effective equations of motion coincide with consistent equations derived in field theory

Equations of motion for massive spin 2 field coupled to gravity

Energy Technology Data Exchange (ETDEWEB)

Buchbinder, I.L. E-mail: ilb@mail.tomsknet.ru; Gitman, D.M. E-mail: gitman@fma.if.usp.br; Krykhtin, V.A. E-mail: krykhtin@phys.dfe.tpu.edu.ru; Pershin, V.D. E-mail: pershin@ic.tsu.ru

2000-09-18

We investigate the problems of consistency and causality for the equations of motion describing massive spin two field in external gravitational and massless scalar dilaton fields in arbitrary spacetime dimension. From the field theoretical point of view we consider a general classical action with non-minimal couplings and find gravitational and dilaton background on which this action describes a theory consistent with the flat space limit. In the case of pure gravitational background all field components propagate causally. We show also that the massive spin two field can be consistently described in arbitrary background by means of the lagrangian representing an infinite series in the inverse mass. Within string theory we obtain equations of motion for the massive spin two field coupled to gravity from the requirement of quantum Weyl invariance of the corresponding two-dimensional sigma-model. In the lowest order in {alpha}' we demonstrate that these effective equations of motion coincide with consistent equations derived in field theory.

Fibonacci-like Differential Equations with a Polynomial Non-Homogeneous Part

NARCIS (Netherlands)

Asveld, P.R.J.

1989-01-01

We investigate non-homogeneous linear differential equations of the form $x''(t) + x'(t) - x(t) = p(t)$ where $p(t)$ is either a polynomial or a factorial polynomial in $t$. We express the solution of these differential equations in terms of the coefficients of $p(t)$, in the initial conditions, and

Non-Gaussian lineshapes and dynamics of time-resolved linear and nonlinear (correlation) spectra.

Science.gov (United States)

Dinpajooh, Mohammadhasan; Matyushov, Dmitry V

2014-07-17

Signatures of nonlinear and non-Gaussian dynamics in time-resolved linear and nonlinear (correlation) 2D spectra are analyzed in a model considering a linear plus quadratic dependence of the spectroscopic transition frequency on a Gaussian nuclear coordinate of the thermal bath (quadratic coupling). This new model is contrasted to the commonly assumed linear dependence of the transition frequency on the medium nuclear coordinates (linear coupling). The linear coupling model predicts equality between the Stokes shift and equilibrium correlation functions of the transition frequency and time-independent spectral width. Both predictions are often violated, and we are asking here the question of whether a nonlinear solvent response and/or non-Gaussian dynamics are required to explain these observations. We find that correlation functions of spectroscopic observables calculated in the quadratic coupling model depend on the chromophore's electronic state and the spectral width gains time dependence, all in violation of the predictions of the linear coupling models. Lineshape functions of 2D spectra are derived assuming Ornstein-Uhlenbeck dynamics of the bath nuclear modes. The model predicts asymmetry of 2D correlation plots and bending of the center line. The latter is often used to extract two-point correlation functions from 2D spectra. The dynamics of the transition frequency are non-Gaussian. However, the effect of non-Gaussian dynamics is limited to the third-order (skewness) time correlation function, without affecting the time correlation functions of higher order. The theory is tested against molecular dynamics simulations of a model polar-polarizable chromophore dissolved in a force field water.

Application of non-linear discretetime feedback regulators with assignable closed-loop dynamics

Directory of Open Access Journals (Sweden)

Dubljević Stevan

2003-01-01

Full Text Available In the present work the application of a new approach is demonstrated to a discrete-time state feedback regulator synthesis with feedback linearization and pole-placement for non-linear discrete-time systems. Under the simultaneous implementation of a non-linear coordinate transformation and a non-linear state feedback law computed through the solution of a system of non-linear functional equations, both the feedback linearization and pole-placement design objectives were accomplished. The non-linear state feedback regulator synthesis method was applied to a continuous stirred tank reactor (CSTR under non-isothermal operating conditions that exhibits steady-state multiplicity. The control objective was to regulate the reactor at the middle unstable steady state by manipulating the rate of input heat in the reactor. Simulation studies were performed to evaluate the performance of the proposed non-linear state feedback regulator, as it was shown a non-linear state feedback regulator clearly outperformed a standard linear one, especially in the presence of adverse disturbance under which linear regulation at the unstable steady state was not feasible.

Reciprocal link for a coupled Camassa–Holm type equation

International Nuclear Information System (INIS)

Li, Nianhua; Zhang, Jinshun; Wu, Lihua

2016-01-01

Highlights: • We construct a reciprocal transformation for a coupled Camassa–Holm type equation proposed by Geng and Xue. • The transformed coupled Camassa–Holm type system is a reduction of the first negative flow in a modified Drinfeld–Sokolov III hierarchy. • The Lax pair and bi-Hamiltonian structure behaviors of the coupled Camassa–Holm type equation under the reciprocal transformation are analyzed. - Abstract: A coupled Camassa–Holm type equation is linked to the first negative flow in a modified Drinfeld–Sokolov III hierarchy by a transformation of reciprocal type. Meanwhile the Lax pair and bi-Hamiltonian structure behaviors of this coupled Camassa–Holm type equation under the reciprocal transformation are analyzed.

Non-linear buckling of an FGM truncated conical shell surrounded by an elastic medium

International Nuclear Information System (INIS)

Sofiyev, A.H.; Kuruoglu, N.

2013-01-01

In this paper, the non-linear buckling of the truncated conical shell made of functionally graded materials (FGMs) surrounded by an elastic medium has been studied using the large deformation theory with von Karman–Donnell-type of kinematic non-linearity. A two-parameter foundation model (Pasternak-type) is used to describe the shell–foundation interaction. The FGM properties are assumed to vary continuously through the thickness direction. The fundamental relations, the modified Donnell type non-linear stability and compatibility equations of the FGM truncated conical shell resting on the Pasternak-type elastic foundation are derived. By using the Superposition and Galerkin methods, the non-linear stability equations for the FGM truncated conical shell is solved. Finally, influences of variations of Winkler foundation stiffness and shear subgrade modulus of the foundation, compositional profiles and shell characteristics on the dimensionless critical non-linear axial load are investigated. The present results are compared with the available data for a special case. -- Highlights: • Nonlinear buckling of FGM conical shell surrounded by elastic medium is studied. • Pasternak foundation model is used to describe the shell–foundation interaction. • Nonlinear basic equations are derived. • Problem is solved by using Superposition and Galerkin methods. • Influences of various parameters on the nonlinear critical load are investigated

Three-point phase correlations: A new measure of non-linear large-scale structure

CERN Document Server

Wolstenhulme, Richard; Obreschkow, Danail

2015-01-01

We derive an analytical expression for a novel large-scale structure observable: the line correlation function. The line correlation function, which is constructed from the three-point correlation function of the phase of the density field, is a robust statistical measure allowing the extraction of information in the non-linear and non-Gaussian regime. We show that, in perturbation theory, the line correlation is sensitive to the coupling kernel F_2, which governs the non-linear gravitational evolution of the density field. We compare our analytical expression with results from numerical simulations and find a very good agreement for separations r>20 Mpc/h. Fitting formulae for the power spectrum and the non-linear coupling kernel at small scales allow us to extend our prediction into the strongly non-linear regime. We discuss the advantages of the line correlation relative to standard statistical measures like the bispectrum. Unlike the latter, the line correlation is independent of the linear bias. Furtherm...

On preconditioning incompressible non-Newtonian flow problems

NARCIS (Netherlands)

He, X.; Neytcheva, M.; Vuik, C.

2013-01-01

This paper deals with fast and reliable numerical solution methods for the incompressible non-Newtonian Navier-Stokes equations. To handle the nonlinearity of the governing equations, the Picard and Newton methods are used to linearize these coupled partial differential equations. For space

On modulated complex non-linear dynamical systems

International Nuclear Information System (INIS)

Mahmoud, G.M.; Mohamed, A.A.; Rauh, A.

1999-01-01

This paper is concerned with the development of an approximate analytical method to investigate periodic solutions and their stability in the case of modulated non-linear dynamical systems whose equation of motion is describe. Such differential equations appear, for example, in problems of colliding particle beams in high-energy accelerators or one-mass systems with two or more degrees of freedom, e.g. rotors. The significance of periodic solutions lies on the fact that all non-periodic responses, if convergent, would approach to periodic solutions at the steady-state conditions. The example shows a good agreement between numerical and analytical results for small values of ε. The effect of the periodic modulation on the stability of the 2π-periodic solutions is discussed

The role of non-equilibrium fluxes in the relaxation processes of the linear chemical master equation.

Science.gov (United States)

de Oliveira, Luciana Renata; Bazzani, Armando; Giampieri, Enrico; Castellani, Gastone C

2014-08-14

We propose a non-equilibrium thermodynamical description in terms of the Chemical Master Equation (CME) to characterize the dynamics of a chemical cycle chain reaction among m different species. These systems can be closed or open for energy and molecules exchange with the environment, which determines how they relax to the stationary state. Closed systems reach an equilibrium state (characterized by the detailed balance condition (D.B.)), while open systems will reach a non-equilibrium steady state (NESS). The principal difference between D.B. and NESS is due to the presence of chemical fluxes. In the D.B. condition the fluxes are absent while for the NESS case, the chemical fluxes are necessary for the state maintaining. All the biological systems are characterized by their "far from equilibrium behavior," hence the NESS is a good candidate for a realistic description of the dynamical and thermodynamical properties of living organisms. In this work we consider a CME written in terms of a discrete Kolmogorov forward equation, which lead us to write explicitly the non-equilibrium chemical fluxes. For systems in NESS, we show that there is a non-conservative "external vector field" whose is linearly proportional to the chemical fluxes. We also demonstrate that the modulation of these external fields does not change their stationary distributions, which ensure us to study the same system and outline the differences in the system's behavior when it switches from the D.B. regime to NESS. We were interested to see how the non-equilibrium fluxes influence the relaxation process during the reaching of the stationary distribution. By performing analytical and numerical analysis, our central result is that the presence of the non-equilibrium chemical fluxes reduces the characteristic relaxation time with respect to the D.B. condition. Within a biochemical and biological perspective, this result can be related to the "plasticity property" of biological systems and to their

The role of non-equilibrium fluxes in the relaxation processes of the linear chemical master equation

International Nuclear Information System (INIS)

Oliveira, Luciana Renata de; Bazzani, Armando; Giampieri, Enrico; Castellani, Gastone C.

2014-01-01

We propose a non-equilibrium thermodynamical description in terms of the Chemical Master Equation (CME) to characterize the dynamics of a chemical cycle chain reaction among m different species. These systems can be closed or open for energy and molecules exchange with the environment, which determines how they relax to the stationary state. Closed systems reach an equilibrium state (characterized by the detailed balance condition (D.B.)), while open systems will reach a non-equilibrium steady state (NESS). The principal difference between D.B. and NESS is due to the presence of chemical fluxes. In the D.B. condition the fluxes are absent while for the NESS case, the chemical fluxes are necessary for the state maintaining. All the biological systems are characterized by their “far from equilibrium behavior,” hence the NESS is a good candidate for a realistic description of the dynamical and thermodynamical properties of living organisms. In this work we consider a CME written in terms of a discrete Kolmogorov forward equation, which lead us to write explicitly the non-equilibrium chemical fluxes. For systems in NESS, we show that there is a non-conservative “external vector field” whose is linearly proportional to the chemical fluxes. We also demonstrate that the modulation of these external fields does not change their stationary distributions, which ensure us to study the same system and outline the differences in the system's behavior when it switches from the D.B. regime to NESS. We were interested to see how the non-equilibrium fluxes influence the relaxation process during the reaching of the stationary distribution. By performing analytical and numerical analysis, our central result is that the presence of the non-equilibrium chemical fluxes reduces the characteristic relaxation time with respect to the D.B. condition. Within a biochemical and biological perspective, this result can be related to the “plasticity property” of biological

Asymptotic stability of linearly evolving non-stationary modes in a ...

Indian Academy of Sciences (India)

attention and it is believed to shed important light on the unresolved .... number assumption and is termed as the triple-deck theory. Having ... to analyse asymptotically the linear and weakly non-linear stability features of the station- ..... A numerical integration of equations (7–9) was implemented first to obtain the basic flow.

Minimally coupled N-particle scattering integral equations

International Nuclear Information System (INIS)

Kowalski, K.L.

1977-01-01

A concise formalism is developed which permits the efficient representation and generalization of several known techniques for deriving connected-kernel N-particle scattering integral equations. The methods of Kouri, Levin, and Tobocman and Bencze and Redish which lead to minimally coupled integral equations are of special interest. The introduction of channel coupling arrays is characterized in a general manner and the common base of this technique and that of the so-called channel coupling scheme is clarified. It is found that in the Bencze-Redish formalism a particular coupling array has a crucial function but one different from that of the arrays employed by Kouri, Levin, and Tobocman. The apparent dependence of the proof of the minimality of the Bencze-Redish integral equations upon the form of the inhomogeneous term in these equations is eliminated. This is achieved by an investigation of the full (nonminimal) Bencze-Redish kernel. It is shown that the second power of this operator is connected, a result which is needed for the full applicability of the Bencze-Redish formalism. This is used to establish the relationship between the existence of solutions to the homogeneous form of the minimal equations and eigenvalues of the full Bencze-Redish kernel

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.

Non-linear calculation of PCRV using dynamic relaxation

International Nuclear Information System (INIS)

Schnellenbach, G.

1979-01-01

A brief review is presented of a numerical method called the dynamic relaxation method for stress analysis of the concrete in prestressed concrete pressure vessels. By this method the three-dimensional elliptic differential equations of the continuum are changed into the four-dimensional hyperbolic differential equations known as wave equations. The boundary value problem of the static system is changed into an initial and boundary value problem for which a solution exists if the physical system is defined at time t=0. The effect of non-linear stress-strain behaviour of the material as well as creep and cracking are considered

Strong-coupling theory of superconductivity

International Nuclear Information System (INIS)

Rainer, D.; Sauls, J.A.

1995-01-01

The electronic properties of correlated metals with a strong electron-phonon coupling may be understood in terms of a combination of Landau''s Fermi liquid theory and the strong-coupling theory of Migdal and Eliashberg. In these lecture notes we discuss the microscopic foundations of this phenomenological Fermi-liquid model of correlated, strong-coupling metals. We formulate the basic equations of the model, which are quasiclassical transport equations that describe both equilibrium and non-equilibrium phenomena for the normal and superconducting states of a metal. Our emphasis is on superconductors close to equilibrium, for which we derive the general linear response theory. As an application we calculate the dynamical conductivity of strong-coupling superconductors. (author)

Similarity transformed coupled cluster response (ST-CCR) theory--a time-dependent similarity transformed equation-of-motion coupled cluster (STEOM-CC) approach.

Science.gov (United States)

Landau, Arie

2013-07-07

This paper presents a new method for calculating spectroscopic properties in the framework of response theory utilizing a sequence of similarity transformations (STs). The STs are preformed using the coupled cluster (CC) and Fock-space coupled cluster operators. The linear and quadratic response functions of the new similarity transformed CC response (ST-CCR) method are derived. The poles of the linear response yield excitation-energy (EE) expressions identical to the ones in the similarity transformed equation-of-motion coupled cluster (STEOM-CC) approach. ST-CCR and STEOM-CC complement each other, in analogy to the complementarity of CC response (CCR) and equation-of-motion coupled cluster (EOM-CC). ST-CCR/STEOM-CC and CCR/EOM-CC yield size-extensive and size-intensive EEs, respectively. Other electronic-properties, e.g., transition dipole strengths, are also size-extensive within ST-CCR, in contrast to STEOM-CC. Moreover, analysis suggests that in comparison with CCR, the ST-CCR expressions may be confined to a smaller subspace, however, the precise scope of the truncation can only be determined numerically. In addition, reformulation of the time-independent STEOM-CC using the same parameterization as in ST-CCR, as well as an efficient truncation scheme, is presented. The shown convergence of the time-dependent and time-independent expressions displays the completeness of the presented formalism.

Half-trek criterion for generic identifiability of linear structural equation models

NARCIS (Netherlands)

Foygel, R.; Draisma, J.; Drton, M.

2012-01-01

A linear structural equation model relates random variables of interest and corresponding Gaussian noise terms via a linear equation system. Each such model can be represented by a mixed graph in which directed edges encode the linear equations, and bidirected edges indicate possible correlations

Half-trek criterion for generic identifiability of linear structural equation models

NARCIS (Netherlands)

Foygel, R.; Draisma, J.; Drton, M.

2011-01-01

A linear structural equation model relates random variables of interest and corresponding Gaussian noise terms via a linear equation system. Each such model can be represented by a mixed graph in which directed edges encode the linear equations, and bidirected edges indicate possible correlations

Introduction to linear systems of differential equations

CERN Document Server

Adrianova, L Ya

1995-01-01

The theory of linear systems of differential equations is one of the cornerstones of the whole theory of differential equations. At its root is the concept of the Lyapunov characteristic exponent. In this book, Adrianova presents introductory material and further detailed discussions of Lyapunov exponents. She also discusses the structure of the space of solutions of linear systems. Classes of linear systems examined are from the narrowest to widest: 1)�autonomous, 2)�periodic, 3)�reducible to autonomous, 4)�nearly reducible to autonomous, 5)�regular. In addition, Adrianova considers the following: stability of linear systems and the influence of perturbations of the coefficients on the stability the criteria of uniform stability and of uniform asymptotic stability in terms of properties of the solutions several estimates of the growth rate of solutions of a linear system in terms of its coefficients How perturbations of the coefficients change all the elements of the spectrum of the system is defin...

Pattern formation due to non-linear vortex diffusion

Science.gov (United States)

Wijngaarden, Rinke J.; Surdeanu, R.; Huijbregtse, J. M.; Rector, J. H.; Dam, B.; Einfeld, J.; Wördenweber, R.; Griessen, R.

Penetration of magnetic flux in YBa 2Cu 3O 7 superconducting thin films in an external magnetic field is visualized using a magneto-optic technique. A variety of flux patterns due to non-linear vortex diffusion is observed: (1) Roughening of the flux front with scaling exponents identical to those observed in burning paper including two distinct regimes where respectively spatial disorder and temporal disorder dominate. In the latter regime Kardar-Parisi-Zhang behavior is found. (2) Fractal penetration of flux with Hausdorff dimension depending on the critical current anisotropy. (3) Penetration as ‘flux-rivers’. (4) The occurrence of commensurate and incommensurate channels in films with anti-dots as predicted in numerical simulations by Reichhardt, Olson and Nori. It is shown that most of the observed behavior is related to the non-linear diffusion of vortices by comparison with simulations of the non-linear diffusion equation appropriate for vortices.

GLOBAL LINEARIZATION OF DIFFERENTIAL EQUATIONS WITH SPECIAL STRUCTURES

Institute of Scientific and Technical Information of China (English)

无

2011-01-01

This paper introduces the global linearization of the differential equations with special structures.The function in the differential equation is unbounded.We prove that the differential equation with unbounded function can be topologically linearlized if it has a special structure.

Geometric Insight into Scalar Combination of Linear Equations

Indian Academy of Sciences (India)

... Journals; Resonance – Journal of Science Education; Volume 14; Issue 11. Geometric Insight into Scalar Combination of Linear Equations. Ranjit Konkar. Classroom Volume 14 Issue 11 November 2009 pp 1092-1097 ... Keywords. Linear algebra; linear dependence; linear combination; family of lines; family of planes.

Quantum non-local charges and absence of particle production in the two-dimensional non-linear sigma-model

International Nuclear Information System (INIS)

Luescher, M.

1977-12-01

Conserved non-local charges are shown to exist in the quantum non-linear sigma-model by a non-perturbative method. They imply the absence of particle production and the 'factorization equations' for the two particle S-matrix, which can then be calculated explicitly. (Auth.)

Local control of globally competing patterns in coupled Swift-Hohenberg equations

Science.gov (United States)

Becker, Maximilian; Frenzel, Thomas; Niedermayer, Thomas; Reichelt, Sina; Mielke, Alexander; Bär, Markus

2018-04-01

We present analytical and numerical investigations of two anti-symmetrically coupled 1D Swift-Hohenberg equations (SHEs) with cubic nonlinearities. The SHE provides a generic formulation for pattern formation at a characteristic length scale. A linear stability analysis of the homogeneous state reveals a wave instability in addition to the usual Turing instability of uncoupled SHEs. We performed weakly nonlinear analysis in the vicinity of the codimension-two point of the Turing-wave instability, resulting in a set of coupled amplitude equations for the Turing pattern as well as left- and right-traveling waves. In particular, these complex Ginzburg-Landau-type equations predict two major things: there exists a parameter regime where multiple different patterns are stable with respect to each other and that the amplitudes of different patterns interact by local mutual suppression. In consequence, different patterns can coexist in distinct spatial regions, separated by localized interfaces. We identified specific mechanisms for controlling the position of these interfaces, which distinguish what kinds of patterns the interface connects and thus allow for global pattern selection. Extensive simulations of the original SHEs confirm our results.

An efficient formulation for linear and geometric non-linear membrane elements

Directory of Open Access Journals (Sweden)

Mohammad Rezaiee-Pajand

Full Text Available Utilizing the straingradient notation process and the free formulation, an efficient way of constructing membrane elements will be proposed. This strategy can be utilized for linear and geometric non-linear problems. In the suggested formulation, the optimization constraints of insensitivity to distortion, rotational invariance and not having parasitic shear error are employed. In addition, the equilibrium equations will be established based on some constraints among the strain states. The authors' technique can easily separate the rigid body motions, and those belong to deformational motions. In this article, a novel triangular element, named SST10, is formulated. This element will be used in several plane problems having irregular mesh and complicated geometry with linear and geometrically nonlinear behavior. The numerical outcomes clearly demonstrate the efficiency of the new formulation.

Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations

Directory of Open Access Journals (Sweden)

Matt Challacombe

2014-03-01

Full Text Available A new approach to solving the Time-Dependent Self-Consistent-Field equations is developed based on the double quotient formulation of Tsiper 2001 (J. Phys. B. Dual channel, quasi-independent non-linear optimization of these quotients is found to yield convergence rates approaching those of the best case (single channel Tamm-Dancoff approximation. This formulation is variational with respect to matrix truncation, admitting linear scaling solution of the matrix-eigenvalue problem, which is demonstrated for bulk excitons in the polyphenylene vinylene oligomer and the (4,3 carbon nanotube segment.

What happens to linear properties as we move from the Klein-Gordon equation to the sine-Gordon equation

International Nuclear Information System (INIS)

Kovalyov, Mikhail

2010-01-01

In this article the sets of solutions of the sine-Gordon equation and its linearization the Klein-Gordon equation are discussed and compared. It is shown that the set of solutions of the sine-Gordon equation possesses a richer structure which partly disappears during linearization. Just like the solutions of the Klein-Gordon equation satisfy the linear superposition principle, the solutions of the sine-Gordon equation satisfy a nonlinear superposition principle.

Response statistics of rotating shaft with non-linear elastic restoring forces by path integration

Science.gov (United States)

Gaidai, Oleg; Naess, Arvid; Dimentberg, Michael

2017-07-01

Extreme statistics of random vibrations is studied for a Jeffcott rotor under uniaxial white noise excitation. Restoring force is modelled as elastic non-linear; comparison is done with linearized restoring force to see the force non-linearity effect on the response statistics. While for the linear model analytical solutions and stability conditions are available, it is not generally the case for non-linear system except for some special cases. The statistics of non-linear case is studied by applying path integration (PI) method, which is based on the Markov property of the coupled dynamic system. The Jeffcott rotor response statistics can be obtained by solving the Fokker-Planck (FP) equation of the 4D dynamic system. An efficient implementation of PI algorithm is applied, namely fast Fourier transform (FFT) is used to simulate dynamic system additive noise. The latter allows significantly reduce computational time, compared to the classical PI. Excitation is modelled as Gaussian white noise, however any kind distributed white noise can be implemented with the same PI technique. Also multidirectional Markov noise can be modelled with PI in the same way as unidirectional. PI is accelerated by using Monte Carlo (MC) estimated joint probability density function (PDF) as initial input. Symmetry of dynamic system was utilized to afford higher mesh resolution. Both internal (rotating) and external damping are included in mechanical model of the rotor. The main advantage of using PI rather than MC is that PI offers high accuracy in the probability distribution tail. The latter is of critical importance for e.g. extreme value statistics, system reliability, and first passage probability.

Stability of braneworlds with non-minimally coupled multi-scalar fields

Energy Technology Data Exchange (ETDEWEB)

Chen, Feng-Wei; Gu, Bao-Min [Lanzhou University, Institute of Theoretical Physics, Lanzhou (China); Lanzhou University, Key Laboratory for Magnetism and Magnetic Materials of the Ministry of Education, Lanzhou (China); Liu, Yu-Xiao [Lanzhou University, Research Center of Gravitation, Lanzhou (China)

2018-02-15

Linear stability of braneworld models constructed with multi-scalar fields is very different from that of single-scalar field models. It is well known that both the tensor and the scalar perturbations of the latter are stable at linear level. However, in general there is no effective method to deal with the stability problem of the scalar perturbations for braneworld models constructed with non-minimally coupled multi-scalar fields. In this work we present a systematic covariant approach to deal with the scalar perturbations. By introducing the orthonormal bases in field space and making the Kaluza-Klein decomposition, we get a set of coupled Schroedinger-like equations of the scalar perturbation modes. Using the nodal theorem, we show that the result is model-dependent. For superpotential derived brane models, the scalar perturbations are stable, but there exist normalizable scalar zero modes, which will result in unacceptable fifth force on the brane. We also use this method to analyze the f(R) braneworld model with an explicit solution and find that the scalar perturbations are stable and the scalar zero modes cannot be localized on the brane, which ensures that there is no extra long-range force and the Newtonian potential on the brane can be recovered. (orig.)

Non Linear signa models probing the string structure

International Nuclear Information System (INIS)

Abdalla, E.

1987-01-01

The introduction of a term depending on the extrinsic curvature to the string action, and related non linear sigma models defined on a symmetric space SO(D)/SO(2) x SO(d-2) is descussed . Coupling to fermions are also treated. (author) [pt

Semi-analog Monte Carlo (SMC) method for time-dependent non-linear three-dimensional heterogeneous radiative transfer problems

International Nuclear Information System (INIS)

Yun, Sung Hwan

2004-02-01

Radiative transfer is a complex phenomenon in which radiation field interacts with material. This thermal radiative transfer phenomenon is composed of two equations which are the balance equation of photons and the material energy balance equation. The two equations involve non-linearity due to the temperature and that makes the radiative transfer equation more difficult to solve. During the last several years, there have been many efforts to solve the non-linear radiative transfer problems by Monte Carlo method. Among them, it is known that Semi-Analog Monte Carlo (SMC) method developed by Ahrens and Larsen is accurate regard-less of the time step size in low temperature region. But their works are limited to one-dimensional, low temperature problems. In this thesis, we suggest some method to remove their limitations in the SMC method and apply to the more realistic problems. An initially cold problem was solved over entire temperature region by using piecewise linear interpolation of the heat capacity, while heat capacity is still fitted as a cubic curve within the lowest temperature region. If we assume the heat capacity to be linear in each temperature region, the non-linearity still remains in the radiative transfer equations. We then introduce the first-order Taylor expansion to linearize the non-linear radiative transfer equations. During the linearization procedure, absorption-reemission phenomena may be described by a conventional reemission time sampling scheme which is similar to the repetitive sampling scheme in particle transport Monte Carlo method. But this scheme causes significant stochastic errors, which necessitates many histories. Thus, we present a new reemission time sampling scheme which reduces stochastic errors by storing the information of absorption times. The results of the comparison of the two schemes show that the new scheme has less stochastic errors. Therefore, the improved SMC method is able to solve more realistic problems with

A toolbox to solve coupled systems of differential and difference equations

International Nuclear Information System (INIS)

Ablinger, Jakob; Schneider, Carsten; Bluemlein, Johannes; Freitas, Abilio de

2016-01-01

We present algorithms to solve coupled systems of linear differential equations, arising in the calculation of massive Feynman diagrams with local operator insertions at 3-loop order, which do not request special choices of bases. Here we assume that the desired solution has a power series representation and we seek for the coefficients in closed form. In particular, if the coefficients depend on a small parameter ε (the dimensional parameter), we assume that the coefficients themselves can be expanded in formal Laurent series w.r.t. ε and we try to compute the first terms in closed form. More precisely, we have a decision algorithm which solves the following problem: if the terms can be represented by an indefinite nested hypergeometric sum expression (covering as special cases the harmonic sums, cyclotomic sums, generalized harmonic sums or nested binomial sums), then we can calculate them. If the algorithm fails, we obtain a proof that the terms cannot be represented by the class of indefinite nested hypergeometric sum expressions. Internally, this problem is reduced by holonomic closure properties to solving a coupled system of linear difference equations. The underlying method in this setting relies on decoupling algorithms, difference ring algorithms and recurrence solving. We demonstrate by a concrete example how this algorithm can be applied with the new Mathematica package SolveCoupledSystem which is based on the packages Sigma, HarmonicSums and OreSys. In all applications the representation in x-space is obtained as an iterated integral representation over general alphabets, generalizing Poincare iterated integrals.

A toolbox to solve coupled systems of differential and difference equations

Energy Technology Data Exchange (ETDEWEB)

Ablinger, Jakob; Schneider, Carsten [Linz Univ. (Austria). Research Inst. for Symbolic Computation (RISC); Bluemlein, Johannes; Freitas, Abilio de [DESY Zeuthen (Germany)

2016-01-15

We present algorithms to solve coupled systems of linear differential equations, arising in the calculation of massive Feynman diagrams with local operator insertions at 3-loop order, which do not request special choices of bases. Here we assume that the desired solution has a power series representation and we seek for the coefficients in closed form. In particular, if the coefficients depend on a small parameter ε (the dimensional parameter), we assume that the coefficients themselves can be expanded in formal Laurent series w.r.t. ε and we try to compute the first terms in closed form. More precisely, we have a decision algorithm which solves the following problem: if the terms can be represented by an indefinite nested hypergeometric sum expression (covering as special cases the harmonic sums, cyclotomic sums, generalized harmonic sums or nested binomial sums), then we can calculate them. If the algorithm fails, we obtain a proof that the terms cannot be represented by the class of indefinite nested hypergeometric sum expressions. Internally, this problem is reduced by holonomic closure properties to solving a coupled system of linear difference equations. The underlying method in this setting relies on decoupling algorithms, difference ring algorithms and recurrence solving. We demonstrate by a concrete example how this algorithm can be applied with the new Mathematica package SolveCoupledSystem which is based on the packages Sigma, HarmonicSums and OreSys. In all applications the representation in x-space is obtained as an iterated integral representation over general alphabets, generalizing Poincare iterated integrals.

Emmy Noether and Linear Evolution Equations

Directory of Open Access Journals (Sweden)

P. G. L. Leach

2013-01-01

Full Text Available Noether’s Theorem relates the Action Integral of a Lagrangian with symmetries which leave it invariant and the first integrals consequent upon the variational principle and the existence of the symmetries. These each have an equivalent in the Schrödinger Equation corresponding to the Lagrangian and by extension to linear evolution equations in general. The implications of these connections are investigated.

Non-Linear Transmission Line (NLTL) Microwave Source Lecture Notes the United States Particle Accelerator School

Energy Technology Data Exchange (ETDEWEB)

Russell, Steven J. [Los Alamos National Laboratory; Carlsten, Bruce E. [Los Alamos National Laboratory

2012-06-26

We will quickly go through the history of the non-linear transmission lines (NLTLs). We will describe how they work, how they are modeled and how they are designed. Note that the field of high power, NLTL microwave sources is still under development, so this is just a snap shot of their current state. Topics discussed are: (1) Introduction to solitons and the KdV equation; (2) The lumped element non-linear transmission line; (3) Solution of the KdV equation; (4) Non-linear transmission lines at microwave frequencies; (5) Numerical methods for NLTL analysis; (6) Unipolar versus bipolar input; (7) High power NLTL pioneers; (8) Resistive versus reactive load; (9) Non-lineaer dielectrics; and (10) Effect of losses.

Simplified Linear Equation Solvers users manual

Energy Technology Data Exchange (ETDEWEB)

Gropp, W. [Argonne National Lab., IL (United States); Smith, B. [California Univ., Los Angeles, CA (United States)

1993-02-01

The solution of large sparse systems of linear equations is at the heart of many algorithms in scientific computing. The SLES package is a set of easy-to-use yet powerful and extensible routines for solving large sparse linear systems. The design of the package allows new techniques to be used in existing applications without any source code changes in the applications.

A linearizing transformation for the Korteweg-de Vries equation; generalizations to higher-dimensional nonlinear partial differential equations

NARCIS (Netherlands)

Dorren, H.J.S.

1998-01-01

It is shown that the Korteweg–de Vries (KdV) equation can be transformed into an ordinary linear partial differential equation in the wave number domain. Explicit solutions of the KdV equation can be obtained by subsequently solving this linear differential equation and by applying a cascade of

Construction of a Roe linearization for the ideal MHD equations

International Nuclear Information System (INIS)

Cargo, P.; Gallice, G.; Raviart, P.A.

1996-01-01

In [3], Munz has constructed a Roe linearization for the equations of gas dynamics in Lagrangian coordinates. We extend this construction to the case of the ideal magnetohydrodynamics equations again in Lagrangian coordinates. As a consequence we obtain a Roe linearization for the MHD equations in Eulerian coordinates. (author)

Non-minimally coupled quintessence dark energy model with a cubic galileon term: a dynamical system analysis

Science.gov (United States)

Bhattacharya, Somnath; Mukherjee, Pradip; Roy, Amit Singha; Saha, Anirban

2018-03-01

We consider a scalar field which is generally non-minimally coupled to gravity and has a characteristic cubic Galilean-like term and a generic self-interaction, as a candidate of a Dark Energy model. The system is dynamically analyzed and novel fixed points with perturbative stability are demonstrated. Evolution of the system is numerically studied near a novel fixed point which owes its existence to the Galileon character of the model. It turns out that demanding the stability of this novel fixed point puts a strong restriction on the allowed non-minimal coupling and the choice of the self-interaction. The evolution of the equation of state parameter is studied, which shows that our model predicts an accelerated universe throughout and the phantom limit is only approached closely but never crossed. Our result thus extends the findings of Coley, Dynamical systems and cosmology. Kluwer Academic Publishers, Boston (2013) for more general NMC than linear and quadratic couplings.

Variational linear algebraic equations method

International Nuclear Information System (INIS)

Moiseiwitsch, B.L.

1982-01-01

A modification of the linear algebraic equations method is described which ensures a variational bound on the phaseshifts for potentials having a definite sign at all points. The method is illustrated by the elastic scattering of s-wave electrons by the static field of atomic hydrogen. (author)

Dynamics of interacting qubits coupled to a common bath: Non-Markovian quantum-state-diffusion approach

International Nuclear Information System (INIS)

Zhao Xinyu; Jing Jun; Corn, Brittany; Yu Ting

2011-01-01

Non-Markovian dynamics is studied for two interacting qubits strongly coupled to a dissipative bosonic environment. We derive a non-Markovian quantum-state-diffusion (QSD) equation for the coupled two-qubit system without any approximations, and in particular, without the Markov approximation. As an application and illustration of our derived time-local QSD equation, we investigate the temporal behavior of quantum coherence dynamics. In particular, we find a strongly non-Markovian regime where entanglement generation is significantly modulated by the environmental memory. Additionally, we study residual entanglement in the steady state by analyzing the steady-state solution of the QSD equation. Finally, we discuss an approximate QSD equation.

Hierarchies of multi-component mKP equations and theirs integrable couplings

International Nuclear Information System (INIS)

Ji Jie; Yao Yuqin; Zhu Fubo; Chen Dengyuan

2008-01-01

First, a new multi-component modified Kadomtsev-Petviashvill (mKP) spectral problem is constructed by k-constraint imposed on a general pseudo-differential operator. Then, two hierarchies of multi-component mKP equations are derived, including positive non-isospectral mKP hierarchy and negative non-isospectral mKP hierarchy. Moreover, new integrable couplings of the resulting mKP soliton hierarchies are constructed by enlarging the associated matrix spectral problem

Rational approximations to solutions of linear differential equations.

Science.gov (United States)

Chudnovsky, D V; Chudnovsky, G V

1983-08-01

Rational approximations of Padé and Padé type to solutions of differential equations are considered. One of the main results is a theorem stating that a simultaneous approximation to arbitrary solutions of linear differential equations over C(x) cannot be "better" than trivial ones implied by the Dirichlet box principle. This constitutes, in particular, the solution in the linear case of Kolchin's problem that the "Roth's theorem" holds for arbitrary solutions of algebraic differential equations. Complete effective proofs for several valuations are presented based on the Wronskian methods and graded subrings of Picard-Vessiot extensions.

Non linear characterisation of optical components of a high power laser chain

International Nuclear Information System (INIS)

Santran, Stephane

2000-01-01

This work concerns the realisation of non linear properties measurement prototypes in glasses in the near infrared and in the visible range. The various devices are time resolved colinear pump probe experiments in which the non linear susceptibility is deduced by the probe beam intensity variations induced by the pump probe coupled in the material. The sensitivity of these experiments allows us to observe unexpected variations, greater than 30%, of several fused silica non linear indexes. As well, this allow us to analyse the origin of the promising oxide glasses non linearity for all optical applications and to understand an d measure non linear processes in the two photons photodiodes. Finally, an original structure for the non linear index measurement in non degenerated configuration by a probe pulse phase measurement approach with a Sagnac interferometer is demonstrated and analysed. (author) [fr

Inflation in non-minimal matter-curvature coupling theories

Energy Technology Data Exchange (ETDEWEB)

Gomes, C.; Bertolami, O. [Departamento de Física e Astronomia and Centro de Física do Porto, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre s/n, 4169-007 Porto (Portugal); Rosa, J.G., E-mail: claudio.gomes@fc.up.pt, E-mail: joao.rosa@ua.pt, E-mail: orfeu.bertolami@fc.up.pt [Departamento de Física da Universidade de Aveiro and CIDMA, Campus de Santiago, 3810-183 Aveiro (Portugal)

2017-06-01

We study inflationary scenarios driven by a scalar field in the presence of a non-minimal coupling between matter and curvature. We show that the Friedmann equation can be significantly modified when the energy density during inflation exceeds a critical value determined by the non-minimal coupling, which in turn may considerably modify the spectrum of primordial perturbations and the inflationary dynamics. In particular, we show that these models are characterised by a consistency relation between the tensor-to-scalar ratio and the tensor spectral index that can differ significantly from the predictions of general relativity. We also give examples of observational predictions for some of the most commonly considered potentials and use the results of the Planck collaboration to set limits on the scale of the non-minimal coupling.

Local Fractional Laplace Variational Iteration Method for Solving Linear Partial Differential Equations with Local Fractional Derivative

Directory of Open Access Journals (Sweden)

Ai-Min Yang

2014-01-01

Full Text Available The local fractional Laplace variational iteration method was applied to solve the linear local fractional partial differential equations. The local fractional Laplace variational iteration method is coupled by the local fractional variational iteration method and Laplace transform. The nondifferentiable approximate solutions are obtained and their graphs are also shown.

New classical r-matrices from integrable non-linear sigma-models

International Nuclear Information System (INIS)

Laartz, J.; Bordemann, M.; Forger, M.; Schaper, U.

1993-01-01

Non-linear sigma models on Riemannian symmetric spaces constitute the most general class of classical non-linear sigma models which are known to be integrable. Using the current algebra structure of these models their canonical structure is analyzed and it is shown that their non-ultralocal fundamental Poisson bracket relation is governed by a field dependent non antisymmetric r-matrix obeying a dynamical Yang Baxter equation. The fundamental Poisson bracket relations and the r-matrix are derived explicitly and a new kind of algebra is found that is supposed to replace the classical Yang Baxter algebra governing the canonical structure of ultralocal models. (Author) 9 refs

Asymptotic properties for half-linear difference equations

Czech Academy of Sciences Publication Activity Database

Cecchi, M.; Došlá, Z.; Marini, M.; Vrkoč, Ivo

2006-01-01

Roč. 131, č. 4 (2006), s. 347-363 ISSN 0862-7959 R&D Projects: GA ČR(CZ) GA201/04/0580 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear second order difference equation * nonoscillatory solutions * Riccati difference equation Subject RIV: BA - General Mathematics

On the structure on non-local conservation laws in the two-dimensional non-linear sigma-model

International Nuclear Information System (INIS)

Zamolodchikov, Al.B.

1978-01-01

The non-local conserved charges are supposed to satisfy a special multiplicative law in the space of asymptotic states of the non-linear sigma-model. This supposition leads to factorization equations for two-particle scattering matrix elements and determines to some extent the action of these charges in the asymptotic space. Their conservation turns out to be consistent with the factorized S-matrix of the non-linear sigma-model. It is shown also that the factorized sine-Gordon S-matrix is consistent with a similar family of conservation laws

Linear and non-linear Modified Gravity forecasts with future surveys

Science.gov (United States)

Casas, Santiago; Kunz, Martin; Martinelli, Matteo; Pettorino, Valeria

2017-12-01

Modified Gravity theories generally affect the Poisson equation and the gravitational slip in an observable way, that can be parameterized by two generic functions (η and μ) of time and space. We bin their time dependence in redshift and present forecasts on each bin for future surveys like Euclid. We consider both Galaxy Clustering and Weak Lensing surveys, showing the impact of the non-linear regime, with two different semi-analytical approximations. In addition to these future observables, we use a prior covariance matrix derived from the Planck observations of the Cosmic Microwave Background. In this work we neglect the information from the cross correlation of these observables, and treat them as independent. Our results show that η and μ in different redshift bins are significantly correlated, but including non-linear scales reduces or even eliminates the correlation, breaking the degeneracy between Modified Gravity parameters and the overall amplitude of the matter power spectrum. We further apply a Zero-phase Component Analysis and identify which combinations of the Modified Gravity parameter amplitudes, in different redshift bins, are best constrained by future surveys. We extend the analysis to two particular parameterizations of μ and η and consider, in addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future surveys will be able to constrain the current values of η and μ at the 2-5% level when using only linear scales (wavevector k < 0 . 15 h/Mpc), depending on the specific time parameterization; sensitivity improves to about 1% when non-linearities are included.

Some mathematical problems in non-linear Physics

International Nuclear Information System (INIS)

1983-01-01

The main results contained in this report are the following: I) A general analysis of non-autonomous conserved densities for simple linear evolution systems. II) Partial differential systems within a wide class are converted into Lagrange an form. III) Rigorous criteria for existence of integrating factor matrices. IV) Isolation of all third-order evolution equations with high order symmetries and conservation laws. (Author) 3 refs

Coupled Higgs field equation and Hamiltonian amplitude equation ...

Indian Academy of Sciences (India)

Home; Journals; Pramana – Journal of Physics; Volume 79; Issue 1. Coupled Higgs field equation and ... School of Mathematics and Computer Applications, Thapar University, Patiala 147 004, India; Department of Mathematics, Jaypee University of Information Technology, Waknaghat, Distt. Solan 173 234, India ...

Stability of generalized Runge-Kutta methods for stiff kinetics coupled differential equations

International Nuclear Information System (INIS)

Aboanber, A E

2006-01-01

A stability and efficiency improved class of generalized Runge-Kutta methods of order 4 are developed for the numerical solution of stiff system kinetics equations for linear and/or nonlinear coupled differential equations. The determination of the coefficients required by the method is precisely obtained from the so-called equations of condition which in turn are derived by an approach based on Butcher series. Since the equations of condition are fewer in number, free parameters can be chosen for optimizing any desired feature of the process. A further related coefficient set with different values of these parameters and the region of absolute stability of the method have been introduced. In addition, the A(α) stability properties of the method are investigated. Implementing the method in a personal computer estimated the accuracy and speed of calculations and verified the good performances of the proposed new schemes for several sample problems of the stiff system point kinetics equations with reactivity feedback

Non-Linear Dynamics of Saturn’s Rings

Science.gov (United States)

Esposito, Larry W.

2015-11-01

Non-linear processes can explain why Saturn’s rings are so active and dynamic. Ring systems differ from simple linear systems in two significant ways: 1. They are systems of granular material: where particle-to-particle collisions dominate; thus a kinetic, not a fluid description needed. We find that stresses are strikingly inhomogeneous and fluctuations are large compared to equilibrium. 2. They are strongly forced by resonances: which drive a non-linear response, pushing the system across thresholds that lead to persistent states.Some of this non-linearity is captured in a simple Predator-Prey Model: Periodic forcing from the moon causes streamline crowding; This damps the relative velocity, and allows aggregates to grow. About a quarter phase later, the aggregates stir the system to higher relative velocity and the limit cycle repeats each orbit.Summary of Halo Results: A predator-prey model for ring dynamics produces transient structures like ‘straw’ that can explain the halo structure and spectroscopy: This requires energetic collisions (v ≈ 10m/sec, with throw distances about 200km, implying objects of scale R ≈ 20km).Transform to Duffing Eqn : With the coordinate transformation, z = M2/3, the Predator-Prey equations can be combined to form a single second-order differential equation with harmonic resonance forcing.Ring dynamics and history implications: Moon-triggered clumping at perturbed regions in Saturn’s rings creates both high velocity dispersion and large aggregates at these distances, explaining both small and large particles observed there. We calculate the stationary size distribution using a cell-to-cell mapping procedure that converts the phase-plane trajectories to a Markov chain. Approximating the Markov chain as an asymmetric random walk with reflecting boundaries allows us to determine the power law index from results of numerical simulations in the tidal environment surrounding Saturn. Aggregates can explain many dynamic aspects

Students' errors in solving linear equation word problems: Case ...

African Journals Online (AJOL)

The study examined errors students make in solving linear equation word problems with a view to expose the nature of these errors and to make suggestions for classroom teaching. A diagnostic test comprising 10 linear equation word problems, was administered to a sample (n=130) of senior high school first year Home ...

Scilab software as an alternative low-cost computing in solving the linear equations problem

Science.gov (United States)

Agus, Fahrul; Haviluddin

2017-02-01

Numerical computation packages are widely used both in teaching and research. These packages consist of license (proprietary) and open source software (non-proprietary). One of the reasons to use the package is a complexity of mathematics function (i.e., linear problems). Also, number of variables in a linear or non-linear function has been increased. The aim of this paper was to reflect on key aspects related to the method, didactics and creative praxis in the teaching of linear equations in higher education. If implemented, it could be contribute to a better learning in mathematics area (i.e., solving simultaneous linear equations) that essential for future engineers. The focus of this study was to introduce an additional numerical computation package of Scilab as an alternative low-cost computing programming. In this paper, Scilab software was proposed some activities that related to the mathematical models. In this experiment, four numerical methods such as Gaussian Elimination, Gauss-Jordan, Inverse Matrix, and Lower-Upper Decomposition (LU) have been implemented. The results of this study showed that a routine or procedure in numerical methods have been created and explored by using Scilab procedures. Then, the routine of numerical method that could be as a teaching material course has exploited.

New Equating Methods and Their Relationships with Levine Observed Score Linear Equating under the Kernel Equating Framework

Science.gov (United States)

Chen, Haiwen; Holland, Paul

2010-01-01

In this paper, we develop a new curvilinear equating for the nonequivalent groups with anchor test (NEAT) design under the assumption of the classical test theory model, that we name curvilinear Levine observed score equating. In fact, by applying both the kernel equating framework and the mean preserving linear transformation of…

Solution of linear and nonlinear matrix systems. Application to a nonlinear diffusion equation

International Nuclear Information System (INIS)

Bonnet, M.; Meurant, G.

1978-01-01

Different methods of solution of linear and nonlinear algebraic systems are applied to the nonlinear system obtained by discretizing a nonlinear diffusion equation. For linear systems, methods in general use of alternating directions type or Gauss Seidel's methods are compared to more recent ones of the type of generalized conjugate gradient; the superiority of the latter is shown by numerical examples. For nonlinear systems, a method on nonlinear conjugate gradient is studied as also Newton's method and some of its variants. It should be noted, however that Newton's method is found to be more efficient when coupled with a good method for solution of the linear system. To conclude, such methods are used to solve a nonlinear diffusion problem and the numerical results obtained are to be compared [fr

Solution of linear and nonlinear matrix systems. Application to a nonlinear diffusion equation

International Nuclear Information System (INIS)

Bonnet, M.; Meurant, G.

1978-01-01

The object of this study is to compare different methods of solving linear and nonlinear algebraic systems and to apply them to the nonlinear system obtained by discretizing a nonlinear diffusion equation. For linear systems the conventional methods of alternating direction type or Gauss Seidel's methods are compared to more recent ones of the type of generalized conjugate gradient; the superiority of the latter is shown by numerical examples. For nonlinear systems, a method of nonlinear conjugate gradient is studied together with Newton's method and some of its variants. It should be noted, however, that Newton's method is found to be more efficient when coupled with a good method for solving the linear system. As a conclusion, these methods are used to solve a nonlinear diffusion problem and the numerical results obtained are compared [fr

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)

Periodic feedback stabilization for linear periodic evolution equations

CERN Document Server

Wang, Gengsheng

2016-01-01

This book introduces a number of recent advances regarding periodic feedback stabilization for linear and time periodic evolution equations. First, it presents selected connections between linear quadratic optimal control theory and feedback stabilization theory for linear periodic evolution equations. Secondly, it identifies several criteria for the periodic feedback stabilization from the perspective of geometry, algebra and analyses respectively. Next, it describes several ways to design periodic feedback laws. Lastly, the book introduces readers to key methods for designing the control machines. Given its coverage and scope, it offers a helpful guide for graduate students and researchers in the areas of control theory and applied mathematics.

General treatment of a non-linear gauge condition

International Nuclear Information System (INIS)

Malleville, C.

1982-06-01

A non linear gauge condition is presented in the frame of a non abelian gauge theory broken with the Higgs mechanism. It is shown that this condition already introduced for the standard SU(2) x U(1) model can be generalized for any gauge model with the same type of simplification, namely the suppression of any coupling of the form: massless gauge boson, massive gauge boson, unphysical Higgs [fr

Coupling of linearized gravity to nonrelativistic test particles: Dynamics in the general laboratory frame

International Nuclear Information System (INIS)

Speliotopoulos, A.D.; Chiao, Raymond Y.

2004-01-01

The coupling of gravity to matter is explored in the linearized gravity limit. The usual derivation of gravity-matter couplings within the quantum-field-theoretic framework is reviewed. A number of inconsistencies between this derivation of the couplings and the known results of tidal effects on test particles according to classical general relativity are pointed out. As a step towards resolving these inconsistencies, a general laboratory frame fixed on the worldline of an observer is constructed. In this frame, the dynamics of nonrelativistic test particles in the linearized gravity limit is studied, and their Hamiltonian dynamics is derived. It is shown that for stationary metrics this Hamiltonian reduces to the usual Hamiltonian for nonrelativistic particles undergoing geodesic motion. For nonstationary metrics with long-wavelength gravitational waves present (GWs), it reduces to the Hamiltonian for a nonrelativistic particle undergoing geodesic deviation motion. Arbitrary-wavelength GWs couple to the test particle through a vector-potential-like field N a , the net result of the tidal forces that the GW induces in the system, namely, a local velocity field on the system induced by tidal effects, as seen by an observer in the general laboratory frame. Effective electric and magnetic fields, which are related to the electric and magnetic parts of the Weyl tensor, are constructed from N a that obey equations of the same form as Maxwell's equations. A gedankin gravitational Aharonov-Bohm-type experiment using N a to measure the interference of quantum test particles is presented

Global non-linear effect of temperature on economic production.

Science.gov (United States)

Burke, Marshall; Hsiang, Solomon M; Miguel, Edward

2015-11-12

Growing evidence demonstrates that climatic conditions can have a profound impact on the functioning of modern human societies, but effects on economic activity appear inconsistent. Fundamental productive elements of modern economies, such as workers and crops, exhibit highly non-linear responses to local temperature even in wealthy countries. In contrast, aggregate macroeconomic productivity of entire wealthy countries is reported not to respond to temperature, while poor countries respond only linearly. Resolving this conflict between micro and macro observations is critical to understanding the role of wealth in coupled human-natural systems and to anticipating the global impact of climate change. Here we unify these seemingly contradictory results by accounting for non-linearity at the macro scale. We show that overall economic productivity is non-linear in temperature for all countries, with productivity peaking at an annual average temperature of 13 °C and declining strongly at higher temperatures. The relationship is globally generalizable, unchanged since 1960, and apparent for agricultural and non-agricultural activity in both rich and poor countries. These results provide the first evidence that economic activity in all regions is coupled to the global climate and establish a new empirical foundation for modelling economic loss in response to climate change, with important implications. If future adaptation mimics past adaptation, unmitigated warming is expected to reshape the global economy by reducing average global incomes roughly 23% by 2100 and widening global income inequality, relative to scenarios without climate change. In contrast to prior estimates, expected global losses are approximately linear in global mean temperature, with median losses many times larger than leading models indicate.

Global non-linear effect of temperature on economic production

Science.gov (United States)

Burke, Marshall; Hsiang, Solomon M.; Miguel, Edward

2015-11-01

Growing evidence demonstrates that climatic conditions can have a profound impact on the functioning of modern human societies, but effects on economic activity appear inconsistent. Fundamental productive elements of modern economies, such as workers and crops, exhibit highly non-linear responses to local temperature even in wealthy countries. In contrast, aggregate macroeconomic productivity of entire wealthy countries is reported not to respond to temperature, while poor countries respond only linearly. Resolving this conflict between micro and macro observations is critical to understanding the role of wealth in coupled human-natural systems and to anticipating the global impact of climate change. Here we unify these seemingly contradictory results by accounting for non-linearity at the macro scale. We show that overall economic productivity is non-linear in temperature for all countries, with productivity peaking at an annual average temperature of 13 °C and declining strongly at higher temperatures. The relationship is globally generalizable, unchanged since 1960, and apparent for agricultural and non-agricultural activity in both rich and poor countries. These results provide the first evidence that economic activity in all regions is coupled to the global climate and establish a new empirical foundation for modelling economic loss in response to climate change, with important implications. If future adaptation mimics past adaptation, unmitigated warming is expected to reshape the global economy by reducing average global incomes roughly 23% by 2100 and widening global income inequality, relative to scenarios without climate change. In contrast to prior estimates, expected global losses are approximately linear in global mean temperature, with median losses many times larger than leading models indicate.

Non-linear numerical studies of the tearing mode

International Nuclear Information System (INIS)

Schnack, D.D. Jr.

1978-01-01

A non-linear, time dependent, hydromagnetic model is developed and applied to the tearing mode, one of a class of instabilities which can occur in a magnetically confined plasma when the constraint of infinite conductivity is relaxed. The model is based on the eight partial differential equations of resistive magnetohydrodynamics (MHD). The equations are expressed as a set of conservation laws which conserves magnetic flux, momentum, mass, and total energy. These equations are then written in general, orthogonal, curvilinear coordinates in two space dimensions, so that the model can readily be applied to a variety of geometries. No assumption about the ordering of terms is made. The resulting equations are then solved by the method of finite differences on an Eulerian mesh. The model is applied to several geometries

A tightly coupled non-equilibrium model for inductively coupled radio-frequency plasmas

International Nuclear Information System (INIS)

Munafò, A.; Alfuhaid, S. A.; Panesi, M.; Cambier, J.-L.

2015-01-01

The objective of the present work is the development of a tightly coupled magneto-hydrodynamic model for inductively coupled radio-frequency plasmas. Non Local Thermodynamic Equilibrium (NLTE) effects are described based on a hybrid State-to-State approach. A multi-temperature formulation is used to account for thermal non-equilibrium between translation of heavy-particles and vibration of molecules. Excited electronic states of atoms are instead treated as separate pseudo-species, allowing for non-Boltzmann distributions of their populations. Free-electrons are assumed Maxwellian at their own temperature. The governing equations for the electro-magnetic field and the gas properties (e.g., chemical composition and temperatures) are written as a coupled system of time-dependent conservation laws. Steady-state solutions are obtained by means of an implicit Finite Volume method. The results obtained in both LTE and NLTE conditions over a broad spectrum of operating conditions demonstrate the robustness of the proposed coupled numerical method. The analysis of chemical composition and temperature distributions along the torch radius shows that: (i) the use of the LTE assumption may lead to an inaccurate prediction of the thermo-chemical state of the gas, and (ii) non-equilibrium phenomena play a significant role close the walls, due to the combined effects of Ohmic heating and macroscopic gradients

Numerical resolution of Navier-Stokes equations coupled to the heat equation

International Nuclear Information System (INIS)

Zenouda, Jean-Claude

1970-08-01

The author proves a uniqueness theorem for the time dependent Navier-Stokes equations coupled with heat flow in the two-dimensional case. He studies stability and convergence of several finite - difference schemes to solve these equations. Numerical experiments are done in the case of a square domain. (author) [fr

A Comparison between Linear IRT Observed-Score Equating and Levine Observed-Score Equating under the Generalized Kernel Equating Framework

Science.gov (United States)

Chen, Haiwen

2012-01-01

In this article, linear item response theory (IRT) observed-score equating is compared under a generalized kernel equating framework with Levine observed-score equating for nonequivalent groups with anchor test design. Interestingly, these two equating methods are closely related despite being based on different methodologies. Specifically, when…

Dual exponential polynomials and linear differential equations

Science.gov (United States)

Wen, Zhi-Tao; Gundersen, Gary G.; Heittokangas, Janne

2018-01-01

We study linear differential equations with exponential polynomial coefficients, where exactly one coefficient is of order greater than all the others. The main result shows that a nontrivial exponential polynomial solution of such an equation has a certain dual relationship with the maximum order coefficient. Several examples illustrate our results and exhibit possibilities that can occur.

New matrix bounds and iterative algorithms for the discrete coupled algebraic Riccati equation

Science.gov (United States)

Liu, Jianzhou; Wang, Li; Zhang, Juan

2017-11-01

The discrete coupled algebraic Riccati equation (DCARE) has wide applications in control theory and linear system. In general, for the DCARE, one discusses every term of the coupled term, respectively. In this paper, we consider the coupled term as a whole, which is different from the recent results. When applying eigenvalue inequalities to discuss the coupled term, our method has less error. In terms of the properties of special matrices and eigenvalue inequalities, we propose several upper and lower matrix bounds for the solution of DCARE. Further, we discuss the iterative algorithms for the solution of the DCARE. In the fixed point iterative algorithms, the scope of Lipschitz factor is wider than the recent results. Finally, we offer corresponding numerical examples to illustrate the effectiveness of the derived results.

Supersymmetric Yang-Mills fields as an integrable system and connections with other non-linear systems

International Nuclear Information System (INIS)

Chau, L.L.

1983-01-01

Integrable properties, i.e., existence of linear systems, infinite number of conservation laws, Reimann-Hilbert transforms, affine Lie algebra of Kac-Moody, and Bianchi-Baecklund transformation, are discussed for the constraint equations of the supersymmetric Yang-Mills fields. For N greater than or equal to 3 these constraint equations give equations of motion of the fields. These equations of motion reduce to the ordinary Yang-Mills equations as the spinor and scalar fields are eliminated. These understandings provide a possible method to solve the full Yang-Mills equations. Connections with other non-linear systems are also discussed. 53 references

Microscopic model for the non-linear fluctuating hydrodynamic of 4 He superfluid helium deduced by maximum entropy method

International Nuclear Information System (INIS)

Alvarez R, J.T.

1998-01-01

This thesis presents a microscopic model for the non-linear fluctuating hydrodynamic of superfluid helium ( 4 He), model developed by means of the Maximum Entropy Method (Maxent). In the chapter 1, it is demonstrated the necessity to developing a microscopic model for the fluctuating hydrodynamic of the superfluid helium, starting from to show a brief overview of the theories and experiments developed in order to explain the behavior of the superfluid helium. On the other hand, it is presented the Morozov heuristic method for the construction of the non-linear hydrodynamic fluctuating of simple fluid. Method that will be generalized for the construction of the non-linear fluctuating hydrodynamic of the superfluid helium. Besides, it is presented a brief summary of the content of the thesis. In the chapter 2, it is reproduced the construction of a Generalized Fokker-Planck equation, (GFP), for a distribution function associated with the coarse grained variables. Function defined with aid of a nonequilibrium statistical operator ρhut FP that is evaluated as Wigneris function through ρ CG obtained by Maxent. Later this equation of GFP is reduced to a non-linear local FP equation from considering a slow and Markov process in the coarse grained variables. In this equation appears a matrix D mn defined with a nonequilibrium coarse grained statistical operator ρhut CG , matrix elements are used in the construction of the non-linear fluctuating hydrodynamics equations of the superfluid helium. In the chapter 3, the Lagrange multipliers are evaluated for to determine ρhut CG by means of the local equilibrium statistical operator ρhut l -tilde with the hypothesis that the system presents small fluctuations. Also are determined the currents associated with the coarse grained variables and furthermore are evaluated the matrix elements D mn but with aid of a quasi equilibrium statistical operator ρhut qe instead of the local equilibrium operator ρhut l -tilde. Matrix

Non-linear elastic thermal stress analysis with phase changes

International Nuclear Information System (INIS)

Amada, S.; Yang, W.H.

1978-01-01

The non-linear elastic, thermal stress analysis with temperature induced phase changes in the materials is presented. An infinite plate (or body) with a circular hole (or tunnel) is subjected to a thermal loading on its inner surface. The peak temperature around the hole reaches beyond the melting point of the material. The non-linear diffusion equation is solved numerically using the finite difference method. The material properties change rapidly at temperatures where the change of crystal structures and solid-liquid transition occur. The elastic stresses induced by the transient non-homogeneous temperature distribution are calculated. The stresses change remarkably when the phase changes occur and there are residual stresses remaining in the plate after one cycle of thermal loading. (Auth.)

Fast simulation of non-linear pulsed ultrasound fields using an angular spectrum approach

DEFF Research Database (Denmark)

Du, Yigang; Jensen, Jørgen Arendt

2013-01-01

A fast non-linear pulsed ultrasound field simulation is presented. It is implemented based on an angular spectrum approach (ASA), which analytically solves the non-linear wave equation. The ASA solution to the Westervelt equation is derived in detail. The calculation speed is significantly...... increased compared to a numerical solution using an operator splitting method (OSM). The ASA has been modified and extended to pulsed non-linear ultrasound fields in combination with Field II, where any array transducer with arbitrary geometry, excitation, focusing and apodization can be simulated...... with a center frequency of 5 MHz. The speed is increased approximately by a factor of 140 and the calculation time is 12 min with a standard PC, when simulating the second harmonic pulse at the focal point. For the second harmonic point spread function the full width error is 1.5% at 6 dB and 6.4% at 12 d...

Dynamics and non-equilibrium steady state in a system of coupled harmonic oscillators

Energy Technology Data Exchange (ETDEWEB)

Ghesquière, Anne, E-mail: Anne.Ghesquiere@nithep.ac.za; Sinayskiy, Ilya, E-mail: sinayskiy@ukzn.ac.za; Petruccione, Francesco, E-mail: petruccione@ukzn.ac.za

2013-10-15

A system of two coupled oscillators, each of them coupled to an independent reservoir, is analysed. The analytical solution of the non-rotating wave master equation is obtained in the high-temperature and weak coupling limits. No thermal entanglement is found in the high-temperature limit. In the weak coupling limit the system converges to an entangled non-equilibrium steady state. A critical temperature for the appearance of quantum correlations is found.

Soliton solutions of coupled nonlinear Klein-Gordon equations

International Nuclear Information System (INIS)

Alagesan, T.; Chung, Y.; Nakkeeran, K.

2004-01-01

The coupled nonlinear Klein-Gordon equations are analyzed for their integrability properties in a systematic manner through Painleve test. From the Painleve test, by truncating the Laurent series at the constant level term, the Hirota bilinear form is identified, from which one-soliton solutions are derived. Then, the results are generalized to the two, three and N-coupled Klein-Gordon equations

Radiation Damping in a Non-Abelian Strongly-Coupled Gauge Theory

OpenAIRE

Chernicoff, Mariano; Garcia, J. Antonio; Guijosa, Alberto

2010-01-01

We study a `dressed' or `composite' quark in strongly-coupled N=4 super-Yang-Mills (SYM), making use of the AdS/CFT correspondence. We show that the standard string dynamics nicely captures the physics of the quark and its surrounding quantum non-Abelian field configuration, making it possible to derive a relativistic equation of motion that incorporates the effects of radiation damping. From this equation one can deduce a non-standard dispersion relation for the composite quark, as well as a...

A linear algebraic approach to electron-molecule collisions

International Nuclear Information System (INIS)

Collins, L.A.; Schnieder, B.I.

1982-01-01

The linear algebraic approach to electron-molecule collisions is examined by firstly deriving the general set of coupled integrodifferential equations that describe electron collisional processes and then describing the linear algebraic approach for obtaining a solution to the coupled equations. Application of the linear algebraic method to static-exchange, separable exchange and effective optical potential, is examined. (U.K.)

Non-linear dynamo waves in an incompressible medium when the turbulence dissipative coefficients depend on temperature

Directory of Open Access Journals (Sweden)

A. D. Pataraya

1997-01-01

Full Text Available Non-linear α-ω; dynamo waves existing in an incompressible medium with the turbulence dissipative coefficients depending on temperature are studied in this paper. We investigate of α-ω solar non-linear dynamo waves when only the first harmonics of magnetic induction components are included. If we ignore the second harmonics in the non-linear equation, the turbulent magnetic diffusion coefficient increases together with the temperature, the coefficient of turbulent viscosity decreases, and for an interval of time the value of dynamo number is greater than 1. In these conditions a stationary solution of the non-linear equation for the dynamo wave's amplitude exists; meaning that the magnetic field is sufficiently excited. The amplitude of the dynamo waves oscillates and becomes stationary. Using these results we can explain the existence of Maunder's minimum.

The collinearly-improved Balitsky–Kovchegov equation

Energy Technology Data Exchange (ETDEWEB)

Iancu, E.; Madrigal, J.D. [Institut de Physique Théorique, CEA Saclay, CNRS UMR 3681, F-91191 Gif-sur-Yvette (France); Mueller, A.H. [Department of Physics, Columbia University, New York, NY 10027 (United States); Soyez, G. [Institut de Physique Théorique, CEA Saclay, CNRS UMR 3681, F-91191 Gif-sur-Yvette (France); Triantafyllopoulos, D.N. [European Centre for Theoretical Studies in Nuclear Physics and Related Areas ECT* and Fondazione Bruno Kessler, Strada delle Tabarelle 286, I-38123 Villazzano (Italy)

2016-12-15

The high-energy evolution in perturbative QCD suffers from a severe lack-of-convergence problem, due to higher order corrections enhanced by double and single transverse logarithms. We resum double logarithms to all orders within the non-linear Balitsky-Kovchegov equation, by taking into account successive soft gluon emissions strongly ordered in lifetime. We further resum single logarithms generated by the first non-singular part of the splitting functions and by the one-loop running of the coupling. The resummed BK equation admits stable solutions, which are used to successfully fit the HERA data at small x for physically acceptable initial conditions and reasonable values of the fit parameters.

Non linear excitation of waves at the vicinity of plasma resonance

International Nuclear Information System (INIS)

Chiron, Arnaud

1992-01-01

This research thesis reports the study of the non linear evolution of ionic acoustic and plasma waves excited by resonant absorption of an electromagnetic wave, in a non collisional plasma, without external magnetic field, and with a parabolic density profile. The plasma resonance occurs about the density profile peak. The numerical resolution of the Zakharov equation system is performed to describe the coupled evolution of the plasma wave electric field envelope, and low frequency density disturbances. Experiments performed in the microwave domain show the existence of a new effect related to the modification of the electromagnetic wave propagation under the influence of plasma density disturbances created by the ponderomotive force. This effect which results in a collisional relaxation of plasma waves trapped in the cavity formed at resonance, cannot be taken into account by a numerical simulation using a capacitive pump field. Measurements showed that plasma waves were trapped and relaxing in a cavity with characteristic dimensions of some thousands of Debye lengths, and that the plasma wave in the cavity was stationary. A new turbulence regime is thus highlighted [fr

Effect of linear and non-linear components in the temperature dependences of thermoelectric properties on the energy conversion efficiency

International Nuclear Information System (INIS)

Yamashita, Osamu

2009-01-01

The new thermal rate equations were built up by taking the linear and non-linear components in the temperature dependences of the Seebeck coefficient α, the electrical resistivity ρ and thermal conductivity κ of a thermoelectric (TE) material into the thermal rate equations on the assumption that their temperature dependences are expressed by a quadratic function of temperature T. The energy conversion efficiency η for a single TE element was formulated using the new thermal rate ones proposed here. By applying it to the high-performance half-Heusler compound, the non-linear component in the temperature dependence of α among those of the TE properties has the greatest effect on η, so that η/η 0 was increased by 11% under the condition of T = 510 K and ΔT = 440 K, where η 0 is a well-known conventional energy conversion efficiency. It was thus found that the temperature dependences of TE properties have a significant influence on η. When one evaluates the accurate achievement rate of η exp obtained experimentally for a TE generator, therefore, η exp should be compared with η the estimated from the theoretical expression proposed here, not with η 0 , particularly when there is a strong non-linearity in the temperature dependence of TE properties.

Non-linear effective Lagrangian treatment of 'Penguin' interaction

International Nuclear Information System (INIS)

Pham, T.N.

1984-01-01

Using the non-linear effective lagrangian technique, we show explicitly that only derivative coupling is allowed for the K - π, K -> 2 π and K -> 3 π transitions induced by the ΔS = 1 Penguin operator of SVZ in agreement with chiral symmetry requirements. From a derivative coupling (3, anti 3) mass term and the SU(3) breaking effect for fsub(K)/fsub(π), we estimate the strength of the Penguin interactions and find it too small to account for the ΔI = 1/2 amplitude. (orig.)

HESS Opinions: Linking Darcy's equation to the linear reservoir

Science.gov (United States)

Savenije, Hubert H. G.

2018-03-01

In groundwater hydrology, two simple linear equations exist describing the relation between groundwater flow and the gradient driving it: Darcy's equation and the linear reservoir. Both equations are empirical and straightforward, but work at different scales: Darcy's equation at the laboratory scale and the linear reservoir at the watershed scale. Although at first sight they appear similar, it is not trivial to upscale Darcy's equation to the watershed scale without detailed knowledge of the structure or shape of the underlying aquifers. This paper shows that these two equations, combined by the water balance, are indeed identical provided there is equal resistance in space for water entering the subsurface network. This implies that groundwater systems make use of an efficient drainage network, a mostly invisible pattern that has evolved over geological timescales. This drainage network provides equally distributed resistance for water to access the system, connecting the active groundwater body to the stream, much like a leaf is organized to provide all stomata access to moisture at equal resistance. As a result, the timescale of the linear reservoir appears to be inversely proportional to Darcy's conductance, the proportionality being the product of the porosity and the resistance to entering the drainage network. The main question remaining is which physical law lies behind pattern formation in groundwater systems, evolving in a way that resistance to drainage is constant in space. But that is a fundamental question that is equally relevant for understanding the hydraulic properties of leaf veins in plants or of blood veins in animals.

Non-resonant oscillations for some third-order differential equations II

International Nuclear Information System (INIS)

Ezeilo, J.O.C.; Omari, P.

1987-11-01

The existence of 2π-periodic solutions to the equation x'''+ax''+g(t,x')+cx=p(t) is proved, under certain non-resonance conditions on the non-linear function g(t,y). Here a,c are constants, but the case where a,c are not necessarily constants is also discussed, subject to some rather special non-resonance conditions on g. The uniqueness of the solutions is also examined. (author). 12 refs

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.

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

Memory loss process and non-Gibbsian equilibrium solutions of master equations

International Nuclear Information System (INIS)

Cataldo, H.M.; Hernandez, E.S.

1988-01-01

The phonon dynamics of a harmonic oscillator coupled to a steady reservoir is studied. In the Markovian limit, the equilibrium is reached through a progressive loss of memory process which involves the moments of the initial distribution. The relationship to the non-Markovian equations of motion and its resolvent poles is settled. As a particular model of the coupling mechanism is adopted, the possibility of non-Gibbsian equilibrium distribution arises, which is analyzed focusing upon the dependence of various parameters of the system on an effective equilibrium temperature

Error bounds on block Gauss-Seidel solutions of coupled multiphysics problems

KAUST Repository

Whiteley, J. P.

2011-05-09

Mathematical models in many fields often consist of coupled sub-models, each of which describes a different physical process. For many applications, the quantity of interest from these models may be written as a linear functional of the solution to the governing equations. Mature numerical solution techniques for the individual sub-models often exist. Rather than derive a numerical solution technique for the full coupled model, it is therefore natural to investigate whether these techniques may be used by coupling in a block Gauss-Seidel fashion. In this study, we derive two a posteriori bounds for such linear functionals. These bounds may be used on each Gauss-Seidel iteration to estimate the error in the linear functional computed using the single physics solvers, without actually solving the full, coupled problem. We demonstrate the use of the bound first by using a model problem from linear algebra, and then a linear ordinary differential equation example. We then investigate the effectiveness of the bound using a non-linear coupled fluid-temperature problem. One of the bounds derived is very sharp for most linear functionals considered, allowing us to predict very accurately when to terminate our block Gauss-Seidel iteration. © 2011 John Wiley & Sons, Ltd.

Error bounds on block Gauss-Seidel solutions of coupled multiphysics problems

KAUST Repository

Whiteley, J. P.; Gillow, K.; Tavener, S. J.; Walter, A. C.

2011-01-01

Mathematical models in many fields often consist of coupled sub-models, each of which describes a different physical process. For many applications, the quantity of interest from these models may be written as a linear functional of the solution to the governing equations. Mature numerical solution techniques for the individual sub-models often exist. Rather than derive a numerical solution technique for the full coupled model, it is therefore natural to investigate whether these techniques may be used by coupling in a block Gauss-Seidel fashion. In this study, we derive two a posteriori bounds for such linear functionals. These bounds may be used on each Gauss-Seidel iteration to estimate the error in the linear functional computed using the single physics solvers, without actually solving the full, coupled problem. We demonstrate the use of the bound first by using a model problem from linear algebra, and then a linear ordinary differential equation example. We then investigate the effectiveness of the bound using a non-linear coupled fluid-temperature problem. One of the bounds derived is very sharp for most linear functionals considered, allowing us to predict very accurately when to terminate our block Gauss-Seidel iteration. © 2011 John Wiley & Sons, Ltd.

Exactly soluble two-state quantum models with linear couplings

International Nuclear Information System (INIS)

Torosov, B T; Vitanov, N V

2008-01-01

A class of exact analytic solutions of the time-dependent Schroedinger equation is presented for a two-state quantum system coherently driven by a nonresonant external field. The coupling is a linear function of time with a finite duration and the detuning is constant. Four special models are considered in detail, namely the shark, double-shark, tent and zigzag models. The exact solution is derived by rotation of the Landau-Zener propagator at an angle of π/4 and is expressed in terms of Weber's parabolic cylinder function. Approximations for the transition probabilities are derived for all four models by using the asymptotics of the Weber function; these approximations demonstrate various effects of physical interest for each model

Linear System of Equations, Matrix Inversion, and Linear Programming Using MS Excel

Science.gov (United States)

El-Gebeily, M.; Yushau, B.

2008-01-01

In this note, we demonstrate with illustrations two different ways that MS Excel can be used to solve Linear Systems of Equation, Linear Programming Problems, and Matrix Inversion Problems. The advantage of using MS Excel is its availability and transparency (the user is responsible for most of the details of how a problem is solved). Further, we…

Contributions of non-intrusive coupling in nonlinear structural mechanics

International Nuclear Information System (INIS)

Duval, Mickael

2016-01-01

This PhD thesis, part of the ANR ICARE project, aims at developing methods for complex analysis of large scale structures. The scientific challenge is to investigate very localised areas, but potentially critical as of mechanical systems resilience. Classically, representation models, discretizations, mechanical behaviour models and numerical tools are used at both global and local scales for simulation needs of graduated complexity. Global problem is handled by a generic code with topology (plate formulation, geometric approximation...) and behaviour (homogenization) simplifications while local analysis needs implementation of specialized tools (routines, dedicated codes) for an accurate representation of the geometry and behaviour. The main goal of this thesis is to develop an efficient non-intrusive coupling tool for multi-scale and multi-model structural analysis. Constraints of non-intrusiveness result in the non-modification of the stiffness operator, connectivity and the global model solver, allowing to work in a closed source software environment. First, we provide a detailed study of global/local non-intrusive coupling algorithm. Making use of several relevant examples (cracking, elastic-plastic behaviour, contact...), we show the efficiency and the flexibility of such coupling method. A comparative analysis of several optimisation tools is also carried on, and the interacting multiple patches situation is handled. Then, non-intrusive coupling is extended to globally non-linear cases, and a domain decomposition method with non-linear re-localization is proposed. Such methods allowed us to run a parallel computation using only sequential software, on a high performance computing cluster. Finally, we apply the coupling algorithm to mesh refinement with patches of finite elements. We develop an explicit residual based error estimator suitable for multi-scale solutions arising from the non-intrusive coupling, and apply it inside an error driven local mesh

A canonical eight-dimensional formalism for linear and non-linear classical spin-orbit motion in storage rings

International Nuclear Information System (INIS)

Barber, D.P.; Heinemann, K.; Ripken, G.

1991-05-01

In the following report we begin to reformulate work by Derbenev on the behaviour of coupled quantized spin-orbit motion. To this end we present a classical symplectic treatment of linear and non-linear spin-orbit motion for storage rings using a fully coupled eight-dimensional formalism which generalizes earlier investigations of coupled synchro-betatron oscillations by introducing two additional canonical spin variables which behave, in a small-angle limit, like those already used in linearised spin theory. Thus in addition to the usual x-z-s couplings, both the spin to orbit and orbit to spin coupling are described canonically. Since the spin Hamiltonian can be expanded in a Taylor series in canonical variables, the formalism is convenient for use in 8-dimensional symplectic tracking calculations with the help, for example, of Lie algebra or differential algebra for the study of chaotic spin motion, for construction of spin normal forms and for the study of the effect of Stern-Gerlach forces. (orig.)

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.

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.

A General Linear Method for Equating with Small Samples

Science.gov (United States)

Albano, Anthony D.

2015-01-01

Research on equating with small samples has shown that methods with stronger assumptions and fewer statistical estimates can lead to decreased error in the estimated equating function. This article introduces a new approach to linear observed-score equating, one which provides flexible control over how form difficulty is assumed versus estimated…

Non-linear dynamic response of reactor containment

International Nuclear Information System (INIS)

Takemori, T.; Sotomura, K.; Yamada, M.

1975-01-01

A computer program was developed to investigate the elasto-plastic behavior of structures. This program is outlined and the problems of non-linear response of structures are discussed. Since the mode superposition method is only valid in an elastic analysis, the direct integration method was adopted here. As the sample model, an actual reactor containment (reactor building) of PWR plant was adopted. This building consists of three components, that is, a concrete internal structure, a steel containment vessel and a concrete outer shield wall. These components are resting on a rigid foundation mat. Therefore they were modeled with a lumped mass model respectively and coupled on the foundation. The following assumptions were employed to establish the properties of dynamic model: rocking and swaying springs of soil can be obtained from an elastic half-space solution, and the hysteretic characteristic of springs is bi-linear; springs connecting each mass are dealt with shear beams so that both bending and shear deflections can be included (Hysteretic characteristics of springs are linear, bi-linear and tri-linear for the internal structure, the containment vessel and the outer shield wall, respectively); generally, each damping coefficient is given for each mode in modal superposition (However, a damping matrix must be made directly in a non-linear response). Therefore the damping matrix of the model was made by combining the damping matrices [C] of each component obtained by Caughy's method and a damping value of the rocking and swaying by the half-space solution. On the basis of above conditions, the non-linear response of the structure was obtained and the difference between elastic and elasto-plastic analysis is presented

Radiation Damping in a Non-Abelian Strongly-Coupled Gauge Theory

International Nuclear Information System (INIS)

Chernicoff, Mariano; Garcia, J. Antonio; Gueijosa, Alberto

2011-01-01

We study the dynamics of a 'composite' or 'dressed' quark in strongly-coupled large-N c N=4 super-Yang-Mills (SYM), making use of the AdS/CFT correspondence. We show that the standard string dynamics nicely captures the physics of the quark and its surrounding non-Abelian field configuration, making it possible to derive a relativistic equation of motion that incorporates the effects of radiation damping. From this equation one can deduce a non-standard dispersion relation for the composite quark, as well as a Lorentz covariant formula for its rate of radiation.

Radiation Damping in a Non-Abelian Strongly-Coupled Gauge Theory

Science.gov (United States)

Chernicoff, Mariano; García, J. Antonio; Güijosa, Alberto

2011-09-01

We study the dynamics of a 'composite` or 'dressed` quark in strongly-coupled large-Nc N=4 super-Yang-Mills (SYM), making use of the AdS/CFT correspondence. We show that the standard string dynamics nicely captures the physics of the quark and its surrounding non-Abelian field configuration, making it possible to derive a relativistic equation of motion that incorporates the effects of radiation damping. From this equation one can deduce a non-standard dispersion relation for the composite quark, as well as a Lorentz covariant formula for its rate of radiation.

A local-global problem for linear differential equations

NARCIS (Netherlands)

Put, Marius van der; Reversat, Marc

An inhomogeneous linear differential equation Ly = f over a global differential field can have a formal solution for each place without having a global solution. The vector space lgl(L) measures this phenomenon. This space is interpreted in terms of cohomology of linear algebraic groups and is

A local-global problem for linear differential equations

NARCIS (Netherlands)

Put, Marius van der; Reversat, Marc

2008-01-01

An inhomogeneous linear differential equation Ly = f over a global differential field can have a formal solution for each place without having a global solution. The vector space lgl(L) measures this phenomenon. This space is interpreted in terms of cohomology of linear algebraic groups and is

Iterative solution of linear equations in ODE codes. [Krylov subspaces

Energy Technology Data Exchange (ETDEWEB)

Gear, C. W.; Saad, Y.

1981-01-01

Each integration step of a stiff equation involves the solution of a nonlinear equation, usually by a quasi-Newton method that leads to a set of linear problems. Iterative methods for these linear equations are studied. Of particular interest are methods that do not require an explicit Jacobian, but can work directly with differences of function values using J congruent to f(x + delta) - f(x). Some numerical experiments using a modification of LSODE are reported. 1 figure, 2 tables.

Exact Solution of Space-Time Fractional Coupled EW and Coupled MEW Equations Using Modified Kudryashov Method

International Nuclear Information System (INIS)

Raslan, K. R.; Ali, Khalid K.; EL-Danaf, Talaat S.

2017-01-01

In the present paper, we established a traveling wave solution by using modified Kudryashov method for the space-time fractional nonlinear partial differential equations. The method is used to obtain the exact solutions for different types of the space-time fractional nonlinear partial differential equations such as, the space-time fractional coupled equal width wave equation (CEWE) and the space-time fractional coupled modified equal width wave equation (CMEW), which are the important soliton equations. Both equations are reduced to ordinary differential equations by the use of fractional complex transform and properties of modified Riemann–Liouville derivative. We plot the exact solutions for these equations at different time levels. (paper)

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.

The strong running coupling from an approximate gluon Dyson-Schwinger equation

International Nuclear Information System (INIS)

Alkofer, R.; Hauck, A.

1996-01-01

Using Mandelstam's approximation to the gluon Dyson-Schwinger equation we calculate the gluon self-energy in a renormalisation group invariant fashion. We obtain a non-perturbative Β function. The scaling behavior near the ultraviolet stable fixed point is in good agreement with perturbative QCD. No further fixed point for positive values of the coupling is found: α S increases without bound in the infrared

A strongly coupled open system with a non-linear bath: fluctuation-dissipation and Langevin dynamics

Science.gov (United States)

Bhadra, Chitrak

2018-03-01

The study of Langevin dynamics and fluctuation-dissipation relation (FDR) for a generic probe system (represented by a mass M ), bilinearly coupled to a bath of harmonic oscillators, has been a standard paradigm for the microscopic theory of stochastic processes for several decades. The question that we probe in this paper is, how robust the structure of the classical FDR is, when one replaces the harmonic bath by an anharmonic one in the limit of strong system-bath coupling? Such a picture carries the signature of the probe system in the zeroth order through a nonlocal time kernel. We observe that the two-time noise correlations hold a rich structure from which the usual FDR emerges only in the leading order of perturbation. Beyond this order, multiple time scales and nontrivial dependence on the temperature starts to manifest. These new aspects conspire to break the time-translational invariance of the noise-correlations. Several other interesting features show up and we discuss them methodically through rigorous calculations order-by-order in perturbation. This formalistic derivation along with a specific example of non-linearity can be easily applied to a huge range of processes and statistical observables that fall under the purview of a system-reservoir theory.

A critical oscillation constant as a variable of time scales for half-linear dynamic equations

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel

2010-01-01

Roč. 60, č. 2 (2010), s. 237-256 ISSN 0139-9918 R&D Projects: GA AV ČR KJB100190701 Institutional research plan: CEZ:AV0Z10190503 Keywords : dynamic equation * time scale * half-linear equation * (non)oscillation criteria * Hille-Nehari criteria * Kneser criteria * critical constant * oscillation constant * Hardy inequality Subject RIV: BA - General Mathematics Impact factor: 0.316, year: 2010 http://link.springer.com/article/10.2478%2Fs12175-010-0009-7

Darboux transformations and linear parabolic partial differential equations

International Nuclear Information System (INIS)

Arrigo, Daniel J.; Hickling, Fred

2002-01-01

Solutions for a class of linear parabolic partial differential equation are provided. These solutions are obtained by first solving a system of (n+1) nonlinear partial differential equations. This system arises as the coefficients of a Darboux transformation and is equivalent to a matrix Burgers' equation. This matrix equation is solved using a generalized Hopf-Cole transformation. The solutions for the original equation are given in terms of solutions of the heat equation. These results are applied to the (1+1)-dimensional Schroedinger equation where all bound state solutions are obtained for a 2n-parameter family of potentials. As a special case, the solutions for integral members of the regular and modified Poeschl-Teller potentials are recovered. (author). Letter-to-the-editor

Piecewise-linear and bilinear approaches to nonlinear differential equations approximation problem of computational structural mechanics

OpenAIRE

Leibov Roman

2017-01-01

This paper presents a bilinear approach to nonlinear differential equations system approximation problem. Sometimes the nonlinear differential equations right-hand sides linearization is extremely difficult or even impossible. Then piecewise-linear approximation of nonlinear differential equations can be used. The bilinear differential equations allow to improve piecewise-linear differential equations behavior and reduce errors on the border of different linear differential equations systems ...

Lax pair and exact solutions of a discrete coupled system related to coupled KdV and coupled mKdV equations

International Nuclear Information System (INIS)

Liu Ping; Jia Man; Lou Senyue

2007-01-01

A modified Korteweg-de Vries (mKdV) lattice is also found to be a discrete Korteweg-de Vries (KdV) equation in this paper. The Lax pair for the discrete equation is found with the help of the Lax pair for a similar discrete equation. A Lax-integrable coupled extension of the lattice is posed, which is a common discrete version of both the coupled KdV and coupled mKdV systems. Some rational expansions of the Jacobian elliptic, trigonometric and hyperbolic functions are used to construct cnoidal waves, negaton and positon solutions of the discrete coupled system

Technological pedagogical content knowledge of junior high school mathematics teachers in teaching linear equation

Science.gov (United States)

Wati, S.; Fitriana, L.; Mardiyana

2018-04-01

Linear equation is one of the topics in mathematics that are considered difficult. Student difficulties of understanding linear equation can be caused by lack of understanding this concept and the way of teachers teach. TPACK is a way to understand the complex relationships between teaching and content taught through the use of specific teaching approaches and supported by the right technology tools. This study aims to identify TPACK of junior high school mathematics teachers in teaching linear equation. The method used in the study was descriptive. In the first phase, a survey using a questionnaire was carried out on 45 junior high school mathematics teachers in teaching linear equation. While in the second phase, the interview involved three teachers. The analysis of data used were quantitative and qualitative technique. The result PCK revealed teachers emphasized developing procedural and conceptual knowledge through reliance on traditional in teaching linear equation. The result of TPK revealed teachers’ lower capacity to deal with the general information and communications technologies goals across the curriculum in teaching linear equation. The result indicated that PowerPoint constitutes TCK modal technological capability in teaching linear equation. The result of TPACK seems to suggest a low standard in teachers’ technological skills across a variety of mathematics education goals in teaching linear equation. This means that the ability of teachers’ TPACK in teaching linear equation still needs to be improved.

Some problems on non-linear semigroups and the blow-up of integral solutions

International Nuclear Information System (INIS)

Pavel, N.H.

1983-07-01

After some introductory remarks, this highly mathematical document considers a unifying approach in the theory of non-linear semigroups. Then a brief survey is given on blow-up of mild solutions from the semilinear case. Finally, the global behavior of solutions to non-linear evolution equations is addressed; it is found that classical results on the behavior of the maximal solution u as t up-arrow tsub(max) hold also for integral solutions

On some perturbation techniques for quasi-linear parabolic equations

Directory of Open Access Journals (Sweden)

Igor Malyshev

1990-01-01

Full Text Available We study a nonhomogeneous quasi-linear parabolic equation and introduce a method that allows us to find the solution of a nonlinear boundary value problem in “explicit” form. This task is accomplished by perturbing the original equation with a source function, which is then found as a solution of some nonlinear operator equation.

Asymptotic method for non-linear magnetosonic waves in an isothermal plasma with a finite conductivity

Energy Technology Data Exchange (ETDEWEB)

Fusco, D [Messina Univ. (Italy). Instituto de Matematica

1979-01-01

The paper is concerned with a three-dimensional theory of non-linear magnetosonic waves in a turbulent plasma. A perturbation method is used that allows a transport equation, like Burgers equation but with a variable coefficient to be obtained.

Non-Linear Structural Dynamics Characterization using a Scanning Laser Vibrometer

Science.gov (United States)

Pai, P. F.; Lee, S.-Y.

2003-01-01

This paper presents the use of a scanning laser vibrometer and a signal decomposition method to characterize non-linear dynamics of highly flexible structures. A Polytec PI PSV-200 scanning laser vibrometer is used to measure transverse velocities of points on a structure subjected to a harmonic excitation. Velocity profiles at different times are constructed using the measured velocities, and then each velocity profile is decomposed using the first four linear mode shapes and a least-squares curve-fitting method. From the variations of the obtained modal \\ielocities with time we search for possible non-linear phenomena. A cantilevered titanium alloy beam subjected to harmonic base-excitations around the second. third, and fourth natural frequencies are examined in detail. Influences of the fixture mass. gravity. mass centers of mode shapes. and non-linearities are evaluated. Geometrically exact equations governing the planar, harmonic large-amplitude vibrations of beams are solved for operational deflection shapes using the multiple shooting method. Experimental results show the existence of 1:3 and 1:2:3 external and internal resonances. energy transfer from high-frequency modes to the first mode. and amplitude- and phase- modulation among several modes. Moreover, the existence of non-linear normal modes is found to be questionable.

A frequency domain linearized Navier-Stokes equations approach to acoustic propagation in flow ducts with sharp edges.

Science.gov (United States)

Kierkegaard, Axel; Boij, Susann; Efraimsson, Gunilla

2010-02-01

Acoustic wave propagation in flow ducts is commonly modeled with time-domain non-linear Navier-Stokes equation methodologies. To reduce computational effort, investigations of a linearized approach in frequency domain are carried out. Calculations of sound wave propagation in a straight duct are presented with an orifice plate and a mean flow present. Results of transmission and reflections at the orifice are presented on a two-port scattering matrix form and are compared to measurements with good agreement. The wave propagation is modeled with a frequency domain linearized Navier-Stokes equation methodology. This methodology is found to be efficient for cases where the acoustic field does not alter the mean flow field, i.e., when whistling does not occur.

Mathematical Methods in Wave Propagation: Part 2--Non-Linear Wave Front Analysis

Science.gov (United States)

Jeffrey, Alan

1971-01-01

The paper presents applications and methods of analysis for non-linear hyperbolic partial differential equations. The paper is concluded by an account of wave front analysis as applied to the piston problem of gas dynamics. (JG)

Linearized gyro-kinetic equation

International Nuclear Information System (INIS)

Catto, P.J.; Tsang, K.T.

1976-01-01

An ordering of the linearized Fokker-Planck equation is performed in which gyroradius corrections are retained to lowest order and the radial dependence appropriate for sheared magnetic fields is treated without resorting to a WKB technique. This description is shown to be necessary to obtain the proper radial dependence when the product of the poloidal wavenumber and the gyroradius is large (k rho much greater than 1). A like particle collision operator valid for arbitrary k rho also has been derived. In addition, neoclassical, drift, finite β (plasma pressure/magnetic pressure), and unperturbed toroidal electric field modifications are treated

Nonlinear dynamics and control strategies: On a energy harvester vibrating system with a linear form to non-ideal motor torquet

Directory of Open Access Journals (Sweden)

de Pontes B. R.

2012-07-01

Full Text Available In this paper, we deal with the research of a vibrating model of an energy harvester device, including the nonlinearities in the model of the piezoelectric coupling and the non-ideal excitation. We show, using numerical simulations, in the analysis of the dynamic responses, that the harvested power is influenced by non-linear vibrations of the structure. Chaotic behavior was also observed, causing of the loss of energy throughout the simulation time. Using a perturbation technique, we find an approximate analytical solution for the non-ideal system. Then, we apply both two control techniques, to keep the considered system, into a stable condition. Both the State Dependent Ricatti Equation (SDRE control as the feedback control by changing the energy of the oscillator, were efficient in controlling of the considered non-ideal system.

Experimental quantum computing to solve systems of linear equations.

Science.gov (United States)

Cai, X-D; Weedbrook, C; Su, Z-E; Chen, M-C; Gu, Mile; Zhu, M-J; Li, Li; Liu, Nai-Le; Lu, Chao-Yang; Pan, Jian-Wei

2013-06-07

Solving linear systems of equations is ubiquitous in all areas of science and engineering. With rapidly growing data sets, such a task can be intractable for classical computers, as the best known classical algorithms require a time proportional to the number of variables N. A recently proposed quantum algorithm shows that quantum computers could solve linear systems in a time scale of order log(N), giving an exponential speedup over classical computers. Here we realize the simplest instance of this algorithm, solving 2×2 linear equations for various input vectors on a quantum computer. We use four quantum bits and four controlled logic gates to implement every subroutine required, demonstrating the working principle of this algorithm.

Solution methods for large systems of linear equations in BACCHUS

International Nuclear Information System (INIS)

Homann, C.; Dorr, B.

1993-05-01

The computer programme BACCHUS is used to describe steady state and transient thermal-hydraulic behaviour of a coolant in a fuel element with intact geometry in a fast breeder reactor. In such computer programmes generally large systems of linear equations with sparse matrices of coefficients, resulting from discretization of coolant conservation equations, must be solved thousands of times giving rise to large demands of main storage and CPU time. Direct and iterative solution methods of the systems of linear equations, available in BACCHUS, are described, giving theoretical details and experience with their use in the programme. Besides use of a method of lines, a Runge-Kutta-method, for solution of the partial differential equation is outlined. (orig.) [de

Perturbations of linear delay differential equations at the verge of instability.

Science.gov (United States)

Lingala, N; Namachchivaya, N Sri

2016-06-01

The characteristic equation for a linear delay differential equation (DDE) has countably infinite roots on the complex plane. This paper considers linear DDEs that are on the verge of instability, i.e., a pair of roots of the characteristic equation lies on the imaginary axis of the complex plane and all other roots have negative real parts. It is shown that when small noise perturbations are present, the probability distribution of the dynamics can be approximated by the probability distribution of a certain one-dimensional stochastic differential equation (SDE) without delay. This is advantageous because equations without delay are easier to simulate and one-dimensional SDEs are analytically tractable. When the perturbations are also linear, it is shown that the stability depends on a specific complex number. The theory is applied to study oscillators with delayed feedback. Some errors in other articles that use multiscale approach are pointed out.