From ab784af17b1a77be9ff13ed05804621c43b6e7d5 Mon Sep 17 00:00:00 2001
From: Paul Dechamps <paul.dechamps@uliege.be>
Date: Mon, 1 Jul 2024 15:34:17 +0200
Subject: [PATCH] (tests) Updated 2D transonic RAE mesh
Increased number of elements in the TE region and updated tests values.
---
blast/validation/raeValidation.py | 16 +++++++---------
1 file changed, 7 insertions(+), 9 deletions(-)
diff --git a/blast/validation/raeValidation.py b/blast/validation/raeValidation.py
index 7e75a65..2ee1444 100644
--- a/blast/validation/raeValidation.py
+++ b/blast/validation/raeValidation.py
@@ -42,7 +42,7 @@ def cfgInviscid(nthrds, verb):
'Verb' : verb, # verbosity
# 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
- '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
'Format' : 'gmsh', # save format (vtk or gmsh)
# Markers
@@ -64,7 +64,7 @@ def cfgInviscid(nthrds, verb):
'y_ref' : 0.0, # reference point for moment computation (y)
'z_ref' : 0.0, # reference point for moment computation (z)
# 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_tol' : 1e-5, # tolerance for GMRES
'G_restart' : 50, # restart for GMRES
@@ -132,14 +132,12 @@ def main():
# Test solution
print(ccolors.ANSI_BLUE + 'PyTesting...' + ccolors.ANSI_RESET)
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 integral', isol.getCd() + vsol.Cdf, 0.0138, 1e-3, forceabs=True))
- tests.add(CTest('Cdf', vsol.Cdf, 0.0069, 1e-3, forceabs=True))
- # if icfg['LSolver'] == 'PARDISO':
- # tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 43, 0, forceabs=True))
- # elif icfg['LSolver'] == 'SparseLu':
- # tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 42, 0, forceabs=True))
+ tests.add(CTest('Cd integral', isol.getCd() + vsol.Cdf, 0.0134, 1e-3, forceabs=True))
+ tests.add(CTest('Cdf', vsol.Cdf, 0.0068, 1e-3, forceabs=True))
+ if icfg['LSolver'] == 'SparseLu':
+ tests.add(CTest('Iterations', len(aeroCoeffs['Cl']), 29, 0, forceabs=True))
tests.run()
expResults = np.loadtxt(os.path.dirname(os.path.dirname(os.path.abspath(__file__))) + '/models/references/rae2822_AR138_case6.dat')
--
GitLab