CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

mapped boundary condition in multi-region solvers

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 13, 2019, 07:09
Default mapped boundary condition in multi-region solvers
  #1
New Member
 
Benedetto Di Paolo
Join Date: Feb 2016
Location: Santander (Spain)
Posts: 4
Rep Power: 10
benedip is on a distinguished road
Hi,


I’m trying to exchange information between two boundaries (mappedWall types) in a multi-region approach. I built separately the two regions (domain1 and domain2) by means of two blockMeshDicts. A similar approach can be found in the following tutorial case:


OpenFOAM-v1812/tutorials/heatTransfer/chtMultiRegionSimpleFoam/heatExchanger


I transferred the velocity field U from patch1 (mappedWall at the outlet of domain1) to patch2 (mappedWall at the inlet of domain2). I achieved it easily this by mapping velocity on patch1 and by fixing it on patch2 in a similar way to the mappedVelocityFluxFixedValue B.C. So B.C. ends like:

Code:
  operator==(patch1UField);
  
 fixedValueFvPatchField<vector>::updateCoeffs();
Thus, I obtained patch1UField = patch2UField. So far no problems.


Now, If I look at the velocity field just a little bit far from patch2 (in domain2), I find completely different values from the ones I fix (patch1UField). It’s like if my value didn’t come out from patch2.


Finally I can obtain the expected values at the neighboring patch cells (patch().faceCells()) by applying this B.C. (by fixing the patch values and the internal cells values), as follows:

Code:
  operator==(patch1UField);
  
 Field<vector>& intFld = const_cast<Field<vector>&>(this->primitiveField());
 UIndirectList<vector>(intFld, this->patch().faceCells()) = patch1UField;

  
fixedValueFvPatchField<vector>::updateCoeffs();
Does anybody know how to fix this problem?
I just want to fix the velocity on the patch without pushing internal cells values.
Thanks you


Best regards,
Benedetto
p { margin-bottom: 0.1in; line-height: 120%; }a:link { } p { margin-bottom: 0.1in; line-height: 120%; }a:link { }

Last edited by benedip; February 13, 2019 at 11:45.
benedip is offline   Reply With Quote

Reply

Tags
fixed value, fixed velocity, mappedboundarycondition, multi regions, openfoam 1806

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
Mapped Boundary Condition louistse OpenFOAM 5 June 1, 2018 08:23
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 08:44
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
[snappyHexMesh] Multi region meshing & recovering the original patch names fluidpath OpenFOAM Meshing & Mesh Conversion 4 May 19, 2013 19:13
External Radiation Boundary Condition (Two sided wall), Grid Interface CFD XUE FLUENT 0 July 8, 2010 06:49


All times are GMT -4. The time now is 17:55.