Solve Ax=b

Rutger · August 8, 2015, 05:19

Hi,

I searched the forum but couldn't find a similar problem respectively an answer.
I'm running a simulation that stores its results in a volScalarField b ∈ R^n and a volVectorField A ∈ R^{nxn}.
The mathematically relation between them is A*x = b, which I want to solve for x. In another project I used a solver from the Eigen libary but now I want to calculate it in OpenFOAM.

It is possible to solve this linear equation or do I need to use a fvMatrix type for A ?

Rutger · August 10, 2015, 04:44

For everyone with the same problem:

The example in application/test/Matrix/Test-Matrix.C helped a lot. Now I'm copying the values from the vol*Field to the corresponding datatypes and solve the system.

Tushar@cfd · August 10, 2015, 09:13

Dear Rutger,

Have a look to the eigen value solver given in Tensor.C file. I think this will help you.

-
Best Luck!

Rutger · August 13, 2015, 07:26

Thanks, I will have a look at it.

August 8, 2015, 05:19	Solve Ax=b	#1
Rutger New Member Join Date: Aug 2015 Posts: 4 Rep Power: 10	Hi, I searched the forum but couldn't find a similar problem respectively an answer. I'm running a simulation that stores its results in a volScalarField b ∈ R^n and a volVectorField A ∈ R^{nxn}. The mathematically relation between them is A*x = b, which I want to solve for x. In another project I used a solver from the Eigen libary but now I want to calculate it in OpenFOAM. It is possible to solve this linear equation or do I need to use a fvMatrix type for A ?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
solve equation	michaelsmit	OpenFOAM Running, Solving & CFD	4	March 24, 2011 05:35
Solve Flow or VOF simultaneously ?	Ramsey	FLUENT	1	February 16, 2011 13:16
Linearized NS euqations: how to solve them?(problem with Matrix operations..)	matteoL	OpenFOAM Running, Solving & CFD	0	November 18, 2009 06:58
Using Compressible Solver (sonicFoam) to solve subsonic flows	ezsoal	OpenFOAM	0	October 27, 2009 09:13
MAC method solve for pressure at boundry	Ron	Main CFD Forum	2	January 28, 2006 02:14

August 10, 2015, 04:44		#2
Rutger New Member Join Date: Aug 2015 Posts: 4 Rep Power: 10	For everyone with the same problem: The example in application/test/Matrix/Test-Matrix.C helped a lot. Now I'm copying the values from the vol*Field to the corresponding datatypes and solve the system.

August 10, 2015, 09:13		#3
Tushar@cfd Senior Member T. Chourushi Join Date: Jul 2009 Posts: 321 Blog Entries: 1 Rep Power: 17	Dear Rutger, Have a look to the eigen value solver given in Tensor.C file. I think this will help you. - Best Luck!

August 13, 2015, 07:26		#4
Rutger New Member Join Date: Aug 2015 Posts: 4 Rep Power: 10	Thanks, I will have a look at it.