Home > atmlab > scattering > Bruggeman_mix_rule.m

Bruggeman_mix_rule

PURPOSE ^

Bruggeman_mix_rule Effective refractive index and effective density of

SYNOPSIS ^

function [n rho]= Bruggeman_mix_rule(n1,n2,rho1,rho2,fraction2)

DESCRIPTION ^

 Bruggeman_mix_rule     Effective refractive index and effective density of 
                        a two-medias mixture.
                         

                        This function uses the Bruggeman formulation
                        to calculate the effective refractive index and
                        effective density of a mixture that includes two
                        medias.

                        
 
 FORMAT     [n rho] = Bruggeman_mix_rule(n1,n2,rho1,rho2,fraction2)
        
 OUT  n            Complex dielectric constant           real or imaginary
      rho          Uniform density of the mixture        [Kg/m3]

 IN   n1           refractive index of Matrix medium     real or imaginary
      n2           refractive index of Inclusion         real or imaginary
      rh1          Density of matrix medium              [Kg/m3]
      rho2         Density of inclusion                  [Kg/m3]
      fraction2    The volum fraction of the inclusion

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

Bruggeman_mix_rule.m

SOURCE CODE ^

0001 % Bruggeman_mix_rule     Effective refractive index and effective density of
0002 %                        a two-medias mixture.
0003 %
0004 %
0005 %                        This function uses the Bruggeman formulation
0006 %                        to calculate the effective refractive index and
0007 %                        effective density of a mixture that includes two
0008 %                        medias.
0009 %
0010 %
0011 %
0012 % FORMAT     [n rho] = Bruggeman_mix_rule(n1,n2,rho1,rho2,fraction2)
0013 %
0014 % OUT  n            Complex dielectric constant           real or imaginary
0015 %      rho          Uniform density of the mixture        [Kg/m3]
0016 %
0017 % IN   n1           refractive index of Matrix medium     real or imaginary
0018 %      n2           refractive index of Inclusion         real or imaginary
0019 %      rh1          Density of matrix medium              [Kg/m3]
0020 %      rho2         Density of inclusion                  [Kg/m3]
0021 %      fraction2    The volum fraction of the inclusion
0022        
0023 % 2013-08-01        Created by Maryam Jamali
0024 
0025 
0026 
0027 
0028 function [n rho]= Bruggeman_mix_rule(n1,n2,rho1,rho2,fraction2)
0029 
0030 % eeff is calculated by solving the following equation
0031 %A=solve('(fraction2.*(einc-x)./(einc+2.*x))+(((1-fraction2).*(emat-x))./(emat+2.*x))==0');
0032 
0033 emat=power(n1,2);
0034 einc=power(n2,2);
0035 
0036 eeff=emat./2 - einc./4 + (3.*einc.*fraction2)./4 - (3.*emat.*fraction2)./4 +...
0037      (9.*einc.^2.*fraction2.^2 - 6.*einc.^2.*fraction2 + einc.^2 - 18.*einc.*emat.*fraction2.^2 +...
0038      18.*einc.*emat.*fraction2 + 4.*einc.*emat + 9.*emat.^2.*fraction2.^2 - 12.*emat.^2.*fraction2 +...
0039      4.*emat.^2).^(1/2)/4;
0040 
0041 n=power(eeff,.5);
0042 
0043 rho=(1-fraction2)*rho1+fraction2*rho2;

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005