NavierStokesElement.cpp File Reference

#include <MoFEM.hpp>
#include <MethodForForceScaling.hpp>
#include <NavierStokesElement.hpp>

Go to the source code of this file.

Functions
VectorDouble3	stokes_flow_velocity (double x, double y, double z)

Function Documentation

stokes_flow_velocity()

VectorDouble3 stokes_flow_velocity	(	double	x,
		double	y,
		double	z
	)

Examples

NavierStokesElement.cpp.

Definition at line 1071 of file NavierStokesElement.cpp.

                                                                 {
  double r = std::sqrt(x * x + y * y + z * z);
  double theta = acos(x / r);
  double phi = atan2(y, z);
  double ur = cos(theta) * (1.0 + 0.5 / (r * r * r) - 1.5 / r);
  double ut = -sin(theta) * (1.0 - 0.25 / (r * r * r) - 0.75 / r);
  VectorDouble3 res(3);
  res[0] = ur * cos(theta) - ut * sin(theta);
  res[1] = ur * sin(theta) * sin(phi) + ut * cos(theta) * sin(phi);
  res[2] = ur * sin(theta) * cos(phi) + ut * cos(theta) * cos(phi);
  return res;
}