olb::MagneticForceFromCylinder3D< T, S > Class Template Reference

Magnetic field that creates magnetization in wire has to be orthogonal to the wire. More...

#include <frameChangeF3D.h>

Inheritance diagram for olb::MagneticForceFromCylinder3D< T, S >:

Collaboration diagram for olb::MagneticForceFromCylinder3D< T, S >:

Public Member Functions
	MagneticForceFromCylinder3D (CartesianToCylinder3D< T, S > &car2cyl, T length, T radWire, T cutoff, T Mw, T Mp=T(1), T radParticle=T(1), T mu0=4 *3.14159265e-7)
	Magnetic field that creates magnetization in wire has to be orthogonal to the wire.
bool	operator() (T output[], const S x[]) override
	operator writes the magnetic force in a point x util::round a cylindrical wire into output field
Public Member Functions inherited from olb::AnalyticalF< D, T, S >
AnalyticalF< D, T, S > &	operator- (AnalyticalF< D, T, S > &rhs)
AnalyticalF< D, T, S > &	operator+ (AnalyticalF< D, T, S > &rhs)
AnalyticalF< D, T, S > &	operator* (AnalyticalF< D, T, S > &rhs)
AnalyticalF< D, T, S > &	operator/ (AnalyticalF< D, T, S > &rhs)
Public Member Functions inherited from olb::GenericF< T, S >
virtual	~GenericF ()=default
int	getSourceDim () const
	read only access to member variable _m
int	getTargetDim () const
	read only access to member variable _n
std::string &	getName ()
	read and write access to name
std::string const &	getName () const
	read only access to name
bool	operator() (T output[])
	wrapper that call the pure virtual operator() (T output[], const S input[]) from above
bool	operator() (T output[], S input0)
bool	operator() (T output[], S input0, S input1)
bool	operator() (T output[], S input0, S input1, S input2)
bool	operator() (T output[], S input0, S input1, S input2, S input3)

Protected Attributes
CartesianToCylinder3D< T, S > &	_car2cyl
T	_length
	length of the wire, from origin to _car2cyl.axisDirection
T	_radWire
	wire radius
T	_cutoff
	maximal distance from wire cutoff/threshold
T	_Mw
	saturation magnetization wire, linear scaling factor
T	_Mp
	magnetization particle, linear scaling factor
T	_radParticle
	particle radius, cubic scaling factor
T	_mu0
	permeability constant, linear scaling factor
T	_factor
	factor = mu0*4/3.*PI*radParticle^3*_Mp*radWire^2/r^3

Additional Inherited Members
Public Types inherited from olb::AnalyticalF< D, T, S >
using	identity_functor_type = AnalyticalIdentity<D,T,S>
Public Types inherited from olb::GenericF< T, S >
using	targetType = T
using	sourceType = S
Public Attributes inherited from olb::GenericF< T, S >
std::shared_ptr< GenericF< T, S > >	_ptrCalcC
	memory management, frees resouces (calcClass)
Static Public Attributes inherited from olb::AnalyticalF< D, T, S >
static constexpr unsigned	dim = D
Protected Member Functions inherited from olb::AnalyticalF< D, T, S >
	AnalyticalF (int n)
Protected Member Functions inherited from olb::GenericF< T, S >
	GenericF (int targetDim, int sourceDim)

Detailed Description

template<typename T, typename S>
class olb::MagneticForceFromCylinder3D< T, S >

Magnetic field that creates magnetization in wire has to be orthogonal to the wire.

To calculate the magnetic force on particles around a cylinder (J. Lindner et al. / Computers and Chemical Engineering 54 (2013) 111-121) Fm = mu0*4/3.*PI*radParticle^3*_Mp*radWire^2/r^3 * [radWire^2/r^2+util::cos(2*theta), util::sin(2*theta), 0]

Definition at line 498 of file frameChangeF3D.h.

Constructor & Destructor Documentation

MagneticForceFromCylinder3D()

template<typename T , typename S >

olb::MagneticForceFromCylinder3D< T, S >::MagneticForceFromCylinder3D	(	CartesianToCylinder3D< T, S > &	car2cyl,
		T	length,
		T	radWire,
		T	cutoff,
		T	Mw,
		T	Mp = T(1),
		T	radParticle = T(1),
		T	mu0 = 4*3.14159265e-7 )

Magnetic field that creates magnetization in wire has to be orthogonal to the wire.

To calculate the magnetic force on particles around a cylinder (J. Lindner et al. / Computers and Chemical Engineering 54 (2013) 111-121) Fm = mu0*4/3.*PI*radParticle^3*_Mp*radWire^2/r^3 * [radWire^2/r^2+util::cos(2*theta), util::sin(2*theta), 0]

Definition at line 1162 of file frameChangeF3D.hh.

  : AnalyticalF3D<T, S>(3),
    _car2cyl(car2cyl),
    _length(length),
    _radWire(radWire),
    _cutoff(cutoff),
    _Mw(Mw),
    _Mp(Mp),
    _radParticle(radParticle),
    _mu0(mu0)
{
  // Field direction H_0 parallel to fluid flow direction, if not: *(-1)
  _factor = -_mu0*4/3*M_PI*util::pow(_radParticle, 3)*_Mp*_Mw*util::pow(_radWire, 2);
}

References olb::MagneticForceFromCylinder3D< T, S >::_factor, olb::MagneticForceFromCylinder3D< T, S >::_Mp, olb::MagneticForceFromCylinder3D< T, S >::_mu0, olb::MagneticForceFromCylinder3D< T, S >::_Mw, olb::MagneticForceFromCylinder3D< T, S >::_radParticle, olb::MagneticForceFromCylinder3D< T, S >::_radWire, M_PI, and olb::util::pow().

Here is the call graph for this function:

Member Function Documentation

operator()()

template<typename T , typename S >

bool olb::MagneticForceFromCylinder3D< T, S >::operator() ( T output[], const S x[] )

overridevirtual

operator writes the magnetic force in a point x util::round a cylindrical wire into output field

Implements olb::GenericF< T, S >.

Definition at line 1180 of file frameChangeF3D.hh.

{
 
  Vector<T, 3> magneticForcePolar;
  T outputTmp[3] = {T(), T(), T()};
  Vector<T, 3> normalAxis(_car2cyl.getAxisDirection());
  normalAxis = normalize(normalAxis);
 
  Vector<T, 3> relPosition;
  relPosition[0] = (x[0] - _car2cyl.getCartesianOrigin()[0]);
  relPosition[1] = (x[1] - _car2cyl.getCartesianOrigin()[1]);
  relPosition[2] = (x[2] - _car2cyl.getCartesianOrigin()[2]);
 
  T tmp[3] = {T(),T(),T()};
  _car2cyl(tmp,x);
  T rad = tmp[0];
  T phi = tmp[1];
 
  T test = relPosition * normalAxis;
 
  if ( (rad > _radWire && rad <= _cutoff*_radWire) &&
       (T(0) <= test && test <= _length) ) {
 
    magneticForcePolar[0] = _factor/util::pow(rad, 3)
                            *(util::pow(_radWire, 2)/util::pow(rad, 2) + util::cos(2*phi));
    magneticForcePolar[1] = _factor/util::pow(rad, 3)*util::sin(2*phi);
 
    // changes radial and angular force components to cartesian components.
    outputTmp[0] = magneticForcePolar[0]*util::cos(phi)
                   - magneticForcePolar[1]*util::sin(phi);
    outputTmp[1] = magneticForcePolar[0]*util::sin(phi)
                   + magneticForcePolar[1]*util::cos(phi);
 
    // if not in standard axis direction
    if ( !( util::nearZero(normalAxis[0]) && util::nearZero(normalAxis[1]) && util::nearZero(normalAxis[2]-1) ) ) {
      Vector<T, 3> e3(T(), T(), T(1));
      Vector<T, 3> orientation =
        crossProduct3D(Vector<T,3>(_car2cyl.getAxisDirection()),e3);
 
      AngleBetweenVectors3D<T,S> angle(util::fromVector3(e3), util::fromVector3(orientation));
      T alpha[1] = {T()};
      T tmp2[3] = {T(),T(),T()};
      for (int i = 0; i<3; ++i) {
        tmp2[i] = _car2cyl.getAxisDirection()[i];
      }
      angle(alpha,tmp2);
 
      std::vector<T> origin(3, T());
      RotationRoundAxis3D<T, S> rotRAxis(origin, util::fromVector3(orientation), alpha[0]);
      rotRAxis(output,outputTmp);
    }
    else {
      output[0] = outputTmp[0];
      output[1] = outputTmp[1];
      output[2] = T();
    }
  }
  else {
    output[0] = outputTmp[0];
    output[1] = outputTmp[1];
    output[2] = outputTmp[1];
  }
  return true;
}

References olb::util::cos(), olb::crossProduct3D(), olb::util::fromVector3(), olb::util::nearZero(), olb::normalize(), olb::util::pow(), and olb::util::sin().

Here is the call graph for this function:

Member Data Documentation

_car2cyl

template<typename T , typename S >

CartesianToCylinder3D<T,S>& olb::MagneticForceFromCylinder3D< T, S >::_car2cyl

protected

Definition at line 500 of file frameChangeF3D.h.

_cutoff

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_cutoff

protected

maximal distance from wire cutoff/threshold

Definition at line 506 of file frameChangeF3D.h.

_factor

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_factor

protected

factor = mu0*4/3.*PI*radParticle^3*_Mp*radWire^2/r^3

Definition at line 516 of file frameChangeF3D.h.

_length

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_length

protected

length of the wire, from origin to _car2cyl.axisDirection

Definition at line 502 of file frameChangeF3D.h.

_Mp

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_Mp

protected

magnetization particle, linear scaling factor

Definition at line 510 of file frameChangeF3D.h.

_mu0

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_mu0

protected

permeability constant, linear scaling factor

Definition at line 514 of file frameChangeF3D.h.

_Mw

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_Mw

protected

saturation magnetization wire, linear scaling factor

Definition at line 508 of file frameChangeF3D.h.

_radParticle

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_radParticle

protected

particle radius, cubic scaling factor

Definition at line 512 of file frameChangeF3D.h.

_radWire

template<typename T , typename S >

T olb::MagneticForceFromCylinder3D< T, S >::_radWire

protected

wire radius

Definition at line 504 of file frameChangeF3D.h.

---

The documentation for this class was generated from the following files:

• src/functors/analytical/frameChangeF3D.h
• src/functors/analytical/frameChangeF3D.hh