Fluids

Fluids

Vapor Pressure calculation

    • anilwakale
      Subscriber

      Dear all,


      Is anyone working on the Vapor pressure model using UNIFAQ method UDF from Properties of GAS and Liquids by B.E. Poling 


      As we all know during combustion the Vapor pressure of liquids will vary with temperature and it should be accurately entered in CFD simulations at least until critical temp and pressure of the liquid. I am working on both Spray simulation using Steady flamelet model and details you can see my previous blogs. I want to know the vapor pressure at every time step or at with increment of 30-50 in both spray simulation as well as diesel unsteady model for Compression ignition engine.


      And another question regarding Droplet injection methods 


      is it DPM and Multicomponent will work similar in Fluent of not      


      Liquid _penetration for Simulation and Exp


      As we can see In both cases of Droplet and the multicomponent method I need to change The value of KH-RT constat CL from 25 to 22 for droplet to multicomponent to match ECN Spray-A validation. 

    • DrAmine
      Ansys Employee

      The only difference which comes to my mind is the interpretation of boiling mode: For droplet the boiling temperature is used whereas for multi-component the sum of saturation pressure of each components is used (compared to surrounding). 

    • anilwakale
      Subscriber

      But then is it valid to change the KH-RT constants as per the methods (in order to match expe. results )


      Because in droplet method, I am simulating the total of 1mole (n-dodecane) and for the same case in Multicomponent we are giving 0.999999mole of (n-dodecane) and 1e-6 mole of (ethanol or any alternative fuel ). In the next cases, we want to study the effect of various fuel blends (methanol, ethanol and n-butanol in 5,10,20,40% blending ) with diesel (n-dodecane). For all these cases we want to use the Multicomponent method and we will not change the KH-RT constants. 

    • DrAmine
      Ansys Employee

      The Levich constant Cl is adjustable and might be related to the cross-sectional area of the contracted liquid jet or just to B1.

    • anilwakale
      Subscriber

      That's Helpful for me but again I will go back to my original question 


      can we find the Vapor pressure data as output file from Fluent? 

    • DrAmine
      Ansys Employee

      There are some default for some materials like for water with 32 points  (data here is from NIST). Is you material a default one? 

    • anilwakale
      Subscriber

      No. Mine material is not default.  Ansys Fluent don't have properties of n-butanol and n-dodecane. I will be happy if fluent database increase to at least to c12-16 higher hydrocarbons and alcohols like (n-butanol, Biodiesel, DME, DEE if possible)


      I am manually entering the value for other fuels (n-dodecane/butanol) but at blends, our mixture fraction should vary and that's why we are giving UDF to saturation vapor pressure data by UNIFAQ method from the book Not only for this but for rest properties to Cp, viscosity, surface tension. Additionally, we want to vary latent heat too but the option was not available. we compared our all data points with NIST std values with less than 6 to 8% we are maintaining all the thermo property values 


      I want to share an animation of spray combustion and want to ask some question too if you permit 

    • DrAmine
      Ansys Employee

      So you are providing the saturation property and other materials via DEFINE_DPM_PROPERTRY. There is not a way to plot that in Fluent. Regarding material properties, it is worth mentioning that we all CHEMKIN power in our portfolio. Saying that one can find for many fuels suitable descriptions.


      As ANSYS Stuff I cannot answer outside of this community..

    • anilwakale
      Subscriber

       Exactly and we modified some code DEFINE_DPM_PROPERTRY . 


      Completely agree with you


      I will share my spray animation soon and the question also   

    • anilwakale
      Subscriber

      Temp cont.


      If you see the Time is 5e-5 and my starting time itself was 1e-6 and I plotted with 2000-2600 temp range


      so in this gif. we are not getting any ignition delay. using Steady diffusion flamelet model. I have to simulate the same case as given below (link)



      Not with LES but with RANS using Fluent. We are getting above lift of length for non-reacting spray but in reacting our spray will start combusting as soon it enters into the chamber. In the link above its clear that spray will come in liquid form enters into the hot environment and with some Ignition delay spray starts combusting. Do you suggest any tutorial for it? Initial conditions of the chamber are mimic of ECN spray A expe. (Pin=61MPA, Tin=900K, fuel Temp=363K and rest all BC)

    • DrAmine
      Ansys Employee

      I recommend using LES and detailed reaction mechanism with EDC or FGM.


      For the ignition delay you can check if one of the methods in Fluent are suitable or not.

Viewing 10 reply threads
  • You must be logged in to reply to this topic.