Method of Bioammonia Production from Wastewater Through Application of Mass-Transfer Reaction Kinetics

Abstract

Systems and methods are disclosed in which ammonia may be recovered efficiently from wastewater through the application of the mass-transfer reaction kinetic model. The recovery process may involve the stripping of NH.sub.3 from wastewater by aeration and the absorption of this stripped NH.sub.3 gas into a recovery solution by utilizing stoichiometric acid-base reactions and the operating parameters optimized by the model which facilitate both processes. Another step may be performed in which saturated solution comprising the recovered NH.sub.3 may be further evaporated and liquefy to produce a pure bioammonia product. Applying the mass-transfer reaction kinetic to the recovery can not only enhance the NH.sub.3 removal and recovery yield but also increase the NH.sub.3 removal and recovery rate. Numerous operating parameters can be adjusted to maximize the recovery yield and recovery rate.

Claims

1. A method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+), comprising the steps of: supplying wastewater, a Brnsted base, a Brnsted acid and water; supplying a stripping column equipped with one or more air diffusers, a recovery column equipped with one or more air diffusers, an evaporator, and, a condenser; the stripping column, recovery column, evaporator and condenser being in fluid communication; mixing, in the stripping column, the wastewater with the Brnsted base to form a stripping mixture, the Brnsted base being operative to react with ammonium ions present in the wastewater to produce ammonia gas and other gasses in the stripping mixture; stripping, in the stripping column, the ammonia from the stripping mixture by aerating the stripping mixture to produce ammonia gas and other gasses, wherein the mass transfer rate per unit volume of the conversion of ammonia in the stripping mixture to ammonia gas is at least 0.034 min.sup.1; recovering, in the recovery column, the ammonia gas by absorbing the ammonia gas into a recovery solution comprising the Brnsted acid and water, wherein the mass transfer rate per unit volume of the ammonia gas to the recovery solution is at least 0.006 min.sup.1, the recovery solution and ammonia gas absorbed therein defining a liquid ammonia product; heating, in the evaporator, the liquid ammonia product to form an evaporated mixture comprising substantially ammonia gas; and condensing, in the condenser, the evaporated mixture into liquid ammonia.

2. A method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+), comprising the steps of: supplying wastewater, a Brnsted base, a Brnsted acid and water; supplying a stripping column equipped with one or more air diffusers, a recovery column equipped with one or more air diffusers, an evaporator, and, a liquefaction device; the stripping column, recovery column, evaporator and liquefaction device being in fluid communication; mixing, in the stripping column, the wastewater with the Brnsted base to form a stripping mixture, the Brnsted base being operative to react with ammonium ions present in the wastewater to produce ammonia gas and other gasses in the stripping mixture; stripping, in the stripping column, the ammonia from the stripping mixture by aerating the stripping mixture to produce ammonia gas and other gasses; and recovering, in the recovery column, the ammonia gas by absorbing the ammonia gas into a recovery solution comprising the Brnsted acid and water, the recovery solution and ammonia gas absorbed therein defining an ammonia product; heating, in the evaporator, the ammonia product to form an evaporated mixture comprising substantially ammonia gas; and liquefying, in the liquefaction device, the evaporated mixture into liquid ammonia (NH.sub.3).

3. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4+) of claim 2, further comprising the step of filtering the wastewater so as to at least partially separate any solid components from the wastewater.

4. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, further comprising repeating the steps stripping ammonia gas from the stripping mixture and, absorbing the ammonia gas in the recovery solution, until the recovery solution is saturated with the ammonia gas.

5. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein the Brnsted base is selected from the group consisting of Ca(OH).sub.2, CaCO.sub.3, Ca(HCO.sub.3).sub.2, NaOH, Na.sub.2CO.sub.3, Na.sub.2(HCO.sub.3).sub.2, KOH, K.sub.2CO.sub.3, and K(HCO.sub.3).sub.2.

6. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein the Brnsted acid is selected from the group consisting of HNO.sub.3, H.sub.2SO.sub.4, H.sub.3PO.sub.4, C.sub.6H.sub.8O.sub.5, formic acid, potassium acetate, trisodium phosphate, and sodium acetate.

7. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein a mole of an anion of the Brnsted base is at least stoichiometric to a mole of the ammonium ions in the wastewater.

8. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein a mole of an anion of the Brnsted base is at least 10% more than stoichiometric to a mole of the ammonium ions in the wastewater.

9. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein a mole of a cation of the Brnsted acid is at least stoichiometric to a mole of the stripped ammonia gas.

10. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein a mole of a cation of the Brnsted acid is at least 10% more than stoichiometric to a mole of the stripped ammonia gas.

11. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein during the step of stripping the ammonia gas from the stripping mixture, the mass transfer rate per unit volume of the conversion of ammonia to ammonia gas is at least 0.034 min.sup.1.

12. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein during the step of recovering the ammonia gas by absorbing the ammonia gas into a recovery solution, the mass transfer rate per unit volume of the ammonia gas to the recovery solution is at least 0.006 min.sup.1.

13. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein the value of a liquid-to-gas mass transfer coefficient (K.sub.L,s) multiplied by an interface area per volume of liquid (a.sub.s) is not be lower than $\frac{Q_{G,s}{H}_s}{V_{L,s}}$ where Q.sub.G,s, V.sub.L,s and H.sub.s are the air flowrate, the volume, and Henry's law constant, respectively.

14. The method of producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, wherein a gas-to-liquid mass transfer coefficient (K.sub.L,r) multiplied by an interface area per volume of liquid (a r) is not be lower than $\frac{Q_{G,r}{H}_r}{V_{L,r}}$ where Q.sub.G,r, V.sub.L,r and H.sub.r are the air flowrate, the volume, and Henry's law constant, respectively.

15. The method for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, further including the step of stripping NH.sub.3 and CO.sub.2 gases from the stripping mixture and absorbing these gasses in the recovery solution to produce ammonium bicarbonate as an organic fertilizer.

16. The method for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2 by applying a pressure to the evaporated ammonia gas.

17. The method for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 2, further comprising the step of reacting the liquid ammonia with a CO.sub.2 supply to produce urea.

18. The method for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 17, wherein the CO.sub.2 supply is selected from the group consisting of CO.sub.2 gas generated by anaerobic digestion of organic wastes, CO.sub.2 gas sequestered from an industrial plant, and CO.sub.2 gas captured from the stripping column.

19. An apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) comprising: a stripping column equipped with one or more air diffusers, a recovery column equipped with one or more air diffusers, an evaporator; and, a liquefaction device; the stripping column, recovery column, evaporator and liquefaction device being in fluid communication; wherein, in the stripping column wastewater is mixed with a Brnsted base to form a stripping mixture, the Brnsted base being operative to react with ammonium ions present in the wastewater to produce ammonia gas and other gasses in the stripping mixture; wherein, the one or more air diffusers are located in the base of the stripping column and function to aerate the stripping mixture to release the ammonia gas and other gasses contained therein; and wherein, in the recovery column, the ammonia gas is absorbing into a recovery solution comprising a Brnsted acid and water, the recovery solution and ammonia gas absorbed therein defining an ammonia product; wherein, in the evaporator, the ammonia product is heated to form an evaporated mixture comprising substantially ammonia gas; and wherein, in the liquefaction device, the evaporated mixture is cooled to form liquid ammonia (NH.sub.3).

20. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein at least one of the one or more air diffusers in the stripping column or recovery column is selected from the group consisting of a sintered stone diffuser, a membrane diffuser, a sparger diffuser, or combinations thereof.

21. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein at least one of the one or more air diffusers has an orifice diameter of 10 mm or less and a surface area per length in the range of 0.5-3 inch.sup.2 inch.sup.1.

22. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the number of the one or more air diffusers in the stripping column or recovery columns is the square of the diameter of the respective column in inches.

23. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein at least one of the one or more air diffusers in the stripping column or the recovery columns is a sparger pipe diffuser, and wherein the number of the one or more air diffusers is equal to the diameter of the respective column in inches.

24. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein at least one of the one or more air diffusers in the stripping column or recovery columns is a sparger pipe diffuser, and wherein the number of the one or more air diffusers is equal to one inch more than the value of the diameter of the respective column diameter in inches.

25. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the stripping column further comprises one or more sub-columns each equipped with one or more air diffusers, the one or more sub-columns having a diameter in the rage of 1 in to 3 in, the one or more sub-columns being operative to increase the efficiency of the step of stripping the ammonia gas from the filtered wastewater.

26. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein an air pump is placed between the stripping column and the recovery column to increase the flow rate of ammonia gas to the recovery column.

27. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the liquefaction device uses a coolant selected from the group consisting of iso-propanol, acetone, 1-, or 2-butanol, 2-butanone, ethanol, diethyl ether, heptane, n-hexane, pentane, 1-propanol, tetrahydrofuran, and triethyl amine.

28. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the coolant is supplied to the liquefaction device from a cryogenic storage tank, the cryogenic storage tank having a cryogenic substance, the cryogenic substance comprising liquid nitrogen.

29. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the stripping column or the recovery column further includes a mesh screen operative to break up air bubbles rising up the column from the one or more air diffusers.

30. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the said ammonia product in the recovery column is a renewable N fertilizer.

31. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the stripping column comprises a first stripping column fluidly coupled to a second stripping column, wherein the process of stripping ammonia gas from wastewater takes place in both the first and the second stripping columns, and wherein a valve is operative to direct the flow of wastewater from the first stripping column to the second stripping column when the first stripping column reaches a predetermined volume.

32. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein a pressure is applied to the evaporated ammonia gas for liquefaction.

33. The apparatus for producing liquid ammonia (NH.sub.3) from wastewater containing ammonium (NH.sub.4.sup.+) of claim 19, wherein the recovery column comprises a first recovery column fluidly coupled to a second recovery column, wherein the process of recovering ammonia gas in the recovery solution takes place in both the first and the second recovery columns, and wherein a valve is operative to direct the flow of ammonia gas from the first recovery column to the second recovery column when the saturation level of ammonia gas in the recovery solution of the first recovery column reaches a predetermined level.

34. A method of producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+), comprising the steps of: supplying wastewater, a Brnsted base, a Brnsted acid and water; supplying a stripping column equipped with one or more air diffusers and a recovery column equipped with one or more air diffusers; the stripping column and recovery column, being in fluid communication; mixing, in the stripping column, the wastewater with the Brnsted base to form a stripping mixture, the Brnsted base being operative to react with ammonium ions present in the wastewater to produce ammonia gas and other gasses in the stripping mixture; stripping, in the stripping column, the ammonia from the stripping mixture by aerating the stripping mixture to produce ammonia gas and other gasses; and recovering, in the recovery column, the ammonia gas by absorbing the ammonia gas into a recovery solution comprising the Brnsted acid and water, the recovery solution and ammonia gas absorbed therein defining an ammonia product.

35. The method of producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 34, wherein the ammonia product is an N fertilizer having a nitrogen concentration of at least 15%.

36. The method of producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 34, wherein the value of a liquid-to-gas mass transfer coefficient (K.sub.L,s) multiplied by an interface area per volume of liquid (a.sub.s) is not be lower than $\frac{Q_{G,s}{H}_s}{V_{L,s}}$ where Q.sub.G,s, V.sub.L,s and H.sub.s are the air flowrate, the volume, and Henry's law constant, respectively.

37. The method of producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 34, wherein a gas-to-liquid mass transfer (K.sub.L,r) multiplied by an interface area per volume of liquid (a.sub.r) is not be lower than $\frac{Q_{G,r}{H}_r}{V_{L,r}}$ where Q.sub.G,r, V.sub.L,r and H.sub.r are the air flowrate, the volume, and Henry's law constant, respectively.

38. An apparatus for producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) comprising: a stripping column equipped with one or more air diffusers and a recovery column equipped with one or more air diffusers; the stripping column and recovery column being in fluid communication; wherein, in the stripping column wastewater is mixed with a Brnsted base to form a stripping mixture, the Brnsted base being operative to react with ammonium ions present in the wastewater to produce ammonia gas and other gasses in the stripping mixture; wherein, the one or more air diffusers are located in the base of the stripping column and function to aerate the stripping mixture to release the ammonia gas and other gasses contained therein; and wherein, in the recovery column, the ammonia gas is absorbing into a recovery solution comprising a Brnsted acid and water, the recovery solution and ammonia gas absorbed therein defining an ammonia product.

39. The apparatus for producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 38, wherein the ammonia product is an N fertilizer having a nitrogen concentration of at least 15%.

40. The apparatus for producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 38, wherein at least one of the one or more air diffusers has an orifice diameter of 10 mm or less and a surface area per length in the range of 0.5-3 inch.sup.2 inch.sup.1.

41. The apparatus for producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 38, wherein the number of the one or more air diffusers in the stripping column or recovery columns is the square of the diameter of the respective column in inches.

42. The apparatus for producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 38, wherein at least one of the one or more air diffusers in the stripping column or the recovery columns is a sparger pipe diffuser, and wherein the number of the one or more air diffusers is equal to the diameter of the respective column in inches.

43. The apparatus for producing an ammonia product from wastewater containing ammonium (NH.sub.4.sup.+) of claim 38, wherein at least one of the one or more air diffusers in the stripping column or recovery columns is a sparger pipe diffuser, and wherein the number of the one of more air diffusers is equal to one inch more than the value of the diameter of the respective column diameter in inches.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0048] The foregoing aspects and the attendant advantages of the present methods and systems will become more readily appreciated by reference to the following detailed description, when taken in conjunction with the accompanying drawings, wherein:

[0049] FIG. 1 is a diagram of the steps of an exemplary NH.sub.3 recovery process, using flushed manure as the feed wastewater;

[0050] FIG. 2 illustrates a schematic of an example of the production process of a saturated solution by using a stripping column and a recovery column;

[0051] FIG. 3 displays a schematic of an example of the process of using the saturated solution from the process of FIG. 2 to produce bioammonia by heating and liquefaction:

[0052] FIG. 4 is a perspective view an embodiment of a stripping column packed with sub columns, each of which is equipped with an air diffuser at the bottom of the column;

[0053] FIG. 5 is a perspective view of one embodiment of a recovery column which incorporates a partition and wall baffles;

[0054] FIG. 6 is an exemplary configuration of a semi-continuous process consisting of a pair of stripping columns and a pair of recovery columns;

[0055] FIG. 7A is a schematic view showing an arrangement of sintered stone air diffusers and sparger pipe diffusers which may be incorporated into the bottom of a stripping or recovery column;

[0056] FIG. 7B is a schematic view showing another arrangement of sintered stone air diffusers and sparger pipe diffusers which may be incorporated into the bottom of a stripping or recovery column;

[0057] FIG. 8 illustrates a setup for a bench-scale stripping experiment testing a prior art method;

[0058] FIG. 9 shows the results of the pH of a carbonate solution according to the step of FIG. 8 as a function of time during aeration;

[0059] FIG. 10 displays the results of an experiment according to the setup of FIG. 8 in which the NH.sub.4.sup.+ concentration is plotted as a function of time during the stripping process.

[0060] FIG. 11 shows the results of an experiment in which the NH.sub.4.sup.+ concentration is plotted as a function of time during the stripping using two columns with differing column diameters;

[0061] FIG. 12 presents the results of an experiment in which the NH.sub.4.sup.+ concentration is plotted as a function of time during the stripping using two types of air diffusers that differ in orifice diameter;

[0062] FIG. 13 illustrates an experimental setup for a bench-scale stripping and recovery experiment according to the methods and systems disclosed herein;

[0063] FIG. 14 shows the experimental data according to the experimental setup of FIG. 13 in which the NH.sub.4.sup.+ concentration is plotted as a function of time during the stripping process under two conditions, one according to a prior art method and one according to the methods and systems disclosed herein;

[0064] FIG. 15 shows the results of an experiment in which the NH.sub.4.sup.+ concentration is plotted as a function of time during the NH.sub.3 absorption process into the recovery solution under two conditions, one according to a prior art method and one according to the methods and systems disclosed herein;

[0065] FIG. 16 plots the experimental results of the dependence of the normalized NH.sub.4.sup.+ concentration on the stripping optimization variables K.sub.L,sa.sub.s and

[00001]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

as a function of time for the stripping process

[0066] FIG. 17 plots the experimental results of the dependence of the normalized NH.sub.4.sup.+ concentration on the recovery optimization variables K.sub.L,ra.sub.r and

[00002]$\frac{Q_{G,r}{H}_r}{V_{L,r}}$

as a function of time for the absorption process;

[0067] FIG. 18 illustrates an exemplary schematic of an apparatus to produce liquid NH.sub.3 from the aerated NH.sub.3 gas according to a prior art method; and

[0068] FIG. 19 illustrates one embodiment of an apparatus to produce liquid NH.sub.3 from the evaporated NH.sub.3 gas from a saturated solution according to the methods and systems disclosed herein.

[0069] Each figure has reference figures that are unique to the figure they are associated with, and as such like reference numerals to not necessarily correspond to like elements from one figure to the next.

DETAILED DESCRIPTION

[0070] The detailed description set forth below is intended as a description of presently preferred embodiments of the methods and systems contemplated herein and is not intended to represent the only forms in which the presently contemplated systems and methods may be constructed and/or utilized. However, it is to be understood that the same or equivalent functions and results may be accomplished by different embodiments that are also intended to be encompassed within the spirit and scope of the presently contemplated subject matter, and additional variations of the present subject matter may be devised without departing from the inventive concept. The description itself is not intended to limit the scope of any patent issuing from this description. Rather, the inventors have contemplated that the claimed subject matter might also be embodied in other ways, to include different elements or combinations of elements similar to the ones described in this document, in conjunction with other present or future technologies.

[0071] An embodiment of the present disclosure includes a method to produce bioammonia by recovering NH.sub.3 from wastewater. Wastewater may be a mixture that contains ammonia and/or ammonium ions, and as such it may include, but is not limited to, industrial wastewater, municipal wastewater, livestock manure liquid, anaerobic manure digestate liquid, non-manure based anaerobic digestate liquid, and combinations thereof.

[0072] A system to produce bioammonia according to the methods disclosure herein may comprise two types of columns: one for the NH.sub.3 stripping where the step of stripping ammonia gas may take place, to be referred to as the stripping column, and another for the NH.sub.3 absorption into a recovery solution where the step of recovering the stripped ammonia gas through absorption may take place, to be referred to as the recovery column,. The product may be a saturated solution which can be used as is, or alternatively processed further into bioammonia. The system may further comprise an evaporator and a liquefaction device which may act upon the saturated solution produced from the recovery column to produce this desirous bioammonia product.

[0073] Before detailing these systems and methods, models for the mass-transfer reaction kinetics of the NH.sub.3 stripping and absorption processes will be developed and discussed first. Then, those models will be used to optimize the various operating parameters so as to describe how the systems and methods herein may be arrived and the preferred embodiments. Acid-base reactions for both NH.sub.3 stripping and absorption processes will expressed as follows, respectively in Equations 1 and 2:

NH.sub.4.sup.++XB.fwdarw.NH.sub.3+BH+X.sup.+(1)

and NH.sub.3+AH.fwdarw.NH.sub.4.sup.++A.sup.(2)

where XB and AH refer to a Brnsted base and acid, respectively. XB is preferably a strong base, including, but is not limited to, X(OH).sub.n, XCO.sub.3, and XHCO.sub.3, where X can be, but is not limited to, Ca, K, and Na, and n=1 or 2 depending on the counter cation (as would be understood by those skilled in the art). AH is preferably a strong acid including, but is not limited to, sulfuric acid, nitric acid, phosphoric acid, and organic acids. As will soon be shown, the Brnsted base can be mixed with the wastewater fed to the stripping column and the Brnsted acid can be mixed with the stripped ammonia gas fed to the recovery column to substantially increase the effectiveness of the processes associated with both of those columns. NH.sub.4.sup.+ in Equations 1 and NH.sub.3 in Equations 2 may act as an acid and a base respectively. Hence, Equations 1 and 2 are acid-base reactions. Acid-base reactions are known to be thermodynamically favorable, exothermic, and kinetically fast, provided sufficient concentrations of the base (XB) and acid (AH). Explicitly incorporating the acid-base reactions to capture NH.sub.3 in the wastewater and absorb the captured NH.sub.3 through the reactions, hence increasing the efficiency of the NH.sub.3 stripping and absorption processes simultaneously, is understood to be novel. The NH.sub.3 gas may be stripped by Equation 1 via aeration of a stripping mixture in the stripping column, and this stripped NH.sub.3 gas may be absorbed into an acidic recovery solution by Equation 2 in the recovery column.
The removal/recovery yield (.sub.s/r) and the removal/recovery rates (.sub.s/r) will be defined as follows in Equations 3-5:

.sub.s=100{1([NH.sub.4.sup.+].sub.s.sup.f)/[NH.sub.4.sup.+].sub.s.sup.i)}(3)

.sub.r=100([NH.sub.4.sup.+].sub.r.sup.f)/[NH.sub.4.sup.+].sub.s.sup.i) (4)

.sub.s/r=a([NH.sub.4.sup.+].sub.s/r.sup.i[NH.sub.4.sup.+].sub.s/r.sup.f)/t.sub.p (5)

where [NH.sub.4.sup.+] represents the NH.sub.4.sup.+ concentration, the subscripts s and r refer to the NH.sub.4.sup.+ concentration in the stripping step and that in the recovery step, respectively. The superscripts i and f refer to the initial and final concentrations of NH.sub.4.sup.+ respectively, and as such [NH.sub.4.sup.+].sub.s.sup.f[NH.sub.4.sup.+].sub.s.sup.i and [NH.sub.4.sup.+].sub.r.sup.f[NH.sub.4.sup.+].sub.s.sup.i. a=1 and 1 for the stripping and recovery processes, respectively. The removal/recovery yield (.sub.s/r) and the removal/recovery rates (.sub.s/r) from these Equations 3 through 5 will be referred to as performance parameters. These performance parameters will ideally be maximized so that the largest amount of ammonia can be extracted from a given amount of wastewater in the quickest amount of time. t.sub.p is the time required for the NH.sub.4.sup.+ concentration to become a plateau, or approximately reach steady state.

[0074] These systems and methods may yield particularly effectual results by optimizing one or more operating parameters. These operating parameters may include, but are not limited to, the concentrations of the Brnsted base in the stripping column, the concentration of the Brnsted acid in the recovery column, the orifice diameter of one or more aeration diffusers in the stripping and/or recovery columns, the diameter of the stripping and/or recovery column for a given column height, the surface area per length of the one or more air diffusers, the number of air diffusers per column, the number of air blowers per column, the ratio of the air flowrate to the volume in each column, and more. The performance parameters, .sub.r and .sub.s/r, may ideally be maximized using the mass-transfer reaction kinetics models and adjusting the operating parameters.

[0075] For both cases, Equations 1 and 2, the chemical reactions may occur very fast; hence, the processes may be primarily controlled by the more physical process of the mass transfer of NH.sub.3 between the liquid state and the air bubbles. Then, Equations 1 and 2 become a liquid-to-gas transfer process: the transfer of NH.sub.3 from the NH.sub.4.sup.+ solution to the air bubble in the stripping column for Equations 1 and the transfer of NH.sub.3 from the air bubbles to an aqueous acid solution for Equations 2.

[0076] Using the two-film theory, the mass balance of NH.sub.3 from liquid to gas may be written as follows in Equation 6:

[00003]$\begin{array}{cc}-V_{L,s}\frac{{d\left[N{H}_3\right]}_L}{dt}=F_{L.fwdarw.G}& \left(6\right)\end{array}$

where [NH.sub.3].sub.L, V.sub.L,s, and F.sub.L.fwdarw.G refer to the NH.sub.3 concentration in the liquid, the volume of stripping solution, and the mass transfer rate of NH.sub.3 from the solution to the gas bubbles generated by air diffusers, respectively, according to Matter-Mller et al. [Matter-Mller, C.; Gujer, W.; Giger, W. Transfer of Volatile Substances from Water to the Atmosphere, Water Research, 15, 1271-1279 (1981).].

[0077] According to Matter-Mller, F.sub.L.fwdarw.G is defined as follows in Equation 7:

[00004]$\begin{array}{cc}{F}_{L.fwdarw.G}=Q_{G,s}{H_s\left[N{H}_3\right]}_L\left\{1-\exp \left(-\frac{K_{L,s}{a}_s{V}_{L,s}}{H_s{Q}_{G,s}}\right)\right\}& \left(7\right)\end{array}$

where Q.sub.G,s, H.sub.s, K.sub.L,s, and a.sub.s represent the air flowrate at the inlet to the stripping column, Henry's law constant for NH.sub.3 in the stripping column, the mass transfer coefficient of NH.sub.3 from liquid to the air bubble, and the interface area per unit volume of liquid in the stripping column, respectively. [Matter et al.]

[0078] It can be seen that in Equation 1, increased [NH.sub.3].sub.L could be achieved by the Brnsted base shifting the reaction of Equation 1 to the right side via the base taking a proton away from NH.sub.4.sup.+, thus raising the value of F.sub.L.fwdarw.G.

[0079] H.sub.s is defined by Matter et al. as follows in Equation 8:

[00005]$\begin{array}{cc}{H}_s=\frac{{\left[N{H}_3\right]}_G}{{\left[N{H}_3\right]}_L}& \left(8\right)\end{array}$

where [NH.sub.3].sub.G and [NH.sub.3].sub.L represent the NH.sub.3 concentrations in the gas and liquid phases, respectively. H.sub.s and K.sub.L,s depend on pH and temperature, the latter of which may ideally be the ambident temperature of the system (Li, S.; Fan, J.; Xu, S.; Li, R.; Luan, J. The Influence of pH on Gas-Liquid Mass Transfer, Chem. Ind. Chem. Eng. Q. 23 (3) 321-327 (2017); Wilson, G. M. A New Correlation of NH.sub.3, CO.sub.2, and H.sub.2S Volatility Data from Aqueous Sour Water Systems, EPA-600/2-80-067, 1980; Kim, E. J.; Kim, H.; Lee, E. Influence of Ammonia Stripping Parameters on the Efficiency and Mass Transfer Rate of Ammonia Removal, Appl. Sci. 11, 441 (2021)).

[0080] [NH.sub.4.sup.+] is more easily measured, especially in a stripping column, when compared to [NH.sub.3].sub.L. Hence, [NH.sub.3].sub.L will be replaced by [NH.sub.4.sup.+] in Equation 6.

[0081] Solving the differential equation of Equation 6 with the boundary conditions of [NH.sub.4.sup.+]=[NH.sub.4.sup.+].sup.t at t=t and [NH.sub.4.sup.+]=[NH.sub.4.sup.+].sup.0 at t=0 gives the following Equation 9:

[00006]$\begin{array}{cc}-\ln \left(\frac{{\left[N{H}_4^{+}\right]}^t}{{\left[N{H}_4^{+}\right]}^0}\right)=\frac{Q_{G,s}{H}_s}{V_{L,s}}\left\{1-\exp \left(-\frac{K_{L,s}{a}_s{V}_{L,s}}{H_s{Q}_{G,s}}\right)\right\}t& \left(9\right)\end{array}$

[0082] Plotting the negative logarithm of the concentration ratio in Equation 9 against t gives a straight line with the following slope in Equation 10:

[00007]$\begin{array}{cc}{}_s=\frac{Q_{G,s}{H}_s}{V_{L,s}}\left\{1-\exp \left(-\frac{K_{L,s}{a}_s{V}_{L,s}}{H_s{Q}_{G,s}}\right)\right\}& \left(10\right)\end{array}$

where .sub.s is the mass transfer rate per unit volume for the stripping process. This value is ideally maximized, as a higher value corresponds to a faster and thus more efficient rate of stripping ammonia gas. Therefore, from experiments under different conditions using different sets of operating parameters, the right-side values of Equation 10 are obtained for a given set of V.sub.L,s, Q.sub.G,s, H.sub.s, a.sub.s, and K.sub.L,s. Equation 10 teaches us that in order to increase .sub.s, the air flowrate and a.sub.s should be higher, and the Brnsted base should have a higher pH to increase the concentration of NH.sub.3 in the liquid, which in turn raises both H.sub.s, and K.sub.L,s.

[0083] There are five variables in Equation 10: V.sub.L,s, Q.sub.G,s, H.sub.s, a.sub.s, and K.sub.L,s. These variables can be consolidated to two combined variables:

[00008]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

and K.sub.L,sa.sub.s. Accordingly, the operating parameters can be optimized to increase the value of these two combined variables. These combined variables may be referred to as stripping optimization variables.

[0084] As would be understood by those skilled in the art, K.sub.L,s. a.sub.s can be increased by generating a larger number of smaller air bubbles per volume. Smaller air bubbles tend to have a longer retention time in a column which is favorable for both stripping and absorption of NH.sub.3. The orifice diameter of air diffuser may primarily determine the size of air bubbles. The number of air bubbles per volume may depend on the number of diffuser orifices across the column they are associated with. Hence, the diameter of a column also matters for a given column height, as a larger diameter means that more air diffusers can be used. In general, a higher ratio of the column height to the diameter tends to give rise to an increased axial dispersion coefficient (Alvar, J.; Al-Dahhan, M. H., Liquid Phase Mixing in Trayed Bubble Column Reactors, Chem. Eng. Scien., 61, 1819 (2006)). a.sub.s also may depend on the surface area per length of the diffuser, provided that the surface area is covered by the orifice pores.

[0085] As to

[00009]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

Q.sub.G,s should ideally be maximized for a given V.sub.L,s. In reality, V.sub.L,s is often determined by the flowrate of wastewater at a plant or a livestock farm. Q.sub.G,s, then, can be chosen to maximize

[00010]$\frac{Q_{G,s}}{V_{L,s}}.$

Still, the ratio of the air flowrate to the volume is determined more or less based primarily on the economy and the amount of wastewater feed available, since a high Q.sub.G,s can be costly.

[0086] That leaves H.sub.s as a more easily adjustable parameter. It is well-known in the art that the equilibrium between [NH.sub.3].sub.G and [NH.sub.3].sub.L shifts towards [NH.sub.3].sub.G at higher values of pH. When enough base is present in the stripping step to shift Equation 1 completely to the right side, a sudden increase in [NH.sub.3].sub.L raises [NH.sub.3].sub.G rapidly, increasing H.sub.s quickly in the beginning of the stripping. The high concentration of the base has a similar effect on K.sub.L,s. This cascade process, along with the increased K.sub.L,s, helps .sub.s increase immediately. The high mass transfer rate, in the beginning of this process, helps reduce the bulk of the NH.sub.4.sup.+ concentration in the stripping column within a short period, promoting the high removal rate in mg L.sup.1 h.sup.1. Accordingly, the optimization variables can be optimized by changing the aforementioned operating parameters to maximize .sub.r and .sub.r.

[0087] Since Equation 10 establishes the relationship between the mass transfer rate and the stripping optimization variables, which are functions of the operating parameters, it can not only provide guidance for the ideal of the operating parameters, but also significantly reduce the potentially lengthy time required to ascertain the ideal operating parameters to result in an increased mass transfer rate. Through rearrangement, Equation 9 becomes Equation 11 as follows:

[NH.sub.4.sup.+].sup.t=[NH.sub.4.sup.+].sup.0e.sup..sup.s.sup.t (11)

[0088] Equation 11 describes the reduction of [NH.sub.4.sup.+] as a function of time by aeration based on Equation 1. Still, Equation 11 is for an ideal situation where all the initial NH.sub.4.sup.+ concentration is removed, which may not always be the case. Then, Equation 11 may be written more realistically as follows in

[00011]$\begin{array}{cc}\mathrm{Equation}12& \\ \left[{\mathrm{NH}}_4^{+}\right]={\left[{\mathrm{NH}}_4^{+}\right]}^0\left\{1-\frac{{}_s}{100}\left(1-e^{-{}_{s}t}\right)\right\}& \left(12\right)\end{array}$

[0089] Equation 12 gives

[00012]${\left[N{H}_4^{+}\right]}^t={\left[N{H}_4^{+}\right]}^0\mathrm{at}t=0{\mathrm{and}\left[{\mathrm{NH}}_4^{+}\right]}^t={\left[N{H}_4^{+}\right]}^0\left(1-\frac{{}_s}{100}\right)$

at t.fwdarw. both of which are consistent with Equation 11 when .sub.s is equal to 100%. .sub.s primarily may depend on the condition of the solution from which NH.sub.4.sup.+ is removed. For example, if a sufficient mass of the Brnsted base is supplied, Equation 1 should go completely to the right side, leaving little NH.sub.4.sup.+ concentration in the stripping step. On the other hand, .sub.s as determines the rate of the stripping. If the experiment follows Equation 12, the information on .sub.s can be obtained for each set of the stripping optimization variables K.sub.L,sa.sub.s and

[00013]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

from the experiments by the same procedure as described above by using the following Equation 13:

[00014]$\begin{array}{cc}\ln \left[\left(\frac{{}_s}{100}\right)\left(\frac{{\left[N{H}_4^{+}\right]}^t}{{\left[N{H}_4^{+}\right]}^0}-1+\frac{{}_s}{100}\right)\right]=-{}_{s}t& \left(13\right)\end{array}$

[0090] Plotting the left hand of Equation 13 against t gives a straight line with a slope, .sub.s. After a number of experiments to cover enough ranges of the stripping operating parameters within reasonable ranges, a relationship between .sub.s and the stripping operating parameters can be established using Equations 5, 10, 12, and 13. From the relationship, the information on the concentration of the Brnsted base, the orifice diameter of the air diffuser, the number of air diffusers in the column, the ratio of the column diameter over the height, and the surface area per length of the air diffusers, and the ratio of the air flowrate over the liquid volume can be obtained for the stripping process. Using those relationships, Equations 3, 5, 10, 12, and 13, it is possible to find a set of stripping operating parameters that gives the best set of .sub.s and .sub.s.

[0091] The same analysis will now be performed on the step of recovering the stripped ammonia gas by absorption. As to Equation 2, the gas-to-liquid mass transfer may be written as follows in Equation 14:

[00015]$\begin{array}{cc}-V_{L,r}\frac{d\left[N{H}_4^{+}\right]}{dt}=F_{G.fwdarw.L}& \left(14\right)\end{array}$

where V.sub.L,r and F.sub.G.fwdarw.L refer to the volume of the recovery solution used and the mass transfer rate of stripped NH.sub.3 in the gas bubbles to NH.sub.4.sup.+ in the recovery solution, respectively. The latter may be written as follows as Equation 15:

[00016]$\begin{array}{cc}\begin{array}{c}{F}_{G.fwdarw.L}=Q_{G,r}{H}_r({\left[N{H}_4^{+}\right]}_{sat}-\\ \left[{\mathrm{NH}}_4^{+}\right])\left\{1-\exp \left(-\frac{K_{L,r}{a}_r{V}_{L,r}}{H_r{Q}_{G,r}}\right)\right\}\\ =Q_{G,r}\left(\frac{{}_r}{100}{\frac{{}_s}{100}\left[N{H}_4^{+}\right]}_s^i-\left[N{H}_4^{+}\right]\right)\\ \left\{1-\exp \left(-\frac{K_{L,r}{a}_r{V}_{L,r}}{H_r{Q}_{G,r}}\right)\right\}\end{array}& \left(15\right)\end{array}$

Where Q.sub.G,r, H.sub.r, K.sub.L,r, and a.sub.r represent the gas flowrate at the inlet to the stripping column, Henry's law constant for NH.sub.3 in the recovery column, the mass transfer coefficient of NH.sub.3 from the air bubble to the liquid state, and the interface area per unit volume of liquid of the recovery solution, respectively.

[0092] At a low pH, [NH.sub.4.sup.+] increases, while [NH.sub.3].sub.G decreases, which in turn accelerates F.sub.G.fwdarw.L. Accordingly, Henry's constant is expressed as follows:

[00017]$\begin{array}{cc}{H}_r=\frac{{\left[N{H}_3\right]}_L}{{\left[N{H}_3\right]}_G}& \left(16\right)\end{array}$

[0093] [NH.sub.4.sup.+].sub.sat is the NH.sub.4.sup.+ concentration at the saturation point or the final NH.sub.4.sup.+ concentration in the recovery column which is replaced by

[00018]$\frac{{}_r}{100}{\frac{{}_a}{100}\left[N{H}_4^{+}\right]}_s^i$

using Equation 4. For Q.sub.G,r and a.sub.r, the same argument can be made as for Q.sub.G,s and a.sub.s, as discussed above.

[0094] As to H.sub.r and K.sub.r, a Brnsted acid with a low pH may be used to increase these variables. What Equation 2 enables us to recognize is that increasing the proton concentration at the interface between the liquid and the air bubbles is beneficial since when the NH.sub.3 gas molecules inside the bubbles have contact with the interface, they can immediately react with the protons outside the bubbles to be dissolved as NH.sub.4.sup.+, hence raising F.sub.G.fwdarw.L.

[0095] Solving the differential equation of Equation 14 with the boundary conditions of

[00019]${\left[N{H}_4^{+}\right]}^t=\frac{{}_r}{100}{\frac{{}_a}{100}\left[N{H}_4^{+}\right]}_s^i$

at t.fwdarw. and [NH.sub.4.sup.+].sup.t=0 at t=0 gives the following Equation 17:

[00020]$\begin{array}{cc}{\left[N{H}_4^{+}\right]}^t=\frac{{}_r}{100}{\frac{{}_a}{100}\left[N{H}_4^{+}\right]}_s^i\left(1-e^{-{}_{r}t}\right)& \left(17\right)\end{array}$

where .sub.r is expressed as follows in Equation 18:

[00021]$\begin{array}{cc}{}_r=\frac{Q_{G,r}{H}_r}{V_{L,r}}\left\{1-\exp \left(-\frac{K_{L,r}{a}_r{V}_{L,r}}{H_r{Q}_{G,r}}\right)\right\}& \left(18\right)\end{array}$

where .sub.r is the mass transfer rate per unit volume for the recovering process. Equation 17 combines the stripping and absorption processes for the reaction kinetics of the multiphase liquid-to-gas and then the gas-to-liquid mass-transfers by one equation. [NH.sub.4.sup.+].sub.r.sup.t thus depends not only the kinetics of the absorption, but also that of the stripping. Q.sub.G,r, H.sub.r, and K.sub.L,r of .sub.r in Equation 18 are also functions of the kinetics in the stripping. Now the NH.sub.3 recovery can be properly described, which is a result of not only the mass transfer in the stripping process, but also that of the absorption process. Thus, from experiments under different conditions using different sets of operating parameters, the values of the right side of Equation 18 may be obtained for a given set of V.sub.L,r, Q.sub.G,r, H.sub.r, a.sub.r, and K.sub.L,r. Equation 18 teaches us that in order to increase .sub.r, the ratio of the gas flowrate to the volume, a.sub.r should be preferably high, and the Brnsted acid should preferably have a low pH in order to increase the concentration of protons in the liquid, which in turn raises both H.sub.r and K.sub.L,r.

[0096] There are five variables in Equation 18: V.sub.L,r, Q.sub.G,r, H.sub.r, a.sub.r, and K.sub.L,r. These variables can be consolidated to two combined variables:

[00022]$\frac{Q_{G,r}{H}_r}{V_{L,r}}$

and K.sub.L,ra.sub.r. Accordingly, the operating parameters may be optimized to increase these two combined variables. This set of combined variables may be referred to as recovery optimization variables.

[0097] As would be understood by those skilled in the art, K.sub.L,r. a.sub.r can be increased by generating a large number of small air bubbles per volume. Small air bubbles may have a longer retention time in the column which is favorable for the absorption of NH.sub.3 into the acidic recovery solution. The orifice diameter of air diffuser may primarily determine the size of air bubbles. The number of air bubbles per volume may depend on the number of air diffuser orifices used across the recovery column. Hence, the diameter of the recovery column matters for a given column height, as a larger diameter gives more space for more air diffusers. a.sub.r also may depends on the surface area per length of the diffuser, provided that the surface area is covered by the orifice pores.

[0098] As to

[00023]$\frac{Q_{G,r}{H}_r}{V_{L,r}},$

Q.sub.G,r may be maximized for a given V.sub.L,r. In reality, V.sub.L,r is often determined by the flowrate of wastewater at a plant or a livestock farm. Q.sub.G,r, then, is chosen to maximize

[00024]$\frac{Q_{G,r}}{V_{L,r}}.$

Still, the ratio of the air flowrate to the volume is determined more or less based primarily on the economy, since a high Q.sub.G,r can be costly.

[0099] That leaves H.sub.r as a more easily adjustable parameter. It is well-known in the art that the equilibrium between [NH.sub.3].sub.G and [NH.sub.3].sub.L shifts towards [NH.sub.3].sub.L at lower values of pH. When enough acid is present in the recovery column to shift Equation 2 completely to the right side, a sudden increase in [NH.sub.3].sub.L may raise Hr rapidly at the beginning of the recovery. The high concentration of the acid has a similar effect on K.sub.L,r. This process, along with the increased, K.sub.L,r, may accelerate the value of .sub.r immediately. The high mass transfer rate in the beginning of the process helps reduce the bulk of the NH.sub.4.sup.+ concentration in the recovery step within a short period, promoting the high removal rate in mg L.sup.1 h.sup.1.

[0100] Accordingly, the recovery optimization variables can be optimized by changing the operating parameters such as the concentration of the Brnsted acid, the orifice diameter of the air diffuser, the number of diffusers in the column, the ratio of the column diameter over the height, and the surface area per length of the diffuser, the ratio of the air flowrate over the liquid volume among others to maximize .sub.r and .sub.r.

[0101] Since Equation 18 establishes the relationship between the mass transfer rate and the recovery optimization variables, which are functions of the operating parameters, it can not only provide guidance for the ideal values of the operating parameters, but also significantly reduce the potentially lengthy time required to ascertain the ideal values of the operating parameters to increase the mass transfer rate.

[0102] Rearrangement of Equation 17 gives the following Equation 19:

[00025]$\begin{array}{cc}\ln \left[1-\frac{{\left[{\mathrm{NH}}_4^{+}\right]}^t}{{\frac{{}_r}{100}\left[{\mathrm{NH}}_4^{+}\right]}_s^i}\right]=-{}_{r}t& \left(19\right)\end{array}$

[0103] Plotting the left side of Equation 19 against t gives a straight line with the slope, .sub.r. As mentioned above, after a number of experiments to cover enough ranges of the operating parameters, a relationship between .sub.r or .sub.r and the operating parameters can be established, using Equations 4, 5, 17, 18 and 19.

[0104] From the relationship, the information on the concentration of the Brnsted acid, the orifice diameter of the air diffuser, the number of air diffusers in the column, the ratio of the column diameter over the height, and the surface area per length of the air diffusers, and the ratio of the air flowrate over the liquid volume can be obtained for the absorption process. Using the relationships, Equations 4, 5, 17, 18 and 19, it is possible to find a set of operating parameters that maximize .sub.r and .sub.r.

[0105] Now that these models have been developed, the present disclosure will turn to the Figures and describe how methods and systems for bioammonia production may be arrived at.

[0106] With reference to FIG. 1, in general terms, in step 130, wastewater 12 is introduced to the ammonia production system 10 of the present invention. In step 132, the wastewater is filtered by passing the wastewater through a solid-liquid separator. In step 134, the wastewater is pumped into a stripping column where the wastewater is aerated and a Brnsted base is added to the solution. In step 136, the aerated wastewater is pumped into a recovery column where stripped ammonia gas is absorbed into a recovery solution with the addition of a Brnsted acid, the Brnsted acid being operative to react with the stripped ammonia gas to produce an ammonia product, the mixing of the stripped ammonia gas with the recovery solution defining a nitrogen solution or the ammonia product. In step 138, the saturation level of NH.sub.3 is tested, if the level of NH.sub.3 is below a predetermined level, the solution is returned to the stripping column. If the level of NH.sub.3 is above a predetermined level, in step 140, the solution is pumped to an evaporator where an evaporation process takes place. Thereafter, in step 142, the evaporated ammonia product may be liquefied.

[0107] Turning now to FIGS. 1, 2, and 3, the steps of an exemplary process and the schematic system 10 for the recovery of NH.sub.3 according to the methods and systems discussed herein are illustrated. NH.sub.3 may be recovered from wastewater 12 and liquid bioammonia 62 may be subsequently produced.

[0108] Though these systems and methods can be applied to any wastewater 12 containing NH.sub.4.sup.+, here flushed manure liquid is used as an example. The wastewater 12 may comprise solid components and liquid components, and if so, the wastewater 12 may be filtered to at least partially separate the liquid components from the solid components. Wastewater filtered in this fashion may be referred to as filtered wastewater. The filtered wastewater can be obtained from a wastewater mixture by a solid-liquid separator 14 which can include, but is not limited to, a screw separator, a screen separator, a centrifuge, a rotary separator, and combinations thereof. The filtered wastewater can be optionally stored in a lagoon or a sedimentation tank 42 before the liquid skimmed from the surface is pumped into a stripping column 24. Alternatively, the filtered wastewater may be directly sent to the stripping column 24. A centrifugal pump 16 may facilitate this process. The total suspended solid (TSS) level of the filtered wastewater from the solid-liquid separator 14 is preferably as low as possible, ideally no more than 10,000 mg/L.

[0109] Once the wastewater is pumped into the stripping column 24, whether it is filtered wastewater or just the original wastewater, air is sent by an air pump or an air blower 22 to the stripping column 24 through one or more air diffusers 28 equipped at the bottom of the stripping column 24 for aeration in the stripping step. The one or more air diffusers (i) may be, but they are not limited to, sintered stone diffusers, spargers, membrane diffusers, and combinations thereof.

[0110] Q.sub.G,s of the air pump 22 should preferably be high.

[00026]$\left(\frac{Q_{G,s}}{V_s}\right)$

is ideally more than 1 t.sup.1, and the higher this value is, the faster the kinetics for the stripping becomes, according to Equation 10.

[0111] When the alkalinity of the wastewater fed to the stripping column 24 is high, as is the case for livestock manure such as dairy manure, the pH of the wastewater is normally high enough to transfer NH.sub.4.sup.+ to NH.sub.3 in the wastewater to some extent. At the standard condition, the shift to NH.sub.3 occurs above pH=7. The high alkalinity is often caused by a high concentration of CaCO.sub.3 or Ca(HCO.sub.3).sub.2 in the wastewater.

[0112] However, the NH.sub.3 stripping may cease once such chemicals are consumed by the following reaction of Equation 20:

2 NH.sub.4.sup.++CaCO.sub.3+H.sub.2O.fwdarw.2 NH.sub.3+Ca(OH).sub.2+CO.sub.2+2H.sup.+(20)

[0113] The invented process may use a Brnsted base in a stored in a hopper 26 to strip NH.sub.3 more completely from the flushed manure liquid by Equation 1. The Brnsted base may be mixed with the filtered wastewater inside the stripping column. When the Brnsted base is mixed with the filtered wastewater, the resulting solution may be referred to as a stripping mixture. The mole of the anion of the Brnsted base used may be stoichiometric to the mole of NH.sub.4.sup.+ in the wastewater 12, which may be measured beforehand with conventional measurement methods such as a UV/vis spectrometer. Preferably, an excess amount of the Brnsted base, at least 10% or more stoichiometric to the mole of NH.sub.4.sup.+ in the original wastewater 12, is desirably used to ensure the rapid mass-transfer kinetics and a more complete removal of NH.sub.4.sup.+. The Brnsted base can also be added as needed, which may be determined via a pH sensor 20 measuring the pH of the stripping mixture.

[0114] The air diffuser(s) 28 for the stripping column 24 may preferably be configured such that they generates small air bubbles for desirable diffusions of bubbles inside the stripping column (g) to ensure that the NH.sub.3 gas produced by Equation 1 is swiftly transported to the bubbles for stripping. The desirable diffusion can be made possible by small air bubbles with a slow rising velocity and a long retention time inside the stripping column 24. The orifice diameter for the air diffuser(s) 28 may be less than 50 m, preferably less than 5 m, to produce these small air bubbles.

[0115] The one or more air diffusers 28 that can be used for the stripping process include, but are not limited to, sintered stone diffusers, membrane diffusers, spargers, and combinations thereof; a higher surface area to orifice diameter ratio, which is often expressed as the surface area per length of a porous sparger, is another ideal property in these air diffusers. The preferred surface area per length of a porous diffuser is in the range of 0.5-3 inch.sup.2/inch.

[0116] The progress of the stripping process can be monitored with an NH.sub.4.sup.+ sensor 18. Once NH.sub.3 has been stripped to a desirable extent, the spent wastewater may be repurposed, such as by discharging it to a lagoon or a storage tank for spraying on croplands later. A new batch of wastewater may then be filtered and pumped into the stripping column and this cycle could be repeated until the product in the recovery column 32 is saturated.

[0117] The stripped NH.sub.3 gas 30, along with the N.sub.2, O.sub.2, and often CO.sub.2 gases, all of which may be generated by aeration in the stripping column 24, may be pumped by a pneumatic pump 40 into the recovery column 32. The piping that may transport these stripped gasses to the recovery column may incorporate a heat exchanger operative to raise the temperature of the stripped NH.sub.3 gas (although this is not depicted in FIG. 2). This may increase the kinetic energy of the stripped NH.sub.3 gas; NH.sub.3 gas molecules with a higher kinetic energy has more collision frequency against the bubble walls, hence increasing the concentration at the interface with the recovery solution, raising the mass transfer to the recovery solution in the recovery column. The stripped NH.sub.3 gas may be bubbled through air diffusers 28 equipped at the bottom of the recovery column 32. The recovery column 32 could be at least partially filled by an acidic recovery solution including water which may be operative to absorb and recover the stripped ammonia gas. The pneumatic pump 40 pumping the stripped ammonia gas can increase the flowrate, Q.sub.G,r, to the recovery column 32, which raises the mass-transfer kinetic rate of the gas-liquid NH.sub.3 transfer, according to Equation 18. Recovery solution 34 containing a Brnsted acid may be supplied as needed, which can be determined by the readings of a pH sensor (e).

[0118] Q.sub.G,r may preferably be such that

[00027]$\left(\frac{Q_{G,r}}{V_r}\right)$

is more than 1 t.sup.1, and the higher the number is, the faster the kinetics for the stripping becomes, according to Equation 18.

[0119] The one or more air diffusers 28 in the recovery column 32 may be configured so as to generate small air bubbles for desirable diffusions of bubbles inside the recovery column 32, ensuring that the air bubbles containing the NH.sub.3 inside have sufficient contact with the acid 34 in the recovery column 32. The desirable diffusion can be made possible by a slow rising velocity and a long retention time of the air bubbles inside the recovery column 32. The one or more air diffusers 28 for bubbling the NH.sub.3 gas into the acidic aqueous solution may have small pores operative to generate small bubbles which have a long retention time to ensure the diffusion of air bubbles throughout the recovery column 32. The orifice diameter for the one or more air diffusers 28 may be less than 50 m, preferably less than 5 m to achieve this result.

[0120] The one or more air diffusers 28 for the recovery process may include, but are not limited to, sintered stone diffusers, membrane diffusers, spargers, and combinations thereof. Having a higher surface area to orifice diameter ratio, which is often expressed as the surface area per length of a porous sparger, is an ideal property in these air diffusers. The preferred surface area per length of a porous diffuser may be in the range of 0.5-3 inch.sup.2/inch.

[0121] Although not depicted, a mesh screen can be used in the stripping column 24 or recovery column 32 to break up the rising air bubbles from the one or more air diffusers 28. It is known that as the air bubbles rise through a liquid column, their sizes grow. A mesh screen, otherwise known as bubble breakers, can reduce the bubble size (Kalbfleisch, A. The Effect of Mesh-Type Bubble Breakers On Two-Phase Vertical Co-Flow, Electronic Thesis and Dissertation Repository, 3946 (2016). https://ir.lib.uwo.ca/etd/3946). The size of the mesh screen can be less than 1 mm. The screen can help maintain the original bubble size through the column.

[0122] The aerated gases 30 coming from the stripping column may include not only NH.sub.3, but N.sub.2, O.sub.2, and possibly CO.sub.2 as a result of the aeration process in the stripping column 24. These latter gases can interfere with the stripped NH.sub.3's absorption into the recovery solution in the recovery column. The acid also may ensure the absorption of the NH.sub.3 gas through the acid-base reaction, as shown by Equation (2), by a sufficient concentration of protons at the interface between the air bubbles and the liquid which facilitates the transfer of the NH.sub.3 gas into the recovery solution.

[0123] The mass of the acid supplied to the recovery column 32 can be determined by the solubility of the acid in water. For example, strong acids such as sulfuric acid and nitric acid are completely miscible in water. Then, the mass of the acid is determined by the maximum solubility of the product in Equation 2, an ammonium salt. The NH.sub.3 gas may continually be sent to the recovery column 32 by replacing the wastewater in the stripping column 24 until the product concentration reaches its maximum in the recovery column 32. For example, the solubility of (NH.sub.4).sub.2SO.sub.4, the product when the acid is sulfuric acid, is 744 g L.sup.1 in water at room temperature. Accordingly, the concentration of NH.sub.4.sup.+ at the maximum solubility is 19%. The use of reaction in Equation 2 is to facilitate the absorption of NH.sub.3 into the water. The progress of the recovery process may be measured via an NH.sub.4.sup.+ sensor 18.

[0124] Once ammonium acid reaches the maximum solubility in the recovery column 32, the air bubbling in the stripping column may be stopped. Then the pump 40 should also be stopped. The resulting solution, which can be referred to as a saturated solution 38, can discharged from the column, and stored under ambient conditions. Depending on the species in the recovery solution, this saturated solution can differ in chemical composition and species produced via reactions in the recovery column 32. The product can be a highly concentrated nitrogen solution that can be used as is, such as a renewable nitrogen fertilizer.

[0125] Alternatively, the saturated solution 38 can be pumped into an evaporator 64 via a pump 66, as is illustrated in FIG. 3. Steam can be circulated through a steam inlet 44 and a steam outlet 46 inside the evaporator 64 to bring the temperature inside the evaporator to around 80 C. or until the saturated solution starts to boil. Alternative heating elements may be utilized as well, including heating coils. This temperature may cause NH.sub.3 in the saturated solution 38 to start to transfer from the liquid phase to the gas phase. This step may be continued until an NH.sub.4.sup.+ sensor 18 shows a desired small value of NH.sub.4.sup.+ concentration inside the evaporator 64.

[0126] The NH.sub.3 gas evaporated in the evaporator 64 may then be sent to the liquefaction device 52. A pipe connecting the evaporator 64 and the liquefaction device 52 may be attached to a coiled pipe 58, preferably made from stainless-steel, inside the liquefaction device 52 whereby the NH.sub.3 gas may be put into thermal contact with a coolant circulating inside the liquefaction device through a coolant inlet 54 and coolant outlet 56. As the NH.sub.3 gas travels through the coiled pipe 58, the NH.sub.3 gas molecules can undergo condensation if the temperature is below 33.3 C., the boiling temperature of NH.sub.3 at one atmospheric pressure. The condensed bioammonia liquid 62 may be collected in the liquid collector 68 and discharged from the liquid outlet 70. The other gases introduced to the liquefaction device 52, which may include N.sub.2 and O.sub.2, can be released from the liquefaction device 52 via a release valve 60. The air gases, along with CO.sub.2 gas, can cause substantial interference with the liquefaction of the NH.sub.3 gas, as will be demonstrated later in this disclosure, and may thus be desirously removed in this fashion.

[0127] Once those extraneous gasses are removed from the liquefaction device 52, the gas being sent from the evaporator 64 may mostly be the NH.sub.3 gas. Hence, there is less interference of the NH.sub.3 liquefaction arising from other gases. The melting point of NH.sub.3 is 77.7 C. at one atmospheric pressure. Hence, the temperature inside the liquefaction device may preferably stay between 70 C. and 40 C. at atmospheric pressure to prevent freezing. Accordingly, the coolant for the liquefaction device should stay as a liquid at temperatures between 70 C. and 40 C. The temperature range can thus preferably be between 60 C. and 50 C. The coolants that can be used include, but are not limited to, iso-propanol, acetone, 1- or 2-butanol, 2-butanone, ethanol, diethyl ether, heptane, n-hexane, pentane, 1-propanol, tetrahydrofuran, and triethyl amine. The liquefaction of the NH.sub.3 gas has been industrially performed for decades and the cryogenic technology associated with the liquefaction is well established. NH.sub.3 gasification can occur at 33.3 C. If the upper temperature inside the liquefaction device is close to 33 C., there is a risk of losing some of the recovered NH.sub.3. Further, temperatures below 70 C. is too close to the melting point of CO.sub.2 under atmospheric pressure, 78.46 C. Hence, there is a risk of contaminating liquid NH.sub.3 with some CO.sub.2 solid at these temperature ranges. Accordingly, the temperature range for the liquefaction of NH.sub.3 should ideally be set between 70 C. and 40 C., preferably between 60 C. and 50 C.

[0128] In another embodiment, the liquefaction of evaporated NH.sub.3 gas described above can be applied to any N solutions recovered from livestock manure or anaerobic digestate liquid by any means including conventional NH.sub.3 stripping/scrubbing processes or membrane recovery processes.

[0129] In another embodiment, the liquefaction of evaporated NH.sub.3 gas described above can be carried out by applying a pressure. For example, the NH.sub.3 gas can be liquefied at a pressure of 7.5 bar at 20 C.

[0130] In yet another embodiment, the concentration of the Brnsted base, the type of the base, the orifice diameter of the one or more air diffusers 28, the number of air diffusers in the column, the ratio of the column diameter to the height, the surface area per length of the one or more air diffusers, and the ratio of Q.sub.G,s to V.sub.L,s for the stripping process are optimized such that the mass transfer rate per unit volume, .sub.s, in the stripping process is at least 3 min.sup.1.

[0131] In yet another embodiment, the concentration of the Brnsted acid, the type of the acid, the orifice diameter of the one or more air diffusers 28, the number of air diffusers in the column, the ratio of the column diameter to the height, the surface area per length of the one or more air diffusers, and the ratio of Q.sub.G,r to V.sub.L,r for the recovery process are optimized such that the mass transfer rate per unit volume, , in the recovery process should be at least 0.006 min.sup.1.

[0132] In yet another embodiment, the concentration of the Brnsted base, the type of the base, the orifice diameter of the one or more air diffusers 28, the number of air diffusers in the column, the ratio of the column diameter to the height, the surface area per length of the diffuser, and the ratio of Q.sub.G,s to V.sub.L,s for the stripping process are optimized such that K.sub.L,s.Math.a.sub.s is larger than

[00028]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

in the mass transfer rate per unit volume for the stripping process.

[0133] In yet another embodiment, the concentration of the Brnsted acid, the type of the acid, the orifice diameter of the air diffuser 28, the number of diffusers in the column, the ratio of the column diameter to the height, the surface area per length of the diffuser, and the ratio of Q.sub.G,r to V.sub.L,r for the recovery process are optimized such that K.sub.L,rs.Math.a.sub.r is larger than

[00029]$\frac{Q_{G,r}{H}_r}{V_{L,r}}$

in the mass transfer rate per unit volume for the absorption process. Looking now to FIG. 4, an illustration of an exemplary stripping column is shown. A stripping column may comprise a number of sub columns 72 (1-2 in diameter) which may be tightly packed together inside a stripping column housing 74. Equation 10 shows that a smaller sub column 72 diameter gives a higher mass-transfer rate for a given number of orifices across a column by increasing the flux of air bubbles, hence raising a.sub.s. This design shown in FIG. 4 allows a high mass-transfer rate at a full scale where a high flowrate of wastewater needs to be treated. The stripping may be performed in each sub stripping column 72; each sub column 72 may be equipped with an air diffuser 28 at the bottom of each sub column 72. The air diffusers may receive the air from a pump 76 connected to the stripping column housing 74 via a pipe 78. Saturated solution from the sub-columns exits the column housing 74 at a common outlet 80. As will be shown under EXAMPLE 3, the stripping rate may rise when the diameter of the column is reduced.

[0134] Turning now to FIG. 5, an illustration of an exemplary recovery column 32 is shown. This recovery column may receive stripped ammonia gas 82 through a sparger diffuser 98. A partition plate 88 may be at the center of the column to promote liquid circulation inside the recovery column and to distribute air bubbles evenly through the recovery column. A number of slant baffles 86 which may be present on the column walls may be operative to slow down the rising velocity of air bubbles from the sparger diffuser 98. These features may thus enhance the capability of the recovery solution to recover the stripped ammonia gas 82 received by the recovery column.

[0135] Bringing our attention now to FIG. 6, another embodiment comprising a pair of stripping columns 24 and a pair of recovery columns 32 is shown. This embodiment may allow for an operation to be carried out in a semi-continuous fashion by using pairs of coupled stripping columns and recovery columns. When one of the paired stripping columns 24 is used for stripping, the other stripping column 24 may discharge the spent wastewater through a first discharge 90 before being washed for subsequent stripping with a new batch of wastewater. This switch between the two stripping columns can be carried out by a valve 92 that can change the direction of the flow of wastewater 12 being fed from a centrifugal pump 16). Likewise, when one of the pair of the recovery columns 32 is used to absorb NH.sub.3, the other recovery column may discharge the saturated solution through a second discharge 94 and be washed. Either the striping columns or the recovery columns or both may comprise a plurality of sub-columns 72. This design allows the system to be operated continuously without the need for a holding tank.

[0136] Looking now to FIGS. 7a and 7b, cross sectional views of two embodiments of a column showcasing their air diffusers 28 are shown. These embodiments can be used in the recovery column 32, the stripping column 24, or both in a certain system. These embodiments present how the number of air diffusers used in a column of a given size can be arrived at. The dots (.circle-solid.) and the straight lines (-) across the column bottoms represent sintered stone air diffusers 96 and sparger pipe diffusers 98, respectively. The diameter of the circles inside the column are all each 1. These circles are not columns, but imaginary circles for the discussion. In certain embodiments, stone air diffusers (.circle-solid.) 96 can be used, and the number of these air diffusers (N.sub.dif) for both columns, the stripping and the recovery columns, can be determined by the following equation:

N.sub.dif=n.sup.2 (21)

where n is an integer ranging from 1 and above representing the diameter of a given column in inches. If the diameter is a real number, the rounded number can be used. FIG. 7A illustrates the relationship between N.sub.diff and n for a 2-diameter column and FIG. 7B for a 3-diameter column. Accordingly, in FIG. 7A there are four sintered stone air diffusers 96 while in FIG. 7B there are nine sintered stone air diffusers 96.

[0137] In other embodiments, just sparger-pipe diffusers (-) 98 may be used for the aeration of both columns, and the number of these diffusers (N.sub.dif) per stripping column can be equal to the number of pipes as follows:

N.sub.dif=n when n is an even number (22)

N.sub.dif=n+1 when n is an odd number (23)

[0138] In FIG. 7a, it can be seen that there are two sparger-pipe diffusers 98, while in FIG. 7B there are four sparger-pipe diffusers 98. The number of air pumps which pump the air/stripped ammonia gas to these air diffusers may also correlate to the geometry of the column they are associated with. As an example, the number of air pumps (N.sub.pump) per column may be chosen by the following equations:

N.sub.pump=n (24)

when n is an even number and

N.sub.pump=n+1 when n is an odd number (25)

where n is the same as in Equation 21.

[0139] Once bioammonia is obtained by the process described above, another step can be performed in which it is mixed with the CO.sub.2 gas under high pressures (110 atm) at temperatures 60 C. to produce urea by the following reactions of Equations 26 and 27:

2NH.sub.3+CO.sub.2[H.sub.2NCO.sub.2][NH.sub.4](26)

[H.sub.2NCO.sub.2][NH.sub.4](NH.sub.2).sub.2CO+H.sub.2O (27)

The CO.sub.2 gas used in Equation 26 can includes, but it is not limited to, the CO.sub.2 gas generated by an anaerobic digestor using livestock manure and/or organic wastes in general as the feed, the CO.sub.2 gas generated by coal-fired power plants, any CO.sub.2-generating industrial plants, the CO.sub.2 gas produced as a by-product in the stripping column when alkali carbonate or alkali bicarbonate is used as the Brnsted base. For every ton of urea produced by this process, not only is 0.73 MT of CO.sub.2 consumed by Equation 21, but also 1.6 MT of CO.sub.2 is reduced by using bioammonia produced by these methods and systems in place of synthetic NH.sub.3 which produces 1.6 MT of CO.sub.2 for every MT of synthetic NH.sub.3 produced. [Ammonia: Zero-Carbon Fertilizer, Fuel, and Energy Store, The Royal Society, February, 2020.] Hence, the total of 2.33 MT of CO.sub.2 emissions per 1 MT of urea can be eliminated by the urea synthesis Equation 27.

[0140] The difference from other NH.sub.3 recovery processes found in the literature may include the use of the stoichiometric chemical reactions or more than stoichiometric chemical reactions shown in Equations 1 and 2 to facilitate the NH.sub.3 stripping and absorption into the water and the optimization of the mass-transfer reaction kinetics of the NH.sub.3 stripping and absorption.

[0141] In another embodiment, when the alkalinity of wastewater is high enough, aeration of such wastewater can strip not only NH.sub.3 gas, but CO.sub.2 gas as well without addition of any base. When the stripped gases are absorbed in water, ammonium bicarbonate may be produced without any acid, according to the following equation:

NH.sub.3+CO.sub.2+H.sub.2O.fwdarw.NH.sub.4HCO.sub.3 (28)

The kinetic model described in this invention can be applied to facilitate the liquid-to-gas mass transfers of the NH.sub.3 and CO.sub.2 gases and the gas-to-liquid mass transfers of the gases to water.

EXAMPLE 1

Effects of Aeration on pH of Stripping Mixture

[0142] Aeration may be an important component in these processes and could be a key factor in driving the NH.sub.3 stripping in the stripping column. By comparing the results of experiments according to the systems and processes disclosed herein with a prior art method, specifically U.S. Pat. No. 11,364,463, the surprising improvement that the former may yield will be more clearly shown. In one of the embodiments in U.S. Pat. No. 11,364,463, it states that the aeration of digestate liquid generates CO.sub.2 gas when the alkalinity of the liquid is high, according to the following formula in Equation 28:

HCO.sub.3.sup.+air.fwdarw.CO.sub.2+OH.sup.(28)

The patent further asserts that the produced OH.sup. raises pH in the liquid, shifting the equilibrium of NH.sub.4.sup.30 to NH.sub.3. This is important to note because Equation 28 shows how NH.sub.3 may be stripped from the wastewater. We shall examine whether Equation 28 holds. Looking now to FIGS. 8, 9, and 10, the experiment setup of a prior art ammonia gas stripping method and the results of this experiment are illustrated. A 1 mM bicarbonate solution was prepared by mixing Ca(HCO.sub.3).sub.2 in water to simulate a high alkalinity solution. The pH value was 8.3. This solution was introduced to a stripping column 32 with a height of 40 and a diameter of 4. Sintered stone diffusers 96 with a 2 m orifice diameter were placed at the bottom of the stripping column 32 which were connected through a pipe to an air blower 22. The air flowrate was 30 L min.sup.1. A pH meter 20 was inserted into the solution to monitor the pH, as well as a NH.sub.4.sup.+18 sensor to measure the NH.sub.4.sup.+ concentration. The aeration continued for 10 hrs. FIG. 9 plots the pH of the solution as a function of time. According to Equation 28, pH should increase by the aeration process. Yet, the pH of the solution stayed flat throughout the aeration. This is due to the following equilibrium shown in Equation 29:

HCO.sub.3CO.sub.2+OH.sup.(29)

[0143] This equilibrium may be determined primarily by pH and temperature to a certain extent. Once pH is determined by the concentration of the bicarbonate, pH does not change, regardless of aeration, supported by our observation. The result demonstrates that aeration alone does not raise pH of the alkaline solution.

EXAMPLE 2

Prior Art Ammonia Gas Stripping

[0144] Next, the NH.sub.3 stripping process was conducted, using ADL separated by a solid-liquid separator as the filtered wastewater. The sample was taken from a dairy farm with 5,000 cows in the Central Valley, California. The ADL had the NH.sub.4.sup.+ and NO.sub.3.sup. concentrations of 2,086 mg L.sup.1, and 5.3 mg L.sup.1, respectively, the total nitrogen content of 2,300 mg L.sup.1, the alkalinity of 8,800 mg L.sup.1, the total solid and TSS of 9,300 mg L.sup.1 and 4,000 mg L.sup.1, respectively, phosphorus of 130 mg L.sup.1, and pH of 8.2. The same experimental setup as EXAMPLE 1 was used for this experiment.

[0145] First, a similar condition described in U.S. Pat. No. 11,364,463 was used for the stripping without applying the acid-base reaction Equation 1: pH in the original ADL was 9 in U.S. Pat. No. 11,364,463. Hence, the pH of the ADL was adjusted by adding NaOH to 9; however, no additional NaOH was supplied to maintain the same pH during the aeration. Though the temperature was set to be at 35 C. in U.S. Pat. No. 11,364,463, the temperature of ADL was ambient for all the examples presented in all the examples. The air flowrate was 30 L min.sup.1, and the volume of the digestate liquid was 2 L. For aeration, sintered stone air diffusers 96 with an orifice diameter of 50 m were used. The number of air diffusers was eight. The dimension of the stripping column was 40 high and 4 wide. The same volume of the digestate liquid, diffusers, and air flow were used for all the examples.

[0146] FIG. 10 shows the decay of [NH.sub.4.sup.+] as a function of time during aeration in the stripping column under the condition described above. Inside the figure, the equation corresponding to Equation 12 is included. .sub.s was determined by Equation 3, while .sub.s was obtained from the slope of the line by plotting the left side of Equation 9 against t. The values of R.sup.2 for the linear regressions in determining .sub.s and .sub.s were 0.996 and 0.997, respectively. The experimental data are well reproduced by the equation, displayed by the solid line, which suggests that the stripping process follows Equation 12. The standard deviations of calculated [NH.sub.4.sup.+] by the equations shown by the solid and the broken lines from the experimental data are 20.0 mg L.sup.1 and 19.8 mg L.sup.1, respectively. The removal yield was 50% which is close to what is described in U.S. Pat. No. 11,364,463. Only half the original NH.sub.3 was removed under this condition. As NH.sub.4.sup.+ donates its proton to OH.sup. to transform itself to NH.sub.3, OH.sup. is consumed. As less OH.sup. is available for the proton donation, NH.sub.4.sup.+ has no choice, but to stay as is without transforming itself to NH.sub.3. Once OH.sup. is consumed, the following reaction stops, and the equilibrium is reached via Equation 30:

NH.sub.4.sup.++OH.sup..fwdarw.>NH.sub.3+H.sub.2O (30)

That is when [NH.sub.4.sup.+] became nearly plateau after 250 min. This process wastes another half of NH.sub.3 potentially available for the recovery. .sub.s, p.sub.s, and .sub.s were 50.2%, 254.2 mg L.sup.1 h.sup.1, and 0.012 t.sup.1, respectively.

EXAMPLE 3

Effects of Adjusting the Diameter of the Stripping Column

[0147] This example examined the effect of one of the operating parameters, the stripping column diameter, on .sub.s and .sub.s. Two columns were tested: the diameters of 4 and 1. The height was the same for both columns: 40. Everything else was the same as EXAMPLE 2 except that pH of ADL was adjusted to be 12 by adding NaOH before aeration, but no additional NaOH was supplied during aeration. The same experimental setup as EXAMPLE 1 was used for this experiment.

[0148] Turning now to FIG. 11, the experimental results of the decay of [NH.sub.4.sup.+] as a function of time in the stripping columns with a 4 diameter () and a 1 diameter () during aeration are displayed. The equations corresponding to Equation 12 are included for each plot. .sub.s and .sub.s were determined as described above in EXAMPLE 2. The values of R.sup.2 for the linear regressions in determining .sub.s were 0.998 and 0.995 for the solid and broken lines, respectively. Both sets of experimental data are well reproduced by the equations, the broken line for the 4 diameter column and the solid line for the 1 diameter column. The standard deviations of calculated [NH.sub.4.sup.+] by the equations shown by the solid and the broken lines from the experimental data are 19.9 mg L.sup.1 and 19.8 mg L.sup.1, respectively. The agreements suggest that both stripping processes follow Equation 12. The column with a 1 diameter had a 34% more removal yield and a 9.6% higher removal rate than the column with a 4 diameter. This may be due to more air bubbles generated across the column per inch.sup.2 of the smaller diameter column, giving rise to a higher interface area per volume for the mass transfer. These results demonstrate the validity of our mass-transfer kinetics model. .sub.s, .sub.s, and .sub.sfor the column with a 1 diameter were 70.1%, 291.2 mg L.sup.1 h.sup.'1, and 0.018 r.sup.1, respectively.

EXAMPLE 4

Effects of Adjusting the Orifice Diameter of the Stripping Column

[0149] This example examined the effect of another operating parameter, namely the orifice diameter of the air diffuser, on .sub.s and .sub.s. Two orifice diameters were tested: diameters of 2 m and 50 m. Everything else was the same as EXAMPLE 3. The column diameter was 1. The same setup as EXAMPLE 1 was used for this experiment.

[0150] Looking now to FIG. 12, the results of the decays of [NH.sub.4.sup.+] as a function of time, using the air diffuser with a 50-m diameter orifice () and a 2-m diameter orifice () during aeration are shown. The equations corresponding to Equation 13 are included for each plot. .sub.s and .sub.s were determined as described above in EXAMPLE 2. The values of R.sup.2 for the linear regressions in determining as were 0.995 and 0.997 for the solid and broken lines, respectively. The equations, the broken line for the diffuser with the 50-m orifice diameter and the solid line for the diffuser with the 2-m orifice diameter, well reproduced both sets of the experimental data. The standard deviations of calculated [NH.sub.4.sup.+] by the equations shown by the solid and the broken lines from the experimental data are 19.9 mg L.sup.1 and 20.0 mg L.sup.1, respectively. The agreements suggest that both stripping processes follow Equation 12. The diffuser with a 2-m diameter orifice produced a 16.2% more removal yield and a 15.8%% higher removal rate than the diffuser with a 50-m diameter. This could be due to more surface area of the air bubbles generated by the diffuser with a 2-m diameter orifice, giving rise to a higher interface area per volume for the mass transfer. The results demonstrate the validity of our previously derived mass-transfer kinetics model. .sub.s, .sub.s, and .sub.s for the diffuser with the 2-m orifice diameter were 81.2%, 480.4 mg L.sup.1 h.sup.1, and 0.025 t.sup.1, respectively.

EXAMPLE 5

Comparison Between Usage and Omission of a Brnsted Acid/Base

[0151] Looking now to FIG. 13, the experimental setup of an ammonia gas stripping and recovery process according to the methods and systems disclosed herein is illustrated. This example examined the effects of the acid-base reactions, Equations 1 and 2, on .sub.s/r and .sub.s/r. Ca(OH).sub.2 and H.sub.2SO.sub.4 were used as the Brnsted base and acid for Equations 1 and 2, respectively. 12 g of Ca(OH).sub.2 was added to 2 L of ADL in the stripping column 24 to form the stripping mixture, while 12 g of 100% H.sub.2SO.sub.4 was mixed with 2 L of distilled water in the recovery column 32 to form the recovery solution. The mole of the former was 30% more than what is stoichiometric to the mole of NH.sub.4.sup.+ in the original ADL. On the other hand, the mole of the latter chemical was stoichiometric to the mole of NH.sub.4.sup.+ in the original ADL, assuming a 100% removal yield. The condition was the same as EXAMPLE 3, and sintered stone diffusers 96 with an orifice diameter of 2 m was used. pH in the stripping column 24 was 12 in the beginning of the experiment and gradually declined to 11 at the end, while pH in the recovery column 32 started at a value of 1.55 and ended with a pH of 5.5. In another experiment for comparison, the pH of ADL in the stripping column 24 was adjusted to be 12 by adding NaOH, but no further NaOH was added during the aeration, which is the same as EXAMPLE 3. For the recovery column 32, the pH was adjusted to be 5 by adding HCl. This condition in the recovery column is similar to U.S. Pat. No. 11,364,463. The stripped NH.sub.3 gas 82 was sent from the stripping column 24 to the recovery column 32 via piping 100 along with other gases including N.sub.2, O.sub.2, and CO.sub.2 and bubbled at the bottom of the recovery column 32 through the sintered stone diffusers 96. The pH and the NH.sub.4.sup.+ concentrations in both the stripping 24 and the recovery column 32 were monitored during the aeration via pH meters 20 and NH.sub.4.sup.+ sensors 18.

[0152] Turning now to FIG. 14, the experimental results of the decays of [NH.sub.4.sup.+] as a function of time in the stripping column, using the base according to Equation 1 () and not using any base, but with the ADL having the initial pH of 12 (), are shown. The equations corresponding to Equation 12 are included for each plot. .sub.s and .sub.s were determined as described above in EXAMPLE 2. The values of R.sup.2 for the linear regressions in determining as were 0.997 and 0.998 for the solid and broken lines, respectively. The equations, the broken line for the stripping without the base and the solid line for the stripping with the base, well reproduce both sets of experimental data. The standard deviations of calculated [NH.sub.4.sup.+] by the equations shown by the solid and the broken lines from the experimental data are 19.1 mg L.sup.1 and 21.0 mg L.sup.1, respectively. The agreements suggest that both stripping processes follow Equation 12. .sub.s and .sub.s were determined as described above in EXAMPLE 3. The difference between with and without the use of the base is clear. Applying Equation 1 not only accelerated the mass transfer of NH.sub.3 from liquid to gas, but also removed NH.sub.3 nearly completely. The use of acid-base reaction kinetics produced a 21% more removal yield and a 165% higher removal rate than the stripping without it. This may be due to the promotion of the removal of NH.sub.3 through the acid-base reaction which is fast and complete. .sub.s, .sub.s, and .sub.s with the use of Equation 1 were 98%, 1274.5 mg L.sup.1 h.sup.1, and 0.034 min.sup.1, respectively. Our experiments conclude that the effect of the acid-base reaction has the largest impact on .sub.s and .sub.s.

[0153] Focusing now on FIG. 15, the experimental results of the growth of [NH.sub.4.sup.+] as a function of time in the recovery column, using H.sub.2SO.sub.4 according to eq 2 () and not using any chemical, but with an initial pH of 5 () set by adding HCl to distilled water, are displayed. The equations corresponding to Equation 17 are included for each plot. .sub.r and .sub.r were determined as described above in EXAMPLE 3. The values of R.sup.2 for the linear regressions in determining .sub.r were 0.998 and 0.997 for the solid and broken lines, respectively. The equations, the broken line for the absorption without H.sub.2SO.sub.4 and the solid line for the absorption with H.sub.2SO.sub.4, well reproduce both sets of the experimental data. The standard deviations of calculated [NH.sub.4.sup.+] by the equations shown by the solid and the broken lines from the experimental data are 29.7 mg L.sup.1 and 6.8 mg L.sup.1, respectively. The agreements suggest that both absorption processes follow Equation 17. The difference between with and without the use of Equation 2 is clear. Applying Equation 2 not only accelerated the mass transfer of NH.sub.3 from gas to liquid, but also absorbed the stripped NH.sub.3 gas into the acidic recovery solution nearly completely. Using the acid-base reaction kinetics produced a 145% more recovery yield and a 199% higher recovery rate than the stripping without it. This could be due to the accelerated NH.sub.3 absorption through the acid-base reaction which is fast and complete. .sub.r and .sub.r with the use of eq 2 are 98% and 200.01 mg L.sup.1 h.sup.1, respectively, significantly higher than those without H.sub.2SO.sub.4, 40% and 66.89 mg L.sup.1 h.sup.1 , respectively. When Equation 2 was not used, once the original acid concentration was consumed, the mass transfer rate significantly dropped, with little transfer of the NH.sub.3 gas to the solution. .sub.r with the acid-base reaction was 0.0061 min.sup.1.

[0154] Table 1 summarizes the removal yields/rates and the recovery yield/rate for EXAMPLE 2 -5. The results demonstrate that the use of acid-base reactions in the NH.sub.3 stripping/absorption significantly increases the removal yield/rate and the recovery yield/rate which translates to a high NH.sub.3 recovery and a high productivity of the NH.sub.3 stripping, which can lead to a low production cost.

TABLE-US-00001 TABLE 1 The Removal Yield and Rate and Recovery Yield and Rate under Different Conditions.sup.a EXAMPLE 2 EXAMPLE 3 EXAMPLE 4 EXAMPLE 5 pH 9 12 12 12 orifice 50 50 2 50 2 diameter (m) column diameter () 4 1 4 1 1 acid-base reaction no no no no no yes No .sub.s (%) 50.2 70.1 52.2 81.2 70.1 98.0 81.2 .sub.s (mg/L/h) 254.2 291.2 265.8 480.4 291.2 1,274.5 480.4 .sub.r (%) 98 40 .sub.r (mg/L/h) 200.1 66.9 Temperature was ambient.

EXAMPLE 6

Continuous Stripping/Recovery of the Brnsted Acid/Base Method

[0155] In this example, the stripping process was continued under the same condition as that for EXAMPLE 5 by replacing ADL in the stripping column once the stripping was completed, while keeping the same solution in the recovery column by absorbing NH.sub.3 continuously into the acidic recovery solution. This cycle was repeated until the solubility of (NH.sub.4).sub.2SO.sub.4 in water reached its limit. The solubility of (NH.sub.4).sub.2SO.sub.4 in water is 744 g L.sup.1, corresponding to 202.9 g L.sup.1 of NH.sub.4.sup.+. The UV-vis spectroscopy measurement of the final solution indicated an NH.sub.4.sup.+ concentration of 19.7%.

[0156] Table 2 compares the results of this EXAMPLE 6 with others published in the literature. After an exhaustive literature search, no publication reporting .sub.s higher than 400 mg L.sup.1 h.sup.1 was found. On sharp contrast, the experiments according to the methods and systems discussed herein produced a three-fold higher .sub.s. This may have arisen due to a rigorous application of the acid-base reaction to facilitate the mass-transfer reaction kinetics of NH.sub.3 from liquid to gas.

TABLE-US-00002 TABLE 2 Comparison of Results by This Invention Against Others in Literature. [NH.sub.4.sup.+].sub.s.sup.0 .sub.s/r .sub.s/r pH T, C. mg/L % mg/L/h Comments Ref. tower stripping 8.5 70 1,100 63 115 air flowrate = 117 L h.sup.1 .sup.a with circulation 7.66 50 1,298 95 329 air flowrate = 2.7 L h.sup.1 .sup.b absorption 20 1,298 >100 .sup.b no aeration; N = 4% .sup.b tower stripping 10.4 42 2,834 98 46.3 air flowrate = 70,000 L h.sup.1 .sup.c with no circulation 9 75 2,244 91 .sup.d absorption N = 4 % .sup.d tower stripping 11.9 45.8 61.04 91 8 .sup.e with no circulation 12 25 2,500 99 178 air flowrate = 720 L min.sup.1 .sup.f 10.5 30 521~681 66.5~74.9 194~246 .sup.g aeration stripping 12 20 2,086 98 1,274 air flowrate = 1,800 L h.sup.1 .sup.h with acid-base reaction absorption with 98 200 N = 20% .sup.h acid-base reaction .sup.aKim, et al., 2021. .sup.bLaureni, M.; Palatsi, J.; Llovera, M.; Bonmat, A. Influence of Pig Slurry Characteristics on Ammonia Stripping Efficiencies and Quality of the Recovered Ammonium-Sulfate Solution, J. Chem. Technol. Biotechnol. 88, 1654-1662 (2013). .sup.cWang, Y.; Pelkonen, M.; Kotro, M. Treatment of High Ammonium-Nitrogen Wastewater from Composting Facilities by Air Stripping and Catalytic Oxidation, Water Air Soil Pollut. 208, 259-273 (2010). .sup.dMenkveld, H. W. H.; Broeders, E. Recovery of Ammonia from Digestate as Fertilizer, Water Practice & Technology 13, 382-387 (2018). .sup.eZangeneh, A.; Sabzalipour, S.; Takdatsan, A.; Yengejeh, R. J.; Khafaie, M. A. Ammonia Removal form Municipal Wastewater by Air Stripping Process: An Experimental Study, South African Journal of Chemical Engineering 36, 134-141 (2021). .sup.fOzyonar, F.; Karagozoglu, B.; Kobya, M. Air Stripping of Ammonia from Coke Wastewater, JESTECH, 15(2), 85-91, (2012). .sup.gZou, M.; Dong, H.; Zhu, Z.; Zhan, Y. Optimization of Ammonia Stripping of Piggery Biogas Slurry by Response Surface Methodology, Int. J. Environ. Res. Public Health 2019, 16, 3819; doi: 10.3390/ijerph16203819. .sup.hThis invention.

EXAMPLE 7

Sensitivity Analysis

[0157] Turning now to FIG. 16 the sensitivity analysis of the normalized NH.sub.4.sup.+ concentration of the stripping process,

[00030]$\frac{{\left[{\mathrm{NH}}_4^{+}\right]}_s^t}{{\left[{\mathrm{NH}}_4^{+}\right]}_s^i},$

to the stripping optimization variables

[00031]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

and K.sub.L,sa.sub.r in the mass transfer rate per unit volume, .sub.s, when

[00032]$K_{L,s}{a}_s\frac{1}{Q_{G,s}{H}_s}=0.5\left(\right),$

1(), and 2 () is displayed. The plots in the inset showing the range of

[00033]$\frac{{\left[{\mathrm{NH}}_4^{+}\right]}_s^t}{{\left[{\mathrm{NH}}_4^{+}\right]}_s^i}$

between 0.45 and 0.55 and the time between 15 min and 25 min depict the difference in the dependencies of

[00034]$\frac{{\left[{\mathrm{NH}}_4^{+}\right]}_s^t}{{\left[{\mathrm{NH}}_4^{+}\right]}_s^i}$

on the two combined variables. The time to reduce

[00035]$\frac{{\left[{\mathrm{NH}}_4^{+}\right]}_s^t}{{\left[{\mathrm{NH}}_4^{+}\right]}_s^i}$

by half is about 38% faster when K.sub.L,sa.sub.s is twice as much as

[00036]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

when the two optimization variables share the same value. On the other hand, the time to reduce

[00037]$\frac{{\left[{\mathrm{NH}}_4^{+}\right]}_s^t}{{\left[{\mathrm{NH}}_4^{+}\right]}_s^i}$

by half is about 22% faster when

[00038]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

is twice as much as K.sub.L,sa.sub.s when compared to when the two optimization variables are the same. Hence, it appears that it is more effective to adjust K.sub.L,sa.sub.s than to adjust

[00039]$\frac{Q_{G,s}{H}_s}{V_{L,s}}$

in to increase .sub.s. It is also important to observe that increasing either one optimization variable relative to the other raises the mass transfer rate compared to when both optimization variables are the same.

[0158] Looking now to FIG. 17, the sensitivity analysis of the normalized NH.sub.4.sup.+ concentration of the recovery process,

[00040]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i},$

to the recovery optimization variables

[00041]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i},$

and K.sub.L,ra.sub.r in the transfer rate per unit volume, .sub.r, when

[00042]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i},$

1(), and 2 () are displayed. The plots in the inset showing the range of

[00043]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i}$

between 0.40 and 0.60 and the time between 70 min and 120 min clearly depict the contrast in the dependencies of

[00044]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i}$

on the two optimization variables. The time to reduce

[00045]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i}$

by half is about 38% faster when K.sub.L,ra.sub.r is twice as much as

[00046]$\frac{Q_{G,r}{H}_r}{V_{L,r}}$

when the two optimization variables are the same. On the other hand, the time to reduce

[00047]$\frac{{\left[N{H}_4^{+}\right]}_r^t}{{\left[N{H}_4^{+}\right]}_s^i}$

by half is about 11% faster when

[00048]$\frac{Q_{G,r}{H}_r}{V_{L,r}}$

is twice as much as K.sub.L,sa.sub.s when the two optimization variables are the same. Hence, it appears that it is more effective to adjust K.sub.L,ra.sub.r than

[00049]$\frac{Q_{G,r}{H}_r}{V_{L,r}}$

to increase .sub.r. It is also important to observe that increasing either one optimization variable relative to the other raises the mass transfer rate compared to when both optimization variables are the same.

EXAMPLE 8

Prior Art Liquefaction of NH.SUB.3 .Gas

[0159] The liquefaction of the NH.sub.3 gas recovered from ADL was performed in this EXAMPLE 8 by two methods: one described in U.S. Pat. No. 11,364,463 and one according to the systems and methods disclosed herein. First, the method described in U.S. Pat. No. 11,364,463 will be examined.

[0160] Brining our attention now to FIG. 18, the experimental apparatus for the NH.sub.3 liquefaction based on the method described in U.S. Pat. No. 11,364,463 is illustrated. 2 L of ADL was placed in the stripping column 24, and pH was adjusted to 9. Aeration was conducted under the same condition as EXAMPLE 2 via air diffusers 28 and an air pump 22. The aerated gases were sent to a stainless-steel coiled pipe 58 connected to a liquid receiver 68 at the bottom inside the liquefaction device 52. Inside the liquefaction device 52, a coolant 102 continuously circulates to maintain the temperature inside below 70 C., but above 50 C. Since water freezes at these temperatures, despite U.S. Pat. No. 11,364,463's recommendation to use such a coolant, iso-propanol was used as the coolant 102 which was pumped by a pneumatic pump 40 from a cryogenic storage 106 to and out of the liquefaction device 52. Dry ice 104 was added to the cryogenic storage 106 to maintain the temperature. The internal temperatures of the liquefaction device 52 and the cryogenic storage 106 were monitored by a thermocouple thermometer 108. The gas release valve 60 was left open to release the other gases, N.sub.2 and O.sub.2. Liquid NH.sub.3 was collected in a liquid receptacle 68 before being discharged through a liquid discharge outlet 70.

[0161] The initial NH.sub.4.sup.+ concentration was 2,086 mg L.sup.1 and the volume was 2 L. With the 50% removal yield, found in EXAMPLE 2, it was expected to recover about 2 g of liquid NH.sub.3. However, as soon as the aeration started, a pungent smell emerged from the gas release valve.

[0162] Without a trace of liquid and the characteristic smell of the NH.sub.3 gas from the gas release valve, it was observed that the stripped NH.sub.3 gas came through the coiled pipe without being liquefied despite the thermometer showing 60 C. inside the liquefaction device 52. It was suspected that the other gases such as N.sub.2 and O.sub.2, and possibly CO.sub.2 interfered with the liquefaction of the stripped NH.sub.3 gas, preventing the contact of the NH.sub.3 gas with the coolant through the coiled pipe 58.

EXAMPLE 9

Liquefaction of NH.SUB.3 .Gas According to the Present Disclosure

[0163] Looking now to FIG. 19, a bench-scale setup for the second liquefaction experiment according to the methods and systems disclosed herein is illustrated. A boiling round flask 110 was connected to two serial condensers, a water vapor condenser 112 and a moisture removal condenser 114 incorporating a desiccant 116, attached to the same liquefaction device 52 as that used for EXAMPLE 8. 500 mL of (NH.sub.4).sub.2SO.sub.4 (19.1% (NH.sub.3)) solution produced by absorbing the stripped NH.sub.3 gas in EXAMPLE 6 was heated at 70 C. inside the round flask 110 via a heating mantle 118 generating bubbles on the solution surface indicating the gas evaporation. The boiling point of NH.sub.3 is33.3 C. Hence, the liquid NH.sub.3 was easily vaporized, but the evaporated NH.sub.3 gas was also pushed towards the liquefaction device by the thermal energy. The temperature was measured with a thermometer 108. The evaporated mixture 120 was sent to the liquefaction device 52. The liquefaction device 52 was operated in the same way as was the case for EXAMPLE 8, incorporating similar features such as the coiled pipe 58, thermometers 108, liquid receptacle 68, liquid discharge outlet 70, gas release valve 60, coolant 102, pneumatic pumps 40, cryogenic storage 106, and dry ice 104.

[0164] After the bubbling of the solution inside the boiling flask 110 ceased, the circulation of the coolant 102 inside the liquefaction device 52 continued for another 30 min. Then, the coolant 102 circulation was stopped, and the liquid outlet valve 70 was released. The liquid NH.sub.3 was discharged, and its volume was measured by a graduate cylinder. The collected volume was 111 mL which was equivalent to 87.6 g. A 500 mL of 19.1% NH.sub.3 solution should yield 95.7 g of liquid NH.sub.3 in principle with a density of 0.73 g mL.sup.1 for liquid NH.sub.3. Hence, the NH.sub.3 recovery efficiency was 91.5%. The difference between this experimental setup and EXAMPLE 8 is that the NH.sub.4.sup.+ concentration in the solution, from which the NH.sub.3 gas is to be recovered, increased by two orders of magnitude in EXAMPLE 9 when compared to EXAMPLE 8. Very little other gases such as N.sub.2 and O.sub.2 were present in the evaporated gas. As a result, a high flux of NH.sub.3 gas enters the liquefaction device with very few interferences by other gases and therefore, it is easier for condensation to occur in EXAMPLE 9.

[0165] While the present invention has been described with regards to particular embodiments, it is recognized that additional variations of the present invention may be devised without departing from the inventive concept.