
[Sponsors] 
Getting inflow conditions for a simulation from another simulation 

LinkBack  Thread Tools  Display Modes 
August 29, 2012, 12:31 
Getting inflow conditions for a simulation from another simulation

#1 
New Member
Join Date: May 2009
Posts: 10
Rep Power: 8 
Hi,
I need to solve the flow in a interior domain B (for example flow around an airfoil inside a nonconstant section channel), and the inlet flow is not known, but must be extracted by the results of a computation already performed, in a domain A which is upstream of B. How do I get the inlet conditions to B, from the results of the computation of A? If the two computations were coupled, this would be easy: I would just mesh A and B so that the outlet of A corresponds to the inlet of B. Then, at each iteration, I would read the flow at A's outlet and interpolate it at B's inlet. However, in this case the solution in A has been already computed. To do this, some kind of outlet BC must have been assigned. These BCs always have some kind of "artificiality": for example, an imposed pressure BC will force uniform pressure at A outlet, which may be false in reality. In order to obtain a realistic flow field inside A, one usually extends the domain at inlet and outlet, so that BCs are not imposed too close to the zone of interest., Thus, the most realistic flow coming out from A may not be the flow exactly at outlet, but maybe the flow in a section upstream of the outlet. The same is true at B inlet. Here again, the flow domain may need to be elongated in order to assign the boundary conditions farther from the flow region of interest. What's the most correct approach in this case? Thanks, Best Regards Sergio Rossi 

August 29, 2012, 14:36 

#2 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,003
Rep Power: 27 
first: are the two simulations done with the same formulation (RANS/URANS/LES/DNS)?
second: are the outlet BCs of domain A physically congruent as the inflow BC of domain B? 

August 31, 2012, 09:24 

#3 
New Member
Join Date: May 2009
Posts: 10
Rep Power: 8 
Ciao, Filippo,
thank you very much for the answer. 1. yes, the formulation is the same, even though the code is different. In both domains, RANS equations are solved. 2. the outlet BC of A is a pressure BC (i.e., static pressure fixed, and all other quantities extrapolated from the interior), while the inlet BC of B is a radial profile for p0, T0 and flow angles. However, as I was saying, since the two computations are uncoupled, A's outlet and B's inlet don't necessarily coincide. In my specific case, A's outlet has been extended with an "unphysical" part just for convergence reasons. Luckily, B's inlet hasn't been modified, so in my case the solution should be easy: I just extract the flow profiles at the section of A which corresponds to the inlet of B, and use them as the BC of B. However, had B's inlet been extended for the same reasons, what should have I done? Some kind of "indirect" BC, e.g., iteratively change the BCs B's "numerical" inlet, so as to minimize the discrepancies between the profiles at B's physical inlet and A's physical inlet? I kind of recall seeing something like that, in a paper on LES for an underexpanded jet... 

August 31, 2012, 10:12 

#4 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,003
Rep Power: 27 
I think that using RANS is much more easier than LES (I don't think that the methods used in LES can help you much ...), you don't need to let your numerical solution to correlate and wait to overcome a numerical transient. However, my doubt remains on the fact the the outlet in A has a fixed static pressure, you should have such congruence at the inlet of B.
Depending on the formulation, you could also try to use the pressure (static and total) profiles of the outlet of A as inlet condition of B, is your case is compressible, you can let the velocity inlet to assess at the energy value fixed by the pressure ... 

August 31, 2012, 12:59 

#5 
New Member
Join Date: May 2009
Posts: 10
Rep Power: 8 
I wasn't thinking about using LES, just copying that type of strategy, i.e., impose a BC somewhere, and iteratively adjust the imposed values to match some values somewhere else.
But let's leave this aside: I don't see why the BCs at A's outlet and B's inlet should agree, since they are not even the same section geometrically. I paste a picture to be sure we're on the same page: A is the red domain, and B is the green. The geometry is not the actual one, which I cannot show, but at least it gives a reasonable idea. A has been extended to the right with a fictitious (red dotted lines) outlet, in order to impose the pressure BC farther downstream, so that the solution in the area around B's inlet would be more realistic. Note that the extension of A to the right is only done for numerical stability reasons, i.e., it doesn't coincide with the real, physical geometry, as it can be seen from the fact that the red dotted lines do not overlap B's geometry. My idea: interpolate the solution of A in the section corresponding to the inlet of B, and use that as the inlet boundary condition of B. Does it seem sensible to you? 

August 31, 2012, 12:59 

#6 
New Member
Join Date: May 2009
Posts: 10
Rep Power: 8 
Whoops! Forgot the picture


August 31, 2012, 15:54 

#7  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,003
Rep Power: 27 
Quote:
Now that I see the sketch the problem appears clear. You can do interpolation according to your idea, however bearing in mind that the solution of problem in domain A, computed at the section you want then to prescribe as inlet of domain B, depends on the specific outflow of A. The problem is not purely hyperbolic, the flow is governed by the (statistically averaged) NS equations and, in principle, the change in the outlet can affect the inlet. Of course you can do some approximations and working on ... 

August 31, 2012, 16:00 

#8 
Senior Member
Joern Beilke
Join Date: Mar 2009
Location: Dresden
Posts: 186
Rep Power: 9 
Just do it this way.
There is one thing to remember. You can not specify velocity and pressure at an inlet at the same time. So you will probably choose velocity, temperature and the turbulence parameters for mapping. 

August 31, 2012, 16:28 

#9 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,003
Rep Power: 27 
Yes, perhaps from my past experience of simulation of compressible flows, I suggest prescribing the total pressure at inlet and the temperature and let the velocity to develop, when I tried the dual setting I always get problems ...


September 6, 2012, 03:54 

#10 
New Member
Join Date: May 2009
Posts: 10
Rep Power: 8 
I agree: as a matter of fact, I'm prescribing p0, T0 and flow angles at B's inlet (interpolated from A's solution) precisely for this reason, i.e., because in compressible flow prescribing velocity creates convergence difficulties. This is expecially true when the inlet is not too far from to an obstacle (an impeller blade in this case, not shown in my sketch for semplicity): I guess the reason is that prescribing velocity in this case would lead to unphysical total pressure gradients in this case, while prescribing total quantities allows the correct velocity profile to develop.


September 6, 2012, 10:28 

#11  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,003
Rep Power: 27 
Quote:
let know what happens 

September 6, 2012, 10:42 

#12 
Member
Ren/Xingyue
Join Date: Jan 2010
Location: Nagoya , Japan
Posts: 44
Rep Power: 8 
I am wondering overlapping grid is a better choice.
Without overlapping, it is difficult to find the interaction between those two domains. 

September 8, 2012, 14:08 

#13 
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 7 
When looking at your picture, there is something that I don't understand. You have one code to perform the computation in domain A and another code for domain B. You want to use outlet of A as inlet of B. As you have computed the solution in domain A, it means you have an inlet BC for this domain, so why don't you solve your pb in the complete domain A+B and thus you won't have to deal with an inlet for B? Can you explain us? 

September 16, 2012, 07:07 

#14 
New Member
Join Date: May 2009
Posts: 10
Rep Power: 8 
Hi leflix,
I don't want to use the outlet of A as inlet of B. As you can see from the picture, the outlet of A is far downstream from the inlet of B, which is "in the middle" of A. This is done to reduce the influence of A's outlet BC on the flow field near B's inlet. The reason why A and B aren't solved together is the following: the sketch I did is only indicative, since I cannot post the real geometry on the forum. But, for people who have experience with multistage centrifugal compressors for industrial application, I can say that A is a sidestream (injection), and as such it is a domain without cyclic symmetry, which must be meshed with a unstructured mesh, and solved with an unstructured CFD code. Domain B is a standard centrifugal compressor stage, so it has cyclic symmetry, and it can be meshed with a multiblock structured mesh, and solved with a multiblock structured CFD code. Thus, two different codes are used to solve the different domains. This is done also to reduce the computational time of the simulations, the number of processors/RAM needed, and to optimize schedules in an industrial environment. By having two different engineers working on two different simulations, some degree of overlapping of their schedules can be allowed. For example, I can modify/fix CAD geometry for B, create meshes, prepare all BCs except for the inlet one, initialize my flow field and prepare postprocessing/final report templates, while the coworker working on A completes his simulation. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Boundary conditions for a fan simulation  billynoe  OpenFOAM  4  July 8, 2013 04:06 
Initial conditions for simulation  lazersos  OpenFOAM PreProcessing  0  June 1, 2011 09:47 
Boundary conditions for pipework simulation  metro  OpenFOAM Running, Solving & CFD  0  August 12, 2010 03:13 
prescribed inflow conditions  Thorsten  Main CFD Forum  2  January 29, 1999 16:03 
Need help about boundary conditions for a francis runner flow simulation  Rodrigo Escobar Moragas  Main CFD Forum  1  October 26, 1998 09:20 