Fdtxjr

What initially awakened the Balrog?

Taylor expansion of ln(1-x)

Morning, Afternoon, Night Kanji

How often does castling occur in grandmaster games?

Can anything be seen from the center of the Boötes void? How dark would it be?

When a candle burns, why does the top of wick glow if bottom of flame is hottest?

Is CEO the "profession" with the most psychopaths?

What is a fractional matching?

Do I really need to have a message in a novel to appeal to readers?

What is the difference between globalisation and imperialism?

Selecting user stories during sprint planning

Can a new player join a group only when a new campaign starts?

What is the appropriate index architecture when forced to implement IsDeleted (soft deletes)?

Illegal assignment from sObject to Id

Time to Settle Down!

Trademark violation for app?

How were pictures turned from film to a big picture in a picture frame before digital scanning?

What is "gratricide"?

What is this clumpy 20-30cm high yellow-flowered plant?

Is a ledger board required if the side of my house is wood?

Is it possible for SQL statements to execute concurrently within a single session in SQL Server?

What was the first language to use conditional keywords?

Denied boarding although I have proper visa and documentation. To whom should I make a complaint?

Disembodied hand growing fangs

Shock in the Inviscid Burgers Equation

0

$begingroup$

I'm trying to implement the Finite-Volume method using local Lax-Friedrichs flux function. I found an example to the Riemann-problem, here is the solution plot:

So from characteristics $$u=begincases
1, x<frac12t \
0, x>frac12t
endcases$$
How can we see from the plot where the shock is?

pde numerical-methods characteristics

share|cite|improve this question

asked Mar 27 at 17:20

dxdydzdxdydz

49110

$endgroup$

add a comment | 

0

$begingroup$

I'm trying to implement the Finite-Volume method using local Lax-Friedrichs flux function. I found an example to the Riemann-problem, here is the solution plot:

So from characteristics $$u=begincases
1, x<frac12t \
0, x>frac12t
endcases$$
How can we see from the plot where the shock is?

pde numerical-methods characteristics

share|cite|improve this question

asked Mar 27 at 17:20

dxdydzdxdydz

49110

$endgroup$

add a comment | 

$begingroup$

I'm trying to implement the Finite-Volume method using local Lax-Friedrichs flux function. I found an example to the Riemann-problem, here is the solution plot:

So from characteristics $$u=begincases
1, x<frac12t \
0, x>frac12t
endcases$$
How can we see from the plot where the shock is?

pde numerical-methods characteristics

share|cite|improve this question

asked Mar 27 at 17:20

dxdydzdxdydz

49110

$endgroup$

share|cite|improve this question

asked Mar 27 at 17:20

dxdydzdxdydz

49110

share|cite|improve this question

asked Mar 27 at 17:20

dxdydzdxdydz

49110

asked Mar 27 at 17:20

dxdydzdxdydz

49110

asked Mar 27 at 17:20

dxdydzdxdydz

49110

1 Answer
1

active

oldest

votes

0

$begingroup$

If you want a numerical way of finding the shock, you could take your array $u[i,j]$ and consider the partial derivative in the x-direction:
dudx = diff(u) / dx

And then if you contour plot dudx, you should get an image where the shock shows up as yellow. This should be along the line $t = 2x$, which you can see in your plot above.

share|cite|improve this answer

answered Mar 27 at 18:50

gdepaulgdepaul

3759

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  I tried to follow your advice, and I got smth like this: imgur.com/FLzCyMo, seems a little bit awkward, I picked the number of nodes $x$ to be $100$ on $[-3,3], dt = 0.01$ Can you show what I should get?
  $endgroup$
  – dxdydz
  Mar 27 at 22:40
  
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  That’s more or less what you should get, with the exception of the loop that juts our from the origin to the left.
  $endgroup$
  – gdepaul
  Mar 27 at 23:10
  

  
  
  
  

add a comment | 

Your Answer

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "69"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3164805%2fshock-in-the-inviscid-burgers-equation%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

0

$begingroup$

If you want a numerical way of finding the shock, you could take your array $u[i,j]$ and consider the partial derivative in the x-direction:
dudx = diff(u) / dx

And then if you contour plot dudx, you should get an image where the shock shows up as yellow. This should be along the line $t = 2x$, which you can see in your plot above.

share|cite|improve this answer

answered Mar 27 at 18:50

gdepaulgdepaul

3759

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  I tried to follow your advice, and I got smth like this: imgur.com/FLzCyMo, seems a little bit awkward, I picked the number of nodes $x$ to be $100$ on $[-3,3], dt = 0.01$ Can you show what I should get?
  $endgroup$
  – dxdydz
  Mar 27 at 22:40
  
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  That’s more or less what you should get, with the exception of the loop that juts our from the origin to the left.
  $endgroup$
  – gdepaul
  Mar 27 at 23:10
  

  
  
  
  

add a comment | 

0

$begingroup$

If you want a numerical way of finding the shock, you could take your array $u[i,j]$ and consider the partial derivative in the x-direction:
dudx = diff(u) / dx

And then if you contour plot dudx, you should get an image where the shock shows up as yellow. This should be along the line $t = 2x$, which you can see in your plot above.

share|cite|improve this answer

answered Mar 27 at 18:50

gdepaulgdepaul

3759

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  I tried to follow your advice, and I got smth like this: imgur.com/FLzCyMo, seems a little bit awkward, I picked the number of nodes $x$ to be $100$ on $[-3,3], dt = 0.01$ Can you show what I should get?
  $endgroup$
  – dxdydz
  Mar 27 at 22:40
  
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  That’s more or less what you should get, with the exception of the loop that juts our from the origin to the left.
  $endgroup$
  – gdepaul
  Mar 27 at 23:10
  

  
  
  
  

add a comment | 

$begingroup$

If you want a numerical way of finding the shock, you could take your array $u[i,j]$ and consider the partial derivative in the x-direction:
dudx = diff(u) / dx

And then if you contour plot dudx, you should get an image where the shock shows up as yellow. This should be along the line $t = 2x$, which you can see in your plot above.

share|cite|improve this answer

answered Mar 27 at 18:50

gdepaulgdepaul

3759

$endgroup$

share|cite|improve this answer

answered Mar 27 at 18:50

gdepaulgdepaul

3759

answered Mar 27 at 18:50

gdepaulgdepaul

3759

answered Mar 27 at 18:50

gdepaulgdepaul

3759