A Comparison of First- and Second Order-in-Time Finite Difference Methods Applied to Nonlinear Reactive Transport

Authors

  • Anastasia Wilson Augusta University
  • Carson Morris Augusta University
  • Kayli Hendricks Augusta University
  • Karen Lawrence Augusta University

DOI:

https://doi.org/10.33043/28xy39zc

Abstract

In this paper, we consider solution methods for the nonlinear reactive transport equation used to model the protein adsorption process. Efficient methods for simulating this process are necessary to aid in the development of novel adsorptive chromatography media to ensure high-volume production of purified product for the purposes of biotherapeutics. Using MATLAB®, we compare four finite difference schemes used to solve the nonlinear reactive transport equation, focusing on the differences of efficacy between implicit and explicit methods. As such, two of the methods are semi-implicit and two are explicit with one of each kind using a first-order temporal scheme and one of each using a second-order temporal scheme. The semi-implicit methods evaluate almost all terms implicitly while lagging the nonlinear coefficient function in time to linearize the equations. We include numerical results that indicate optimal convergence of the schemes, and we compare the effectiveness of the schemes in matching experimental data using two different boundary conditions.

Author Biographies

Anastasia Wilson, Augusta University

Anastasia Wilson is an Associate Professor of Mathematics at Augusta University with a Ph.D. in Computational Mathematics from Clemson University and an M.S. in Applied Mathematics from NC State. Formerly at NASA Langley Research Center, her research focuses on mathematical modeling, computational methods, numerical analysis, and optimization, particularly in interdisciplinary.

Carson Morris, Augusta University

Carson Morris graduated with a bachelors degree in mathematics at Augusta University in 2022, and he is scheduled to receive his masters degree in mathematics in May of 2024. He has a wide variety of research interests, including but not limited to probability, statistics, numerical analysis, and differential equations.

Kayli Hendricks, Augusta University

Kayli Hendricks contributed to this paper as an undergraduate at Augusta University where, in 2023, she graduated with a bachelors degree in mathematics and a Minor in Accounting. Since graduating, she has been promoted to Office Manager at a local tax and accounting firm in Augusta, Georgia.

Karen Lawrence, Augusta University

Karen Lawrence graduated from Augusta University in 2022 with a Bachelors Degree of Science in Mathematics, having worked on this paper from her Junior to Senior year. Karen currently works full-time in the legal field and plans to attend law school in 2025.

References

Davies, N. The future of biologics. Thepharmaletter. (2017)

Research, T. Global Biological Drugs Market to be Worth US 287,139.7 Million by 2020. (http://globenewswire.com/newsrelease/2014/10/20/674317/10103285/en/Global-Biological-Drugs-Market-to-be-Worth-US-287-139-7-Million-by-2020-Transparency-Market-Research.html,2014)

Hiller, A. Fast Growth Foreseen for Protein Therapeutics. Genet. Eng. Biotechn.. 29 pp. 153-155 (2009), http://www.genengnews.com/genarticles/fast-growth-foreseen-for-protein-therapeutics/2722/

Langer, E. Focus on Efficiency: Single-use, analytical methods and downstream processing at the forefront. Pharm. Manuf.. March pp. 3-11 (2013)

Bhut, B., Wickramasinghe, S. & Husson, S. Preparation of high-capacity, weak anion-exchange membranes for protein separations using surfaceinitiated atom transfer radical polymerization. J. Membr. Sci.. 325 pp. 176-183 (2008)

Bhut, B. & Husson, S. Dramatic performance improvement of weak anionexchange membranes for chromatographic bioseparations. J. Membr. Sci.. 337 pp. 215-233 (2009)

Bhut, B., Christensen, K. & Husson, S. Membrane chromatography: Protein purification from E.colilysate using newly designed and commercial anionexchange stationary phases. J. Chromatogr. A. 1217, 4946-4957 (2010)

Bhut, B., Christensen, K. & Husson, S. Membrane chromatography: Protein purification from E. Colilysate using newly designed and commercial anionexchange stationary phases. J. Chromatogr. A. 1217, 4946-3957 (2010)

Chenette, H., Robinson, J., Hobley, E. & Husson, S. Development of highproductivity, strong cation-exchange adsorbers for protein capture by graft polymerization from membranes with different pore sizes. J. Membrane Sci.. 423424 pp. 43-52 (2012)

Wang, J., Sproul, R., Anderson, L. & Husson, S. Development of multimodal membrane adsorbers for antibody purification using atom transfer radical polymerization. Polymer. 55, 1404-1411 (2014)

Wang, J., Wilson, A., Robinson, J., Jenkins, E. & Husson, S. A new multimodal membrane adsorber for monoclonal antibody purifications. J. Membr. Sci.. 492 pp. 137-146 (2015)

Marsily, G. Quantitative Hydrogeology: Groundwater Hydrology for Engineers. (Academic Press,1986)

Nfor, B., Noverraz, M., Chilamkurthi, S., Verhaert, P., Wielen, L. & Ottens, M. High-throughput Isotherm Determination and Thermodynamic Modeling of Protein Adsorption on Mixed Mode Adsorbents. J. Chromatogr. A. 1217 pp. 6829-6850 (2010)

Singh, N., Husson, S., Zdyrko, B. & Luzinov, I. Surface modification of microporous PVDF membranes by ATRP. J. Membrane Sci.. 262 pp. 81-90 (2005)

Singh, N., Wang, J., Ulbrict, M., Wickramasinghe, S. & Husson, S. Surfaceinitiated atom transfer radical polymerization: A new method for the preparation of polymeric membrane adsorbers. J. Membrane Sci.. 309 pp. 64-72 (2008)

Boyer, T., Miller, C. & Singer, P. Modeling the Removal of Dissolved Organic Carbon by Ion Exchange in a Completely Mixed Flow Reactor. Water Res.. 42 pp. 1897-1906 (2008)

Farthing, M., Kees, C., Russell, T. & Miller, C. An ELLAM Approximation for Advective-Dispersive Transport with Nonlinear Sorption. Adv. Water Resour.. 29 pp. 657-675 (2006)

Kaur, J., Malengier, B. & Remesíková, M. Convergence of an operator splitting method on a bounded domain for a convection-diffusion-reaction system. J. Math. Anal. Appl.. 348 pp. 894-914 (2008)

Poulain, C. & Finlayson, B. A Comparison of Numerical Methods Applied to Non-linear Adsorption Columns. Int. J. Numer. Meth. Fl.. 17 pp. 839-859 (1993)

Remesíková, M. Solution of Convection-Diffusion Problems with Nonequilibrium Adsorption. J. Comput. Appl. Math.. 169 pp. 101-116 (2004)

Suen, S. & Etzel, M. A mathematical analysis of affinity membrane bioseparations. Chem. Eng. Sci.. 47, 1355-1364 (1992)

Arbogast, T., Wheeler, M. & Zhang, N. A nonlinear mixed finite element method for a degenerate parabolic equation arising in flow in porous media. SIAM J. Numer. Anal.. 33 pp. 1669-1687 (1996)

Wilson, A. & Jenkins, E. Numerical Simulation of Solid Phase Adsorption Models Using Time-Integrated, Up-winded Finite Element Strategies. Comput. Sci. Eng.. (2019)

Wilson, A. & Jenkins, E. Analysis of a fully implicit SUPG scheme for a filtration and separation model. Comp. Appl. Math.. (2020)

Wilson, A. & Jenkins, E. Towards Higher Order Methods for Nonlinear Adsorption Problems.

Tarafder, A. Modeling and Multi-Objective Optimization of a Chromatographic System. Multi-Objective Optimization In Chemical Engineering: Developments And Applications. (2013)

Yang, H., Bitzer, M. & Etzel, M. Analysis of Protein Purification Using Ion-Exchange Membranes. Ind. Eng. Chem. Res.. 38 pp. 4044-4050 (1999)

Scopes, R. Protein Purification: Principles and Practice. (Springer-Verlag,1994)

Mott, H. & Green, Z. On Danckwerts’ Boundary Conditions for the Plug-Flow with Disperson/Reaction Model. Chem. Eng. Comm.. 202, 739-745 (2015)

Mollerup, J. A review of the thermodynamics of protein association to ligands, protein adsorption and adsorption isotherm. Chem. Eng. Technol.. 31, 864-874 (2008)

Strikwerda, J. Finite Difference Schemes and Partial Differential Equations. (SIAM,2004)

Cheney, E. & Kincaid, D. Numerical Mathematics and Computing. (Brooks Cole,2007)

Suen, S. & Etzel, M. A mathematical analysis of affinity membrane bioseparations. Chem. Eng. Sci.. 47, 1355-1364 (1992)

Agarwal, N., Semmens, M., Novak, P. & Hozalski, R. Zone of influence of a gas permeable membrane system for delivery of gases to groundwater. Water Resour. Res.. 41 (2005) Comparing Finite Difference Methods to Nonlinear Transport 25

Xu, M. & Eckstein, Y. Statistical analysis of the relationships between dispersivity and other physical properties of porous media. Hydrogeol. J.. 5, 4-20 (1997)

Dimartino, S., Boi, C. & Sarti, G. A validated model for the simulation of protein purification through affinity membrane chromatography. J Chrom A. 1218 pp. 1677-1690 (2022)

Riske, F. & Ransohoff, T. Development of Continuous Capture Steps in Bioprocess Applications. Preparative Chromatography For Separation Of Proteins. (2017)

Published

2025-03-28

How to Cite

Wilson, A., Morris, C., Hendricks, K., & Lawrence, K. (2025). A Comparison of First- and Second Order-in-Time Finite Difference Methods Applied to Nonlinear Reactive Transport. Mathematics Exchange, 18(1), 2–25. https://doi.org/10.33043/28xy39zc

Funding data