From 1cbb1e0a088ce8a5473aecaef7abb98cababb4e6 Mon Sep 17 00:00:00 2001 From: Guillaume Giudicelli Date: Fri, 4 Dec 2020 14:53:52 -0700 Subject: [PATCH] Add test for location dependent material properties in FV, evaluated at boundary by FVDiffusion kernel, #15194 --- .../fv_material_quadrature.i | 78 ++++++++++++++++++ .../gold/fv_material_quadrature_out.e | Bin 0 -> 32464 bytes test/tests/materials/boundary_material/tests | 8 ++ 3 files changed, 86 insertions(+) create mode 100644 test/tests/materials/boundary_material/fv_material_quadrature.i create mode 100644 test/tests/materials/boundary_material/gold/fv_material_quadrature_out.e diff --git a/test/tests/materials/boundary_material/fv_material_quadrature.i b/test/tests/materials/boundary_material/fv_material_quadrature.i new file mode 100644 index 000000000000..c85946bb068b --- /dev/null +++ b/test/tests/materials/boundary_material/fv_material_quadrature.i @@ -0,0 +1,78 @@ +# +# Parsed material properties depend on the physical location of the element +# This requires the initialization of the quadrature in the FVFlux loop + +[Mesh] + type = GeneratedMesh + dim = 2 + xmin = 0 + xmax = 1 + ymin = 0 + ymax = 1 + nx = 4 + ny = 4 + elem_type = QUAD9 +[] + +[Functions] + [linear_x] + type = ParsedFunction + value = 'x' + [] +[] + +[Variables] + [u] + family = MONOMIAL + order = CONSTANT + fv = true + [] +[] + +[FVKernels] + [diff] + type = FVDiffusion + variable = u + coeff = k1 + [] + + [r] + type = FVReaction + variable = u + [] +[] + +[FVBCs] + [all] + type = FVDirichletBC + variable = u + boundary = 'left right bottom top' + value = 1 + [] +[] + +[Materials] + [k1] + type = ADGenericFunctionMaterial + prop_names = 'k1' + prop_values = linear_x + block = 0 + [] +[] + +[Executioner] + type = Steady + solve_type = NEWTON + + petsc_options_iname = '-pc_type' + petsc_options_value = 'lu' +[] + +[Problem] + kernel_coverage_check = false +[] + +[Outputs] + execute_on = 'timestep_end' + exodus = true +[] diff --git a/test/tests/materials/boundary_material/gold/fv_material_quadrature_out.e b/test/tests/materials/boundary_material/gold/fv_material_quadrature_out.e new file mode 100644 index 0000000000000000000000000000000000000000..37c9df8fee3eff355e5d163eddd858735498ca27 GIT binary patch literal 32464 zcmeHQ&y(F&R_p5j?4r{JuFPKi_}G9}kUzJLF@2Z;5_C7}ya*63j)0*`7_N zFM*CFexvW*CySKvIrW|q*>UKJAfQh;m@ZSu1}a%i9_rcQXFGshxn#^u!QjA z3{P+(1`<0TwpJ2$wj^vl4SQC7$Cq#m{6;Q}aK$jf6~hQ8!=x-sBjdni`crQ*-CW-Z z;JpF=anLiMXU&mo9N1yx20pwCUj#jI`FqEDxIR-6{*vPbW(>_6qrmlJ<3SJ(kRA`I zN2Z1R9mjX^2f1f_ejG2JUxz`M@t7)7@8Ndgi?JJf1cRd;Y{YgHM`wm+Y=^Gt zMQ8fK*dLhTk)eXlI0xCi@q1%)5CS%|jbI!vCqO5L=Ak>JJ;B}#80ucUPki-Vf%|7c zydIOb`d~jug`m7&04eu0ZWa^Ix8K5b$?i`g0>6rX!jgQ0-+^0t1wV`6ABDk)wxf8d z58gkO<_E_63)i;rJN^D^@Gtq%%=c(B%JEyYU$~{O;b#&2`R_>^M)mr!ofcNBOT;X0^!OwBYKmC1{Bg%n);d0_RLHQOXeHb% z_Z2P+Y2v2w<|C=^)PIs!gF`sE5!XELogA-< z+ldNZ=2C^Lju}K&;A!QAFxur$gD8V?xx9Gt0|dT!fPz1JZ`JP8py9~J zvYB?R9LtM8NIc)Zg==E7s0)qQ?8DD@nD_Zj<@oLN`zJrbJ?bYXFhbixNfg1)cfJgw z?vz8A`@4zfiW+W zeihCbz-_22Q=!bS!aq)*t1C+zxQ(J;g{SdJxWJeCwdJpqw@$v#<(bw~g}#)Y)HTyr zQI8)a9NBZ_Pp`?-iM}(x;;XWQhTqBX8uD*^n7ol)l|89&s z=rm{Pl+6CPN-2%M>YLbSm>s30k{yz>J%5UCGO| z<#Uz3O#6HqQwL0E>XF}2=}lc0)11%w4W>Q$sMAv*>VUk|1E2F7)Cc*fU!@oFQa5}~ zT~kNoqpn{EQCH-pzNmA4L+OpWFY1oZ`3>rje3V1!lDyO>pYt2kEBR<2O26c#j`^Ij zya6H~a70^|XvH?ElH+DdmL9{7l6Xc~`@HxNn4G{TgXK#aO zC*-BA@Hu7K1(A=or)-D3?|_=v&oS^(A9=f>&Cni-@6%r-=IbR8?TGeKB<6Ku6Yt3a z@tx;E>U*@EJgM*V8;JR)qd%Ma&M9#*nCXxz%%)bz$x7;Ov@Xt@PK;l?yA*z+PK1FFHlYTibVTTeEONoLyc$yUdDL)7`$l zeXCQztK%?)8H+oZcbLQ%Nx2icV#~IKSQo47m(HzRT01Y3uTJ3YL)#i-0Ku~*CPuD@ zBMCeWxO(Bz>cvKYbszk{f6k0|%^!_pv4fHPPg}NQ`rP|g$EhTQ>zh{W9@vvyH+9z& zula#7upM*k#cb`=I0CvcrY#I742FiB65;#;r@|y~Ph8);vwizc4NmsB3=m^m6Cr!z z+V0)W>*?e{_Fm{fptHc0I56uD5roJJ43OpU?x|6K#s~HHyaXh0;KMG zV)NZQY%Vl(LR0{?2b5#i)u{Rsa!g){Q*!NzH#gtCeY;_VKY^Um_{&Wh?P|cC zBwGt%j^hARQ5ePb2grnhk4ej1TJ{8+`OV!YQ8+|5`|o3ZW)o!&>hYq}m$24sZpXGc zIBFVcq@iWmBjh|EGfs!b$hAz&_PUsgmCYQ{6YPpLcd1g$!GY;pWd^Qad{N76bXi7Z znOJ|2S1}S)9}rC}RC>0dO774`?jEBE@7qp*X=PXxRuC9{D{2&d6tF2^DjyPaD(d=L zhC#=ursQf1tY$DnYu_+EFL+=N49^_)2d06>3dLh5sMb-!SwJriywo7uq@f>&frlxV zy^6SVGuE41n)n0`tS8fI_C0yZ1Q8NcX^pevom!aftdLnHbHp5$(0yIXZb>G~ZGM%0n zVYQOg5Dp!1!E*Z0!SbK+Xn?_ooN)5CC)U?4u3tQVVeMk0cmvRJ21a$o(i0Cxjl~0_ z5Po7eYy+3&J^h@L0+!`t##!GE0XT3E%7jbPaMgZ&RmTniNWC~?AGOA;(>LrCuqvzc zghpYgl)|#GQ>9cDu1tU=>**A*Tz0Y-x`VO=OaOJvZ~?u5ag2(-0xrkYJ>VX{TIIoE}*+Yc`9yND31r;cP|K6JEp7(&U2xOj?9h^?-iXO99jR zew6f*RfJn9C6N$KHw#OrP7_f4Uk3eE`2|<}70EaOuM;$03EVb<2(l7#|3qJBpOO8Gl z>O4kScci_h3*JU4_F){F=v!bf&f&~xODhQ_-~om*usXdseA5Qd9ikSnFlSlj+C;e` zvI%@oL`Ueu3`P#9C!p3rKuFYx6G z{U*3Q;UgoVsx3*cvLKd+%V!?B6r4am7Pa2ncIew)mo?uc0Xn#omDIp>oUZgj+xUM` zNjc6~?RTq5XGU30i}I$-wJ8s+z;>K6FoyL+x&yIBGWD46)Yf{brK)%4UI{pNoVZQ# z_35vqlFwRW~be=Z+e7b6 zCLY$5oI^xk!L=EXlz_78E2*on#3*e_m73j5Amxc{?4hi3ZM58RKjkx6jbf6lyDExy zs3zc&-KR6k*3T((>pT_?x+{X_8Vvz;`^&Cxd1IXCg^gF(U4Zik?(9UnJJ#&n+JiVZAV$qr0axfKf{j8@dWSC8it$I8YH;<@qHe7q|*`_lCmB7Me zwEDN$l7l5~kz3MSSAbf03jnMx7zl}7NF}h#KdCqkh}V z=?T428w;YV)mnHW#33v@lV|ourz}jAayaBV?n~O&|JRwKak#-~Wcvf7Y~L-5ApO!N2WYT6U`{Wi)^9n;edBUb zRcV>1FQFB9#rMLamP@O6=$=huKxo-x$T1}tQ@%(%h>9$eECX{rlXeQQrGIENGX((0W1+%+R^(Bn< zg9j>pC0Bj3L6s*LoE!TX5DvmzCrZM%knqUI(EvPZu0%c-U%CT9v2Ar*9v;N^OGQfq z^`m?^FuMjU?Rcw?w2*LERklh@UBJhnTX8E4a9+T$!!~yJwcfmqvTfQnpq&S-qzZ7|3(+<}6Fl zos$Z&qaPbXLnA?sm)n%s@{hDCMk5jPp7^czetWK= zOh3dNeRA$;VV*imMBU=+6BxR;!w}n+al8i36$vmw<>SN`|Dhf7&wMxmW~pNfn<_%TH_%}VAZHOiu26ICPI;P z9M~KYvVChR8TA3>>T`dLLnk7f636}OSSZ77tVts3N4XkHBI*NX(;s~^3|$noXnIB_ zj!{Y4+xmdDNBT6&!>SMyB76w_;*eege+F5|c8qgRyr*~@UDH(r<{&v-7E5@P#2Ei$ zOH@C~a#;e7c1zAMDB5)!KpyEM6CjZ!@f#Z8~`vLk5P68$; zwWA%DGApchv%08VoyyfV5=z4?FY0;L+%T}Dd_wYmRZ-vZ02?> zdP3<#V;R4@XP3_`Kw$kUpr#+tCd$sd$$&P1HLLD*niH-@G)3|p8Ypy z@*hh7l>Ud(ne;!DzM1}q(w}eq@Lx`w-}&B!jUT`D{lEI?#JTi8m_ADXgXviM{Cweg I