Scribd Logo
Upload

Welcome to Scribd!

  • 3 views

    Solve-a-nonlinear-equation-1

    Uploaded by

    akhileshpandey

    Copyright:

    © All Rights Reserved
    Download as XLSX, PDF, TXT or read online from Scribd

    Solve-a-nonlinear-equation-1

    Uploaded by

    akhileshpandey
    3 views2 pages

    Copyright

    © © All Rights Reserved

    Available Formats

    XLSX, PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    Copyright:

    © All Rights Reserved
    Download as XLSX, PDF, TXT or read online from Scribd
    Download as xlsx, pdf, or txt
    3 views2 pages

    Solve-a-nonlinear-equation-1

    Uploaded by

    akhileshpandey

    Copyright:

    © All Rights Reserved
    Download as XLSX, PDF, TXT or read online from Scribd
    Download as xlsx, pdf, or txt
    You are on page 1of 2

    Solve volume in Peng-Robinson equation of state

    Demonstrate three ways to solve for V


    1. from graph Solve for V so that f(V) = 0
    2. use sliders to vary V until f(V) = 0 a 396195 f(V) 1617.511394
    3. use Solver to find V b 26.65 1st term -1256.073
    k 0.71 2nd term 361.3455
    a 1.138
    Use these sliders to change volume P 0.093 MPa
    (cell C11) to minimize f(V) T 250 K
    V 25.00 cm3
    Rg 8.314 cm3 MPa/(mol K)
    TC 304.2 K
    Solutions (cm3) for P = 0.093 MPa PC 7.382 MPa
    41.54, 123.1, 22,158 w 0.228

    Peng-Robinson EOS for CO2 at 250 K


    1.0E+02

    1.0E+01
    Pressure (MPa)

    1.0E+00

    1.0E-01

    1.0E-02
    1.0E+01 1.0E+02 1.0E+03 1.0E+04 1.0E+05
    Volume (cm3)
    University of Colorado Boulder
    www.LearnChemE.com

    V P
    cm3 MPa

    The constants (a, b, 𝛼, 𝜅, 𝜔) are functions of


    34.0 83.4 Peng-Robinson equation of state

    criticial pressure and temperature (𝑃_𝐶, 𝑇_𝐶), the


    36.0 42.5

    (𝑇_𝑟=𝑇/𝑇_𝐶 )
    38.0 19.9 ideal gas constant (R) and reduced temperature

    𝑓(𝑉)= P-
    40.0 6.61

    𝑅𝑇/(𝑉−𝑏)+𝑎𝛼/(𝑉^2+2𝑏𝑉−𝑏^2 )
    41.0 2.10
    41.5 0.223

    𝑃=𝑅𝑇/(𝑉−𝑏)−𝑎𝛼/(𝑉^2+2𝑏𝑉−𝑏^2 )
    41.56 0.014

    𝑎=(0.45724𝑅^2 𝑇_𝐶^2)/𝑃_𝐶
    41.561 0.010
    122.0 0.003
    122.1 0.012 𝑏=(0.07780𝑅𝑇_𝐶)/𝑃_𝐶
    122.3 0.029 𝛼=(1+𝜅(1−𝑇_𝑟^0.5 ))^2
    𝜅=0.37464+1.54226𝜔 −0.26992𝜔^2
    122.6 0.054
    123.0 0.087
    124 0.169
    130 0.618
    150 1.72
    190 2.82
    300 3.32
    500 2.76
    1500 1.22
    5000 0.400
    22160 0.093
    90000 0.023

    You might also like

    pFad - Phonifier reborn

    Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

    Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


    Alternative Proxies:

    Alternative Proxy

    pFad Proxy

    pFad v3 Proxy

    pFad v4 Proxy