MHB How Can Sources and Sinks be Incorporated in a Parabolic PDE Algorithm?

  • Thread starter Thread starter ognik
  • Start date Start date
  • Tags Tags
    Algorithm Pde
AI Thread Summary
The discussion revolves around incorporating sources and sinks into a parabolic PDE algorithm, specifically using a Fortran program based on a given equation. The user expresses confusion about the program's output, which shows similar results across various parameters, indicating a potential validation issue. They seek clarification on the concept of sources and sinks, suggesting that these represent inflow and outflow in the context of the diffusion-like equation. The user is particularly interested in how to implement these concepts within the existing program structure. Overall, they are looking for guidance on enhancing the program to study the effects of sources and sinks on the solutions.
ognik
Messages
626
Reaction score
2
Hi - on the last chapter of this course and was feeling much better about it all, but I now confess to being back in my normal state - confused. I am given a simple fortran program (code attached in the zip file) and asked to investigate its accuracy and stability, for various values of $$\Delta$$t and lattice spacings. The program is an implementation of:
$ {\phi}^{n+1} = \frac{1}{1 + H\Delta t}\left[{\phi}^{n} + S^n \Delta t \right] $ (H is hermitian)

I have run this program for various sets of values - and the output all looks so similar that I can't see anything to discuss. The errors are similar magnitude. Some combinations of input don't produce any output - but that should be just a validation issue, as I say it is a simple program with no frills.
Someone give me a clue or 2 please...
 

Attachments

Mathematics news on Phys.org
I see a few people have looked without responding - I think I should be finding limitations of this approach, but so far haven't, so I might just be missing something obvious to someone else...either way, hopefully talking about the 2nd part of the exercise will help.

The second (and maybe prime) part of the exercise says 'Incorporate sources or sinks along the lattice and study the solutions that arise when ∅ vanishes everywhere at t = 0

The text hasn't used the concept of sources/sinks before, but I think the PDE in question is similar to a diffusion equation so - please correct me - sources and sinks would be where there is inflow/outflow from the volume under study? And they are related to the Sn term in the equation? Earlier in this chapter I did some exercises on discretization, so I am familiar with that and lattices, but I am clueless otherwise (this course is about computational physics and as it happens I won't do equations like this until next year, C'est la vie)
 
Hi, I have edged a bit further along with this:

Please correct me carefully here - sources and sinks would be where there is inflow/outflow from the volume/area under study? Therefore should I be looking at something like $ \frac{\partial \phi}{\partial t} < 0 \: $ for a sink? $\: > 0 $ for a source?

If so, how does one incorporate them, along the lattice, into the attached program? I really am just blank about this...an example would be very useful! Thanks.
 
Seemingly by some mathematical coincidence, a hexagon of sides 2,2,7,7, 11, and 11 can be inscribed in a circle of radius 7. The other day I saw a math problem on line, which they said came from a Polish Olympiad, where you compute the length x of the 3rd side which is the same as the radius, so that the sides of length 2,x, and 11 are inscribed on the arc of a semi-circle. The law of cosines applied twice gives the answer for x of exactly 7, but the arithmetic is so complex that the...
Is it possible to arrange six pencils such that each one touches the other five? If so, how? This is an adaption of a Martin Gardner puzzle only I changed it from cigarettes to pencils and left out the clues because PF folks don’t need clues. From the book “My Best Mathematical and Logic Puzzles”. Dover, 1994.
Thread 'Imaginary Pythagoras'
I posted this in the Lame Math thread, but it's got me thinking. Is there any validity to this? Or is it really just a mathematical trick? Naively, I see that i2 + plus 12 does equal zero2. But does this have a meaning? I know one can treat the imaginary number line as just another axis like the reals, but does that mean this does represent a triangle in the complex plane with a hypotenuse of length zero? Ibix offered a rendering of the diagram using what I assume is matrix* notation...
Back
Top