Synthetic Turbulence Inflow Generator - Revision history

Jrwrigh: /* plotBias_TC.m */

2021-10-12T17:27:45Z

‎plotBias_TC.m

Jrwrigh: /* Creation of STGInflow.dat */ : Updated to the "modern" file format

2021-09-16T14:24:03Z

‎Creation of STGInflow.dat: : Updated to the "modern" file format

Jopa6460: /* Creation of STGRand.dat */

2021-03-04T04:03:16Z

‎Creation of STGRand.dat

Jopa6460: /* Creation of STGSpectra.dat */

2021-03-04T04:02:28Z

‎Creation of STGSpectra.dat

Jopa6460: /* Creation of STGSpectra.dat */

2021-03-04T04:02:06Z

‎Creation of STGSpectra.dat

Jopa6460: /* Creation of STGSpectra.dat */

2021-03-04T04:01:12Z

‎Creation of STGSpectra.dat

Jopa6460: /* Solver.inp */

2021-03-02T05:45:28Z

‎Solver.inp

Jopa6460: /* Solver.inp */

2021-03-02T05:40:10Z

‎Solver.inp

Jopa6460: /* Running the code */

2021-03-02T05:36:27Z

‎Running the code

Jopa6460: /* PlotSTG_meanFeilds.m */

2021-03-02T05:29:22Z

‎PlotSTG_meanFeilds.m

@@ Line 459: / Line 459: @@
 ==== plotBias_TC.m ====
- This script output plots that are meant to demonstrate bias and the expected time convergence relative to the second order stresses. Better visualization of the time convergence can be visualized via analyzeTC.m. That script requires output from this script. In a matlab file save the SpanWise_rms4k8k and SpanWise_rms4k8kMax with the loadedCase string appending it. Do the same in the same matlab file with other cases you wish to compare against.
+This script output plots that are meant to demonstrate bias and the expected time convergence relative to the second order stresses. Better visualization of the time convergence can be visualized via analyzeTC.m. That script requires output from this script. In a matlab file save the SpanWise_rms4k8k and SpanWise_rms4k8kMax with the loadedCase string appending it. Do the same in the same matlab file with other cases you wish to compare against.
 After this script has executed enter the following in the command line to save the required data for analyzeTC.m:

@@ Line 73: / Line 73: @@
 ===Creation of STGInflow.dat===
-This auxiliary file will be required not only for the running of PHASTA, but also, will be required for the random number optimization. The user will need mean flow field data at the position of the inflow from some  preliminary simulation such as RANS data or data repository. The code requires this file to have the following format:
+This auxiliary file will be required not only for the running of PHASTA, but also, will be required for the random number optimization. The user will need mean flow field data at the position of the inflow from some  preliminary simulation such as RANS data or data repository. The file format itself can take a few different forms, but the primary is given as the 14 columns version:
 '''Where:'''
-*N is the number of rows that this file contains. It is the number of nodes in the specially varying direction.
+* <code>N</code> is the number of rows that this file contains. It is the number of nodes in the specially varying direction.
-*d2wall is the distance to the wall
+* <code>d2wall</code> is the distance to the wall
-*U_i are the three components of the average velocity feild
+* <code>U_i</code> are the three components of the average velocity field
-*uu, vv, ww, uv, uw, vw are the averages of the Reynolds' stress tensor
+* <code>uu, vv, ww, uv, uw, vw</code> are the Reynolds' stress tensor components
-*sclr1 and sclr2 are the scalars that represent the turbulence model.
+* <code>sclr1</code> and <code>sclr2</code> are the scalars that represent the turbulence model.
 The input profiles are varied by the type of auxiliary simulation. So the post-processing of that preliminary simulation will be varied and so the generation of this data is left to the user. However for creating the STGInflow.dat file it is useful to know the format that the fortran file reader is setup for. The below chunk of MATLAB code is useful once the profiles are known:
@@ Line 109: / Line 100: @@
    for j=1:nPoints
        i=index(j);
-       tmp=[d(i) u(i,:) R(i,:) nut(i)];
+       tmp=[d(i) u(i,:) R(i,:) nut(i) lt(i) eps(i)];
        if (j~=nPoints)
-           fprintf(fid,'%.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E\n',tmp);
+           fprintf(fid,'%.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E\n',tmp);
        elseif (j==nPoints)
-           fprintf(fid,'%.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E',tmp);
+           fprintf(fid,'%.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E %.12E',tmp);
        end
    end

← Older revision		Revision as of 04:03, 4 March 2021
Line 322:		Line 322:

	The result of this file will output npool number of valid STGRand.dat number sets named according to the STG_basename and the process that created it. Additionally, the method will output in the command window the Dmag for each process. The lowest of these numbers is the best STGRand.dat file that was created from the MATLAB proccesses.		The result of this file will output npool number of valid STGRand.dat number sets named according to the STG_basename and the process that created it. Additionally, the method will output in the command window the Dmag for each process. The lowest of these numbers is the best STGRand.dat file that was created from the MATLAB proccesses.
		+	[[File:Command_postChooseRand.PNG]]

	=== Pre-Simulation-Analyzation ===		=== Pre-Simulation-Analyzation ===

@@ Line 208: / Line 208: @@
 * alphaWave: This is the logarithmic growth of the spacing between one wavenumber and the next. It may be set from 0.05 to 0.01. Normally we set this to 0.01.
 * setEps: This is set to the dissipation rate of the desired inflow spectrum.
 === Creation of STGRand.dat ===

@@ Line 209: / Line 209: @@
 * setEps: This is set to the dissipation rate of the desired inflow spectrum.
-[[command:Command_postChooseRand.PNG]]
+[[Command_postChooseRand.PNG:Command_postChooseRand.PNG]]
 === Creation of STGRand.dat ===

← Older revision		Revision as of 05:45, 2 March 2021
Line 529:		Line 529:
	568 Number of Points for Eng. Calculation: 0.0		568 Number of Points for Eng. Calculation: 0.0
	569 }		569 }
		+
		+	Many of these are known from the preceding steps with the exception of a few:
		+	* STG Dump Spectra: When this is 1 a PHASTA run will dump spectra files that getSTGspectra.m can read-in.
		+	* STG Start Step: this is the timestep that PHASTA will turn on the STG boundary condition
		+	* STG Turn Off Fluctuations: Turning this to True will provide the mean flowfield at the inflow without the STG fluctuations.
		+	* STG SurfID: This will be the SurfId that was assigned at the problem definition step of the workflow.

@@ Line 496: / Line 496: @@
 === Solver.inp===
-This is the final remaining step in the workflow to include the STG inflow boundary condition.
+This is the final remaining step in the workflow to include the STG inflow boundary condition. Below are the default conditions pertaining to the STG boundary and initial condition
 BOUNDARY AND INITIAL CONDITIONS
 {
       Ramp Down Nut in STG : False
       x Position Nut STG Ramp 0 : 0.123
       x Position Nut STG Ramp 1 : 0.456
       Load and set STG IC: False
 }

← Older revision		Revision as of 05:36, 2 March 2021
Line 491:		Line 491:

	== Running the code ==		== Running the code ==
		+
		+	In addition to the case set-up the additional auxiliary files within the n-procs case directory that are needed for PHASTA to use the STG boundary condition are STGInflow.dat, STGRand.dat, and STGSpectra.dat. The only other specifications needed are within the Solver.inp file.
		+
	=== Solver.inp===		=== Solver.inp===
		+
		+	This is the final remaining step in the workflow to include the STG inflow boundary condition.
		+
		+	137 BOUNDARY AND INITIAL CONDITIONS
		+	138 {
		+	139 Remove No Penetration Vel BC: False
		+	140 Add Random Fluctuations to Velocity: False
		+	141 Dump Coordinates for Shepard Interpolation: False
		+	142 Ramp Down TTS : False
		+	143 Ramp Down Nut Near Wall : False
		+	144 dwal Position Nut Wall Ramp 0 : 0.123 # values only read if True
		+	145 dwal Position Nut Wall Ramp 1 : 0.456
		+	146 Ramp Down Nut in STG : False
		+	147 x Position Nut STG Ramp 0 : 0.123
		+	148 x Position Nut STG Ramp 1 : 0.456
		+	149 Load and set 2D IC: False
		+	150 Load and set 3D IC: False
		+	151 Load and set STG IC: False
		+	152 Modify Top Surface of Bump Domain: False #iMoveTopWall
		+	153 }
		+
		+	547 #STG input Parameters
		+	548 {
		+	549
		+	550 Use STGinflow BCs: False
		+	551 STG Dump Spectra:0
		+	552 STG Start Step: 0
		+	553 STG Turn Off Fluctuations: False
		+	554 STG SurfID: NODEFAULT
		+	555 STG BL Height: 1.0
		+	556 STG U_0: 1.0
		+	557 STG Number of Modes: 1.0
		+	558 STG Mode Growth: 1.0
		+	559 STG Mesh Growth: 1.0
		+	560 STG Dissip. Rate: 1.0e6
		+	561 STG Channel: False
		+	562 STG SST Method: 1
		+	563 STG SST uTau: 1.0
		+	564 STG SST Eddy Viscosity Factor: 1.0
		+	565 Collect STG Energy Spectrum Data: False
		+	566 STG Energy Distance Tolerance: 1.0
		+	567 Energy Spectrum at Distances: 0.0 0.0 0.0
		+	568 Number of Points for Eng. Calculation: 0.0
		+	569 }

@@ Line 487: / Line 487: @@
 ==== PlotSTG_meanFeilds.m ====
 == Running the code ==
 === Solver.inp===