Good morning,
I'm sorry for
Good morning,
I'm sorry for disturbing during coffee break (in France, time is coming). But I'm still a little bit confused. In the articles I've read about LES on OpenFOAM, to get real turbulent inlet, we used to define a "mapping surface". I agree. Here comes my first question, to be sure I didn't make any mistakes here, I think that it is defined in the constant/polymesh/boundary files in the inlet by changing the offset (if the inlet is directMappedPatch) isn't it ? so, doing this, we get a periodic volume between inlet patch and the Mapping surface right ? My mask function should be 1 in this volume and 0 elsewhere. by writing these 2 lines I will never get the volume but only the cells which are in my patch "inlet": label patchINLET = mesh.boundaryMesh().findPatchID("inlet"); toto.boundaryField()[patchINLET]=1.0; So my question is is there a way to take all these cells and to impose "toto" of all these cell =1 (not only boundary patch) Thank you again, Regards, Cedric |
hi to all,
I m a beginner i
hi to all,
I m a beginner in OpenFoam and I want to program fluctuations with tensor... for example I want the following program to do with fluctuations. magSqr(fvc::grad(U))+ sqr(gradU.component(tensor::YX) + gradU.component(tensor::XY)) does someone can help me? best regards |
hi to all,
I m a beginner
hi to all,
I m a beginner in OpenFoam and I want to program fluctuations with tensor... for example I want the following program to do with fluctuations. magSqr(fvc::grad(U))+ sqr(gradU.component(tensor::YX) + gradU.component(tensor::XY)) does someone can help me? best regards |
Hello,
this is a question r
Hello,
this is a question regarding the DirectMappedFixedValue Boundary Patch. I would like to use the patch for generating a nice developed flow as an inlet condition (LES). Unfortunally I am working with the OpenFOAM Version 1.3. Is there a version of this b.c. for OF 1.3? Trying to figure out the structure of this boundary condition brought me fast to the directMapped file that allows to transfer the needed inlet values to the boundary patch. This class is refering to PStream. I think, in this file, there are some big changes from 1.3 to 1.5. How deep are the changes going? Is it possible to include this boundary condition to 1.3? Perhaps somebody can help me with an other idea to generate an inlet profile with a small number of cells as it can be done by this b.c. I belive that a pipelength of 4d (2d for recirculation and 2d to take into account the interaction with the pipe exit) should be enough. Thank you so fare Wilko |
Hi,
I'm working on Klein i
Hi,
I'm working on Klein inflow condition which is an effective algebraic Methods. It is used to target the statistics of turbulence by performing operation on them. Now, I'm implementing it directly on a LES solver. I've to define the Cholesky decomposotion "aij" of the Reynolds stress tensor. The expression of aij is not so hard to implement but, while compiling, I still have a error: error: no matching function for call to âFoam::dimensioned<foam::tensor<double> >::dimensioned(const char [5], Foam::dimensionSet)â /home/users/duprat/OpenFOAM/OpenFOAM-1.4.1/src/OpenFOAM/lnInclude/dimensionedTyp e.C:95: note: candidates are: ... I create aij as a volTensorField: volTensorField aij ( IOobject ( "aij", runTime.constant(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh, dimensionedTensor("null",dimensionSet(0, 1, -1, 0, 0, 0, 0)) ); And in these last lines is the mistake. Then, I write : aij[cellI] = (....) where cellI is the mesh at the inlet patch Can someone give me a hint to solve this problem ? I'll use this tensor every time step but it is constant in time that's why it is defined in constant directory Thank you for helping, Cedric PS: aij tensor is not a symetric tensor and I need it only in the inlet patch |
I got it
dimensionedTenso
I got it
dimensionedTensor("null",dimensionSet(0, 1, -1, 0, 0, 0, 0)) should be dimensionedTensor("null",dimensionSet(0, 1, -1, 0, 0, 0, 0),tensor::zero) ..... initial value, everything is there !! I hope I can give to the community this solver but .... I've still some work to do :-) Cedric |
Hi,
for the implementation
Hi,
for the implementation of a synthetic inlet of L. Davidson one needs random numbers with the probability distribution of p(phi)=1/2*Pi with 0<=phi<=2*Pi. Does anyone know, how to achieve this? So far I know of: Random perturbation(1234567); scalar phi_n=phiValues.GaussNormal(); Would be nice, if you can give me a hint! Fabian |
have you tried this?
some p
have you tried this?
some pdf's can be found in thermophysicalModels Random rndGen(0); autoPtr<pdf> phiPDF ( pdf::New ( <dictionary> rndGen ) ); scalar phi = phiPDF->sample(); the pdf is defined in the the <dictionary> of which you can find examples in sprayProperties where I use them. N |
Hi Niklas,
thansk for the h
Hi Niklas,
thansk for the hint! Fabian |
Hello Cedric
I am working on the break up study of liquid sheets. I am using lesInterfoam for this purpose. I am looking into start working on the implementation of Kleins initial conditions for turbulence. I saw your post and thought as you have already worked on it you could give me some hints on how to start along with it or just explain me the steps you followed to implement it in the lessolver. Any help is appreciated. bye with regards K.Suresh kumar |
1 Attachment(s)
Quote:
Dear eugene, It seems that OpenFOAM only can simply map the internal value back to inlet. If I want to use mapping position value multiply an function and then feed back to inlet, is it possible to realise ? < > is time-averaged value, is boundary layer thickness recy is mapping plane position |
I see no reason why modifying the mapped value should not be possible. In fact several people have already done this. Try to contact Gavin Tabor at Exeter Uni. I think he has some code for this. (I do not unfortunately.)
|
Dear all,
It seems that pitzDailyDirectMapped case can't be run in version-1.6. When I run, I have the following mistakes: patch type 'patch' not type 'directMappedPatchBase' for patch inlet of field U in file "/home/panda60/OpenFOAM/run/tutorials/incompressible/pisoFoam/les/pitzDailyDirectMapped/0/U" From function directMappedFixedValueFvPatchField<Type>::directMa ppedFixedValueFvPatchField ( const fvPatch& p, const DimensionedField<Type, volMesh>& iF, const dictionary& dict ) in file fields/fvPatchFields/derived/directMappedFixedValue/directMappedFixedValueFvPatchField.C at line 109. FOAM exiting who knows how to use DirectMapped condition crectly ? Thank you very much. |
Inflow boundary question
hey there foamers. i have been using openfoam for about 6months now, and have a question. if anyone can assist please do so, it would be greatly appreciated. i am trying to implement the turbulentinlet boundary for my simulation, but i am not able to find any documentation on how the following quantities are chosen. can anyone explain to me? :confused:
boundaryField { inlet { type turbulentInlet; referenceField uniform (10 0 0); fluctuationScale (0.02 0.01 0.01); value uniform (10 0 0); |
Quote:
1. blockMesh 2. changeDictionary 3. pisoFoam. This should do the trick. Regards, Francois. |
Quote:
Regards, Francois. |
Inlet velocity Profile
Hi,
I am using parabolic inlet velocity profile (solved analytically) for shear-driven multiphase flow using InterFoam solver. I could observe numerical diffusion due to the analytical velocity profile which is imposed at inlet. Please do let me know how to get interpolated discretized velocity profile from internal field on to inlet boundary. Thanks, Jagannath |
directMapped boundary conditions
Good morning to the OpenFOAM community!
I would like to use the directMapped method in the LES of the pitzDaily3D testcase without. In order to do that after having previously read the suggestions contained in the former threads, I'have followed these steps:
Code:
[CODE]
Code:
Anyway I'm not sure to know the correct meaning of the keyword sampleMode contained in the dictionary. Could someone give any further information about the sampleMode nearestCell keywords ? Regards, Andrea Aprovitola |
Hello Andrea,
A little late, but better late than never! To answer your questions:
I hope this helps! Kind regards, Francois. |
Dear Francois,
thank you very much for your help. Now I'm confident about the meaning of the offset value. Anyway, only the sample mode "nearestCell" seems working for my application. NEARESTFACE and NEARESTPATCHFACE model are not working. I'm not sure of what I'm doing wrong, Regards Andrea Code:
Quote:
|
All times are GMT -4. The time now is 06:19. |