Κεφάλαιο 1 - Nemertes

ΑΝΑΠΤΥΞΗ ΚΩ∆ΙΚΑ ΜΟΡΙΑΚΗΣ ∆ΥΝΑΜΙΚΗΣ
242
}
else if (x > xBoundaryEnd)
{
result = xBoundaryStart + ((Math.abs(x - xBoundaryEnd)) % xBoundarySize);
}
return result;
}
public double yBoundary(double y)
{
double result = y;
if (y < yBoundaryStart)
{
result = yBoundaryEnd - ((Math.abs(yBoundaryStart - y)) % yBoundarySize);
}
else if (y > yBoundaryEnd)
{
result = yBoundaryStart + ((Math.abs(y - yBoundaryEnd)) % yBoundarySize);
}
return result;
}
public double zBoundary(double z)
{
//special conditions apply here:
// z bottom is reflecting
double result = z;
if (z > zBoundaryEnd)
{
result = zBoundaryEnd - Math.abs(z - zBoundaryEnd);
}
// z top is free
if (z < 0 && MDCoordinator.DEBUG)
System.out.println("Z NEGATIVE");
return result;
}
Για τους άξονες x και y οι συνοριακές συνθήκες είναι περιοδικές, ενώ για τον άξονα z στην
κάτω πλευρά της υπολογιστικής κυψελίδας αναφοράς, ανακλαστικά ενώ στην πάνω πλευρά
της ελεύθερα. Κλείνοντας, µια ακόµα συνάρτηση που υλοποιείται είναι η ανάθεση ταχύτητας
σύµφωνα µε την κατανοµή Maxwell-Boltzmann:
/**
* Used in MaxwellBoltzmann velocity to give UD[0,1]
*/
Random r = new Random(System.currentTimeMillis());
public double getNextMaxwellBoltzmannVelocity()
{
double Xn = 0;