Particle source problem on windows v4

Hello, I have a basic simulation with a spacecraft and two different sources which works fine in windows v3.6 and v3.7.

When I run it on windows v4 with exactly the same input, I have no exceptions raised but the simulation runs with absurd time steps (changing but allways around E-70). When I turn off the particle sources (sourceFlags=0), it runs normally, but as soon as I turn on sourceFlag1 and sourceFlag2, the problem appears. Does someone understand what is happening? Thanks a lot

This is the beginning of the log:

Time 0.0 s => POTENTIALS: 0.0 V on SC ground (node0), 0.0 V on node1 ground, 0.0 V on node2 ground, 0.0 V on top of node0, 0.0 V on top of node1, 0.0 V on top of node2, 0.0 Vmin onTopOfNode0, 0.0 Vmin onTopOfNode1, 0.0 Vmin onTopOfNode2, 0.0 Vmax onTopOfNode0, 0.0 Vmax onTopOfNode1, 0.0 Vmax onTopOfNode2,

Monitoring numerics of simulation at time 0.0 ThreeDUnstructVolMesh: 93033 tetrahedra, 188545 triangles, 111306 edges, 15796 nodes Mesh quality (best = 1): min = 3.0169587E-4, max = 0.9892087, average = 0.5915012 histogram of mesh quality: number per interval between min = 3.0169587E-4 and max = 0.9892087 (unit -): 1122.0, 2097.0, 4283.0, 6749.0, 9802.0, 15510.0, 23212.0, 20358.0, 8762.0, 1138.0 PICVolDistrib ions1: 318393 particles, last numerical time step used = 0.0s (max allowed = 6.3390494E-6s), numerical speed up = 1.0 PICVolDistrib elec1: 318341 particles, last numerical time step used = 0.0s (max allowed = 1.4847186E-7s), numerical speed up = 1.0 PICVolDistrib source1: 0 particles, last numerical time step used = 0.0s (max allowed = 4.1936855E-6s), numerical speed up = 1.0 PICVolDistrib source2: 0 particles, last numerical time step used = 0.0s (max allowed = 8.572027E-9s), numerical speed up = 1.0 Last local volume potential change: 0.0V maximum change, and 0.0V average change (last plasma time step = 0.0s)

Time: 0.0, Dt = 4.286327E-7

Number of particles reflected on symmetry boundaries : 0 Number of particles reflected on symmetry boundaries : 0 Number of particles reflected on symmetry boundaries : 0 Number of particles reflected on symmetry boundaries : 0 Variation on node 0 : -0.0011697058798745275 Ratio of validity : 1.0000000150474662E30 Validity ratio too big : 1.0000000150474662E30 New time step : 2.1431633617487145E-37 Integration duration performed so far : 0.0 Validity Table : 0.0 Variation on node 0 : -1.1700966277001054E-33 Ratio of validity : 1.0000000150474662E30 Validity ratio too big : 1.0000000150474662E30 New time step : 2.1431633294995365E-67 Integration duration performed so far : 0.0 Validity Table : 0.0 Variation on node 0 : -0.0 Ratio of validity : 0.0 Variation on node 0 : -0.0 Ratio of validity : 0.0 Variation on node 0 : -0.0 … Variation on node 0 : -1.4013132620366003E-36 Ratio of validity : 1.0000000150474662E30 Validity ratio too big : 1.0000000150474662E30 New time step : 2.143163297250355E-70 Integration duration performed so far : 2.381292588332814E-41 Variation on node 0 : -0.0 Ratio of validity : 0.0 Variation on node 0 : -0.0 … and so on (I have not been able to finish a simulation)

Matias

posted by Matias Wartelski at Feb 18, 2010 9:25 AM