Exemplo n.º 1
0
 def test_plot2dcontourf(self):
     xc, yc = self.mesh2d.centers()
     f = field.fdata(self.model, self.mesh2d, [
         1.,
         euler.datavector(0., 0.),
         self.fn(np.sqrt(np.square(xc - 8) + np.square(yc - 4)))
     ])
     fig, ax = plt.subplots(1, 1)
     f.contourf('pressure', axes=ax, style='r-')
     return fig
Exemplo n.º 2
0
def fuv(x,y):
    vmag = .01 ; k = 10.
    return euler.datavector(0.*x+.4, 0.*x+.2)
Exemplo n.º 3
0
 def fuv(x, y):
     return euler.datavector(0. * x + .4, 0. * x + .2)
Exemplo n.º 4
0
def fuv(x, y):
    r = np.sqrt((x - lx / 2)**2 + (y - ly / 2)**2)
    return euler.datavector(
        0.5 - b / (2 * np.pi) * np.exp(0.5 * (1 - r**2)) * (y - ly / 2),
        0.5 + b / (2 * np.pi) * np.exp(0.5 * (1 - r**2)) * (x - lx / 2))