org.jscience.physics.fluids.dynamics.characteristics
Class QuadraticCharacteristics

java.lang.Object
  extended by org.jscience.physics.fluids.dynamics.characteristics.QuadraticCharacteristics

public class QuadraticCharacteristics
extends java.lang.Object

This class calculates the values of the different magnitudes cx, cy, cc y cd at the root of the caracteristics for a 2D meshing with a speed field. Tipically these magnitudes correspond with the two speed components x and y and with the concentration of a contaminant or the turbulent mu.


Constructor Summary
QuadraticCharacteristics(KernelADFC kadfc, NavierStokesMesh m, double deltaT, boolean contam)
          Constructor.
 
Method Summary
 void antidifusorPoisson(double[] a, int coord)
          experimental Method
 void calculate(double[] cxa, double[] cya, double[] cca, double[] cda)
          principal method of the characteristics.
 void characteristicForwardCorrection()
          experimental Method
 void correctAntidifusive(double[] c)
          experimental Method
 void correctionForwardCharacteristics2()
          experimental Method
 void correctMinMax(double[] orig, double[] trans)
          experimental Method
 void initializeCaches()
          Call this method before calculate().
 void obtainTransformationCurve(double[] x, double[] y)
          experimental Method
 void setSpeedFields(double[] vxa, double[] vya, double[] vxb, double[] vyb)
          Indicatestes the values of the two speed components in each node at the instants n and (n-1)
 
Methods inherited from class java.lang.Object
clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait
 

Constructor Detail

QuadraticCharacteristics

public QuadraticCharacteristics(KernelADFC kadfc,
                                NavierStokesMesh m,
                                double deltaT,
                                boolean contam)
Constructor. Takes as arguments the mesh and the time step duration

Parameters:
kadfc - kernel parent of the current process
m - mesh to consider
deltaT - time step
contam - indicates if the contaminant should be transported or not
Method Detail

calculate

public void calculate(double[] cxa,
                      double[] cya,
                      double[] cca,
                      double[] cda)
principal method of the characteristics. Takes the values of the magnitudes to be estudied in each of the nodes (cx, cy, cc, cd) and returns a solution in those given vectors.

Parameters:
cxa - value component X of the magnitude in the nodes
cya - value component Y of the magnitude in the nodes
cca - value auxiliary magnitude 1 in the nodes
cda - value auxiliary magnitude 2 in the nodes

characteristicForwardCorrection

public void characteristicForwardCorrection()
experimental Method


correctionForwardCharacteristics2

public void correctionForwardCharacteristics2()
experimental Method


antidifusorPoisson

public void antidifusorPoisson(double[] a,
                               int coord)
experimental Method

Parameters:
a - DOCUMENT ME!
coord - DOCUMENT ME!

obtainTransformationCurve

public void obtainTransformationCurve(double[] x,
                                      double[] y)
experimental Method

Parameters:
x - DOCUMENT ME!
y - DOCUMENT ME!

correctAntidifusive

public void correctAntidifusive(double[] c)
experimental Method

Parameters:
c - DOCUMENT ME!

correctMinMax

public void correctMinMax(double[] orig,
                          double[] trans)
experimental Method

Parameters:
orig - DOCUMENT ME!
trans - DOCUMENT ME!

initializeCaches

public void initializeCaches()
Call this method before calculate(). The invocation is done by default if needed.


setSpeedFields

public void setSpeedFields(double[] vxa,
                           double[] vya,
                           double[] vxb,
                           double[] vyb)
Indicatestes the values of the two speed components in each node at the instants n and (n-1)

Parameters:
vxa - component x of the speed at instant n
vya - component y of the speed at instant n
vxb - component x of the speed at instant n-1
vyb - component y of the speed at instant n-1