Initialisation Dipoles: both dipoles: velocity: U = 2; maximum initial vorticity 22.1 radius: a = 1 dipole 1: position: (-3.5,0) angle of dipole axis with X-axis: 0 degrees dipole 2: position: (0,+3.5) angle of dipole axis with X-axis: -90 degrees Computation: domain: X = (-5, ..., 3), Y = (-5, ..., 5) no. of grid points from wall to wall: (129, 129) boundary conditions: stress-free along all walls standard Arakawa scheme used no rotation effects included no bottom topography viscosity: nu = 10^-3 (i.e. Re=1000 at T=0) time step: dt = 0.01 final time: T = 4 Tracers: along vorticity contours (i.e. along streamlines) of both Lamb dipoles at vorticity levels: 0.1, 6.0, 12.0, 18.0; and a single tracers at each extremum Note The tracer lines are not always smooth. The reason for this is that the representation of the lines by points is not always accurate enough. But since these tracers only passively follow the flow, this has no effect on the computations. It needs repair, though.
<=== Collision under 90 degrees of two Lamb dipoles
Jos van Geffen --
Home |
Site Map |
Contact Me