heaviside?

ginette · February 25, 2011, 14:41

Dear Open Foam Community,

I am a brand new to OpenFoam.

I would like to create a field like

Code:

volScalarField toto = toto1 * heaviside(phi) + toto2 * heaviside(-phi);

where toto1 and toto2 are two "constants" and phi is a "volScalarField". Is there a way to do this?

Many thanks in advance,
Martin.

ngj · February 25, 2011, 15:00

Hi Martine

You could use either of the functions Foam::neg or Foam:

os. These functions return 1 if the argument is negative (positive) and 0 if the argument is positive (negative). I do not know how the functions handles 0, however.

Best regards,

Niels

ginette · February 25, 2011, 15:24

Indeed!

Thanks so much, Niels!

Martin.

PS: Is there a way to close the thread and tag is "solved"?

February 25, 2011, 14:41	heaviside?	#1
ginette New Member Martine Ginette Join Date: Feb 2011 Posts: 5 Rep Power: 0	Dear Open Foam Community, I am a brand new to OpenFoam. I would like to create a field like Code: volScalarField toto = toto1 * heaviside(phi) + toto2 * heaviside(-phi); where toto1 and toto2 are two "constants" and phi is a "volScalarField". Is there a way to do this? Many thanks in advance, Martin.

February 25, 2011, 15:00		#2
ngj Senior Member Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Copenhagen, Denmark Posts: 1,903 Rep Power: 37	Hi Martine You could use either of the functions Foam::neg or Foam:os. These functions return 1 if the argument is negative (positive) and 0 if the argument is positive (negative). I do not know how the functions handles 0, however. Best regards, Niels

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Heaviside Function	MCDG	CFX	1	December 3, 2007 04:32
smoothed interface capturing method: compressible.	jinwon	Main CFD Forum	0	May 14, 2007 01:20
Heaviside Function?	Zhou Hua	Main CFD Forum	2	December 8, 2003 06:10

February 25, 2011, 15:24		#3
ginette New Member Martine Ginette Join Date: Feb 2011 Posts: 5 Rep Power: 0	Indeed! Thanks so much, Niels! Martin. PS: Is there a way to close the thread and tag is "solved"?