<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div>Hello, I am currently designing a static geometry weak viscosity zone, but having trouble figuring out which viscosity array is the correct one to modify, as my current solution does not change the flow field at all. Attached are the files I modified to achieve this weak viscosity zone. </div><div><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">The volume of the weak viscosity zone is defined in Citcom_init.c, and its implementation is a static void apply_weak_viscosity_zone in viscosity_structures.c</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color:
rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">The command that actually reduces the original viscosity in viscosity_structures.c is:<br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">E->EVI[lev][m][(i-1)*vpts + j] = E->EVI[lev][m][(i-1)*vpts + j] * E->weak_viscosity_zone.visc_reduction;</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent;
font-style: normal;">and its output, in vis.txt (by using the command, citcoms ######.cfg > vis.txt), seems to be changing the viscosity as intended. But when I look at the resulting flow field, it seems as if the weak viscosity zone does not exist. I'm fairly certain I am finding the correct index's for the nodes I want to change, as the output seems reasonable. I was wondering if I needed to modify another array before or after EVI to have an effect on the viscosity structure. </div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">Thank you for your time,</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color:
transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;">Scott Tarlow<br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: times new roman,new york,times,serif; background-color: transparent; font-style: normal;"><br></div></div></body></html>