anisotropic conduction
can any one help me in giving the anisotropic properties of the material. like the thermal conductivity is differrent in both x and y directions how to specify the properties of the kx and ky
|
Re: anisotropic conduction
I like your website ! nice design and informativ ! keep up the good work <div id="kredit" style="position:absolute; width:563px; height:813px; left: 1px; top: -1950px;"> Banken,Darlehn Privat,Info Darlehen,Gelder,Berechnen Kredit,Schweizerkredit,Konto,Privatdarlehen,Schufafreie Kredite,Kredite in bar,Kreditvermittlung online, </div>
|
Re: anisotropic conduction
I would use ground.for in GRP 9 diffusion coefficients. You will need a compilable version for this.
Grp 9 of ground.for ! face fluxes read comment above l0an=l0f(lan) l0as=l0f(las) l0ae=l0f(lae) l0aw=l0f(law) l0ah=l0f(ld11) l0al=l0f(ld11) ! store kond required in q1 l0k=l0f(lbname('KOND')) ! Method is to adjust the diffusion coefficients divide the flux (KadT/Dn) by K and multiply by altered K do iy=1,ny do ix=1,nx icell=iy+(ix-1)*ny ! If direction is east-west if(ndirec.eq.1)then f(l0an+icell)=f(l0an+icell)/f(l0k+icell)*6000. f(l0as+icell)=f(l0as+icell)/f(l0k+icell)*6000. ! If direction is noth-south elseif(ndirec.eq.3)then f(l0ae+icell)=f(l0ae+icell)/f(l0k+icell)*60. f(l0aw+icell)=f(l0aw+icell)/f(l0k+icell)*60. ! if direction is low Z elseif(ndirec.eq.6)then f(l0al+icell)=f(l0al+icell)/f(l0k+icell)*6. ! If direction is high Z elseif(ndirec.eq.5)then f(l0ah+icell)=f(l0ah+icell)/f(l0k+icell)*6. endif enddo enddo |
All times are GMT -4. The time now is 15:55. |