更新算例

P01_Flate_Plate_1_Part/bin/GenerateGrid.hypara

@@ -6,10 +6,10 @@
 //!                    2 -- Flat plate boundary layer with 1 part.
 //!                    3 -- Flat plate boundart later with 2 parts.
 //!                    4 -- Poiseuille flow. 
-//int gridDataBaseIndex = 1;
-//string gridDataBaseFileName = "./bin/GridDataBase_TaylorGreenVortex.hypara";
-int gridDataBaseIndex = 2;
-string gridDataBaseFileName = "./bin/GridDataBase_Flat_Plate_Boundary_Layer_1_part.hypara";
+int gridDataBaseIndex = 1;
+string gridDataBaseFileName = "./bin/GridDataBase_TaylorGreenVortex.hypara";
+//int gridDataBaseIndex = 2;
+//string gridDataBaseFileName = "./bin/GridDataBase_Flat_Plate_Boundary_Layer_1_part.hypara";
 
 //! gridOutputFile: The name of output grid.
 string gridOutputFile = "./grid/test_grid.fts";

P01_Flate_Plate_1_Part/bin/GridDataBase_TaylorGreenVortex.hypara

@@ -1,7 +1,7 @@
 int nBlocks = 1;
 
-int nNodeX = 32;
-int nNodeY = 32;
+int nNodeX = 64;
+int nNodeY = 64;
 int nNodeZ = 1;
 
 double lengthX = 6.283185307179586476925286766559;

P01_Flate_Plate_1_Part/bin/cfd_para.hypara

@@ -921,4 +921,8 @@ int     numberOfMovingBodies       = 2;
 #########################################################################
 int isPlotVolumeField = 0;
 
-
+#########################################################################
+#                       user define parameter                           #
+#########################################################################
+int UDFInitFluid = 1;
+string UDFInitFluidType = "TGV";

P01_Flate_Plate_1_Part/bin/cfd_para_subsonic.hypara

@@ -7,9 +7,9 @@
 // intervalStepForce:   the step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
 // intervalStepRes:    the step intervals for residual 'res.dat' saved.
 
-int maxSimuStep       = 15000;
+int maxSimuStep       = 5000;
 
-int intervalStepFlow  = 1;
+int intervalStepFlow  = 5;
 int intervalStepPlot  = 100;
 int intervalStepForce = 100;
 int intervalStepRes   = 10;
@@ -17,7 +17,7 @@ int intervalStepRes   = 10;
 // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
 //                   0 -- no precondition process. (default, mach > 0.3)
 //                   1 -- carry out precondition process. (mach number <= 0.3)
-int ifLowSpeedPrecon = 0;
+int ifLowSpeedPrecon = 1;
 
 #########################################################################
 #                       Inflow Parameter                                #
@@ -44,12 +44,12 @@ int ifLowSpeedPrecon = 0;
 // forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
 // TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
 
-double refMachNumber = 0.1;
+double refMachNumber = 0.3;
 double attackd       = 0.00;
 double angleSlide    = 0.00;
 
 int inflowParaType = 0;
-double refReNumber = 2.0e5;
+double refReNumber = 100;
 double refDimensionalTemperature = 288.15;
 
 //int inflowParaType = 1;

P01_Flate_Plate_1_Part/bin/particlepoint_para.hypara

@@ -62,17 +62,17 @@ int useSlipVelocityType =0;
 int forceDragType = 1;
 
 //! The lift force
-int forceMagnusType = 1;
-int forceSaffmanType = 1;
+int forceMagnusType = 0;
+int forceSaffmanType = 0;
 
 //! The stress-gradient force
-int forceFluidAccelerationType = 1;
+int forceFluidAccelerationType = 0;
 
 //! The added-mass force
-int forceAddedMassType = 1;
+int forceAddedMassType = 0;
 
 //! The Basset history force;
-int forceBassetType = 1;
+int forceBassetType = 0;
 
 //! #########################################################################
 //! # Output Variable                                                       #