Skip to content
Snippets Groups Projects
Verified Commit ab784af1 authored by Paul Dechamps's avatar Paul Dechamps :speech_balloon:
Browse files

(tests) Updated 2D transonic RAE mesh

Increased number of elements in the TE region and updated tests values.
parent 6b71dd47
No related branches found
No related tags found
1 merge request!1BLASTER v1.0
Pipeline #37686 failed
...@@ -42,7 +42,7 @@ def cfgInviscid(nthrds, verb): ...@@ -42,7 +42,7 @@ def cfgInviscid(nthrds, verb):
'Verb' : verb, # verbosity 'Verb' : verb, # verbosity
# Model (geometry or mesh) # Model (geometry or mesh)
'File' : os.path.dirname(os.path.dirname(os.path.abspath(__file__))) + '/models/dart/rae_2.geo', # Input file containing the model 'File' : os.path.dirname(os.path.dirname(os.path.abspath(__file__))) + '/models/dart/rae_2.geo', # Input file containing the model
'Pars' : {'xLgt' : 50, 'yLgt' : 50, 'msF': 10, 'msTe' : 0.01, 'msLe' : 0.001}, # parameters for input file model 'Pars' : {'xLgt' : 50, 'yLgt' : 50, 'msF': 10, 'msTe' : 0.0075, 'msLe' : 0.001}, # parameters for input file model
'Dim' : 2, # problem dimension 'Dim' : 2, # problem dimension
'Format' : 'gmsh', # save format (vtk or gmsh) 'Format' : 'gmsh', # save format (vtk or gmsh)
# Markers # Markers
...@@ -64,7 +64,7 @@ def cfgInviscid(nthrds, verb): ...@@ -64,7 +64,7 @@ def cfgInviscid(nthrds, verb):
'y_ref' : 0.0, # reference point for moment computation (y) 'y_ref' : 0.0, # reference point for moment computation (y)
'z_ref' : 0.0, # reference point for moment computation (z) 'z_ref' : 0.0, # reference point for moment computation (z)
# Numerical # Numerical
'LSolver' : 'PARDISO', # inner solver (Pardiso, MUMPS or GMRES) 'LSolver' : 'SparseLu', # inner solver (Pardiso, MUMPS or GMRES)
'G_fill' : 2, # fill-in factor for GMRES preconditioner 'G_fill' : 2, # fill-in factor for GMRES preconditioner
'G_tol' : 1e-5, # tolerance for GMRES 'G_tol' : 1e-5, # tolerance for GMRES
'G_restart' : 50, # restart for GMRES 'G_restart' : 50, # restart for GMRES
...@@ -132,14 +132,12 @@ def main(): ...@@ -132,14 +132,12 @@ def main():
# Test solution # Test solution
print(ccolors.ANSI_BLUE + 'PyTesting...' + ccolors.ANSI_RESET) print(ccolors.ANSI_BLUE + 'PyTesting...' + ccolors.ANSI_RESET)
tests = CTests() tests = CTests()
tests.add(CTest('Cl', isol.getCl(), 0.766, 5e-2)) tests.add(CTest('Cl', isol.getCl(), 0.759, 5e-2))
tests.add(CTest('Cd wake', vsol.Cdt, 0.0093, 1e-3, forceabs=True)) tests.add(CTest('Cd wake', vsol.Cdt, 0.0093, 1e-3, forceabs=True))
tests.add(CTest('Cd integral', isol.getCd() + vsol.Cdf, 0.0138, 1e-3, forceabs=True)) tests.add(CTest('Cd integral', isol.getCd() + vsol.Cdf, 0.0134, 1e-3, forceabs=True))
tests.add(CTest('Cdf', vsol.Cdf, 0.0069, 1e-3, forceabs=True)) tests.add(CTest('Cdf', vsol.Cdf, 0.0068, 1e-3, forceabs=True))
# if icfg['LSolver'] == 'PARDISO': if icfg['LSolver'] == 'SparseLu':
# tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 43, 0, forceabs=True)) tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 29, 0, forceabs=True))
# elif icfg['LSolver'] == 'SparseLu':
# tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 42, 0, forceabs=True))
tests.run() tests.run()
expResults = np.loadtxt(os.path.dirname(os.path.dirname(os.path.abspath(__file__))) + '/models/references/rae2822_AR138_case6.dat') expResults = np.loadtxt(os.path.dirname(os.path.dirname(os.path.abspath(__file__))) + '/models/references/rae2822_AR138_case6.dat')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment