incompressible_flow/Inc_Combustion/1__premixed_hydrogen/H2_burner.cfg

@@ -0,0 +1,144 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% SU2 configuration file                                                       %
+% Case description: Laminar premixed flame stabilized on isothermal burner     %
+% Author: Evert Bunschoten                                                     %
+% Institution: Delft University of Technology                                  %
+% Date: 16/06/2023                                                             %
+% File Version 7.5.1 "Blackbird"                                               %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER = INC_NAVIER_STOKES
+KIND_TURB_MODEL= NONE
+%
+% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
+%
+INC_DENSITY_MODEL= VARIABLE
+INC_ENERGY_EQUATION = YES
+INC_DENSITY_INIT= 1.00
+INC_VELOCITY_INIT= (0.575, 0.0, 0.0)
+INC_TEMPERATURE_INIT= 300.0
+INC_NONDIM= DIMENSIONAL
+%
+CUSTOM_OUTPUTS= 'outlet_YNO : AreaAvg{SPECIES[3]}[outlet];\
+                 outlet_T : AreaAvg{TEMPERATURE}[outlet];\
+                 inlet_P : AreaAvg{PRESSURE}[inlet];'
+OBJECTIVE_FUNCTION= CUSTOM_OBJFUNC
+CUSTOM_OBJFUNC= outlet_YNO
+% -------------------- FLUID MODEL --------------------------------------- %
+%
+FLUID_MODEL= FLUID_FLAMELET
+
+% List the MLP files describing the respective architectures.
+INTERPOLATION_METHOD= MLP
+FILENAMES_INTERPOLATOR= (MLP_TD1.mlp, MLP_TD2.mlp, MLP_PD.mlp, MLP_SPV.mlp, MLP_PNO.mlp, MLP_null.mlp)
+
+% FGM controlling variable names (must include mixture fraction for partial 
+% premixed cases)
+PREFERENTIAL_DIFFUSION=YES
+CONTROLLING_VARIABLE_NAMES=(ProgressVariable,EnthalpyTot,MixtureFraction)
+
+CONTROLLING_VARIABLE_SOURCE_NAMES=(ProdRateTot_PV,NULL,NULL)
+
+USER_SCALAR_NAMES= (Y_NO)
+USER_SOURCE_NAMES = ( \
+    Y_dot_net-NO,  NULL \
+)
+
+LOOKUP_NAMES=(Heat_Release)
+% -------------------- SCALAR TRANSPORT ---------------------------------------%
+%
+KIND_SCALAR_MODEL= FLAMELET      
+DIFFUSIVITY_MODEL= FLAMELET
+VISCOSITY_MODEL= FLAMELET
+CONDUCTIVITY_MODEL= FLAMELET
+
+FLAME_INIT_METHOD= SPARK
+SPARK_INIT= (0.001, 0.0004, 0.00, 1e-4, 100, 5)
+SPARK_REACTION_RATES=(1000, 0, 0)
+
+CONV_NUM_METHOD_SPECIES= BOUNDED_SCALAR
+
+MUSCL_SPECIES= YES
+SLOPE_LIMITER_SPECIES= NONE
+TIME_DISCRE_SPECIES= EULER_IMPLICIT
+% SCALAR CLIPPING
+SPECIES_CLIPPING= YES
+SPECIES_CLIPPING_MIN= -1.5064702033996582e+00 -3.4190245000000000e+06 +5.0012799911201000e-03   0
+SPECIES_CLIPPING_MAX= +4.6642374992370605e-01 +2.1569957500000000e+06 +1.4977467060089111e-01   1
+%
+% Inlet and initial condition corresponding to mixture at T=300K, eq ratio 0.5
+SPECIES_INIT=(-0.575, 2.227e3, 1.447e-2, 0)
+MARKER_INLET_SPECIES = (inlet, -0.575, 2.227e3, 1.447e-2, 0)
+
+INC_INLET_TYPE= VELOCITY_INLET
+MARKER_INLET=(inlet, 300.0, 0.575,1,0,0)
+CFL_REDUCTION_SPECIES= 1.0
+
+MARKER_ISOTHERMAL=(burner_wall, 350, \
+                        cylinder_wall, 400)
+MARKER_SPECIES_STRONG_BC=(burner_wall,cylinder_wall)
+MARKER_INTERNAL=(fluid)
+% ---------------------- REFERENCE VALUE DEFINITION ---------------------------%
+%
+REF_ORIGIN_MOMENT_X = 0.25
+REF_ORIGIN_MOMENT_Y = 0.00
+REF_ORIGIN_MOMENT_Z = 0.00
+REF_LENGTH= 1.0
+REF_AREA= 1.0
+%
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+
+MARKER_SYM= (sides)
+INC_OUTLET_TYPE= PRESSURE_OUTLET
+INC_INLET_DAMPING = 0.05
+INC_OUTLET_DAMPING = 0.05
+MARKER_OUTLET= (outlet, 0.0)
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+%
+NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
+CFL_NUMBER= 150
+CFL_ADAPT= NO
+ITER= 5
+%
+% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
+%
+LINEAR_SOLVER= FGMRES
+LINEAR_SOLVER_PREC= LU_SGS
+LINEAR_SOLVER_ERROR= 1E-04
+LINEAR_SOLVER_ITER= 5
+%
+% -------------------------- MULTIGRID PARAMETERS -----------------------------%
+%
+MGLEVEL= 0
+%
+% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
+%
+CONV_NUM_METHOD_FLOW= FDS
+MUSCL_FLOW= YES
+SLOPE_LIMITER_FLOW = NONE
+TIME_DISCRE_FLOW= EULER_IMPLICIT
+%
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+%
+CONV_RESIDUAL_MINVAL= -15
+CONV_STARTITER= 10
+CONV_CAUCHY_ELEMS= 100
+CONV_CAUCHY_EPS= 1E-6
+SCREEN_OUTPUT = INNER_ITER RMS_VELOCITY-X RMS_PRESSURE RMS_ProgressVariable RMS_EnthalpyTot RMS_MixtureFraction
+HISTORY_OUTPUT = RMS_RES 
+VOLUME_OUTPUT = SOLUTION
+CONV_FIELD = RMS_PRESSURE, RMS_VELOCITY-Z, RMS_VELOCITY-Y, RMS_MixtureFraction
+WRT_ZONE_HIST= YES
+%
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+RESTART_SOL=YES
+OUTPUT_PRECISION= 16
+TABULAR_FORMAT = CSV
+MESH_FILENAME=H2_burner.su2
+OUTPUT_FILES=(RESTART, PARAVIEW_MULTIBLOCK)

incompressible_flow/Inc_Combustion/1__premixed_hydrogen/H2_burner.geo

@@ -0,0 +1,97 @@
+domain_width=0.0008;
+hex_center=0.002;
+hex_radius=0.0003;
+
+//+
+Point(1) = {-0.003, 0, 0, 1.0};
+//+
+Point(2) = {-0.003, domain_width, 0, 1.0};
+//+
+Point(3) = {-0.001, domain_width, 0, 1.0};
+//+
+Point(4) = {-0.001, 0.5*domain_width, 0, 1.0};
+//+
+Point(5) = {-0.00, 0.5*domain_width, 0, 1.0};
+//+
+Point(6) = {-0.00, domain_width, 0, 1.0};
+//+
+Point(7) = {0.005, domain_width, 0, 1.0};
+//+
+Point(8) = {0.005, 0.000, 0, 1.0};
+//+
+Point(9) = {hex_center, 0.000, 0, 1.0};
+//+
+Point(10) = {hex_center+hex_radius, 0.000, 0, 1.0};
+//+
+Point(11) = {hex_center-hex_radius, 0.000, 0, 1.0};
+//+
+Line(1) = {2, 1};
+//+
+Line(2) = {2, 3};
+//+
+Line(3) = {3, 4};
+//+
+Line(4) = {4, 5};
+//+
+Line(5) = {5, 6};
+//+
+Line(6) = {6, 7};
+//+
+Line(7) = {7, 8};
+//+
+Line(8) = {10, 8};
+//+
+Line(9) = {1, 11};
+//+
+Circle(10) = {10, 9, 11};
+//+
+Curve Loop(1) = {1, 9, -10, 8, -7, -6, -5, -4, -3, -2};
+//+
+Plane Surface(1) = {1};
+//+
+//Physical Curve("inlet", 11) = {1};
+//+
+//Physical Curve("outlet", 12) = {7};
+//+
+//Physical Curve("burner_wall", 13) = {3, 4, 5};
+//+
+//Physical Curve("cylinder_wall", 14) = {10};
+//+
+//Physical Curve("sides", 15) = {2, 9, 6, 8};
+//+
+//Physical Surface("fluid", 16) = {1};
+//+
+Transfinite Curve {1, 7} = 30 Using Progression 1;
+//+
+Transfinite Curve {3, 5} = 30 Using Progression 1;
+//+
+Transfinite Curve {4} = 80 Using Progression 1;
+//+
+Transfinite Curve {2} = 100 Using Progression 1;
+//+
+Transfinite Curve {6} = 300 Using Progression 1;
+//+
+Transfinite Curve {8} = 200 Using Progression 1;
+//+
+Transfinite Curve {9} = 350 Using Progression 1;
+//+
+Transfinite Curve {10} = 100 Using Progression 1;
+
+//Mesh 2;
+
+//Save "H2_burner.su2";//+
+Extrude {{1, 0, 0}, {0, -0.001, 0}, Pi/8} {
+  Surface{1}; Layers{10}; Recombine;
+}
+//+
+Physical Surface("inlet", 63) = {25};
+//+
+Physical Surface("outlet", 64) = {41};
+//+
+Physical Surface("sides", 65) = {62, 29, 61, 1, 45, 37};
+//+
+Physical Surface("burner_wall", 66) = {57, 53, 49};
+//+
+Physical Surface("cylinder_wall", 67) = {33};
+//+
+Physical Volume("fluid", 68) = {1};

incompressible_flow/Inc_Combustion/1__premixed_hydrogen/MLP_PD.mlp

@@ -0,0 +1,105 @@
+<header>
+
+[number of layers]
+5
+
+[neurons per layer]
+3
+12
+14
+17
+4
+
+[activation function]
+linear
+gelu
+gelu
+gelu
+linear
+
+[input names]
+ProgressVariable
+EnthalpyTot
+MixtureFraction
+
+[input normalization]
+-1.5064702033996582e+00	+4.6642374992370605e-01
+-3.4190245000000000e+06	+2.1581340000000000e+06
++5.0033628940582275e-03	+1.4977467060089111e-01
+
+[output names]
+Beta_ProgVar
+Beta_Enth
+Beta_Enth_Thermal
+Beta_MixFrac
+
+[output normalization]
+-3.1124658584594727e+00	+5.7480305433273315e-01
+-4.2545765000000000e+06	+5.4751185000000000e+06
+-2.2948374023437500e+03	+1.2485579833984375e+03
+-2.0932110026478767e-02	+3.7706804275512695e-01
+
+</header>
+
+[weights per layer]
+<layer>
++5.5936765670776367e-01	-2.2262737751007080e+00	+8.2052564620971680e-01	-4.1941456794738770e+00	+2.1595144271850586e+00	-3.1018166542053223e+00	+1.5609126091003418e+00	+8.6116069555282593e-01	+1.5900988578796387e+00	-1.7861762046813965e+00	+1.7673474550247192e+00	-8.4477740526199341e-01
++1.7452495098114014e+00	-3.4178844094276428e-01	+3.9943918585777283e-01	-1.3026353716850281e-01	-7.3038059473037720e-01	-2.8568980693817139e+00	-3.0391344428062439e-01	-3.8516503572463989e-01	+1.5720043182373047e+00	-7.9144310951232910e-01	+1.1181222200393677e+00	+1.2395139932632446e+00
+-1.4442356824874878e+00	-2.5616261959075928e+00	+1.8109363317489624e+00	-1.0304316282272339e+00	+1.1644048690795898e+00	-1.2046636343002319e+00	-2.9849336147308350e+00	-2.1314477920532227e+00	-4.4052114486694336e+00	+1.1522486209869385e+00	-5.7982987165451050e-01	-1.7774406075477600e-01
+</layer>
+<layer>
+-7.8683698177337646e-01	+1.2348909378051758e+00	+5.6039118766784668e-01	+6.0650748014450073e-01	+2.8504414558410645e+00	+1.5738261938095093e+00	+5.7276231050491333e-01	+5.8169031143188477e-01	+3.2371595501899719e-01	-7.0329976081848145e-01	-1.1573293209075928e+00	+4.0381139516830444e-01	-1.1950639486312866e+00	-2.4807786941528320e-01
++1.1483484506607056e+00	-1.2595731019973755e-01	+1.2920708954334259e-01	-1.1038926839828491e+00	+4.1350812911987305e+00	-5.4770808219909668e+00	-2.6621893048286438e-01	+9.1883033514022827e-01	+5.1413297653198242e-01	-3.0619984865188599e-01	-1.7892551422119141e+00	-1.1669331789016724e+00	+4.3047985434532166e-01	+1.4167983531951904e+00
+-5.9545308351516724e-01	-1.0836304426193237e+00	-8.0444717407226562e-01	+1.0556594133377075e+00	+3.4371836185455322e+00	+1.1640481948852539e+00	-2.0730528831481934e+00	+1.4779643714427948e-01	-1.5686173439025879e+00	+9.2952281236648560e-01	+4.9133530259132385e-01	-2.6363868713378906e+00	-5.7484567165374756e-01	+3.5405933856964111e-01
+-5.3732520341873169e-01	-8.2899725437164307e-01	+2.1532373130321503e-01	-7.1387672424316406e-01	+8.9794740676879883e+00	+4.7302904129028320e+00	+1.5216027498245239e+00	-3.8849669694900513e-01	-1.4209500513970852e-02	-6.4040400087833405e-02	-4.3968173861503601e-01	+9.7932510375976562e+00	+6.8080127239227295e-01	+6.2318420410156250e-01
++8.5792666673660278e-01	+1.8092916905879974e-01	-3.7652462720870972e-02	+7.6494835317134857e-02	+5.7821393013000488e+00	+1.8568019866943359e+00	-4.4975766539573669e-01	-1.3800190687179565e+00	-5.4718071222305298e-01	+1.3050122559070587e-01	+5.5000835657119751e-01	-1.3087571859359741e+00	+8.2167816162109375e-01	-1.2398092746734619e+00
++1.0533299446105957e+00	+5.1667013168334961e+00	-5.8652287721633911e-01	+1.8700404167175293e+00	+8.8745584487915039e+00	+9.0429029464721680e+00	+4.2462283372879028e-01	-1.6792943701148033e-02	-1.1622345447540283e+00	-2.9977554082870483e-01	-1.5730124711990356e+00	+1.6794612407684326e+00	-1.7160067558288574e+00	+3.5444295406341553e-01
+-6.6769257187843323e-02	+4.7257164120674133e-01	-4.5346045494079590e-01	-1.2471233606338501e+00	-2.9763574600219727e+00	-1.9769278168678284e-01	+3.1241457462310791e+00	+2.9946759343147278e-01	+4.8229956626892090e-01	+4.1772599220275879e+00	+1.1854103207588196e-01	-4.2974877357482910e+00	+1.0468765497207642e+00	+1.0986258983612061e+00
++4.1323339939117432e-01	-1.0169169902801514e+00	+8.7525016069412231e-01	-1.6123955249786377e+00	-1.2063846588134766e+00	+1.4181346893310547e+00	+2.4689452648162842e+00	-1.5278127789497375e-01	+8.0006361007690430e-01	+4.5333325862884521e-01	+1.2695608139038086e+00	-3.2444796562194824e+00	+4.8076021671295166e-01	+5.5755352973937988e-01
++1.7266492545604706e-01	+8.4896576404571533e-01	-7.4060529470443726e-02	+4.1923028230667114e-01	+1.2362111806869507e+00	+4.9337596893310547e+00	+1.9498051404953003e+00	-2.7439159154891968e-01	+1.2809288501739502e-01	+1.2340784817934036e-01	+1.2390412092208862e+00	+5.3083837032318115e-01	-3.0664911866188049e-01	-9.8217773437500000e-01
++7.7484333515167236e-01	+9.7000545263290405e-01	-6.7400932312011719e-01	-6.5156424045562744e-01	-2.4923715591430664e+00	+1.4137476682662964e-01	-6.2123656272888184e-01	-4.3161913752555847e-01	+5.1852053403854370e-01	+7.9075646400451660e-01	-3.9763647317886353e-01	-2.4046406149864197e-01	+9.3279755115509033e-01	-1.2537041306495667e-01
++4.5890122652053833e-01	-5.4653114080429077e-01	+6.7716372013092041e-01	+1.0643934011459351e+00	+2.8232173919677734e+00	+3.4084556102752686e+00	-6.9353675842285156e-01	+2.5891508907079697e-02	-1.7301708459854126e+00	-7.0677620172500610e-01	-7.4925404787063599e-01	-1.5857355594635010e+00	-4.5951095223426819e-01	-1.2846924066543579e+00
+-2.0459318161010742e+00	-1.3967022299766541e-01	-4.8723414540290833e-01	+4.9151417613029480e-01	-2.4336813390254974e-01	+2.6867616176605225e-01	-1.3156063556671143e+00	+1.2456988096237183e+00	-6.5143465995788574e-01	+2.1156504750251770e-01	-7.7381742000579834e-01	+2.5252652168273926e+00	-1.6944633722305298e+00	-3.7577667832374573e-01
+</layer>
+<layer>
+-2.0184963941574097e-01	+1.1482220888137817e+00	-1.4430931210517883e-01	-1.4198860526084900e-01	-2.7880823612213135e-01	-9.3023681640625000e-01	+3.9403364062309265e-01	-1.9451752901077271e+00	-7.9280453920364380e-01	-2.2843687236309052e-01	-6.9571918249130249e-01	-2.8054302092641592e-03	+4.3372470140457153e-01	+1.4820166826248169e+00	-1.3737711906433105e+00	-1.1117722988128662e+00	-3.0378806591033936e+00
+-5.5540651082992554e-01	+6.4094573259353638e-01	+2.4161060154438019e-01	+2.1091027557849884e-01	-9.4434881210327148e-01	-2.1586446762084961e+00	-3.9282519370317459e-02	-1.6424154043197632e+00	-2.6327818632125854e-01	+1.1203126609325409e-01	+4.8172667622566223e-01	+3.7274917960166931e-01	-1.6386254131793976e-01	+5.1941469311714172e-02	-2.7635736465454102e+00	+8.2758963108062744e-01	-2.9272327423095703e+00
++4.5503288507461548e-01	-7.2760522365570068e-01	-1.1357721090316772e+00	-1.9960451126098633e+00	-7.9151540994644165e-01	-6.5633213520050049e-01	-3.6358233541250229e-02	+1.8928986042737961e-02	-2.1545618772506714e-01	-1.0437804460525513e+00	+4.1001714766025543e-02	+3.4175619482994080e-02	-5.9232139587402344e-01	-3.4576010704040527e-01	-7.5590562820434570e-01	+1.2019655108451843e-01	+9.4450451433658600e-02
+-1.8047140538692474e-01	+2.9581779241561890e-01	+3.5067749023437500e-01	-2.8636634349822998e+00	+3.2111722230911255e-01	+8.1838601827621460e-01	+8.1139966845512390e-02	-7.7054816484451294e-01	-1.8204875290393829e-01	-3.8663303852081299e-01	+2.4034704267978668e-01	-5.4463660717010498e-01	+9.2897361516952515e-01	+1.6271905899047852e+00	-3.5175859928131104e-01	-5.0068655014038086e+00	-4.7041401863098145e+00
+-1.9534336403012276e-02	+4.7664877772331238e-01	-6.4567965269088745e-01	-4.2800779342651367e+00	-9.3973153829574585e-01	-4.7965177893638611e-01	-9.2063266783952713e-03	-2.0184078216552734e+00	-7.4120527505874634e-01	+1.8926399946212769e-01	+5.5009335279464722e-02	+3.4929367899894714e-01	-6.3203439116477966e-02	-8.8821917772293091e-01	-1.2712147235870361e+00	+1.1286169290542603e+00	-6.3628115653991699e+00
++9.5517247915267944e-02	-3.8509157299995422e-01	-1.6671176999807358e-02	-9.3428063392639160e-01	+2.0996043086051941e-01	+4.6733841300010681e-01	+6.2699921429157257e-02	-3.0154204368591309e+00	+7.8265428543090820e-02	-3.5962693393230438e-02	-2.6601490378379822e-01	-1.4848266541957855e-01	-6.1616800725460052e-02	-2.3090261220932007e-01	+3.9950132369995117e-01	-2.2599265575408936e+00	-6.1456360816955566e+00
++1.8454858660697937e-01	+6.4392393827438354e-01	-9.2412662506103516e-01	-1.0717321634292603e+00	-7.1740353107452393e-01	+3.0721921920776367e+00	+1.7082008719444275e-01	-1.5314382314682007e+00	-1.7115883827209473e+00	-4.5960146188735962e-01	-9.8167255520820618e-02	-2.8839845210313797e-02	-4.2424803972244263e-01	-2.2923085689544678e+00	+2.3634102940559387e-01	-5.4382983595132828e-02	-2.0812418460845947e+00
++7.2276979684829712e-02	-2.0215651988983154e+00	+3.9073342084884644e-01	+1.9962520599365234e+00	+1.4543913602828979e+00	+1.1602535247802734e+00	-6.9512927532196045e-01	+2.5073411464691162e+00	+7.0365089178085327e-01	-9.3733049929141998e-02	+1.4115136861801147e-01	-5.5563193559646606e-01	+2.7879530191421509e-01	-4.3389448523521423e-01	+8.8280761241912842e-01	+1.4637989997863770e+00	+1.3671472072601318e+00
+-1.5433320403099060e-01	+9.0719491243362427e-01	+1.6209267079830170e-01	-2.1964335441589355e+00	-2.2867610454559326e+00	+1.8663008213043213e+00	-5.4520970582962036e-01	-2.0296618938446045e+00	-3.4272500872612000e-01	+2.1870354190468788e-02	-1.6796220541000366e+00	+2.3241455852985382e-01	-5.4697149991989136e-01	-1.4633115530014038e+00	+2.6578347682952881e+00	-3.6186618804931641e+00	-2.6407463550567627e+00
+-4.0505900979042053e-01	-5.5341029167175293e-01	+1.5711200237274170e-01	-7.8734397888183594e-01	+4.0377050638198853e-01	-1.5427069365978241e-01	-1.6064138710498810e-01	-7.5215941667556763e-01	-9.4907665252685547e-01	+6.0539180040359497e-01	-6.1539363116025925e-02	-2.6970055699348450e-01	+2.5097161531448364e-01	-6.4160132408142090e-01	+3.1192746758460999e-01	-8.1207108497619629e-01	-1.5614933967590332e+00
++3.3514145016670227e-01	-7.0093733072280884e-01	-4.7175073623657227e-01	-7.1444272994995117e-01	-8.0870187282562256e-01	+4.5800599455833435e-01	+4.3255224823951721e-01	+4.2898580431938171e-01	+3.7949404120445251e-01	+9.8860666155815125e-02	-5.5204015970230103e-01	-4.0922561287879944e-01	-2.4125184118747711e-01	-6.4242467284202576e-02	-6.4699250459671021e-01	-1.2401467561721802e+00	-2.5669851303100586e+00
++1.9164021359756589e-03	+3.5010153055191040e-01	+4.9808942712843418e-03	+1.8590971231460571e+00	-7.3865538835525513e-01	+1.2189323902130127e+00	+6.6397942602634430e-02	+1.3344639539718628e+00	+4.5933112502098083e-01	+8.8727630674839020e-02	-2.4055258929729462e-01	-1.2325777858495712e-01	-1.9092649221420288e-01	+1.3201860189437866e+00	+2.1976089477539062e-01	-8.9640945196151733e-01	+1.7600762844085693e+00
++3.4221750497817993e-01	+1.0432732105255127e+00	+5.7717007398605347e-01	-1.3118093013763428e+00	+1.3004060983657837e+00	+3.8324832916259766e-01	+2.1446179598569870e-03	-9.6294808387756348e-01	-6.8076103925704956e-01	+1.1276344954967499e-01	-4.8275017738342285e-01	+5.8240431547164917e-01	-6.1639112234115601e-01	-1.1057349443435669e+00	+1.5214335918426514e+00	+1.7708278894424438e+00	+5.0165188312530518e-01
++3.9782714843750000e-01	-1.4638626575469971e+00	+7.1814291179180145e-02	+2.1027655601501465e+00	-3.7042555809020996e+00	+2.0856988430023193e+00	+9.2289514839649200e-02	-3.8187915086746216e-01	-1.8466535210609436e-01	-4.8500988632440567e-02	+5.8131277561187744e-01	-3.8134214282035828e-01	-1.6914245486259460e-01	-1.9203429222106934e+00	-2.1438066959381104e+00	+2.3994569778442383e+00	-1.3890436887741089e+00
+</layer>
+<layer>
+-6.6255033016204834e-02	-4.0674057602882385e-01	-5.9483653306961060e-01	-2.9391551017761230e-01
+-1.6143990680575371e-02	-2.4502272903919220e-01	+6.6644951701164246e-02	-2.9729947447776794e-02
+-1.8170315027236938e-01	-2.1033358573913574e-01	+1.9836166501045227e-01	-1.0021093301475048e-02
++1.1217525228857994e-02	-6.6836856305599213e-02	-5.1224365830421448e-02	+3.5130698233842850e-02
+-1.5285339206457138e-02	+3.1436040997505188e-02	-9.9899478256702423e-02	+4.2289614677429199e-02
+-6.9087431766092777e-03	+6.0401394963264465e-02	+1.7236797139048576e-03	-5.5530946701765060e-03
++4.1897991299629211e-01	-1.9636897742748260e-01	+1.0367606580257416e-01	+1.3625358045101166e-01
++1.6224026679992676e-01	-1.1469457149505615e+00	+1.9311760365962982e-01	-5.8333348482847214e-02
++1.3593722879886627e-01	-1.4048127830028534e-01	+5.7687763124704361e-02	-1.0419285297393799e-01
+-4.3920856714248657e-01	+1.2700842320919037e-01	-6.2523370981216431e-01	+6.3199234008789062e-01
+-8.8329017162322998e-02	+3.8248756527900696e-01	+1.1002863943576813e-01	-2.2133979946374893e-02
++2.0028780400753021e-01	+1.6378790140151978e-01	-9.7495861351490021e-02	+2.7645299211144447e-02
++5.3663062863051891e-03	+6.3815879821777344e-01	-2.5993975996971130e-01	+3.1326732039451599e-01
+-3.6526303738355637e-02	+2.8247556090354919e-01	-6.5268911421298981e-02	+3.3453144133090973e-02
++1.3270338065922260e-02	-5.9688974171876907e-02	-1.2011376023292542e-01	-3.8207445293664932e-02
++1.1055923998355865e-02	-8.4312260150909424e-02	+8.2623641937971115e-03	+2.2568788845092058e-03
++7.1002934128046036e-03	-4.9927394837141037e-02	+5.6600286625325680e-03	+5.8617728063836694e-04
+</layer>
+
+[biases per layer]
++0.0000000000000000e+00	+0.0000000000000000e+00	+0.0000000000000000e+00
+-1.0984963178634644e+00	+1.8675320148468018e+00	-1.6051063537597656e+00	+1.1760630607604980e+00	-1.6234400272369385e+00	+2.4758839607238770e+00	-1.1829069852828979e+00	+1.4792251586914062e-01	-1.4212688207626343e+00	+4.1564255952835083e-01	-1.4366428852081299e+00	-2.8742429614067078e-01
++1.0597276501357555e-02	+5.7052063941955566e-01	-4.7234252095222473e-02	-2.9167464375495911e-01	-1.0065306425094604e+00	-1.3775391578674316e+00	-5.6155574321746826e-01	-2.9035386443138123e-01	-1.5129996836185455e-01	-1.5230740420520306e-02	-5.3327292203903198e-01	-3.0793693661689758e-01	-2.5004431605339050e-01	-2.4458266794681549e-02
+-7.0153272151947021e-01	-2.6219487190246582e-01	-1.1805304884910583e-01	-3.3423715829849243e-01	-7.4985355138778687e-01	+4.2404353618621826e-01	+6.7827850580215454e-02	-2.0516588687896729e+00	-5.0613462924957275e-01	+4.0178808569908142e-01	-2.8277575969696045e-01	-3.0988002195954323e-02	-1.4533357322216034e-01	-1.4993039369583130e+00	-8.7631237506866455e-01	-3.3392030000686646e-01	-1.2971979379653931e+00
++6.8486714363098145e-01	+3.1399080157279968e-01	+7.6813256740570068e-01	+1.0490987449884415e-01

incompressible_flow/Inc_Combustion/1__premixed_hydrogen/MLP_null.mlp

@@ -0,0 +1,48 @@
+<header>
+
+[number of layers]
+3
+
+[neurons per layer]
+3
+2
+1
+
+[activation function]
+linear
+linear
+linear
+
+[input names]
+ProgressVariable
+EnthalpyTot
+MixtureFraction
+
+[input normalization]
+-6.8059708053507162e-01	+4.5623852432084366e-01
+-4.9308262569627967e+06	+8.6731375409855954e+06
++2.5577753524576785e-03	+3.6840647718850406e-01
+
+[output names]
+NULL
+
+[output normalization]
+0.0 1.0
+
+</header>
+
+[weights per layer]
+<layer>
+0.0 0.0
+0.0 0.0
+0.0 0.0
+</layer>
+<layer>
+0.0
+0.0
+</layer>
+
+[biases per layer]
+0.0 0.0 0.0
+0.0 0.0
+0.0

incompressible_flow/Inc_Species_Transport/3__gradient-validation/postprocess.py

@@ -49,7 +49,7 @@ def printGradVal(FDgrad, DAgrad):
   nDV = DAgrad_specVar.size
 
   # Load primal values and create FD gradient
-  FDvals = pd.read_csv("doe.his")
+  FDvals = pd.read_csv("doe.csv")
 
   FDstring_specVar = '  specVar'
 

multiphysics/contact_resistance_cht/fluid.cfg

@@ -0,0 +1,55 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer  %
+% with contact resistance.                                                     %
+% Author: E.C.Bunschoten                                                       %
+% Date: August 6, 2024                                                         %
+% File Version 8.0.1 "Harrier"                                                 %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+
+SOLVER= INC_NAVIER_STOKES
+INC_ENERGY_EQUATION= YES
+
+% ------------------- INCOMPRESSIBLE FREE-STREAM DEFINITION -------------------%
+
+FREESTREAM_TEMPERATURE= 300
+
+INC_VELOCITY_INIT=(0.1, 0.0, 0.0)
+FLUID_MODEL=INC_IDEAL_GAS
+
+VISCOSITY_MODEL=CONSTANT_VISCOSITY
+MU_CONSTANT=1e-5
+
+REF_DIMENSIONALIZATION= DIMENSIONAL
+INC_NONDIM= DIMENSIONAL
+
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+INC_INLET_TYPE=VELOCITY_INLET
+INC_OUTLET_TYPE=PRESSURE_OUTLET
+MARKER_SYM= ( side_3 )
+MARKER_INLET= ( inlet, 300, 0.1, 1.0, 0.0, 0.0 )
+MARKER_OUTLET= ( outlet, 0.0)
+MARKER_INTERNAL= (fluid)
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+
+CFL_NUMBER= 100
+
+% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
+%
+CONV_NUM_METHOD_FLOW= FDS
+MUSCL_FLOW= YES
+
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+
+CONV_RESIDUAL_MINVAL= -12
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+MESH_FILENAME= fluid_mesh.su2
+CONV_FILENAME= history_fluid_3
+TABULAR_FORMAT= CSV

multiphysics/contact_resistance_cht/fluid_3.geo

@@ -0,0 +1,60 @@
+//+
+Point(1) = {0, 0, 0, 1.0};
+//+
+Point(2) = {-0.01, 0, 0, 1.0};
+//+
+Point(3) = {0.01, 0, 0, 1.0};
+//+
+Point(4) = {0.01, -0.005, 0, 1.0};
+//+
+Point(5) = {0.0, -0.005, 0, 1.0};
+//+
+Point(6) = {-0.01, -0.005, 0, 1.0};
+//+
+Line(1) = {6, 2};
+//+
+Line(2) = {5, 1};
+//+
+Line(3) = {4, 3};
+//+
+Line(4) = {2, 1};
+//+
+Line(5) = {1, 3};
+//+
+Line(6) = {6, 5};
+//+
+Line(7) = {5, 4};
+//+
+Curve Loop(1) = {1, 4, -2, -6};
+//+
+Plane Surface(1) = {1};
+//+
+Curve Loop(2) = {5, -3, -7, 2};
+//+
+Plane Surface(2) = {2};
+//+
+Physical Curve("inlet", 8) = {1};
+//+
+Physical Curve("outlet", 9) = {3};
+//+
+Physical Curve("side_3", 10) = {6, 7};
+//+
+Physical Curve("cht_interface_3_1", 11) = {4};
+//+
+Physical Curve("cht_interface_3_2", 12) = {5};
+//+
+Physical Surface("fluid", 13) = {1, 2};
+//+
+Transfinite Curve {1, 2, 3} = 20 Using Progression 0.8;
+//+
+Transfinite Curve {4, 6, 7, 5} = 80 Using Progression 1;
+//+
+Transfinite Surface {1} = {6, 5, 1, 2};
+//+
+Transfinite Surface {2} = {5, 4, 3, 1};
+//+
+Recombine Surface {1, 2};
+
+Mesh 2;
+
+Save "fluid_mesh.su2";

multiphysics/contact_resistance_cht/master.cfg

@@ -0,0 +1,32 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer  %
+% with contact resistance.                                                     %
+% Author: E.C.Bunschoten                                                       %
+% Date: August 6, 2024                                                         %
+% File Version 8.0.1 "Harrier"                                                 %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= MULTIPHYSICS
+MATH_PROBLEM= DIRECT
+MULTIZONE_MESH=NO
+CONFIG_LIST = (solid_1.cfg, solid_2.cfg, fluid.cfg)
+MARKER_ZONE_INTERFACE= (cht_interface_1_2, cht_interface_2_1, \
+                        cht_interface_1_3, cht_interface_3_1, \
+                        cht_interface_2_3, cht_interface_3_2)
+MARKER_CHT_INTERFACE= (cht_interface_1_2, cht_interface_2_1, \
+                        cht_interface_1_3, cht_interface_3_1, \
+                        cht_interface_2_3, cht_interface_3_2)
+CHT_INTERFACE_CONTACT_RESISTANCE = (1e-5, 0, 0)
+
+CHT_COUPLING_METHOD= DIRECT_TEMPERATURE_ROBIN_HEATFLUX
+TIME_DOMAIN = NO
+OUTER_ITER = 2000
+WRT_ZONE_HIST=YES
+
+OUTPUT_FILES=(RESTART, PARAVIEW_MULTIBLOCK)
+SCREEN_OUTPUT= OUTER_ITER, RMS_RES[0], RMS_RES[1], RMS_RES[2]

multiphysics/contact_resistance_cht/solid_1.cfg

@@ -0,0 +1,47 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer  %
+% with contact resistance.                                                     %
+% Author: E.C.Bunschoten                                                       %
+% Date: August 6, 2024                                                         %
+% File Version 8.0.1 "Harrier"                                                 %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= HEAT_EQUATION
+
+
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_ISOTHERMAL= ( isothermal_wall_1, 500.0 )
+MARKER_SYM=(side_1)
+MARKER_MONITORING= ( NONE )
+MARKER_INTERNAL=(solid_1)
+
+% ---------------- (SOLIDS) CONDUCTION CONDITION DEFINITION -------------------%
+
+% We should keep the dimensionalization of the coupled flow solver
+INC_NONDIM= DIMENSIONAL
+FREESTREAM_TEMPERATURE= 300.0
+
+% Properties of stainless steel
+MATERIAL_DENSITY= 8935
+SPECIFIC_HEAT_CP= 3850
+THERMAL_CONDUCTIVITY_CONSTANT=26
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+
+CFL_NUMBER= 100.0
+
+
+% -------------------- HEAT NUMERICAL METHOD DEFINITION -----------------------%
+
+TIME_DISCRE_HEAT= EULER_IMPLICIT
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+TABULAR_FORMAT= CSV
+CONV_FILENAME= history_solid_1
+MESH_FILENAME=solid_mesh_1.su2

multiphysics/contact_resistance_cht/solid_1.geo

@@ -0,0 +1,42 @@
+//+
+Point(1) = {0, 0, 0, 1.0};
+//+
+Point(2) = {-0.01, 0, 0, 1.0};
+//+
+Point(3) = {-0.01, 0.01, 0, 1.0};
+//+
+Point(4) = {0, 0.01, 0, 1.0};
+//+
+Line(1) = {2, 3};
+//+
+Line(2) = {3, 4};
+//+
+Line(3) = {4, 1};
+//+
+Line(4) = {1, 2};
+//+
+Curve Loop(1) = {1, 2, 3, 4};
+//+
+Plane Surface(1) = {1};
+//+
+Physical Curve("isothermal_wall_1", 5) = {1};
+//+
+Physical Curve("cht_interface_1_2", 6) = {3};
+//+
+Physical Curve("side_1", 7) = {2};
+//+
+Physical Curve("cht_interface_1_3", 8) = {4};
+//+
+Physical Surface("solid_1", 9) = {1};
+//+
+Transfinite Curve {1, 2, 3, 4} = 40 Using Progression 1;
+//+
+Transfinite Surface {1} = {2, 1, 4, 3};
+//+
+Recombine Surface {1};
+//+
+Recombine Surface {1};
+
+Mesh 2;
+
+Save "solid_mesh_1.su2";

multiphysics/contact_resistance_cht/solid_2.cfg

@@ -0,0 +1,41 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: solid-to-solid and solid-to-fluid conjugate heat transfer  %
+% with contact resistance.                                                     %
+% Author: E.C.Bunschoten                                                       %
+% Date: August 6, 2024                                                         %
+% File Version 8.0.1 "Harrier"                                                 %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= HEAT_EQUATION
+
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_ISOTHERMAL= ( isothermal_wall_2, 300.0 )
+MARKER_SYM=(side_2)
+MARKER_INTERNAL=(solid_2)
+
+% ---------------- (SOLIDS) CONDUCTION CONDITION DEFINITION -------------------%
+
+% We should keep the dimensionalization of the coupled flow solver
+INC_NONDIM= DIMENSIONAL
+FREESTREAM_TEMPERATURE= 300.0
+
+MATERIAL_DENSITY= 8000
+SPECIFIC_HEAT_CP= 4420
+THERMAL_CONDUCTIVITY_CONSTANT= 61
+
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+
+CFL_NUMBER= 100.0
+
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+
+MESH_FILENAME= solid_mesh_2.su2
+TABULAR_FORMAT= CSV
+CONV_FILENAME= history_solid_2
+

multiphysics/contact_resistance_cht/solid_2.geo

@@ -0,0 +1,40 @@
+//+
+Point(1) = {0, 0, 0, 1.0};
+//+
+Point(2) = {0.01, 0, 0, 1.0};
+//+
+Point(3) = {0.01, 0.01, 0, 1.0};
+//+
+Point(4) = {0.0, 0.01, 0, 1.0};
+//+
+Line(1) = {1, 4};
+//+
+Line(2) = {4, 3};
+//+
+Line(3) = {3, 2};
+//+
+Line(4) = {2, 1};
+//+
+Curve Loop(1) = {1, 2, 3, 4};
+//+
+Surface(1) = {1};
+//+
+Physical Curve("cht_interface_2_1", 5) = {1};
+//+
+Physical Curve("cht_interface_2_3", 6) = {4};
+//+
+Physical Curve("isothermal_wall_2", 7) = {3};
+//+
+Physical Curve("side_2", 8) = {2};
+//+
+Physical Surface("solid_2", 9) = {1};
+//+
+Transfinite Curve {1, 4, 3, 2} = 40 Using Progression 1;
+//+
+Transfinite Surface {1} = {1, 2, 3, 4};
+//+
+Recombine Surface {1};
+
+Mesh 2;
+
+Save "solid_mesh_2.su2";