Commit 81645e33 authored by Lambert Theisen's avatar Lambert Theisen 🔥

Add 3 more GLS tests for R13

- Change tau-params by intuition ¯\_(ツ)_/¯
parent f10ff94b
Pipeline #320866 passed with stages
in 28 minutes and 51 seconds
......@@ -59,11 +59,11 @@ stabilization:
delta_p: 0.01
gls:
enable: True
tau_energy: 0.001
tau_heatflux: 0.001
tau_mass: 0.01
tau_momentum: 0.01
tau_stress: 0.01
tau_energy: 0.0001
tau_heatflux: 0.0001
tau_mass: 0.005
tau_momentum: 0.005
tau_stress: 0.005
# Formulation Parameters
# ======================
......
# General
# =======
# - output_folder: Used as output folder
output_folder: r13_1_coeffs_nosources_norot_inflow_p2p2p2p2p2_gls
# Meshes
# ======
# - meshes: List of input meshes in h5 format to run simulations on
meshes:
- ../mesh/ring0.h5
- ../mesh/ring1.h5
- ../mesh/ring2.h5
- ../mesh/ring3.h5
# - ../mesh/ring4.h5
# - ../mesh/ring5.h5
# Numerical Parameters
# ====================
# - elements: Must contain the fields: theta, s, p, u, sigma
# - fields: List of FEM parameters (shape, degree)
# - shape: Element shape, e.g. Lagrange
# - degree: Element degree, e.g. 2
# - stabilization: Must contain cip and gls
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
# - gls: Collection of Garlerkin Least Squares (GLS) parameters
# - enable: Enable GLS stabilization
# - tau_energy: Stabilization with energy eq. residual
# - tau_heatflux: Stabilization with heatflux eq. residual
# - tau_mass: Stabilization with mass eq. residual
# - tau_momentum: Stabilization with momentum eq. residual
# - tau_stress: Stabilization with stress eq. residual
elements:
theta:
shape: Lagrange
degree: 2
s:
shape: Lagrange
degree: 2
p:
shape: Lagrange
degree: 2
u:
shape: Lagrange
degree: 2
sigma:
shape: Lagrange
degree: 2
stabilization:
cip:
enable: False
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
gls:
enable: True
tau_energy: 0.0001
tau_heatflux: 0.0001
tau_mass: 0.005
tau_momentum: 0.005
tau_stress: 0.005
# Formulation Parameters
# ======================
# - nsd: Number of spatial dimensions == 2
# - mode: Formulation mode, one of heat, stress, r13
# - heat_source: Heat source function for mode==heat||r13
# - mass_source: Mass source function for mode==stress||r13
# - body_force: Body force for mode==stress||r13
nsd: 2
mode: r13
heat_source: 0
mass_source: 0
body_force: [0,0]
# Region Parameters
# =================
# - regs: Dictionary of all mesh regions
# - reg_id: Must contain the following parameters:
# - kn: Knudsen number
regs:
4000:
kn: 1.0
# Boundary Conditions
# ===================
# - bcs: Dictionary of all boundary IDs from mesh
# - bc_id: must contain the following parameters
# - chi_tilde: Refaction coefficient in Maxwell accomodation model
# - theta_w: Value for temperature at wall
# - u_t_w: Value for tangential velocity at wall
# - u_n_w: Value for normal velocity at wall
# - p_w: Value for pressure at wall
# - epsilon_w: Inflow-model parameter <=> Weight of pressure prescription
bcs:
3000:
chi_tilde: 1.0
theta_w: 1
u_t_w: 0
u_n_w: 0
p_w: 0
epsilon_w: pow(10,-3)
3100:
chi_tilde: 1.0
theta_w: 2
u_t_w: -1.00 * sin(phi)
u_n_w: +1.00 * cos(phi)
p_w: -0.27 * cos(phi)
epsilon_w: pow(10,+3)
# Convergence Study
# =================
# - enable: Enable convergence study on given meshes
# - exact_solution: Path to exact solution in cpp-format to compare errors
# - plot: Show errors in matplotlib window. PDF output is always per default.
# - write_systemmatrix: Writes out systemmatrix (LHS) to use for analysis
# - rescale_pressure: Rescale numerical pressure result to have zero mean
# - relative_errors: Use relative errors. If exact sol. is zero, use absolute.
convergence_study:
enable: True
exact_solution: esols/1_coeffs_nosources_norot_inflow_positive.cpp
plot: False
write_systemmatrix: False
rescale_pressure: False
relative_error: True
# Postprocessing
# ==============
# - write_pdfs: Write all solution fields as PDF plot
# - write_vecs: Write all solution fields as vectors
# - massflow: List of BC IDs to compute massflow J=int_bc dot(u,n) ds
postprocessing:
write_pdfs: False
write_vecs: False
massflow: []
# Parameter Study
# ==============
# - enable: Repeat simulation with different p. values (study)
# - parameter_key: Key as list, e.g. ["elemenets", "p", "degree"]
# - parameter_values: List of value for parameter, e.g. [0.01,0.1,1,10]
parameter_study:
enable: False
parameter_key: []
parameter_values: []
# General
# =======
# - output_folder: Used as output folder
output_folder: r13_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_gls
# Meshes
# ======
# - meshes: List of input meshes in h5 format to run simulations on
meshes:
- ../mesh/ring0.h5
- ../mesh/ring1.h5
- ../mesh/ring2.h5
- ../mesh/ring3.h5
- ../mesh/ring4.h5
# - ../mesh/ring5.h5
# - ../mesh/ring6.h5
# - ../mesh/ring7.h5
# Numerical Parameters
# ====================
# - elements: Must contain the fields: theta, s, p, u, sigma
# - fields: List of FEM parameters (shape, degree)
# - shape: Element shape, e.g. Lagrange
# - degree: Element degree, e.g. 2
# - stabilization: Must contain cip and gls
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
# - gls: Collection of Garlerkin Least Squares (GLS) parameters
# - enable: Enable GLS stabilization
# - tau_energy: Stabilization with energy eq. residual
# - tau_heatflux: Stabilization with heatflux eq. residual
# - tau_mass: Stabilization with mass eq. residual
# - tau_momentum: Stabilization with momentum eq. residual
# - tau_stress: Stabilization with stress eq. residual
elements:
theta:
shape: Lagrange
degree: 1
s:
shape: Lagrange
degree: 1
p:
shape: Lagrange
degree: 1
u:
shape: Lagrange
degree: 1
sigma:
shape: Lagrange
degree: 1
stabilization:
cip:
enable: False
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
gls:
enable: True
tau_energy: 0.0001
tau_heatflux: 0.0001
tau_mass: 0.005
tau_momentum: 0.005
tau_stress: 0.005
# Formulation Parameters
# ======================
# - nsd: Number of spatial dimensions == 2
# - mode: Formulation mode, one of heat, stress, r13
# - heat_source: Heat source function for mode==heat||r13
# - mass_source: Mass source function for mode==stress||r13
# - body_force: Body force for mode==stress||r13
nsd: 2
mode: r13
heat_source: 1.0 * (1.0 - (5.0*pow(R,2))/(18.0*pow(1.0,2))) * cos(phi)
mass_source: 1.0 * (1.0 - (5.0*pow(R,2))/(18.0*pow(1.0,2))) * cos(phi)
body_force: [0,0]
# Region Parameters
# =================
# - regs: Dictionary of all mesh regions
# - reg_id: Must contain the following parameters:
# - kn: Knudsen number
regs:
4000:
kn: 1.0
# Boundary Conditions
# ===================
# - bcs: Dictionary of all boundary IDs from mesh
# - bc_id: must contain the following parameters
# - chi_tilde: Refaction coefficient in Maxwell accomodation model
# - theta_w: Value for temperature at wall
# - u_t_w: Value for tangential velocity at wall
# - u_n_w: Value for normal velocity at wall
# - p_w: Value for pressure at wall
# - epsilon_w: Inflow-model parameter <=> Weight of pressure prescription
bcs:
3000:
chi_tilde: 1.0
theta_w: 1.0
u_t_w: -10
u_n_w: 0
p_w: 0
epsilon_w: 0
3100:
chi_tilde: 1.0
theta_w: 0.5
u_t_w: 0
u_n_w: 0
p_w: 0
epsilon_w: 0
# Convergence Study
# =================
# - enable: Enable convergence study on given meshes
# - exact_solution: Path to exact solution in cpp-format to compare errors
# - plot: Show errors in matplotlib window. PDF output is always per default.
# - write_systemmatrix: Writes out systemmatrix (LHS) to use for analysis
# - rescale_pressure: Rescale numerical pressure result to have zero mean
# - relative_errors: Use relative errors. If exact sol. is zero, use absolute.
convergence_study:
enable: True
exact_solution: esols/1_coeffs_sources_rot_noinflow.cpp
plot: False
write_systemmatrix: False
rescale_pressure: True
relative_error: True
# Postprocessing
# ==============
# - write_pdfs: Write all solution fields as PDF plot
# - write_vecs: Write all solution fields as vectors
# - massflow: List of BC IDs to compute massflow J=int_bc dot(u,n) ds
postprocessing:
write_pdfs: False
write_vecs: False
massflow: []
# Parameter Study
# ==============
# - enable: Repeat simulation with different p. values (study)
# - parameter_key: Key as list, e.g. ["elemenets", "p", "degree"]
# - parameter_values: List of value for parameter, e.g. [0.01,0.1,1,10]
parameter_study:
enable: False
parameter_key: []
parameter_values: []
# General
# =======
# - output_folder: Used as output folder
output_folder: r13_1_coeffs_sources_rot_noinflow_p2p2p2p2p2_gls
# Meshes
# ======
# - meshes: List of input meshes in h5 format to run simulations on
meshes:
- ../mesh/ring0.h5
- ../mesh/ring1.h5
- ../mesh/ring2.h5
- ../mesh/ring3.h5
# - ../mesh/ring4.h5
# - ../mesh/ring5.h5
# Numerical Parameters
# ====================
# - elements: Must contain the fields: theta, s, p, u, sigma
# - fields: List of FEM parameters (shape, degree)
# - shape: Element shape, e.g. Lagrange
# - degree: Element degree, e.g. 2
# - stabilization: Must contain cip and gls
# - cip: Collection of Continous Interior Penalty (CIP) parameters
# - enable: Enable CIP stabilization
# - delta_theta: Stabilization of grad(T)*grad(T_test) over edge
# - delta_u: Stabilization of grad(u)*grad(u_test) over edge
# - delta_p: Stabilization of grad(p)*grad(p_test) over edge
# - gls: Collection of Garlerkin Least Squares (GLS) parameters
# - enable: Enable GLS stabilization
# - tau_energy: Stabilization with energy eq. residual
# - tau_heatflux: Stabilization with heatflux eq. residual
# - tau_mass: Stabilization with mass eq. residual
# - tau_momentum: Stabilization with momentum eq. residual
# - tau_stress: Stabilization with stress eq. residual
elements:
theta:
shape: Lagrange
degree: 2
s:
shape: Lagrange
degree: 2
p:
shape: Lagrange
degree: 2
u:
shape: Lagrange
degree: 2
sigma:
shape: Lagrange
degree: 2
stabilization:
cip:
enable: False
delta_theta: 1.0
delta_u: 1.0
delta_p: 0.01
gls:
enable: True
tau_energy: 0.0001
tau_heatflux: 0.0001
tau_mass: 0.005
tau_momentum: 0.005
tau_stress: 0.005
# Formulation Parameters
# ======================
# - nsd: Number of spatial dimensions == 2
# - mode: Formulation mode, one of heat, stress, r13
# - heat_source: Heat source function for mode==heat||r13
# - mass_source: Mass source function for mode==stress||r13
# - body_force: Body force for mode==stress||r13
nsd: 2
mode: r13
heat_source: 1.0 * (1.0 - (5.0*pow(R,2))/(18.0*pow(1.0,2))) * cos(phi)
mass_source: 1.0 * (1.0 - (5.0*pow(R,2))/(18.0*pow(1.0,2))) * cos(phi)
body_force: [0,0]
# Region Parameters
# =================
# - regs: Dictionary of all mesh regions
# - reg_id: Must contain the following parameters:
# - kn: Knudsen number
regs:
4000:
kn: 1.0
# Boundary Conditions
# ===================
# - bcs: Dictionary of all boundary IDs from mesh
# - bc_id: must contain the following parameters
# - chi_tilde: Refaction coefficient in Maxwell accomodation model
# - theta_w: Value for temperature at wall
# - u_t_w: Value for tangential velocity at wall
# - u_n_w: Value for normal velocity at wall
# - p_w: Value for pressure at wall
# - epsilon_w: Inflow-model parameter <=> Weight of pressure prescription
bcs:
3000:
chi_tilde: 1.0
theta_w: 1.0
u_t_w: -10
u_n_w: 0
p_w: 0
epsilon_w: 0
3100:
chi_tilde: 1.0
theta_w: 0.5
u_t_w: 0
u_n_w: 0
p_w: 0
epsilon_w: 0
# Convergence Study
# =================
# - enable: Enable convergence study on given meshes
# - exact_solution: Path to exact solution in cpp-format to compare errors
# - plot: Show errors in matplotlib window. PDF output is always per default.
# - write_systemmatrix: Writes out systemmatrix (LHS) to use for analysis
# - rescale_pressure: Rescale numerical pressure result to have zero mean
# - relative_errors: Use relative errors. If exact sol. is zero, use absolute.
convergence_study:
enable: True
exact_solution: esols/1_coeffs_sources_rot_noinflow.cpp
plot: False
write_systemmatrix: False
rescale_pressure: True
relative_error: True
# Postprocessing
# ==============
# - write_pdfs: Write all solution fields as PDF plot
# - write_vecs: Write all solution fields as vectors
# - massflow: List of BC IDs to compute massflow J=int_bc dot(u,n) ds
postprocessing:
write_pdfs: False
write_vecs: False
massflow: []
# Parameter Study
# ==============
# - enable: Repeat simulation with different p. values (study)
# - parameter_key: Key as list, e.g. ["elemenets", "p", "degree"]
# - parameter_values: List of value for parameter, e.g. [0.01,0.1,1,10]
parameter_study:
enable: False
parameter_key: []
parameter_values: []
h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
0.9886573325052778,0.180605313071,0.261114230569,0.377807012621,0.306646244746,0.413002680559,0.317306187036,0.316168943416,0.470672089665,1.91599181943,1.06969671752,0.648739970729,0.698682742711,0.808664933249,0.851958111861,1.09942707914,0.595159375448,1.26705484475,1.13814783912
0.6340332990709842,0.15528892383,0.222440558532,0.249169406581,0.220462494934,0.263565998481,0.21145098483,0.197515549887,0.301427043379,0.796882046201,0.427583880196,0.279701836461,0.284743047308,0.52333139814,0.495456150171,0.647837181536,0.406716508026,0.76698418419,0.627667942693
0.32904683851469807,0.0609243854636,0.169733939726,0.0840591868686,0.077815069942,0.0946018825755,0.104775999231,0.0713854110081,0.104789584129,0.326036736288,0.196051991843,0.112227133269,0.176450314621,0.206512959909,0.174139380121,0.286476100683,0.191625302388,0.283413370874,0.215322914926
0.16754966839339377,0.0209231699323,0.0753524519726,0.026555408525,0.032704825765,0.0250525678851,0.0247859112886,0.0223401437921,0.0307390323959,0.0985192883511,0.0903003117193,0.0347113677067,0.0829711137882,0.066948805771,0.0550913992393,0.093050002277,0.0655105695883,0.0929858751815,0.066300113921
0.08734460120995041,0.00498118391959,0.0256461854468,0.00651868965854,0.00782373825314,0.00627629744909,0.00687484128807,0.00650793370863,0.00974935692896,0.0248769021146,0.0209030082005,0.00964008648834,0.0263561301934,0.0199636758501,0.0175289574347,0.0278166117552,0.0177819509853,0.0279462682097,0.0198214469802
0.9886573325052778,0.631895483532,0.794571098104,0.451069153247,0.351489674899,0.466628176439,0.330194326007,0.262671692186,0.349632690359,2.43094246603,1.39047918469,0.880943606335,0.874899308216,0.820157053283,1.06385466202,1.10024053576,0.715613598921,1.2509096967,1.39965048044
0.6340332990709842,0.714885662006,0.721987497921,0.302003110011,0.259397079019,0.341485194644,0.260134600446,0.161270421776,0.234238111224,0.959576427235,0.511507575345,0.363865582412,0.458968879292,0.507841788485,0.580979016578,0.619500094695,0.470491201045,0.724067505631,0.730183108398
0.32904683851469807,0.197307106092,0.566515789674,0.0925165262998,0.0873639080016,0.109586527909,0.134364488311,0.061047713412,0.0850629104858,0.36731893409,0.280448343774,0.143930350041,0.352619896467,0.188519245648,0.201038433501,0.258641042771,0.217405558061,0.25287811288,0.24502466752
0.16754966839339377,0.074215516999,0.259026674219,0.0283802867764,0.0390318787274,0.0277000207399,0.0292397632362,0.0184570684518,0.0227790911019,0.110842719084,0.203971304836,0.0439866065842,0.16209837418,0.0583617807846,0.0579320652124,0.0799738707375,0.0706208463938,0.0792644158495,0.0731364769453
0.08734460120995041,0.0181163400216,0.103214165496,0.00625251871827,0.00838231680321,0.00644588608632,0.00865685122678,0.00538112391082,0.00756020612376,0.0278972146457,0.0652490761938,0.0122141330829,0.0694265880961,0.0171146534153,0.0185933368402,0.0235807741295,0.0190358024455,0.0234288345416,0.0207604598285
h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
0.9886573325052778,0.0419312837545,0.0328328120826,0.144674503692,0.162306245642,0.173442348824,0.172719508294,0.0731497530603,0.150480400412,0.541140373362,0.208956056846,0.157744020403,0.076390812032,0.195429308309,0.314705261904,0.290917731238,0.351751889486,0.270570071284,0.346086653859
0.6340332990709842,0.0246216101651,0.0197137113697,0.0910616245238,0.119936710874,0.0952751487421,0.107426913214,0.0739393376607,0.0580249908928,0.417810373533,0.168653709284,0.115306837484,0.0837804931443,0.265828714087,0.289982594106,0.318189452031,0.269360569977,0.372539966099,0.291245615572
0.32904683851469807,0.00929031609849,0.0114809024355,0.0394833645198,0.0476762211902,0.0342570544196,0.0358694404627,0.0553673903251,0.0302305931084,0.258403272018,0.121325386931,0.0813115383241,0.0754129069435,0.177808871754,0.120879107348,0.243761231226,0.150719655144,0.239372251144,0.135167113868
0.16754966839339377,0.00198949570048,0.00515068947524,0.0138193905016,0.0156576114023,0.0103379709474,0.010099776496,0.0217717124283,0.0152610019103,0.0905317461933,0.0472297664167,0.0303374203608,0.0385043065158,0.0691324781456,0.0454774244466,0.0945025216645,0.0468613622429,0.093836917621,0.0505751562826
h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
0.9886573325052778,1.04667739136,1.3236931735,0.509493758786,0.440378031329,0.407766022771,0.353495081217,0.17841265789,0.211576746968,0.998444665591,1.27174732292,1.74439722193,2.43261153283,0.443531575782,0.48923085346,0.382442869467,0.424163815361,0.431619708192,0.465421364149
0.6340332990709842,1.17247979977,1.2049688644,0.308871241364,0.286664040443,0.3356685019,0.298736187205,0.122402683749,0.123151497476,0.771443446822,1.44315302905,1.45452156932,2.67244589146,0.272985755916,0.299210676355,0.265426580794,0.29643900693,0.261794630315,0.28192322485
0.32904683851469807,0.320521575238,0.899902601665,0.0954179057242,0.0963722315889,0.103980939483,0.118547287111,0.101822034712,0.26311446658,0.479049357393,1.42983069371,0.833215676678,2.88898079926,0.102783603451,0.136607909466,0.0983426724122,0.147874365338,0.0991345921008,0.127054621869
0.16754966839339377,0.118890316498,0.402982707095,0.0283421386184,0.0429023926559,0.0263630925327,0.0438763335455,0.0334531542975,0.144480241714,0.195101918326,1.09695379458,0.332436914629,1.53302270993,0.0318678923591,0.0510184630634,0.029732706955,0.049046675786,0.0302403689227,0.046423241413
0.08734460120995041,0.029803563456,0.190328461345,0.00599268819319,0.00919521794433,0.006179679532,0.00985574000632,0.00595819875196,0.0262777120116,0.0573190600689,0.454953195708,0.104386747054,0.701371474951,0.00748010876933,0.0144619139351,0.00707557660666,0.0126964297618,0.00716747263195,0.0126276780572
h,theta_L_2,theta_l_inf,sx_L_2,sx_l_inf,sy_L_2,sy_l_inf,p_L_2,p_l_inf,ux_L_2,ux_l_inf,uy_L_2,uy_l_inf,sigmaxx_L_2,sigmaxx_l_inf,sigmaxy_L_2,sigmaxy_l_inf,sigmayy_L_2,sigmayy_l_inf
0.9886573325052778,0.0943069472667,0.0571446999058,0.180615492245,0.207188909887,0.160769822122,0.186054895037,0.231510396231,0.121999373687,0.852357639651,0.562583249179,1.48483762941,0.79244905714,0.248109552988,0.560798988145,0.218656167482,0.536339142705,0.239194710907,0.545803094584
0.6340332990709842,0.0507129517382,0.0303699131681,0.0984222234074,0.128351779888,0.101833485293,0.118699162979,0.0793635199031,0.0847286139736,0.267125945482,0.217931950275,0.497060647796,0.367636006433,0.134448743475,0.311980754866,0.130407360575,0.310657084,0.129113250449,0.30066683757
0.32904683851469807,0.0162277081392,0.0268156980091,0.0326931585621,0.0437056245149,0.0331143887793,0.0373019985474,0.0261815331548,0.0815418445045,0.0851422607568,0.235147387511,0.146885354636,0.579636942312,0.0471283911398,0.101483857904,0.0458338281165,0.0937048540073,0.0456625824847,0.0948541424278
0.16754966839339377,0.00270766556467,0.00810914229264,0.00940928548063,0.0120247703165,0.00951479276388,0.0113424424208,0.00808586791811,0.0411827695883,0.0231087451499,0.118693847541,0.0433461136792,0.230025740621,0.0135701357227,0.0266842392829,0.0135516965874,0.0238484084672,0.0129850130713,0.0257158976604
......@@ -167,6 +167,43 @@ class TestR13Convergence(object):
referrors = "referrors/" + name + "/errors.csv"
self.compare_errors(errors, referrors)
def test_r13_1_coeffs_nosources_norot_inflow_p2p2p2p2p2_gls(self):
r"""
Execute full linear R13 system test and check with reference errors.
Use second order equal elements.
Test case is similar to [TOR2017]_.
.. [TOR2017] Torrilhon, M. et al. (2017). “Hierarchical Boltzmann
simulations and model error estimation”. In: Journal of
Computational Physics 342 (2017), pp. 66–84.
==================== ===================================================
Parameter Value
==================== ===================================================
:math:`Kn` :math:`1.0`
:math:`\dot{m}` :math:`0`
:math:`r` :math:`0`
:math:`\theta_w^1` :math:`1.0`
:math:`v_t^1` :math:`0`
:math:`v_n^1` :math:`0`
:math:`p_w^1` :math:`0`
:math:`\epsilon_w^1` :math:`10^{-3}`
:math:`\theta_w^2` :math:`2.0`
:math:`v_t^2` :math:`-1.00 \sin(\phi)`
:math:`v_n^2` :math:`+1.00 \cos(\phi)`
:math:`p_w^2` :math:`-0.27 \cos(\phi)`
:math:`\epsilon_w^2` :math:`10^{3}`
Elements :math:`P_2P_2P_2P_2P_2`
Stabilization GLS
==================== ===================================================
"""
name = "r13_1_coeffs_nosources_norot_inflow_p2p2p2p2p2_gls"
self.run_solver("inputs/" + name + ".yml")
errors = name + "/" + "errors.csv"
referrors = "referrors/" + name + "/errors.csv"
self.compare_errors(errors, referrors)
def test_r13_1_coeffs_nosources_norot_inflow_p2p2p2p2p2_stab(self):
r"""
Execute full linear R13 system test and check with reference errors.
......@@ -195,7 +232,7 @@ class TestR13Convergence(object):
:math:`p_w^2` :math:`-0.27 \cos(\phi)`
:math:`\epsilon_w^2` :math:`10^{3}`
Elements :math:`P_2P_2P_2P_2P_2`
Stabilization Off
Stabilization CIP
==================== ===================================================
"""
name = "r13_1_coeffs_nosources_norot_inflow_p2p2p2p2p2_stab"
......@@ -204,6 +241,43 @@ class TestR13Convergence(object):
referrors = "referrors/" + name + "/errors.csv"
self.compare_errors(errors, referrors)
def test_r13_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_gls(self):
r"""
Execute full linear R13 system test and check with reference errors.
Test case is similar to [WES2019]_.
.. [WES2019] A. Westerkamp and M. Torrilhon. “Finite Element Methods
for the Linear Regularized 13-Moment Equations Describing Slow
Rarefied Gas Flows”. In: Journal of Computational Physics 389
(2019).
==================== ===================================================
Parameter Value
==================== ===================================================
:math:`Kn` :math:`1.0`
:math:`\dot{m}` :math:`(1-\frac{5R^2}{18{Kn}^2})\cos(\phi)`
:math:`r` :math:`0`
:math:`\theta_w^1` :math:`1.0`
:math:`v_t^1` :math:`10.0`
:math:`v_n^1` :math:`0`
:math:`p_w^1` :math:`0`
:math:`\epsilon_w^1` :math:`0`
:math:`\theta_w^2` :math:`0.5`
:math:`v_t^2` :math:`0.0`
:math:`v_n^2` :math:`0`
:math:`p_w^2` :math:`0`
:math:`\epsilon_w^2` :math:`0`
Elements :math:`P_1P_1P_1P_1P_1`
Stabilization GLS
==================== ===================================================
"""
name = "r13_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_gls"
self.run_solver("inputs/" + name + ".yml")
errors = name + "/" + "errors.csv"
referrors = "referrors/" + name + "/errors.csv"
self.compare_errors(errors, referrors)
def test_r13_1_coeffs_sources_rot_noinflow_p1p1p1p1p1_stab(self):
r"""