How can I solve a system of BVPs with only boundary conditions given?

Click For Summary

Discussion Overview

The discussion revolves around solving a system of boundary value problems (BVPs) consisting of four differential equations with specified boundary conditions for the variables x and y. The focus is on finding numerical or analytical methods to approach the solution given the complexity of the equations.

Discussion Character

  • Exploratory
  • Technical explanation
  • Mathematical reasoning

Main Points Raised

  • One participant presents a system of four equations and requests recommendations for solving them given boundary conditions x(0) = y(0) = -1 and x(π/2) = y(π/2) = 1.
  • Another participant suggests converting the equations to polar coordinates, proposing transformations for x and y, and reformulating the equations accordingly.
  • The same participant elaborates on the transformations and provides the new forms of the equations in polar coordinates, indicating how to manipulate them further.
  • A later reply acknowledges the assistance and corrects a term in one of the equations, indicating a potential oversight in the formulation.
  • The participant expresses uncertainty about how to proceed towards finding a solution after the transformations have been made.

Areas of Agreement / Disagreement

Participants do not reach a consensus on a specific method to solve the system, and multiple approaches are discussed without resolution. The discussion remains exploratory with ongoing questions about the next steps.

Contextual Notes

The discussion highlights the complexity of the equations and the challenges in finding solutions, particularly given the boundary conditions and the transformations proposed. There are unresolved aspects regarding the mathematical steps and the implications of the transformations.

amr07
Messages
5
Reaction score
0
Dear all,

I have system(4 equations) of BVPs. Could anybody recommend me, how to solve this system(whatever numericaly or analytical):

x'=-y/sqrt(x^2+y^2) + u
y'=x/sqrt(x^2+y^2) + v
u' = -xy/(x^2+y^2)^3/2 u - [1/sqrt(x^2+y^2) - x^2 /(x^2+y^2)^3/2] v
v' = xy/(x^2+y^2)^3/2 v - [-1/sqrt(x^2+y^2) + y^2 /(x^2+y^2)^3/2] u

and we have only boundary condations for x and y...x(0)=y(0)=-1, x(pi/2)=y(pi/2)=1

thanks
 
Last edited:
Physics news on Phys.org
With all those "[itex]x^2+ y^2[/itex]" in there, the first thing I would try is to convert x and y to polar coordinates: [itex]x= r cos(\theta)[/itex] so that [itex]x^2+ y^2= r^2[/itex], [itex]x'= r' cos(\theta)- r sin(\theta)\theta'[/itex], and [itex]y= r sin(\theta)[/itex] so that [itex]y'= r' sin(\theta)+ r cos(\theta)\theta'[/itex].
Your first equation becomes [itex]r' cos(\theta)- r sin(\theta)\theta'= - sin(\theta)+ u[/itex] and your second equation [itex]r' sin(\theta)+ r cos(\theta)\theta'= cos(\theta)+ v[/itex]. If you multiply the first equation by [itex]cos(\theta)[/itex] and the second by [itex]sin(\theta)[/itex] and add you get [itex]r'= u cos(\theta)+ v\sin(\theta)[/itex]. If you multiply the first equation by [itex]sin(\theta)[/itex], the second by [itex]cos(\theta)[/itex] and subtract first from second you get [itex]r \theta'= 1+ v cos(\theta)- u sin(\theta)[/itex].

Of course, the third and fourth equations become
[tex]u'= -\fra{1}{r}sin(\theta)cos(\theta)u- \left[1- cos^2(\theta)\right]\frac{v}{r}[/tex]
and
[tex]v'= \frac{1}{r}sin(\theta)cos(\theta)u+ \left[1+ sin^2(\theta)\right]\frac{u}{r}[/tex]
 
thanks
 
HallsofIvy said:
With all those "[itex]x^2+ y^2[/itex]" in there, the first thing I would try is to convert x and y to polar coordinates: [itex]x= r cos(\theta)[/itex] so that [itex]x^2+ y^2= r^2[/itex], [itex]x'= r' cos(\theta)- r sin(\theta)\theta'[/itex], and [itex]y= r sin(\theta)[/itex] so that [itex]y'= r' sin(\theta)+ r cos(\theta)\theta'[/itex].
Your first equation becomes [itex]r' cos(\theta)- r sin(\theta)\theta'= - sin(\theta)+ u[/itex] and your second equation [itex]r' sin(\theta)+ r cos(\theta)\theta'= cos(\theta)+ v[/itex]. If you multiply the first equation by [itex]cos(\theta)[/itex] and the second by [itex]sin(\theta)[/itex] and add you get [itex]r'= u cos(\theta)+ v\sin(\theta)[/itex]. If you multiply the first equation by [itex]sin(\theta)[/itex], the second by [itex]cos(\theta)[/itex] and subtract first from second you get [itex]r \theta'= 1+ v cos(\theta)- u sin(\theta)[/itex].

Of course, the third and fourth equations become
[tex]u'= -\fra{1}{r}sin(\theta)cos(\theta)u- \left[1- cos^2(\theta)\right]\frac{v}{r}[/tex]
and
[tex]v'= \frac{1}{r}sin(\theta)cos(\theta)u+ \left[1+ sin^2(\theta)\right]\frac{u}{r}[/tex]


Dear HallsofIvy,

many thanks for your helping. the second term in the fourth equation(for v') should be [1-sin^2]. so we have no also four differential equations. it seems easily but how i can get the solution?

thanks again
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 17 ·
Replies
17
Views
3K
  • · Replies 22 ·
Replies
22
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 8 ·
Replies
8
Views
4K
Replies
3
Views
1K
Replies
8
Views
2K