[FlyBob91]

Hello to all,
I'm going to add the sensitivity equation to my createFields.H

the equation is

-A*nu*(n*nabla)u*(n*nabla)u

I know how define all the terms because i included them in the outlet boundary files, however i think i missed some header line to be included.

At the end of the file i added

 
    scalar nu = readScalar(db().lookupObject<IOdictionary>("transportProperties").lookup  
    ("nu"));

volScalarField sens
(
    IOobject 
    (
        "sensitivity"
        "alpha",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    -magSf()*nu*snGrad(Ua)*snGrad(U)
);

and this is the error

In file included from /home/roby/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/postProcess.H:129:0,
                 from myAdjointShapeOptimizationFoam.C:72:
./createFields.H: In function ‘int main(int, char**)’:
./createFields.H:112:34: error: ‘patch’ was not declared in this scope
     vectorField U_n = (U & patch().nf())*patch().nf();  // Normale
                                  ^
./createFields.H:118:31: error: ‘db’ was not declared in this scope
     scalar nu = readScalar(db().lookupObject<IOdictionary>("transportProperties
                               ^
./createFields.H:118:58: error: expected primary-expression before ‘>’ token
     scalar nu = readScalar(db().lookupObject<IOdictionary>("transportProperties
                                                          ^
./createFields.H:118:83: error: request for member ‘lookup’ in ‘("transportProperties")’, which is of non-class type ‘const char [20]’
 nu = readScalar(db().lookupObject<IOdictionary>("transportProperties").lookup  
                                                                        ^
./createFields.H:132:12: error: ‘magSf’ was not declared in this scope
     -magSf()*nu*snGrad(Ua_t)*snGrad(U_t)
            ^
./createFields.H:132:28: error: ‘snGrad’ was not declared in this scope
     -magSf()*nu*snGrad(Ua_t)*snGrad(U_t)
                            ^
./createFields.H:132:28: note: suggested alternative:
In file included from /home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvcSnGrad.H:88:0,
                 from /home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvc.H:51,
                 from /home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvCFD.H:8,
                 from myAdjointShapeOptimizationFoam.C:48:
/home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvcSnGrad.C:88:1: note:   ‘Foam::fvc::snGrad’
 snGrad
 ^
In file included from myAdjointShapeOptimizationFoam.C:78:0:
createFields.H:112:34: error: ‘patch’ was not declared in this scope
     vectorField U_n = (U & patch().nf())*patch().nf();  // Normale
                                  ^
createFields.H:118:31: error: ‘db’ was not declared in this scope
     scalar nu = readScalar(db().lookupObject<IOdictionary>("transportProperties
                               ^
createFields.H:118:58: error: expected primary-expression before ‘>’ token
     scalar nu = readScalar(db().lookupObject<IOdictionary>("transportProperties
                                                          ^
createFields.H:118:83: error: request for member ‘lookup’ in ‘("transportProperties")’, which is of non-class type ‘const char [20]’
 nu = readScalar(db().lookupObject<IOdictionary>("transportProperties").lookup  
                                                                        ^
createFields.H:132:12: error: ‘magSf’ was not declared in this scope
     -magSf()*nu*snGrad(Ua_t)*snGrad(U_t)
            ^
createFields.H:132:28: error: ‘snGrad’ was not declared in this scope
     -magSf()*nu*snGrad(Ua_t)*snGrad(U_t)
                            ^
createFields.H:132:28: note: suggested alternative:
In file included from /home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvcSnGrad.H:88:0,
                 from /home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvc.H:51,
                 from /home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvCFD.H:8,
                 from myAdjointShapeOptimizationFoam.C:48:
/home/roby/OpenFOAM/OpenFOAM-4.1/src/finiteVolume/lnInclude/fvcSnGrad.C:88:1: note:   ‘Foam::fvc::snGrad’
 snGrad
 ^
/home/roby/OpenFOAM/OpenFOAM-4.1/wmake/rules/General/transform:8: set di istruzioni per l'obiettivo "Make/linux64GccDPInt32Opt/myAdjointShapeOptimizationFoam.o" non riuscito
make: *** [Make/linux64GccDPInt32Opt/myAdjointShapeOptimizationFoam.o] Errore 1

[FlyBob91]

Ok i added these line above runTime.write(); in the shapeOptimization.C file

// Sensitivity calculation
word patchName = "wing";
label patchID = mesh.boundary().findPatchID(patchName);
scalarField temp = (U.boundaryField()[patchID].snGrad() & 
Ua.boundaryField()[patchID].snGrad());
sensitivity.boundaryField()[patchID] =
mesh.boundary()[patchID].magSf() * (temp);

but now the problem is how correctly to initialize it in the createFields.H
I added this lines at the end of the file

volScalarField sensitivity
(
    IOobject
    (
        "sensitivity",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    ******************
)

but i don't know what to write where i used the stars.
Any idea?

Best Regards,
Roberto

[FlyBob91]

Sorry for the bump, i'm a bit disperate

[verboomj]

Check chapter 3.4: http://www.tfd.chalmers.se/~hani/kur...ortAdjoint.pdf

[FlyBob91]

Originally Posted by verboomj

Check chapter 3.4: http://www.tfd.chalmers.se/~hani/kur...ortAdjoint.pdf

Hi verboomj,
I have already read that article but it was a very simple case compared to mine.
Here the solution of what i did few days ago, maybe it can help someone.

I created a file where i wrote the sensitivity equation and then i initializied it in this manner in the createFields.H

// PRIMAL NORMAL GRADIENT
volVectorField wallGradU
(
	IOobject
	(
	"wallGradU",
	runTime.timeName(),
	mesh,
	IOobject::NO_READ,
	IOobject::NO_WRITE
	),
	mesh,
	dimensionedVector
 	(
  	"wallGradU",
  	U.dimensions()/dimLength,
  	vector::zero
  	)
);


// ADJOINT NORMAL GRADIENT
volVectorField wallGradUa
(
 	IOobject
 	(
  	"wallGradUa",
  	runTime.timeName(),
  	mesh,
  	IOobject::NO_READ,
  	IOobject::NO_WRITE
  	),
 	mesh,
	dimensionedVector
 	(
  	"wallGradUa",
  	Ua.dimensions()/dimLength,
  	vector::zero
  	)
);


// FACE AREA
volScalarField areaMag
(
	IOobject
	(
	"areaMag",
	runTime.timeName(),
	mesh,
	IOobject::NO_READ,
	IOobject::AUTO_WRITE
	),
	mesh,
	dimensionedScalar("areaMag", dimLength*dimLength, 0.0)
);



// SENSITIVITY
volScalarField Sens
(
	IOobject
	(
	"Sens",
	runTime.timeName(),
	mesh,
	IOobject::NO_READ,
	IOobject::AUTO_WRITE
	),	
	mesh,
	dimensionedScalar("Sens", dimLength*dimLength*dimLength*dimLength/(dimTime*dimTime*dimTime), 0.0)
);

Best Regards,
Roberto