Message Boards Message Boards

0
|
1501 Views
|
8 Replies
|
0 Total Likes
View groups...
Share
Share this post:

How to get a plot for an equation

Cesar Cardenas Ruiz
Posted 6 months ago

Hello, I'm trying to code the equation shown below and get a curve or possible solution as shown below as well (not all of them). This is my attempt but not sure if it is correct and how to obtain a similar plot?. Are the units right? Any help will be appreciated. Thanks

ClearAll["Global`*"];
T = Quantity[1.2 * 10^6, "Kelvins"];
kb = Quantity[1, "BoltzmannConstant"];
Ms = Quantity[1.99 *10^30, "Kilograms"];
G = Quantity[6.67 *10^-11, "Meters^2*Newton/Kilograms^2"]; 
mH = Quantity[1.99 * 10^30, "Kilograms"];
r = Quantity[1*10^9, "Meters"]
u = Quantity[17, "Meters"/"Second"]

    uc = Sqrt[2*kb*T/mH]
    rc = G*Ms*mH/4*kb*T
    eq = (u/uc)^2 - 2 * Log[u/uc] == 4 *Log[r/rc] + 4 *(rc/r) - 3

enter image description here

enter image description here
POSTED BY: Cesar Cardenas Ruiz
8 Replies
Sort By:

Agreed, there must be something wrong with your code for uc and rc. Can you post the original equations for those?
POSTED BY: Gareth Russell

These are the equations for rc and uc. Any help will be appreciated. Thank you

enter image description here
POSTED BY: Cesar Cardenas Ruiz

So as Gianluca points out, u and uc have the same units (m/s), but r and rc do not. r is in meters, but r_c based on your code is in newtons * seconds^2 / kilograms. Could there be something wrong with your definitions of Ms and mH? You have the both as (identical) masses. But the different capitalization in your source suggests they might not be the same thing?
POSTED BY: Gareth Russell

Ms ans mH are different things. If rc, r, u and uc are assumed to be dimensionless, is there any way to get a curve like the ones shown in the previous plot? Thanks
POSTED BY: Cesar Cardenas Ruiz

Well that's odd because they are the same in your code. Perhaps you can post the reference for these equations? Something to do with gas particle velocity I assume...
POSTED BY: Gareth Russell

Thanks for your reply. This is the reference.

DYNAMICS OF THE INTERPLANETARY GAS AND MAGNETIC FIELDS
POSTED BY: Cesar Cardenas Ruiz

I see that Figure 1 is what you are trying to recreate, but the relevant equations (16 and 17) have a different form from the ones you posted earlier and, alas, I don't have the time to try to see how your version relates to those in the article. Especially as I'm not a physicist!
POSTED BY: Gareth Russell

If you take the logarithm, as in Log[r/rc], the argument must be dimensionless. In your case r is a length but rc is not.
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract