Engineering Why is A(1-y²) Negative When y² < 1 in Van der Pol Equation?

AI Thread Summary
In the Van der Pol equation, the term A(1-y²) is negative when y² < 1, indicating that small amplitude oscillations experience negative damping, which leads to an increase in amplitude. This behavior contrasts with large amplitude oscillations (y² > 1), where A(1-y²) becomes positive, resulting in positive damping and a decay in amplitude. The analysis shows that for small values of y, the system feeds energy into oscillations, while for larger values, it dissipates energy. This dynamic illustrates the unique nature of the Van der Pol oscillator, where oscillations are amplified at low amplitudes and damped at high amplitudes. Understanding this relationship is crucial for analyzing the stability and behavior of non-linear oscillatory systems.
shivajikobardan
Messages
637
Reaction score
54
Homework Statement
positive and negative damping in van der pol equation-simulation and modeling
Relevant Equations
Van der pol equation
The van der pol non-linear equation is given as-:

y''=A(1-y²)y'-By

y=amplitude

The analysis given by the book is this-:

When y²<1
i.e when y is small
A(1-y²) is negative.
A(1-y²) is called damping term.

I don't understand how is it negative? It obviously becomes positive in this case.

Sth similar is done for when y²>1.
The final conclusion is-:
Thus, small amplitude oscillations will build up and large amplitude oscillations will be damped out.
i.e when y is small=>there should be no damping as shown in what I'm confused with. (If I did according to me, then I'd be getting positive value i.e more damping).
 
Physics news on Phys.org
shivajikobardan said:
The analysis given by the book
which book?
shivajikobardan said:
Homework Statement:: positive and negative damping in van der pol equation-simulation and modeling
Relevant Equations:: Van der pol equation

When y²<1
i.e when y is small
A(1-y²) is negative.
depends on sign of A of course
 
drmalawi said:
which book?
i just googled van der pol equation in google books. almost all of them say the same(any resource says the same).
drmalawi said:
depends on sign of A of course
 
Consider ## \ddot y -\mu (1-y^2 ) \dot y + y = 0## where ##y = y(t)## and ##\mu \geq 0##
Compare with the damped oscillator equation
## \ddot y +b \dot y + y = 0## where ##b## is the amount of "damping".
We have that ##b## plays the role of ## -\mu (1-y^2 ) ##
Now, if ##b>0## then the solution ##y## will decay, and if ##b<0## the solution will "blow up"
We can make this analysis now, if ##|y| > 1## then ##|y|^2 > 1## and your "damping" "b" will be positive and hence solution ##y## will decay.
If ##|y|< 1## then ##|y|^2 < 1## and your "damping" "b" will be negative and hence solution ##y## will "blow up".

This is the idea of the Van Der Pol oscillator, oscilations with small amplitude will get their amplitudes increased, and vice versa. Energy is leaked out from the system when it grows to big, and energy is fed in when it grows too small.
 
  • Like
Likes shivajikobardan
Thread 'PLL - How to find all the gains of a PI corrector and fix Ki ? MATLAB'
Hello everyone, it's not really for my homework (because I'm not at school anymore) but I'm new and I don't know how to start a new forum. I've never done automatic before and I need help. I have a 2nd ordre PLL on MATLAB, with an NCO and a PI corrector; the PLL works via the Mid/End algorithm (that's part is ok). I need to find the link between the theoretical formulas and their values in the code to find the value of the gain Ki, and possibly that of Kp if it's not correct. Unfortunately...
Thread 'Have I solved this structural engineering equation correctly?'
Hi all, I have a structural engineering book from 1979. I am trying to follow it as best as I can. I have come to a formula that calculates the rotations in radians at the rigid joint that requires an iterative procedure. This equation comes in the form of: $$ x_i = \frac {Q_ih_i + Q_{i+1}h_{i+1}}{4K} + \frac {C}{K}x_{i-1} + \frac {C}{K}x_{i+1} $$ Where: ## Q ## is the horizontal storey shear ## h ## is the storey height ## K = (6G_i + C_i + C_{i+1}) ## ## G = \frac {I_g}{h} ## ## C...
Back
Top