Skip to content
Snippets Groups Projects
Commit 6a87b110 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

pp

parent 51bae575
No related branches found
No related tags found
No related merge requests found
Pipeline #
...@@ -17,8 +17,8 @@ ...@@ -17,8 +17,8 @@
File > Open > magnus.pro File > Open > magnus.pro
Run (button at the bottom of the left panel) Run (button at the bottom of the left panel)
------------------------------------------------------------------- */ ------------------------------------------------------------------- */
/*
This model solves a 2D potential flow around a cylinder or a naca airfoil /* This model solves a 2D potential flow around a cylinder or a naca airfoil
placed in a uniform "V_infinity" flow. Potential flows are defined by placed in a uniform "V_infinity" flow. Potential flows are defined by
V = grad phi => curl V = 0 V = grad phi => curl V = 0
...@@ -244,12 +244,11 @@ Resolution{ ...@@ -244,12 +244,11 @@ Resolution{
// While[]{} : iteration // While[]{} : iteration
// Print[{}, Format ..., File ...] : formatted display // Print[{}, Format ..., File ...] : formatted display
// A pseudo-Newton iteration is implemented here // A pseudo-Newton iteration is implemented here to determine the value
// to determine the value of Dmdt (in the Airfoil case) // of Dmdt (in the Airfoil case) that verifies Kutta's condition. The
// that verifies Kutta's condition. // run-time variable $newDmdt is used in Generate[A] whereas $circ,
// The run-time variable $newDmdt is used in Generate[A] // $dmdt, $argV and $phiTrailing are evaluated by the
// whereas $circ, $dmdt, $argV and $phiTrailing are evaluated // PostOperation[Trailing].
// by the PostOperation[Trailing].
DeleteFile["KJiter.txt"]; DeleteFile["KJiter.txt"];
...@@ -289,29 +288,40 @@ PostProcessing{ ...@@ -289,29 +288,40 @@ PostProcessing{
{ Name PotentialFlow; NameOfFormulation PotentialFlow; { Name PotentialFlow; NameOfFormulation PotentialFlow;
Quantity{ Quantity{
{Name phi; Value { {Name phi; Value {
Local{ [ {phi} ] ; In Dom_Vh; Jacobian Vol; } } } Term{ [ {phi} ] ; In Dom_Vh; Jacobian Vol; }
}
}
{ Name phiCont ; Value { { Name phiCont ; Value {
Local { [ { phiCont } ] ; In Dom_Vh ; Jacobian Vol ; } } } Term { [ { phiCont } ] ; In Dom_Vh ; Jacobian Vol ; }
}
}
{ Name phiDisc ; Value { { Name phiDisc ; Value {
Local { [ { phiDisc } ] ; In Dom_Vh ; Jacobian Vol ; } } } Term { [ { phiDisc } ] ; In Dom_Vh ; Jacobian Vol ; }
}
}
{Name velocity; Value { {Name velocity; Value {
Local { [ {d phi} ]; In Dom_Vh; Jacobian Vol; } } } Term { [ {d phi} ]; In Dom_Vh; Jacobian Vol; }
}
}
{Name normVelocity; Value { {Name normVelocity; Value {
Local { [ Norm[{d phi}] ]; In Dom_Vh; Jacobian Vol; } } } Term { [ Norm[{d phi}] ]; In Dom_Vh; Jacobian Vol; }
}
}
{Name pressure; Value { {Name pressure; Value {
Local { [-0.5*rho[]*SquNorm[ {d phi} ]]; Term { [-0.5*rho[]*SquNorm[ {d phi} ]]; In Dom_Vh; Jacobian Vol; }
In Dom_Vh; Jacobian Vol; } } } }
}
{Name Angle; Value { {Name Angle; Value {
Local{ [ argVTrail[{d phi}] ]; Term{ [ argVTrail[{d phi}] ]; In Dom_Vh; Jacobian Vol; }
In Dom_Vh; Jacobian Vol; } } } }
}
{ Name Circ; Value { Local { [ {Circ} ]; In Sur_Cut; } } } { Name Circ; Value { Term { [ {Circ} ]; In Sur_Cut; } } }
{ Name Dmdt; Value { Local { [ {Dmdt} ]; In Sur_Cut; } } } { Name Dmdt; Value { Term { [ {Dmdt} ]; In Sur_Cut; } } }
// Kutta-Jukowski approximation for Lift // Kutta-Jukowski approximation for Lift
{ Name LiftKJ; Value { Local { [ -rho[]*{Circ}*Velocity ]; { Name LiftKJ; Value { Term { [ -rho[]*{Circ}*Velocity ]; In Sur_Cut; } } }
In Sur_Cut; } } }
// Lift computed with the real pressure field // Lift computed with the real pressure field
{ Name Lift; { Name Lift;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment