Solving 2-D partial integro-differential equation

Click For Summary

Discussion Overview

The discussion revolves around solving a 2-D partial integro-differential equation involving the Hilbert transform and Fourier transform, with specific initial and boundary conditions. Participants seek assistance in coding a solution using software like Matlab, Mathematica, or Maple.

Discussion Character

  • Homework-related
  • Technical explanation
  • Exploratory

Main Points Raised

  • One participant presents the integro-differential equation and initial conditions, expressing difficulty in writing code to solve it.
  • Another participant offers to look for comparable examples that may assist in solving the equation.
  • A third participant shares links to external resources that may provide insights or similar problems, although one of the links is noted to contain a potentially incorrect solution.
  • There is a recognition of the specialized knowledge required to address such questions, reflecting on past experiences in academia.

Areas of Agreement / Disagreement

Participants do not reach a consensus on the solution to the equation, and there is acknowledgment of uncertainty regarding the correctness of shared resources.

Contextual Notes

Participants express limitations in their ability to find relevant examples and solutions, indicating a dependency on specialized knowledge and existing resources.

semivermous
Messages
2
Reaction score
0
While reproducing a research paper, I came across the following equation,
∂f/∂t−(H(f)(∂f/∂x)=0
where [H(f)] is hilbert transform of 'f.'
and f=f(x,t) and initial condition is f(x,0)=cos(x) and also has periodic boundary conditions given by
F{H{f(x′,t)}}=i⋅sgn(k)F{f(x,t)},
where F(f(x,t) is Fourier transform of f(x,t).
and here ''t'' runs from 0 to 1.3 seconds

so I think we have to use iterations on basis of 't' while solving this equation.
Please help me in solving this integro differential equation(PDE). I am unable write a code for this(Matlab/Mathematica/Maple)
And suggestions are highly appreciated.
Link for research paper: https://journals.aps.org/prl/pdf/10.1103/PhysRevLett.75.4614
After solving the given equation, we have to get the figure-1 of the paper

Code I have written is:
[CODE title="Mathematica code"]L = Pi; tmax = 1.2; sys = {D[u[x, t], t] -
1/(Pi)*int[u[x, t], x, t]*D[u[x, t], x] == 0, u[0, t] == 1,
u[x, 0] == Cos[x]};
int[u_, x_?NumericQ, t_ /; t == 0] :=
NIntegrate[Cos[xp]/(x - xp), {xp, 0, x, x + 2 L},
Method -> {"InterpolationPointsSubdivision",
Method -> "PrincipalValue"}, MaxRecursion -> 20];
PrintTemporary@Dynamic@{foo, Clock[Infinity]};
Internal`InheritedBlock[{MapThread}, {state} =
NDSolve`ProcessEquations[sys, u, {x, 0, 2 L}, {t, 0, tmax},
StepMonitor :> (foo = t)];
NDSolve`Iterate[state, {0, tmax}];
sol = NDSolve`ProcessSolutions[state]] // AbsoluteTiming
{Plot3D[u[x, t] /. sol, {x, 0, 2 Pi}, {t, 0., 1.}, Mesh -> None,
ColorFunction -> Hue, AxesLabel -> Automatic] // Quiet,
Plot[Evaluate[Table[u[x, t] /. sol, {t, 0., 1., .2}]], {x, -Pi,
Pi}] // Quiet}
[/CODE]
 
Physics news on Phys.org
Hopefully someone here knows mathematica sufficiently to help.

In the meantime, I will look for some comparable examples using the same functions.
 
  • Like
Likes   Reactions: semivermous
The question in first link was posted by myself but the solution is wrong.
 
Ahh okay. Well it may still hope others who read this thread.

One problem with this kind of question is the specialized knowledge needed to answer it. This would happen to me often in grad school before the internet. There was simply no place to turn to get help.
 

Similar threads

Replies
8
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
6
Views
3K
  • · Replies 3 ·
Replies
3
Views
6K
  • · Replies 7 ·
Replies
7
Views
3K