CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Rapid increase of velocity close to the boundary (MHD simulation)

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   May 11, 2017, 07:12
Default Rapid increase of velocity close to the boundary (MHD simulation)
New Member
Roman G.
Join Date: Apr 2017
Posts: 9
Rep Power: 2
Novel is on a distinguished road
I'm trying to simulate a Taylor instability by using the electric potential formulation of the MHD equations (similar to ).

The domain is a cuboid (96x96x120 mm) where the patches inlet,outlet are the sides with the size 96x96 and all other sider are defined as wall. All boundaries are defined as closed domain.
U -> noSlip
p -> zeroGradient
Epot -> zerogradient at the wall, and fixedGradient 0.5... at the inlet at fixedValue 0 at the outlet. This corresponds to a current of around 18kA.
B -> fixedValue where the value is defined with the analytic solution.

The lorentz force is defined as:
lorentz = (J0 ^  B);
the force is added in the momentum equation the following way:
tmp<fvVectorMatrix> tUEqn
    fvm::ddt(U) + fvm::div(phi, U)
  + MRF.DDt(U)
  + turbulence->divDevReff(U)
  - (1.0/rho) * (lorentz)
The simulation runs and the velocity field looks just fine but close to the boundary the velocity increases rapidly and than falls down to zero due to noSlip. I'm wondering if I should add the force term not in the momentum equation or if I need to use other boundary conditions but which one?

In the link bellow I uploaded some pictures where I measured the velocity, pressure and lorentz force along the y-axis.

Maybe one of you had a similar problem and has some ideas what I could try.


Novel is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
3D Windturbine simulation in SU2 k.vimalakanthan SU2 13 June 20, 2016 05:27
Radiation interface hinca CFX 15 January 26, 2014 18:11
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
Low Mixing time Problem Mavier CFX 5 April 29, 2013 00:00
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45

All times are GMT -4. The time now is 19:18.