CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Source terms with gradient and curl

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 31, 2018, 19:06
Unhappy Source terms with gradient and curl
New Member
Join Date: Dec 2011
Posts: 29
Rep Power: 10
msaravia is on a distinguished road
Hi, I am trying to solve an equation of the form:

\nabla^2A = - b \cdot (\nabla A-\nabla A^T)

being A and b vector fields; so, \nabla A is the Jacobian and \nabla A-\nabla A^T is a skew symmetric tensor. The equation can also be written as:

\nabla^2A = -b \times (\nabla \times A)

I am wondering wich is the best way to treat the RHS. I have been thinking of using simply

fvm::laplacian(A) == b ^ fvc::curl(A)
fvm::laplacian(A) == b & ( fvc::grad(A)-fvc::grad(A).T)
I know that positive sources should be treated implicitly, but the RHS is no actually a source. Also, the curl does not have an implicit implementation in OpenFOAM.

At last... I'm lost. Any suggestions to implement this equation would be welcome!


Last edited by msaravia; July 31, 2018 at 20:34.
msaravia is offline   Reply With Quote

Old   February 28, 2020, 17:00
Default any progress
Thomas Flint
Join Date: Jan 2016
Posts: 39
Rep Power: 6
tom_flint2012 is on a distinguished road

I'm implementing a similar equation. I thought I had some success converting one of my curls using identities but cant get rid of another.

Have you made any progress with this in openfoam?


tom_flint2012 is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
access to velocity gradient for Lagrangian particles jiejie OpenFOAM 31 December 2, 2016 04:56
I can not find voricity or curl variables in post processing raminostadi FLUENT 2 September 14, 2016 03:56
Using Curl in CFX11 Cardero CFX 0 July 21, 2009 10:29
Curl Velocity Cardero CFX 0 March 5, 2009 10:24
vortex stretching saba Main CFD Forum 10 December 2, 2002 05:03

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