public void BasicRectangularFlowTest2()
{
CanalStretch canal = new CanalStretch("stretch1", 602, 2.5486, new RectangularSection(4, 0.014, 0.0004));
double rh = canal.CanalSection.GetHydraulicRadius(0.694);
double mv = canal.CanalSection.GetManningVelocity(0.694);
double fr2 = Math.Pow(canal.GetFroudeNumber(0.694), 2.0);
double flow = canal.CanalSection.GetManningFlow(0.694);
double criticalDepth = canal.CanalSection.GetCriticalWaterLevel(flow);
double criticalSlope = canal.CanalSection.GetCriticalSlope(flow);
double normalDepth = canal.CanalSection.GetNormalWaterLevel(flow);
RungeKutta solver = new RungeKutta(1, canal.FlowEquation());
// Act
double result = solver.Solve(0, 0.694);
double actualRh = (4.0 * 0.694) / (4.0 + 2.0 * 0.694);
Assert.IsTrue(rh - actualRh <= double.Epsilon);
Assert.IsTrue(mv - 0.91811 <= 0.0001);
Assert.IsTrue(fr2 - 1.90669 <= 0.0001);
Assert.IsTrue(result - 0.694 <= 0.00001);
Assert.IsTrue(criticalSlope > 0.0004);
Assert.IsTrue(normalDepth > criticalDepth);
Assert.IsTrue(Math.Abs(criticalDepth - 0.34593708152888353) <= double.Epsilon);
}
public void BasicCanalTest2()
{
CanalManager canalManager = new CanalManager();
Canal canal = new Canal();
canal.CanelEdges.Add(new CanalEdge()
{
Id = "SourceCanal",
WaterLevel = 2.9
});
canal.CanelEdges.Add(new CanalEdge()
{
Id = "EndCanal"
});
CanalStretch canalStretch = new CanalStretch("stretch1", 602, 20.32, new RectangularSection(5, 0.028, 0));
canalStretch.FromNode = canal.CanelEdges.First();
canalStretch.ToNode = canal.CanelEdges.Last();
canal.CanalStretches.Add(canalStretch);
canalManager.SetCanal(canal);
var result = canalManager.ExecuteCanalSimulation();
Assert.IsTrue(result != null);
}
public void BasicRectangularFlowTest()
{
CanalStretch canal = new CanalStretch("stretch1", 602, 20.32, new RectangularSection(5, 0.028, 0));
double rh = canal.CanalSection.GetHydraulicRadius(2.9);
double mv = canal.CanalSection.GetManningVelocity(2.9);
double fr2 = Math.Pow(canal.GetFroudeNumber(2.9), 2.0);
double criticalDepth = canal.CanalSection.GetCriticalWaterLevel(20.32);
double criticalSlope = canal.CanalSection.GetCriticalSlope(20.32);
double flow = canal.CanalSection.GetManningFlow(canal.CanalSection.GetCriticalWaterLevel(20.32));
double normalDepth = canal.CanalSection.GetNormalWaterLevel(20.32);
RungeKutta solver = new RungeKutta(1, canal.FlowEquation());
// Act
double result = solver.Solve(0, 2.9);
double actualRh = (5.0 * 2.9) / (5.0 + 2.0 * 2.9);
Assert.IsTrue(rh - actualRh <= double.Epsilon);
Assert.IsTrue(mv - 0 <= double.Epsilon);
Assert.IsTrue(fr2 - 0.0690 <= 0.0001);
Assert.IsTrue(result - 2.89888 <= 0.00001);
Assert.IsTrue(criticalSlope > 0);
Assert.IsTrue(Math.Abs(flow - 0) < Double.Epsilon);
Assert.IsTrue(Math.Abs(criticalDepth - 1.1897678822813287) <= double.Epsilon);
}
public void GetCanalStretchForAbsoluteXTest()
{
CanalEdge initialNode = new CanalEdge()
{
Id = "InitNode"
};
CanalEdge midNode = new CanalEdge()
{
Id = "MidNode"
};
CanalEdge endNode = new CanalEdge()
{
Id = "EndNode"
};
CanalStretch canal = new CanalStretch("stretch1", 600, 20.32, new RectangularSection(5, 0.028, 0));
canal.FromNode = initialNode;
canal.ToNode = midNode;
CanalStretch canal2 = new CanalStretch("stretch2", 400, 20.32, new RectangularSection(5, 0.028, 0.004));
canal2.FromNode = midNode;
canal2.ToNode = endNode;
Canal globalCanal = new Canal();
globalCanal.CanalStretches = new List <ICanalStretchModel> {
canal, canal2
};
globalCanal.CanelEdges = new List <ICanalEdge> {
initialNode, midNode, endNode
};
ICanalStretchModel result = globalCanal.GetCanalStretchForAbsoluteX(0);
Assert.AreEqual(canal, result);
result = globalCanal.GetCanalStretchForAbsoluteX(300);
Assert.AreEqual(canal, result);
result = globalCanal.GetCanalStretchForAbsoluteX(600);
Assert.AreEqual(canal2, result);
result = globalCanal.GetCanalStretchForAbsoluteX(900);
Assert.AreEqual(canal2, result);
result = globalCanal.GetCanalStretchForAbsoluteX(1000);
Assert.AreEqual(canal2, result);
}
public void SaveAndLoadConfiguration()
{
// Arrange
CanalManager canalManager = new CanalManager();
ConfigurationManager.AppSettings[ServerConfigConsts.RelativeModelsLocation] = "ScienceAndMaths";
CanalConfigurationLoader loader = new CanalConfigurationLoader();
loader.CanalManager = canalManager;
Canal canal = new Canal
{
Id = "TestConfig"
};
canal.CanelEdges.Add(new CanalEdge()
{
Id = "SourceCanal",
WaterLevel = 2.9
});
canal.CanelEdges.Add(new CanalEdge()
{
Id = "EndCanal"
});
CanalStretch canalStretch = new CanalStretch("stretch1", 602, 20.32, new RectangularSection(5, 0.028, 0));
canalStretch.FromNode = canal.CanelEdges.First();
canalStretch.ToNode = canal.CanelEdges.Last();
canal.CanalStretches.Add(canalStretch);
CanalConfiguration canalConfiguration = new CanalConfiguration();
canalConfiguration.ConfigId = "TestConfig";
loader.SaveCanalConfiguration(canal);
loader.LoadCanalConfiguration();
Assert.IsNotNull(canalManager.Canal);
Assert.AreEqual(2, canalManager.Canal.CanelEdges.Count);
Assert.AreEqual(1, canalManager.Canal.CanalStretches.Count);
}
