Ggtkuk

Can a creature tell when it has been affected by a Divination wizard's Portent?

How to back up a running remote server?

Why do Ichisongas hate elephants and hippos?

Will tsunami waves travel forever if there was no land?

How to stop co-workers from teasing me because I know Russian?

Were there two appearances of Stan Lee?

Single Colour Mastermind Problem

Can solid acids and bases have pH values? If not, how are they classified as acids or bases?

How deep to place a deadman anchor for a slackline?

Why does processed meat contain preservatives, while canned fish needs not?

Can fracking help reduce CO2?

Confused by notation of atomic number Z and mass number A on periodic table of elements

Pythonic way to find the last position in a string not matching a regex

Is thermodynamics only applicable to systems in equilibrium?

Any examples of headwear for races with animal ears?

Phrase for the opposite of "foolproof"

Is GOCE a satellite or aircraft?

What's the metal clinking sound at the end of credits in Avengers: Endgame?

What are the spoon bit of a spoon and fork bit of a fork called?

Examples of non trivial equivalence relations , I mean equivalence relations without the expression " same ... as" in their definition?

Does a creature that is immune to a condition still make a saving throw?

Toggle Overlays shortcut?

What is a Recurrent Neural Network?

Will a top journal at least read my introduction?

ContourPlot — How do I color by contour curvature?

Custom contour labels in ContourPlotListContourPlot is blocking my geometryHow to plot the contour of f[x,y]==0 if always f[x,y]>=0Contour coloring and (List)ContourPlot projectionMore stream lines in a ListStreamPlotContourPlot - unequal contour spacingContourPlot color problems3D Stack of Disks with dedicated height plotsHow to color Contours in ContourPlot with custom ColorFunctionChanging the color of a specific curve in ContourPlot

11

2

$begingroup$

I'm plotting the stream lines of fluid flow past a cylinder, and I would like the colors to increase with contour curvature (i.e. increase as the velocity of the flow increases. Here's a MWE that seems to color it based on the the y-axis value:

ψ[r_, θ_] := U (r - a^2/r) Sin[θ]
r = Sqrt[x^2 + y^2];
θ = ArcSin[y/r];

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, 
 x, -5,5, y, -5, 5, 
 Contours -> 10 Table[i, i, -10, 10, 0.025]
 ];

cyl = Graphics[Disk[0, 0, 1]];

Show[stream, cyl]

plotting color

share|improve this question

edited Mar 20 at 21:11

m_goldberg

89.3k873200

asked Mar 20 at 19:05

dpholmesdpholmes

350111

$endgroup$

add a comment | 

11

2

$begingroup$

I'm plotting the stream lines of fluid flow past a cylinder, and I would like the colors to increase with contour curvature (i.e. increase as the velocity of the flow increases. Here's a MWE that seems to color it based on the the y-axis value:

ψ[r_, θ_] := U (r - a^2/r) Sin[θ]
r = Sqrt[x^2 + y^2];
θ = ArcSin[y/r];

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, 
 x, -5,5, y, -5, 5, 
 Contours -> 10 Table[i, i, -10, 10, 0.025]
 ];

cyl = Graphics[Disk[0, 0, 1]];

Show[stream, cyl]

plotting color

share|improve this question

edited Mar 20 at 21:11

m_goldberg

89.3k873200

asked Mar 20 at 19:05

dpholmesdpholmes

350111

$endgroup$

add a comment | 

$begingroup$

I'm plotting the stream lines of fluid flow past a cylinder, and I would like the colors to increase with contour curvature (i.e. increase as the velocity of the flow increases. Here's a MWE that seems to color it based on the the y-axis value:

ψ[r_, θ_] := U (r - a^2/r) Sin[θ]
r = Sqrt[x^2 + y^2];
θ = ArcSin[y/r];

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, 
 x, -5,5, y, -5, 5, 
 Contours -> 10 Table[i, i, -10, 10, 0.025]
 ];

cyl = Graphics[Disk[0, 0, 1]];

Show[stream, cyl]

plotting color

share|improve this question

edited Mar 20 at 21:11

m_goldberg

89.3k873200

asked Mar 20 at 19:05

dpholmesdpholmes

350111

$endgroup$

ψ[r_, θ_] := U (r - a^2/r) Sin[θ]
r = Sqrt[x^2 + y^2];
θ = ArcSin[y/r];

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, 
 x, -5,5, y, -5, 5, 
 Contours -> 10 Table[i, i, -10, 10, 0.025]
 ];

cyl = Graphics[Disk[0, 0, 1]];

Show[stream, cyl]

share|improve this question

edited Mar 20 at 21:11

m_goldberg

89.3k873200

asked Mar 20 at 19:05

dpholmesdpholmes

350111

share|improve this question

edited Mar 20 at 21:11

m_goldberg

89.3k873200

asked Mar 20 at 19:05

dpholmesdpholmes

350111

edited Mar 20 at 21:11

m_goldberg

89.3k873200

edited Mar 20 at 21:11

m_goldberg

89.3k873200

asked Mar 20 at 19:05

dpholmesdpholmes

350111

asked Mar 20 at 19:05

dpholmesdpholmes

350111

1 Answer
1

active

oldest

votes

14

$begingroup$

f = ψ[r, θ] /. U -> 10, a -> 1;
gradf = D[f, x, y, 1];
Hessf = D[f, x, y, 2];
normal = gradf[[1]]/Sqrt[gradf[[1]].gradf[[1]]] // Simplify;
secondfundamentalform = -PseudoInverse[gradf].Hessf // ComplexExpand // Simplify;
tangent = RotationMatrix[Pi/2].normal;
curvaturevector = (secondfundamentalform.tangent).tangent // Simplify;
signedcurvature = curvaturevector.normal // Simplify;

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, x, -5, 5, y, -5, 5,
 Contours -> 10 Table[i, i, -10, 10, 0.2],
 ContourShading -> None
 ];
curvatureplot = DensityPlot[signedcurvature, x, -5, 5, y, -5, 5,
 ColorFunction -> "RedGreenSplit",
 ClippingStyle -> Automatic,
 PlotPoints -> 100,
 PlotRange -> -1, 1
 ];
Show[
 curvatureplot,
 stream,
 cyl
 ]

share|improve this answer

edited Mar 21 at 8:40

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171

$endgroup$

add a comment | 

Your Answer

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "387"
;
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: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
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
,
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%2fmathematica.stackexchange.com%2fquestions%2f193665%2fcontourplot-how-do-i-color-by-contour-curvature%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

14

$begingroup$

f = ψ[r, θ] /. U -> 10, a -> 1;
gradf = D[f, x, y, 1];
Hessf = D[f, x, y, 2];
normal = gradf[[1]]/Sqrt[gradf[[1]].gradf[[1]]] // Simplify;
secondfundamentalform = -PseudoInverse[gradf].Hessf // ComplexExpand // Simplify;
tangent = RotationMatrix[Pi/2].normal;
curvaturevector = (secondfundamentalform.tangent).tangent // Simplify;
signedcurvature = curvaturevector.normal // Simplify;

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, x, -5, 5, y, -5, 5,
 Contours -> 10 Table[i, i, -10, 10, 0.2],
 ContourShading -> None
 ];
curvatureplot = DensityPlot[signedcurvature, x, -5, 5, y, -5, 5,
 ColorFunction -> "RedGreenSplit",
 ClippingStyle -> Automatic,
 PlotPoints -> 100,
 PlotRange -> -1, 1
 ];
Show[
 curvatureplot,
 stream,
 cyl
 ]

share|improve this answer

edited Mar 21 at 8:40

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171

$endgroup$

add a comment | 

14

$begingroup$

f = ψ[r, θ] /. U -> 10, a -> 1;
gradf = D[f, x, y, 1];
Hessf = D[f, x, y, 2];
normal = gradf[[1]]/Sqrt[gradf[[1]].gradf[[1]]] // Simplify;
secondfundamentalform = -PseudoInverse[gradf].Hessf // ComplexExpand // Simplify;
tangent = RotationMatrix[Pi/2].normal;
curvaturevector = (secondfundamentalform.tangent).tangent // Simplify;
signedcurvature = curvaturevector.normal // Simplify;

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, x, -5, 5, y, -5, 5,
 Contours -> 10 Table[i, i, -10, 10, 0.2],
 ContourShading -> None
 ];
curvatureplot = DensityPlot[signedcurvature, x, -5, 5, y, -5, 5,
 ColorFunction -> "RedGreenSplit",
 ClippingStyle -> Automatic,
 PlotPoints -> 100,
 PlotRange -> -1, 1
 ];
Show[
 curvatureplot,
 stream,
 cyl
 ]

share|improve this answer

edited Mar 21 at 8:40

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171

$endgroup$

add a comment | 

$begingroup$

f = ψ[r, θ] /. U -> 10, a -> 1;
gradf = D[f, x, y, 1];
Hessf = D[f, x, y, 2];
normal = gradf[[1]]/Sqrt[gradf[[1]].gradf[[1]]] // Simplify;
secondfundamentalform = -PseudoInverse[gradf].Hessf // ComplexExpand // Simplify;
tangent = RotationMatrix[Pi/2].normal;
curvaturevector = (secondfundamentalform.tangent).tangent // Simplify;
signedcurvature = curvaturevector.normal // Simplify;

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, x, -5, 5, y, -5, 5,
 Contours -> 10 Table[i, i, -10, 10, 0.2],
 ContourShading -> None
 ];
curvatureplot = DensityPlot[signedcurvature, x, -5, 5, y, -5, 5,
 ColorFunction -> "RedGreenSplit",
 ClippingStyle -> Automatic,
 PlotPoints -> 100,
 PlotRange -> -1, 1
 ];
Show[
 curvatureplot,
 stream,
 cyl
 ]

share|improve this answer

edited Mar 21 at 8:40

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171

$endgroup$

f = ψ[r, θ] /. U -> 10, a -> 1;
gradf = D[f, x, y, 1];
Hessf = D[f, x, y, 2];
normal = gradf[[1]]/Sqrt[gradf[[1]].gradf[[1]]] // Simplify;
secondfundamentalform = -PseudoInverse[gradf].Hessf // ComplexExpand // Simplify;
tangent = RotationMatrix[Pi/2].normal;
curvaturevector = (secondfundamentalform.tangent).tangent // Simplify;
signedcurvature = curvaturevector.normal // Simplify;

stream = ContourPlot[
 ψ[r, θ] /. U -> 10, a -> 1, x, -5, 5, y, -5, 5,
 Contours -> 10 Table[i, i, -10, 10, 0.2],
 ContourShading -> None
 ];
curvatureplot = DensityPlot[signedcurvature, x, -5, 5, y, -5, 5,
 ColorFunction -> "RedGreenSplit",
 ClippingStyle -> Automatic,
 PlotPoints -> 100,
 PlotRange -> -1, 1
 ];
Show[
 curvatureplot,
 stream,
 cyl
 ]

share|improve this answer

edited Mar 21 at 8:40

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171

answered Mar 20 at 19:28

Henrik SchumacherHenrik Schumacher

61.3k585171