|
[Sponsors] |
icoUncoupledKinematicParcelFoam problem with injection |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 14, 2017, 07:53 |
icoUncoupledKinematicParcelFoam problem with injection
|
#1 |
New Member
Adam Katz
Join Date: Feb 2017
Posts: 23
Rep Power: 9 |
Hi all,
This is my first LPT and I can't seem to figure it out. It is quite frustrating, help will be much appreciated. So far I ran simpleFoam on my model, then I run icoUncoupledKinematicParcelFoam but it doesn't time-loop (and doesn't create lagrangian particles). The output I get is: Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 402 Reading g Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type laminar Constructing kinematicCloud kinematicCloud Constructing particle forces Selecting particle force sphereDrag Constructing cloud functions none Constructing particle injection models Creating injector: model1 Selecting injection model patchInjection Constructing 3-D injection Choosing nParticle to be a fixed value, massTotal variable now does not determine anything. Selecting distribution model fixedValue Selecting dispersion model none Selecting patch interaction model none Selecting stochastic collision model none Selecting surface film model none Selecting U integration scheme Euler Selecting collision model none Reading non-inertial frame fields Starting time loop End controlDict Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application icoUncoupledKinematicParcelFoam; startFrom latestTime; startTime 0; stopAt endTime; endTime 1; deltaT 1e-3; writeControl timeStep; writeInterval 1; purgeWrite 0; writeFormat binary; writePrecision 8; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable true; functions { #includeFunc residuals } // ************************************************************************* // Code:
subModels { particleForces { sphereDrag; } injectionModels { model1 { type patchInjection; parcelBasisType fixed; patchName Inlet; nParticle 1; parcelsPerSecond 300000; massTotal 0; U0 (0 4.55 0); flowRateProfile constant 1; SOI 0; duration 1; // NOTE: set to 1 for steady state sizeDistribution { type fixedValue; fixedValueDistribution { value 15.5e-6; } } } } dispersionModel none; patchInteractionModel none; surfaceFilmModel none; stochasticCollisionModel none; collisionModel none; } Thanks a lot |
|
March 14, 2017, 13:26 |
|
#2 |
Member
Joshua
Join Date: Dec 2016
Location: St. Louis, Missouri
Posts: 91
Rep Power: 10 |
Adamkat,
It does appear to be both files. Your controlDict says that the simulation ends at time = 1. When creating your mesh (and taking your most recent time) the simulation pulls time = 402 (from your log). As 402 > 1, the code reads this as your simulation is over. Change your end time from 1 to 403, and your SOI for your injector from 0 to 402. This should have your Lagrangian solver run from 402 to 403. Hope this helps, Joshua |
|
March 15, 2017, 05:17 |
|
#3 |
New Member
Adam Katz
Join Date: Feb 2017
Posts: 23
Rep Power: 9 |
It worked!
Thanks Joshua! I had to do some modifications to make it really work like adding wall rebound (see code below). It is still missing one part, rotation. The particles must have rotation but for some reason they do not rotate. What am I doing wrong? I see this when I try to plot angular momentum (all zero) . It can also be seen in the runtime output: Code:
Time = 402.00998 Evolving kinematicCloud Solving 3-D cloud kinematicCloud Cloud: kinematicCloud Current number of parcels = 28 Current mass in system = 5.2629444e-11 Linear momentum = (-8.1501672e-11 -1.4583216e-11 -8.5672026e-13) |Linear momentum| = 8.2800524e-11 Linear kinetic energy = 2.8511479e-10 model1: number of parcels added = 297 mass introduced = 5.5824803e-10 Parcel fate (number, mass) : patch Walls - escape = 0, 0 - stick = 0, 0 Parcel fate (number, mass) : patch Inlet - escape = 0, 0 - stick = 0, 0 Parcel fate (number, mass) : patch Outlet - escape = 269, 5.0561818e-10 - stick = 0, 0 Rotational kinetic energy = 0 ExecutionTime = 2.65 s ClockTime = 3 s Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solution { active true; coupled false; transient yes; cellValueSourceCorrection off; maxCo 0.3; sourceTerms { schemes { } } interpolationSchemes { rho cell; U cellPoint; mu cell; curlUcDt cellPoint; } integrationSchemes { U Euler; } } constantProperties { rho0 964; youngsModulus 6e8; poissonsRatio 0.35; } subModels { particleForces { sphereDrag; SaffmanMeiLiftForce { U U; } } injectionModels { model1 { type patchInjection; parcelBasisType fixed; patchName Inlet; nParticle 1; parcelsPerSecond 30000; massTotal 0; U0 (0 4.55 0); flowRateProfile constant 1; SOI 402; duration 1; // NOTE: set to 1 for steady state sizeDistribution { type fixedValue; fixedValueDistribution { value 15.5e-6; } } } } dispersionModel none; patchInteractionModel none; surfaceFilmModel none; stochasticCollisionModel none; collisionModel none; patchInteractionModel localInteraction; localInteractionCoeffs { patches ( Walls { type rebound; e 0.97; mu 0.09; } Inlet { type rebound; e 0.97; mu 0.09; } Outlet { type escape; } ); } StandardWallInteractionCoeffs { type rebound; e 0.97; mu 0.09; } } cloudFunctions {} // ************************************************************************* // |
|
Tags |
particle tracing, particle tracking |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
problem with Injection | Mark_89 | STAR-CCM+ | 1 | December 19, 2013 06:20 |
problem with DPM injection | Hayder Mohammed | Main CFD Forum | 0 | August 12, 2013 03:01 |
jou file command for the injection! | asal | FLUENT | 0 | June 4, 2013 05:15 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 20:13 |
DPM injection problem | akr | FLUENT | 5 | November 6, 2006 11:34 |