|June 27, 2014, 04:37||
[Lagrangian] Injection Model Start of Injection timestep Mistake?
Join Date: Dec 2013
Posts: 13Rep Power: 5
Looking at the lagrangian injection Model : InjectionModel.C,
I find in the membre function PrepareForTheNextTimestep something I don't Understand :
// Make times relative to SOI
scalar t0 = timeStep0_ - SOI_;
scalar t1 = time - SOI_;
In my mind, because timestep0_ is the last time step, so the first delta Time between current time and SOI is miss, indeed function parcelsToinject require t0 > 0
so in my mind,
if SOI = 0.5
at first step time =1 (dt=1) timeStep_0=0
so t0=-0.5 so no parcel to inject ... [But we want to inject parcel, no?]
at second step time time =2 timestep_0 = 1
so t0 = 0.5 and t1 = 1.5 so inject parcel give us t1-t0 = dt =1
But we have forget the first par of injection from SOI = 0.5 to time =1
Is it right? Where is my missunderstanding?
If someone can enlighten me, it will be great!
Have a nice day!
|October 29, 2014, 17:51||
Join Date: Oct 2014
Posts: 4Rep Power: 4
Has anybody looked into this? InjectionModel seems good, but the derived injection model's parcelsToInject disallows time0<0 as you noted Emeline.
e.g. I am using coneNozzleInjection with SOI = 10.010. See log file, the first injection time step 10.0101 is missed...!? May be a bug? (albeit minor I suppose!)
|November 5, 2014, 21:53||
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 320Rep Power: 10
The answer to your question should lie in the call to parcelsToInject function in one of the InjectionModel derived classes as opposed to looking at the base class.
For example, take a look at the PatchFlowRateInjection derived class.
You'll see that the parcesToInject() function won't necessarily return a positive value even if t>SOI. Here it introduces some randomness to the injection process while still meeting the concentration criteria you set in the kinematicCloud dictionary. My understanding is the purpose is to remove a timestep dependency to the injection rates. I actually think its a pretty clever implementation.
I hope that helps!
|Thread||Thread Starter||Forum||Replies||Last Post|
|problem with solving lagrange reaction cloud||Polli||OpenFOAM Running, Solving & CFD||0||April 30, 2014 07:53|
|Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam)||cfdonline2mohsen||OpenFOAM||3||October 21, 2013 09:28|
|An error has occurred in cfx5solve:||volo87||CFX||5||June 14, 2013 17:44|
|How to start a room model with heat object||steveleung||Main CFD Forum||0||February 24, 2011 04:38|
|about group injection of DPM model||zwdi||FLUENT||5||April 12, 2004 19:31|