I am trying to test the FORTRAN routine pt_injection.F from the /examples folder of the CFX installation user routine for injecting particles into a domain using a user FORTRAN routine to set the particle location, velocity, diameter, temperature, and mass flow rate. I keep getting errors when I try to run this case and would like help getting this routine to run.
|
There are a number of things to consider if you use this routine as is. 1. Make sure that the order of the return variables listed in the CCL is the same as that assumed in the FORTRAN. The FORTRAN assumes the following order: Particle Position, Velocity, Diameter, Temperature, and Particle Mass Flow Rate. If the CCL order differs from this, you can change it by editing the Particle Injection region object in the command editor. 2. Even though the Particle Mass Flow Rate is defined by the user FORTRAN routine, you still need to set a particle mass flow rate in the CCL. Tests show that setting a value of 0 kg/s for the mass flow rate in the GUI will work. The actual mass flow rate will then be set by the routine. 3. The pt_injection.F routine included in the /examples folder assumes that the injected particles are uniformly distributed from 4 injection locations. When you set the number of injection points, you should make sure that the number is divisible by 4, or you should modify the FORTRAN accordingly. The CCL for the Particle Injection Routine should look something like what's shown below (note the order of the return variables list matches the FORTRAN, the number of locations is evenly divisible by 4, and the particle mass flow rate is set to 0 kg/s): PARTICLE INJECTION REGION: Particle Injection Region 1 FLUID: Water INJECTION CONDITIONS: INJECTION METHOD: Option = User Defined Particle User Routine = injectthem Particle Volume Injection Return Variables List = Particle Position,Velocity,Mean Particle Diameter,Temperature,Particle Mass Flow Rate NUMBER OF POSITIONS: Number = 40 Option = Direct Specification END END PARTICLE MASS FLOW RATE: Mass Flow Rate = 0.0 [kg s^-1] END END END END END A sample definition file and FORTRAN file are attached. |
||
|