tension.c File Reference

Define Documentation

#define PHI_R_LINEAR	(	in,
		i,
		t1	)

Value:

(((t1 - (in).t[i])*(in).p[i+1] + ((in).t[i+1] - t1)*(in).p[i])/ \
          ((in).t[i+1] - (in).t[i]))

Definition at line 9 of file tension.c.

Referenced by interface_tensionu(), and interface_tensionv().

Function Documentation

void interface_tensionu	(	interface	in,
		interface_cut	cut,
		real2D	u,
		real2D	v,
		real2D	cu,
		real2D	ap,
		real2D	syu,
		real2D	p,
		real2D	pi,
		real2D	ni,
		int	nx,
		int	ny
	)

Definition at line 87 of file tension.c.

References listinter::coord, fill0(), INP, INU, listinter::ip, interface_cut::list, interface_cut::ncut, PHI_R_LINEAR, real, real2D, splint, interface::spx, interface::spy, listinter::t, listinter::type, VER_TYPE, interface_cut::xo, and interface_cut::yo.

Referenced by timestep().

00092 {
00093   int i, j, l;
00094   real t, t1, t2;
00095   polynom3 px, py;
00096 
00097   fill0(pi, nx, ny);
00098   fill0(ni, nx, ny);
00099 
00100   /* For each intersection */
00101   for (l = 0; l < cut.ncut - 1; l++) {
00102     px = in.spx[cut.list[l].ip];
00103     py = in.spy[cut.list[l].ip];
00104 
00105     t1 = cut.list[l].t; t2 = cut.list[l+1].t;
00106     t = 0.5*(t1 + t2);
00107     i = splint(px, t) - cut.xo;
00108     j = splint(py, t) - cut.yo;
00109 
00110     /* Intersection with the verticals */
00111     if (cut.list[l].type == VER_TYPE) {
00112       if (cut.list[l].coord == i) {
00113         if (!INP(i-1,j)) {
00114           pi[i-1][j] += PHI_R_LINEAR(in, cut.list[l].ip, t1)/
00115             (splint(py, t1) - 2.);
00116           ni[i-1][j] += 1.0;
00117         }
00118       }
00119       else if (!INP(i,j)) {
00120         pi[i][j] += PHI_R_LINEAR(in, cut.list[l].ip, t1)/
00121           (splint(py, t1) - 2.);
00122         ni[i][j] += 1.0;
00123       }
00124     }
00125     
00126     if (INU(i,j)) u[i][j] -= phi_r_dr(in, cut.list[l].ip, t1, t2)/cu[i][j];
00127   }
00128   for (i = 2; i <= nx; i++)
00129     for (j = 2; j < ny; j++) 
00130       if (INU(i,j)) {
00131         if (pi[i-1][j] != 0.0)
00132           u[i][j] += syu[i][j]*pi[i-1][j]/(cu[i][j]*ni[i-1][j]);
00133         if (pi[i][j] != 0.0)
00134           u[i][j] -= syu[i+1][j]*pi[i][j]/(cu[i][j]*ni[i][j]);
00135       }
00136 }

void interface_tensionv	(	interface	in,
		interface_cut	cut,
		real2D	u,
		real2D	v,
		real2D	cv,
		real2D	ap,
		real2D	sxv,
		real2D	av,
		real2D	p,
		real2D	pi,
		real2D	ni,
		int	nx,
		int	ny
	)

Definition at line 139 of file tension.c.

References listinter::coord, fill0(), HOR_TYPE, INP, INV, listinter::ip, interface_cut::list, interface_cut::ncut, PHI_R_LINEAR, real, real2D, splint, interface::spx, interface::spy, listinter::t, listinter::type, interface_cut::xo, and interface_cut::yo.

Referenced by timestep().

00144 {
00145   int i, j, l;
00146   real t, t1, t2;
00147   polynom3 px, py;
00148 
00149   fill0(pi, nx, ny);
00150   fill0(ni, nx, ny);
00151 
00152   /* For each intersection */
00153   for (l = 0; l < cut.ncut - 1; l++) {
00154     px = in.spx[cut.list[l].ip];
00155     py = in.spy[cut.list[l].ip];
00156 
00157     t1 = cut.list[l].t; t2 = cut.list[l+1].t;
00158     t = 0.5*(t1 + t2);
00159     i = splint(px, t) - cut.xo;
00160     j = splint(py, t) - cut.yo;
00161 
00162     /* Intersection with the horizontals */
00163     if (cut.list[l].type == HOR_TYPE) {
00164       if (cut.list[l].coord == j) {
00165         if (!INP(i,j-1)) {
00166           pi[i][j-1] += PHI_R_LINEAR(in, cut.list[l].ip, t1)
00167             /(splint(py, t1) - 2.);
00168           ni[i][j-1] += 1.0;
00169         }
00170       }
00171       else if (!INP(i,j)) {
00172         pi[i][j] += PHI_R_LINEAR(in, cut.list[l].ip, t1)
00173           /(splint(py, t1) - 2.);
00174         ni[i][j] += 1.0;
00175       }
00176     }
00177 
00178     if (INV(i,j)) v[i][j] += phi_r_dz(in, cut.list[l].ip, t1, t2)/cv[i][j];
00179   }
00180   for (i = 2; i < nx; i++)
00181     for (j = 3; j <= ny; j++) 
00182       if (INV(i,j)) {
00183         if (pi[i][j-1] != 0.0)
00184           v[i][j] += (sxv[i][j] + av[i][j]/2.)*pi[i][j-1]
00185             /(cv[i][j]*ni[i][j-1]);
00186         else if (!INP(i,j-1))
00187           v[i][j] += av[i][j]/2.*p[i][j]/cv[i][j];
00188         if (pi[i][j] != 0.0)
00189           v[i][j] -= (sxv[i][j+1] - av[i][j]/2.)*pi[i][j]
00190             /(cv[i][j]*ni[i][j]);
00191         else if (!INP(i,j))
00192           v[i][j] += av[i][j]/2.*p[i][j-1]/cv[i][j];
00193       }
00194 }

void interfacial_pressure	(	interface	in,
		real2D	u,
		real2D	v,
		real2D	ap,
		real	sigma,
		real	mu,
		real	pg,
		real	tau,
		int	nx,
		int	ny
	)

Definition at line 197 of file tension.c.

References interface::n, real, real2D, interface::spx, interface::spy, sq, interface::t, and interface::y.

Referenced by timestep().

00200 {
00201   real h = 1.0/(nx - 2.);
00202   int i;
00203 
00204   sigma *= sq(tau)/(h*h*h);
00205   mu *= tau/sq(h);
00206   pg *= sq(tau/h);
00207 
00208   for (i = 0; i < in.n - 1; i++)
00209     in.p[i] = pg*(in.y[i] - 2.) + 
00210       phi_r(in.spx[i], in.spy[i], in.t[i], sigma, mu, 
00211                          u, v, ap, nx, ny);
00212   in.p[in.n-1] = pg*(in.y[in.n-1] - 2.) +
00213     phi_r(in.spx[in.n-2], in.spy[in.n-2], in.t[in.n-1], 
00214           sigma, mu, u, v, ap, nx, ny);
00215 }

Generated on Wed Feb 19 22:28:45 2003 for Markers by

1.2.18


Defines
#define	PHI_R_LINEAR(in, i, t1)
Functions
void	interface_tensionu (interface in, interface_cut cut, real2D u, real2D v, real2D cu, real2D ap, real2D syu, real2D p, real2D pi, real2D ni, int nx, int ny)
void	interface_tensionv (interface in, interface_cut cut, real2D u, real2D v, real2D cv, real2D ap, real2D sxv, real2D av, real2D p, real2D pi, real2D ni, int nx, int ny)
void	interfacial_pressure (interface in, real2D u, real2D v, real2D ap, real sigma, real mu, real pg, real tau, int nx, int ny)

tension.c File Reference

Defines

Functions

Define Documentation

Function Documentation