A nonlinear state-variable method is presented and used to solve the pseudo-2D (P2D) Li-ion cell model under high-frequency input current and temperature signals. The physics-based governing equations are formulated into a nonlinear state variable method (NSVM), in which the mass transfer variables are evaluated using a 1 st order exponential integrator approach at each discrete time point and the electrochemical kinetics (Butler-Volmer) equations are solved by either an iterative or an explicit method. This procedure provides an accurate, computationally efficient method to develop physics-based simulations of the performance of a dual-foil Li-ion cell during practical drive cycles. Physics-based Li-ion cell models are useful tools for analyzing the battery cell performances, characterizing material properties, supporting cell design, and developing battery management systems (BMS). 5 in their work, a dual-foil Li-ion cell was represented by a pseudo-two-dimensional domain; and therefore, this model is usually referred as the "pseudo-2D model (P2D)".3 This model contains several coupled transport phenomena described by nonlinear partial differential equations (PDEs). Due to the numerical complexity of these models, various reduced-order models (ROMs) have been proposed as substitutes for either the single particle (SP) or the P2D full-order model (FOM) in practical online simulation applications.6-16 Also, ROMs have been published to simulate capacity fading 17,18 and for use in parameter estimations. presented an approximate solution to the spherical diffusion problem with a fixed flux which provides high accuracy with only a few terms in a series when compared to the analytic solution with at least 100 terms. This approximate solution provides significant savings in computation time.In order to predict the performance of lithium ion cells operating at higher rates, more complicated models were developed. In 1982, West et al. 4 presented a Nernst-Planck model for porous insertion electrodes for a Li/TiS 2 cell. They pointed out the importance of coupling the transport of lithium ions between the insertion electrode and the electrolyte. In 1993, Mao and White 35 presented a Nernst-Planck model of the Li/TiS 2 cell which included transport of lithium ions through the separator and demonstrated quantitatively how the utilization of the TiS 2 electrode could be improved by decreasing the thickness of the separator. Also in 1993, Doyle et al. 36 published a concentrated solution theory model for a lithium/polymer/TiS 2 cell. Shortly after that in 1994 Fuller et al.5 published a concentrated solution theory model for a dual lithium ion insertion cell ("dual foil"), which is known as "pseudo-2D model (P2D)". In 1996, Doyle et al. 37 published a comparison of their model predictions to experimental data from a Bellcore plastic lithium ion cell. Unfortunately, the computation burden associated with using the program dual foil is too great to use their program for extensive simulations needed for parameter estimation an...