Open access peer-reviewed chapter

Equation of State

By Eman Mohamed Mansour

Submitted: February 5th 2019Reviewed: September 25th 2019Published: March 4th 2020

DOI: 10.5772/intechopen.89919

Downloaded: 278

Abstract

An equation of state (EOS) is a thermodynamic expression that relates pressure (P), temperature (T), and volume (V). This equation is used to describe the state of reservoir fluids at given conditions. The cubic equations of state (CEOS) such as Van der Waals, Redlich-Kwong, Soave, and Peng-Robinson are simple models that have been widely used in the oil industry. This chapter expressed literature for EOS that varies from simple expressions to multiple constant and convoluted types of equations. Many attempts have been made to describe the thermodynamic behavior of fluids to predict their physical properties at given conditions. So, several forms of the equation of state have been presented to the oil industry in order to calculate reservoir fluid properties. The heat exchanger is important in wildly fields as in aerospace, petrochemical industry, refrigeration, and other fields. The optimization design of the heat exchanger is a great significance to industry process to reduce production cost, realize energy conservation, and reduce energy consumption.

Keywords

  • EOS types
  • EOS importance
  • EOS history
  • cubic equation of state
  • heat transfer

1. Introduction

In the last six decades, petroleum engineers realized the importance of using EOS for PVT modeling in addition to the following [1]:

  1. Predicting physical properties of the black oil in the laboratories is very expensive.

  2. Difficulty to obtain a representative sample due to reservoir nature or operation issues.

  3. Insufficient sample volume in order to obtain a complete analysis.

  4. Error in lab analyses.

  5. Samples are nonrepresentative (not monophasic sample, OBM contamination, etc.).

  6. Lab report quality check.

  7. Estimating the hydrocarbon reserves.

1.1 Classification of equation of state

There are different types of EOS which fall into three categories:

1.1.1 First class of EOS

These equations are basically cubic equation of state. The cubic equations of state such as the Van der Waals [2], Redlicha and Kwong [3], Soave-Redlich-Kwong [4], and Peng-Robinson [5] equations give reasonable results for the thermodynamic behavior of real fluids.

1.1.2 Second class of EOS

These EOS are non-cubic in form. They are providing accurate results for both vapor and liquid phases. The Benedict et al. [6] equation is a good example for this class equation.

1.1.3 Third class of EOS

These are nonanalytical EOS that are highly constrained for some specific fluids [7]. Even though they are constrained, they are capable of expressing real fluid thermodynamic properties precisely.

Among all these EOS, the first-class EOS is more useful because it provides an analytical solution than the more complex and complicated non-cubic second type and nonanalytical third type that require time-consuming iterative calculations. In general, the overall performance in fluid properties prediction is somewhat better using the Soave-Redlich-Kwong (SRK) equation than using the Redlich-Kwong (RK) and Van der Waals EOS [8].

2. Development history of the equation of state

Several forms of EOS have been presented to the petroleum industry to estimate hydrocarbon reservoir fluid properties and sought to a better representation of the PVT relationship for fluids [9].

In 1662, Robert Boyle (Boyle’s law) discovered that for a constant temperature, there is an inverse relationship between volume of gas and its pressure (P ∝ V−1). In 1780, Jacques Charles (Charles’s Law) showed that the volume of gas is proportional to the absolute temperature at a constant pressure (V ∝ T). In 1834, Clapeyron combined these two results into the ideal gas law, PV = RT [10], assuming that the molecules are very far and have no attractive or repulsive forces between them and elastic collisions between these molecules. This equation is known as the ideal gas law and/or the general gas law. It is expressed mathematically as [11].

PV=nRTE1

where P: absolute pressure, psia; V: volume, ft3; T: absolute temperature, °R; R: the universal gas constant (10.73159 ft3 psia °R−1 lb-mole−1; n: number of moles of gas, lb-mole.

For gases at low pressures, the ideal gas law is a convenient satisfactory tool. The application of the ideal gas law at higher pressures may lead to errors up to 500%, compared to 2–3% at atmospheric pressure. Real gases behave differently than ideal gases, the reason for this deviation is that the ideal gas law was derived under the assumption that the volume of molecules is very small and neither molecular attraction nor repulsion exists between them, and this is not the real case. In order to write an equation of state for a real gas, a correction factor has to be inserted into the ideal gas equation [12]:

PV=ZnRTE2

where Z: the correction factor which is known as the compressibility factor.

The equation has different names, such as the compressibility equation and/or the real gas equation [13]. A review of recent advances in the empirical cubic EOS field is presented next [11]. Van der Waals [2] is one of the earliest attempts to represent the behavior of real gases by an equation, where the two assumptions were made for the ideal gas EOS:

  1. The gas molecule volume is very small compared to the volume of the container.

  2. There are no attractive or repulsive forces between the gas molecules or the walls of the container.

Van der Waals attempted to eliminate these assumptions in the development of an empirical EOS for the real gases.

First assumption elimination: the gas molecules occupy a considerable fraction of the volume at higher pressures, and the volume of the molecules (b) is subtracted from the actual molar volume (V) to give the following expression:

p=RTvbE3

Second assumption elimination: he added corrective term (a), denoted by (a/V2), in order to account for the attractive forces between molecules.

Van der Waals introduced the following equation (Eq. (4)):

p+aVM2VMb=RTE4

where a: attraction parameter; b: repulsion parameter.

The symbol “a” is considered a measure of the intermolecular attractive forces between the molecules. “b” is known as the co-volume and considered to reflect the volume of molecules [2]. The “a” and “b” values can be obtained from the critical properties of the fluid [14], where the repulsion pressure, prepulsion, is represented by the term RT/(Vm – b), and the attraction pressure, pattraction, is described by a/Vm 2. The Van der Waals equation of state despite its simplicity, while it is provide a correct description and qualitative of the PVT substances behavior in the liquid and gaseous phases. Yet, it is not accurate enough to be suitable for design purposes. The equation of state approach for calculating physical properties and phase equilibrium proved to be a powerful tool, and much energy was devoted to the development of new and accurate equations of state [11]. Other researchers began attempts to improve Van der Waals equation of state for over 100 years. Usually a change of the molecular attraction term (a/VmM2) was proposed. Clausius in 1880 [15] proposed that the molecular attraction term was inversely proportional to temperature [16]:

p+aTVM+c2VMb=RTE5

The addition of a fourth constant (c) enabled better agreement with data. However, mathematical manipulations required in thermodynamic calculations were more difficult. So Berthelot in 1899 [17] removed the constant (c), resulting in the following equation:

p+aTVM2VMb=RTE6

Dieterici in 1899 [18] handled the temperature dependence of the molecular attraction term in a different manner [6]:

PEXPaVMRTVMb=RTE7

Lorentz in 1881 [19] addressed the molecular volume term [20]:

p+aVM2VMbVMVM+b=RT.E8

Wohl in 1927 [21] considered the effect of temperature on the molecular attraction term:

P+aTVMVMbcT2VM3VMb=RTE9

The constants a, b, and c in the equations above have different values for different substances. Several investigations proposed virial-type of EOS. Kammerlingh-Onnes in 1901 [22] proposed the virial equation of state as follows [23]:

PVM=RT1+BVM+CVM2+..E10

where B and C are not constants which are functions of temperature and are called the second and third virial coefficients. Beattie and Bridgeman in 1927 published a five-constant equation that gives a satisfactory representation of volumetric properties except in the critical region [24]:

P=RTVM21cVMT3VM+Bo1bVMAo1a/VMVME11

Benedict et al. [6] suggested a multiparameter equation of state known as the Benedict-Webb-Rubin (BWR) equation [6]:

P=RTVM+BoRTAoCo/T2VM2+bRTaVM3+VM6+cT2VM31+γVM2EXPγVM2E12

This equation may be considered a modification of the Beattie-Bridgeman equation of state where A0, B0, C0, a, b, c, α, and γ are eight adjustable parameters. The BWR equation could treat critical components and was able to work in the critical area. However, the BWR equation suffers from some disadvantages [25]. Perhaps, the most important model for the modification of the Van der Waals equation of state is the Redlich-Kwong (RK) (1949) which is demonstrated by an adjustment of the Van der Waals’s attraction term (a/Vm 2) and includes the system temperature explicitly. They could improve the prediction of the physical and volumetric properties of the vapor phase. In RK EOS, the attraction pressure term was replaced with a generalized temperature-dependent term (Eq. (13)) [3]:

ρ=RTVbαVV+bTE13

For pure substances, the equation parameters a and b are usually expressed as.

b=ΩbRTc/PcE14
a=ΩαR2Tc2.5/PcE15

where Ωa = 0.42747 and Ωb = 0.08664.

Replacing the molar volume (V) in Eq. (13) with (ZRT/P) and rearranging give.

Z3Z2+ABB2ZAB=0E16

where

B=bpRTE17
A=apR2T2.5E18

Three real roots in the two-phase region are yielded. The largest root corresponds to the compressibility factor of the gas phase, Zv, while the smallest positive root corresponded to that of the liquid, ZL [11].

For mixtures, the equation parameters a and b are usually expressed as am and bm for a hydrocarbon liquid mixture with a composition of xi:

am=i=1nXiai2E19
bm=i=1nXibiE20

am and bm for a hydrocarbon gas mixture with a composition of yi:

αm=i=1nyiai2E21
bm=i=1nyibiE22

where n: number of components in the mixture; ai: Redlich-Kwong a parameter for the i component; bi: Redlich-Kwong b parameter for the i component; bm: parameter b for mixture; xi: mole fraction of component i in the liquid phase; yi: mole fraction of component i in the gas phase.

Replacing the molar volume (V) in Eq. (13) with (ZRT/P) and rearranging give.

Z3Z2+ABB2ZAB=0E23

where

B=bmpRTE24
A=amPR2T2.5E25

Then the compressibility factor of the gas phase or the liquid can be calculated.

Joffe and Zudkevitch [26] showed that a substantial improvement in the representation of fugacity of gas mixtures could be obtained by treating interaction parameters as empirical parameters [26]. Spear et al. [27] also states that the RK equation of state could be used to calculate the vapor-liquid critical properties of binary mixtures [28]. Chueh and Prausnitz [29] showed that the RK equation can be adapted to predict both vapor and liquid properties. Spear et al. [28] gave seven examples of systems for which the vapor-liquid critical properties of hydrocarbon mixtures could be calculated by using the RK equation of state. Carnahan and Starling [30] used the Redlich-Kwong equation of state to calculate the gas-phase enthalpies for a variety of substances [30]. Their results showed that the Redlich-Kwong equation was a significant improvement over the Van der Waals equation. Other workers applied the Redlich-Kwong equation to the critical properties and the high-pressure phase equilibria of binary mixtures. The results showed that the accuracy of the Redlich-Kwong equation of state calculations for ternary systems was only slightly less than that for the constituent binaries [31].

The success of the Redlich-Kwong equation has been the impetus for many further empirical improvements. One of the milestones in developing of CEOS was reported by Soave [4]. His development in the evaluation of the parameter in the attraction pressure term for the RK equation is shown in (Eq. (22)). Soave replaced the term (a/T0.5) in Eq. (22) with a more general temperature-dependent term, denoted by a α (T), to give

ρ=RTVbaαTVV+bE26

where α(T) is a dimensionless factor. Soave used vapor pressures of pure components to introduce an expression for the temperature correction parameter α(T). At temperatures other than the critical temperature, the correction parameter α(T) was defined by the following equation:

αT=1+m1Tr2E27

Soave correlated the parameter “m” with the centric factor (ω) to give.

m=0.480+1.57,0.172E28

where Tr: reduced temperature, °R; ω: a centric factor of the substance; T: system temperature, °R.

For pure substances the equation parameters a and b are usually expressed as.

b=ΩbRTc/PcE29
a=ΩaR2Tc2/PcE30

In general, most EOS inputs are only the critical properties, and a centric factor of each components is shown in Table 1.

CompoundFormulaMolecular weightCritical pressure (psla)Critical temperature (°F)Critical volume (ft3/lb)Liquid specific gravity (water = 1)Gas specific (air = 1)Acentric factor
MethaneCH416.042667.0−116.660.0985(0.3)0.554000.0115
EthaneC2H630.069706.689.920.07750.356431.038300.0994
PropaneC3H844.096615.5205.920.07280.507381.522700.1529
IsobutaneC4H1058.122527.9274.410.07150.562952.007100.1865
n-ButaneC4H1058.122550.9305.550.07030.584082.007100.2003
IsopentaneC5H1272.149490.43690.06850.624602.491400.2284
n-PentaneC5H1272.149488.8385.80.06760.631132.491400.2515
n-HexaneC6H1486.175436.9453.30.06880.664042.975800.2993
n-HeptaneC7H16100.202396.8512.90.06820.688193.460200.3483
n-OctaneC8H18114.229360.7564.20.06730.706983.944500.3977
n-NonaneC9H20128.255330.7610.80.6930.721864.428900.4421
n-DecaneC10H22142.282304.6652.20.07030.734064.913300.4875
Carbon monoxideCO28.01506.7−220.630.05270.792650.967200.0510
Carbon dioxideCO244.011070.087.760.03430.822031.519700.2239
Hydrogen sulfideH2S34.0821306.5212.810.04620.802691.176900.1010
Air28.9586551.9−220.970.04580.876031.00000
HydrogenH22.0159190.7−399.90.53190.070870.06961−0.2140
OxygenO231.9988731.4−181.430.03671.142301.105000.0222
NitrogenN228.0135492.5−232.530.05110.806870.967400.0372
WaterH2O18.01533200.1705.10.049751.000000.622100.3443

Table 1.

Physical properties of each components.

where Ωa and Ωb are the SRK dimensionless pure component parameters:

Ωa = 0.42747.

Ωb = 0.08664.

Replacing the molar volume (V) in the equation with (ZRT/p) and rearranging give the compressibility factor Z:

Z3Z2+ABB2ZAB=0E31

where

B=bmpRTE32
A=amPRT2E33

For mixtures, the equation parameters a and b are usually expressed as am and bm for a hydrocarbon liquid mixture with a composition of xi:

am=ij[xixjaiajαiαj1kijE34
bm=iXibiE35

The following is the calculation for am and bm for a gas mixture with a composition of yi:

am=ij[yiyjaiajαiαj1kijE36
bm=iyibiE37

A binary interaction parameter (BI), classically noted as kij, is usually involved in the “a” parameter expression to provide more flexibility to the EOS and designed to characterize any binary system formed by components i and j in the hydrocarbon mixture [32]. Vidal and Daubert [33], Graboski and Daubert [34], and Slot-Petersen [35] suggested that no BIs were required for hydrocarbon systems. However, with no hydrocarbons present, binary interaction parameters can improve the phase in volumetric behavior predictions of the mixture by the SRK EOS for compressibility factor calculations of the gas or the liquid phases [34, 36, 37]. The equilibrium ratio, Ki, that is, Ki = yi /xi, can be redefined in terms of the fugacity of component:

Ki=fiL/XiPfiV/yiP=ΦiLΦivE38

where fi v = fugacity of component “i” in the gas phase; fi L = fugacity of component “i” in the liquid phase; Φiv= fugacity coefficient of component “i” in the vapor phase; ΦiL= fugacity coefficient of component “i” in the liquid phase.

Soave proposed the following expression for the fugacity coefficient of component i in the liquid phase:

lnfiLXiP=lnΦiL=biZl1bmlnZLBAP2ψiambibmln1+BZLE39

where

ψj=jxjaiajαiαj1kijE40
am=ij[xixjaiajαiαj1kijE41

Fugacity coefficient of component i in the gas phase:

lnϕj=biZi1bmlnZiBAB2ψiambibmln1+BZiE42

where:

ψj=jyjaiajαiαj1kijE43
am=ij[yiyjaiajαiαj1kijE44

3. Heat exchanger

Heat exchanger is an energy (heat) exchange equipment, where it transfers the heat from a working medium to another working medium. Knowing heat exchanger is important in wildly fields as in aerospace, petrochemical industry, refrigeration, and other fields. The optimization design of the heat exchanger is a great significance to industry process to reduce production cost, realize energy conservation, and reduce energy consumption [38]. The development technique for different types of the heat exchanger has been reviewed by many researchers. The development method can be by two ways: passive method and active method. The passive method is to generate swirling flow and disturb the thermal boundary layer by installing vortex generator or tabulators such as baffle, rib, winglet, wing, etc. The active method is to add the external power to increase efficiency and heat transfer rate such as vibration. So the use of the active method must consider both benefit of the system and additional power cost [39].

3.1 Optimization techniques of heat exchanger

The optimization techniques of heat exchanger can be shown at three different stages as the following [40]:

  1. Identification of the lowest initial cost of a heat exchanger design that meets the process specifications.

  2. Identification of a heat exchanger design that will work most acceptably over the plant lifetime.

  3. Identification of the minimum total cost of the process by choosing heat exchangers system and auxiliary components that will make the best plant process specifications.

3.2 Criteria of the ideal heat exchanger

The dissimilarities of optimization techniques levels can be understood if we list the required criteria of the ideal heat exchanger as follows [41]:

  1. The heat exchangers should be flexible enough to meet specifications process such as normal fouling transients and seasonal and diurnal changes in service stream temperatures.

  2. Special requirements as weight, length, or inventory standards are important for heat exchangers especially in retrofit applications.

  3. The heat exchanger must endure operation under standard and foreseeable operating conditions, maintain the mechanical stresses of manufacturing transport, and minimize the effects of fouling and corrosion.

  4. The heat exchanger must be maintainable, cleaning, repair or replacement and its components as gaskets and tubes with minimum downtime.

  5. The exchanger must achieve process specifications, i.e., achieve any changes in the thermal conditions by allowable pressure drops.

  6. There are other requirements, such as experience, capability of operating, maintenance personnel, and manufacturing time.

  7. The exchanger should cost as little as possible.

4. Fluid flow and heat transfer mechanism

The heat transfer and fluid flow mechanism in porous media is important in engineering and industrial fields such as petroleum and chemical engineering [42]. This mechanism occurs for many studies such as in microchannels and nanofluids.

Example 1: In case of study characteristics of fluid flow and heat transfer in the (100) silicon microchannel heat sink, the heat convection capabilities in the phase changes as well as in a single-phase flow and the mechanism of bubble nucleation. In the heat transfer characteristics, the results illustrate that changing in the phase process in the microchannels reduces environment working temperature and absorbs the heat. Six different microchannel geometries are selected for the heat transfer experiment as shown in Table 2.

Chip nameWidth (μm)
Wc
Depth (μm)
Hc
Hydraulic diameter (μm)Number of channels
Chip 140026022110
Chip 230013015013
Chip 325018413415
Chip 420014810919
Chip 51501138325
Chip 6100785738

Table 2.

Specification of the sink.

Figure 1 shows that the decreasing wall temperature phenomenon during the phase change is the same as Peng and Wang [43].

Figure 1.

The heat flux and channel wall temperature.

On the aspect of fluid flow characteristics, the effects of the viscosity and friction coefficient of the fluid in the microchannels are much significant than the macros. Where the specifications of the sink are registered in Table 2, Chip 1–4 are prepared for fluid flow experiment. The friction factor is decreasing with the power of Reynolds number as shown in Figure 2 [44].

  1. Nanoparticle types and concentration in the base fluids

  2. Changing in thermophysical properties such as viscosity, surface tension, thermal conductivity, density, and heat capacity

  3. The operation condition especially the mass and heat fluxes

  4. The modification of wettability and capillary wicking force surface roughness

Figure 2.

Exponential relations between the friction factor and the Reynolds number. Example 1: in the case of study characteristics of fluid flow and heat transfer, nanofluid is widely utilized in numerous industrial applications such as boiler tubes, evaporators, and cooling of reactors in a nuclear power plant. The main parameters that directly influence on the heat transfer performance are listed as follows [45].

During the last 2 years, there were some review papers which outlined the subject of boiling heat transfer using nanofluids as a new category in thermal fluids. Available results reported that the effect of nanoparticles on the flow boiling HTC is conflicting, but the CHF could enhance on 50%. During the boiling process, parameters such as flow pattern and pressure drop were affected by the deposition of nanoparticles [46]. Authors concluded that using nanofluids might intensify the boiling HT and CHF, depending on many parameters related to additives, nanoparticles, geometry concentration, and fluid properties [47].

Their work shows how the nanofluids can achieve high heat flux with small temperature differences throughout the boiling process, which modify the critical heat flux [48]. All article reviewers said that nanofluids is a complicated phenomenon and it is not fully understood from mechanism of boiling heat transfer and two-phase flow. Collected studies show enhancement in CHF, and its improvement could achieve more than 50% [49]. HTC behavior could increase or decrease during flow boiling and pool, and it depends on several parameters related to surface texture on thermophysical properties [50]. Wang et al. shows the mechanism of depositing nanoparticles on the heating surface and creating the porous layer as shown in Figure 3. The results show that the CHF of flow boiling is enhanced up to 18% as compared to conventional fluid. This enhancement increases with increasing some parameters, for example, the pressure system and the channel diameter [45].

Figure 3.

Illustration of the mechanism of flow boiling CHF using nanofluid introduced.

Example 2: in the case of study fluid flow and heat transfer characteristics using nanofluid in a single-phase turbulent flow by using helically corrugated tubes, pitch-to-diameter ratio (P/DH = 0.18, 0.22 and 0.27) and rib-height-to-diameter ratio (e/DH = 0.02, 0.04 and 0.06) of helically corrugated tubes on the heat transfer enhancement, isothermal friction, and thermal performance factor in a concentric tube heat exchanger are examined. Results illustrate that the thermal performance of the corrugated tube and heat transfer are increased as compared to those of the smooth tube. The rate increase in heat transfer rate is between 123 and 232%, depending on the rib height/pitch ratios and Reynolds number. The friction factor (average) of the corrugated tube is between 1.46 and 1.93 over the smooth tube [51].

Abbreviations

EOSequation of state
BWRBenedict-Webb-Rubin
RKRedlich-Kwong
SRKSoave-Redlich-Kwong
CEOScubic equations of state
MSRKEmodified Soave-Redlich-Kwong equation of state
psystem pressure, psia
Pccritical pressure, psia
Prreduced pressure, psia
piinitial pressure, psia
Tsystem temperature, oF
Tccritical temperature, oR
Trreduced temperature, oR
Tiinitial temperature, oF
aequation of state attraction parameter
bequation of state repulsion parameter
Bsecond virial coefficient
Cthird virial coefficient
a, b, c, A0, B0, C0constant in Benedict-Webb-Rubin equation
a, b, c, A0, B0constant in Beattie and Bridgeman equation of state
fv ifugacity of component i in the gas phase
fL jfugacity of component j in the liquid phase
Øv ifugacity coefficient of component i in the vapor phase
ØL ifugacity coefficient of component i in the liquid phase
xmole fraction of gas phase
ymole fraction of liquid phase
kequilibrium ratio for the component
Kijinteraction coefficient parameter
A, Bparameter in Soave-Redlich-Kwong equation of state
aTtemperature-dependent coefficient in Soave-Redlich-Kwong equation of state
acconstant coefficient in Soave-Redlich-Kwong equation of state
aTitemperature-dependent coefficient of component i
mparameter in Soave-Redlich-Kwong equation of state
Zcompressibility factor
Zithe mole fraction of component in the mixture
ZLcompressibility factor of the liquid phase
Zvcompressibility factor of the gas phase
ωa centric factor of the substance
nlnumber of moles in liquid phase
nvnumber of moles in gas phase
nttotal number of moles of mixture
Flone mole liquid of total fraction
Fvone mole vapor of total fraction
Ximole fraction in the liquid phase
yimole fraction in the vapor phase
γconstant in Benedict-Webb-Rubin equation of state

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Eman Mohamed Mansour (March 4th 2020). Equation of State, Inverse Heat Conduction and Heat Exchangers, Suvanjan Bhattacharya, Mohammad Moghimi Ardekani, Ranjib Biswas and R. C. Mehta, IntechOpen, DOI: 10.5772/intechopen.89919. Available from:

chapter statistics

278total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Influence of Input Parameters on the Solution of Inverse Heat Conduction Problem

By Rakhab C. Mehta

Related Book

First chapter

Introduction to Infrared Spectroscopy

By Theophile Theophanides

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us