the invention comprises an electrochemical cell with a cathode and anode on either side of a microfluidic channel . the channel is preferably of a height between 100 and 500 microns , and optionally between a height of 50 microns and 1 millimeter . several smaller channels lead into this channel , so that fuel , oxidant , and electrolyte flow can be controlled . on the anode side , fuel is provided either in liquid form via one of these smaller channels or in gaseous form by making the anode out of a porous material that allows gaseous fuel to diffuse through and reach the surface of the anode facing the electrolyte . on the cathode , a liquid oxidant such as bromine is provided via one of the smaller channels leading into the larger channel . the non - turbulent nature of the fluid flow within the device ensures that the mixing zones between fuel , electrolyte , and oxidant can be determined analytically , and kept small enough to eliminate any efficiency loss due to crossover of reactants . with reference now to fig1 a and 1 b , in fig1 a an electrochemical cell 10 includes a porous anode 12 and a cathode 14 . a fuel 16 such as hydrogen is flowed along the anode 12 and will diffuse to the surface of the anode facing the electrolyte . a smaller channel 18 introduces an electrolyte such as hydrobromic acid , sulfuric acid , or potassium hydroxide , and a smaller channel 20 serves to introduce a concentrated liquid oxidant such as bromine , either in concentrated faun , in water , or mixed with the electrolyte . the electrolyte is diluted in water to a concentration between 0 . 1 and 10 molar . with respect to fig1 b , the anode 12 is non - porous . in this embodiment a fuel passage 22 serves to introduce fuel into the cell 10 . it should be noted that in fig1 b , a liquid fuel such as formic acid , methanol , ethanol or hydrazine is suitable and may be supplied to flow along the anode 12 . depletion zones 24 will be determined by the operating current of the cell as will be discussed below and can be analyzed to determine the limiting current of the cell . the non - turbulent nature of the fluid flow within the cell 10 ensures that the mixing zones between fuel , electrolyte , and oxidant can be determined analytically and kept small enough to eliminate any efficiency loss due to crossover of reactants . because non - turbulent flows can be described by relatively simple forms of the navier - stokes equation , and all the products and reactants ( apart from hydrogen , which is kept distinct from the channel ) are liquids , steady state modeling and optimization of this cell is possible without resorting to two phase flow calculations . one additional feature of the device shown in fig1 a is that as the system approaches peak power , the edge of the depletion region 24 will approach the oxidant electrolyte interface , and the bromine concentration in the effluent will asymptotically approach zero , resulting in high fuel utilization . in such a configuration , no recirculation of oxidant would be necessary , and recovery of bromine from the hydrobromic acid can be done electrolytically by reversing the flow of direction to the system and supplying current . likewise , because no product gases or liquids are produced at the anode , there is no need to recirculate or humidify the hydrogen gas being supplied to the system . a more detailed description of the electrochemistry of the system is provided below . the fluid mechanics of the system are such that the system can be easily scaled in one dimension without affecting the underlying physics . in order to fabricate large scale systems , however , two dimensional scaling will be necessary . in order to achieve this , we implement sequential injections of oxidant into the system using an array of electrodes as shown in fig2 a , b and c . such a system could be constructed in a disc - like geometry as shown in fig3 a . as shown in fig3 b , oxidant and electrolyte may be injected at the center of the cell and flow outwardly . a lower , oxidant - only channel is used to provide periodic injections of oxidant to the cathode in order to replenish the depletion zone . the disc - like geometry in fig3 a is able to accommodate the increased flow . in order to evaluate the capabilities of a system based on this technology , a simple model for the performance of a hydrogen - bromine laminar flow fuel cell , as illustrated in fig1 a is presented here . for a first order estimation , we will examine a one dimensional system with planar electrodes . we will assume that the equilibrium voltage is described by the nernst equation , and that the reaction kinetics are accurately described by the butler - volmer equation with symmetric ( α = ½ ) reactions . concentration polarizations will be introduced based on our results shown below . v ( i )= φ c − φ a − ir el = v eq − η c − η a − ir el it is worth noting that because all the reactants and products ( excluding the hydrogen gas , which has a well defined interface with the rest of the system ) are liquid , there will be no two - phase flow , and it is reasonable to assume locally homogeneous concentrations . the equilibrium voltage v eq can be determined from the nernst equation using existing thermodynamic data at reference conditions of 1 atm , 298 . 15 k , and unit molality . [ 18 ] the activity coefficient for hbr over a broad range of molalities is ˜ 0 . 8 , so an ideal dilute solution theory for the liquid reactants is valid . likewise , the hydrogen can be treated as an ideal gas . in this case , the nernst equation reads : under a typical operating condition of m br2 = 3 mol / kg and m hbr , = 2 mol / kg with everything else at standard conditions , veq = 1 . 083 v . to determine the electrolyte resistance , we choose a characteristic channel width of 500 μm , and employ existing data to determine the specific resistance of the cell to be rel = 89 mω · cm 2 for 2 molar hydrobromic acid . [ 19 ] to complete the model , we need to describe the reaction kinetics at both the anode and cathode . in general , the butler - volmer equation can be written as i = i i 0 ( e ( 1 − α ) fη i / rt − e − αfη i / rt ) where η i is the activation overpotential and i i 0 is the exchange current density at either the anode or cathode . if we make the symmetric reaction assumption , butler - volmer can be simplified and written as i i 0 = k i 0 √ { square root over ( a reactants a products )} although concentrated solution theory would be necessary to obtain optimal accuracy , we are justified in approximating the system as ideally dilute because the activity coefficients for all components in the system are close to unity . previous studies have determined the anode and cathode reactions in the presence of platinum to have exchange current densities of order by k a 0 ˜ 1 ma / cm2 and k c 0 ˜ 50 ma / cm2 respectively . [ 8 , 18 ] if we further observe that the anode in this system will be porous with a high surface area catalyst , we are justified in estimating a 100 to 1 ratio for the anode surface area to electrolyte surface area . [ 20 ] these elements can be combined to determine a voltage current relationship . as a final correction , we now include the effects of concentration polarization on the equilibrium voltage , again assuming ideal solution behavior of the reaction kinetics . we know that the bromine concentration mbr2 at the electrode surface can be written in terms of the bulk concentration as so we can easily include concentration polarizations in the nernst equation . there will also be an enrichment effect on the hydrobromic acid , which can be written for simplicity , we consider the case where pe ˜ l / h . because of the weak ( ⅓ power ) dependence on limiting current density that was derived below , variations from this condition will not strongly affect limiting current . the limiting current density can then be written where w = 50 μm is the channel width , d ˜ 10 5 cm 2 is the diffusivity of bromine in an aqueous solution , and ρ = 3 . 10 × 10 − 3 kg / cm 3 is the density of bromine . for 3 molal bromine and 2 molal hydrobromic acid , we obtain a limiting current under our typical operating conditions of ilim ˜ 7 . 179 a / cm 2 . when all of these effects are taken into account , we obtain power and current characteristics ( neglecting external contact losses ) with a peak power density of over 5 w / cm 2 occurring at 6 . 895 a / cm 2 . at 1 . 5 a / cm 2 , a typical high current operating point , the cell produces 1 . 383 w / cm 2 at an efficiency of 84 . 9 %. these efficiencies rival those of batteries while still maintaining the energy density advantages that fuel cells typically exhibit . it worth mentioning that these results do not consider overall depletion of the oxidant in the bulk of the channel , which is significant at high fuel utilization . in order to further verify the predicted results , a two - dimensional model of the system was implemented in comsol multiphysics , a finite element software package . unlike the 1d model , a finite flow rate is assumed in this case , and the polarization losses due to depletion of the oxidant are considered . the model geometry is based on the system shown in fig8 , and the governing equations are shown in fig4 . the results obtained are qualitatively consistent with the one - dimensional boundary layer model shown in fig7 . as expected , the diffusion and depletion zones both grow gradually over the length of the electrode , so that the bromine concentration at the outlet of the cell is significantly lower than at the inlet of the cell . a representative solution is shown in fig5 . by running the model under a number of operating conditions , voltage - current relationships can be obtained . fig6 shows predicted performance over a range of concentrations at a fixed oxidant flow rate . under these conditions , a maximum power density of over 2 w / cm2 is predicted . understanding the characteristics of the depletion region above the electrode in a laminar flow cell as shown in fig7 is critical to estimating the performance characteristics of this type of system . for low reynolds number flow in a channel with height and depth h & lt ;& lt ; d and flow rate q , the unidirectional flow profile can be described as we wish to characterize the oxidant concentration profile above an electrode operating at constant current i for the case where the inlet concentration is c 0 and the oxidant diffusion constant is d . the governing advection diffusion equation for this system is if we treat the concentration boundary layer as thin , the boundary conditions on this system can be written : for thin boundary layers , the quadratic velocity profile in y can be approximated as linear , and the governing equation can be expressed in dimensionless quantities as for distances sufficiently far from the reaction zone , the concentration will not be affected by the reaction , which means that all of the change in concentration must occur within a layer of some thickness δ . this layer can be determined by considering the dominant balance for the governing equation . putting the dimensions back in , the boundary layer is of thickness : this result shows us that as the depletion region develops along the length of the electrode , the solute will have to diffuse further to reach the electrode . this will result in a reduction in the attainable current , and will ultimately allow us to predict a limiting current . we can calculate this limiting current by observing that the maximum oxidant flux occurs when the oxidant concentration at the electrode reaches zero . for a purely diffusion based system , the limiting current can be written as note that for the diffusion only system , the dimensionless sherwood number sh is equal to one . in the proposed system , then the limiting current can be determined by we can see from our boundary condition that so our limiting current is this result demonstrates that in order to avoid being limited by the depletion zone , flow velocity must be sufficiently high such that pe ˜ l / h . in order for the proposed system to exhibit the necessary fluid flow and transport properties for optimal performance , the fluid flow must be of a laminar nature . for this to be the case , the reynolds number of the system must be well below 2100 . for a fluid channel of height and depth h & lt ;& lt ; d , and a fluid of viscosity , μ and density ρ with a volumetric flow rate of q , the reynolds number for the system can be expressed as in a preferred configuration , d = 381 μm , ρ = 1000 kg / m 3 , ρ = 0 . 001 pa s , and q = 1 ml / min . in this system , the reynolds number is approximately 8 , which is well within the laminar flow regime . although the viscosity and density of the fluid will tend not to vary widely , the flow rate and channel depth can be varied freely within this constraint . for example , the flow rate could be increased by a factor often to facilitate transport of oxidant to the cathode without introducing any turbulence into the velocity profile in the channel . an additional transport requirement for this system is that the diffusive mixing zone at the interface between the streams must not reach either the cathode or anode . as discussed above , the diffusive mixing zone scan be estimated to have thickness in the same preferred configuration , d ≈ 10 − 5 cm 2 / s , so pe ≈ 800 . in that case , in order for the mixing zone not to exceed half the channel width , the electrode length must be less than 3 . 8 cm . in a preferred embodiment , a 381 μm layer of 75 durometer viton sheet is sandwiched between a 3 mm hastelloy - c cathode current collector and a 10 mm polyvinylidene fluoride ( pvdf ) porting plate . a channel 2 mm wide is cut into the viton layer to form the fluid channel . the cathode current collector is prepared by applying a 1 mg / cm 2 loading of platinum black to the active area underneath the channel to serve as the cathode catalyst . the pvdf porting plate is prepared by milling a 2 mm wide , 24 mm long , 315 μm deep pocket into the underside of the plate , perpendicular to the fluid channel . this pocket is filled with a piece of porous gas diffusion media of the same size that serves as the anode current collector . the face of the diffusion media facing the fluid channel is coated with a 2 μm microporous layer consisting of high surface area carbon and polytetrafluoroethylene ( ptfe ), followed by a layer of 1 mg / cm 2 platinum black . three ports for oxidant , electrolyte , and effluent are drilled and tapped into the pvdf porting plate along the length of the fluid channel , and two gas ports are drilled and tapped into the plate along the length of the gas diffusion media channel . in order to facilitate sealing but maintain chemical compatibility , a compression fitting employing an ethylene - tetrafluoroethylene ( etfe ) ferrule is used . these components are shown in fig8 . copper wires connect the cathode current collector to an external electronic load initially programmed to maintain zero current . separately , two aqueous solutions are prepared in glass syringes : one is a mixture of one molar hydrobromic acid and one molar bromine , which will serve as the oxidant stream , and one is one molar hydrobromic acid , which will serve as the electrolyte . the syringes are loaded into syringe pumps . the electrolyte solution is fed via a 1 / 16 ″ outer diameter ptfe tube into the electrolyte port on the pvdf plate at a rate of 0 . 8 ml per minute , and it is collected out of the effluent port at the far end of the channel . hydrogen gas at a pressure of 10 kpa is fed into one of the gas ports and allowed to flow through the second port for two minutes to purge out any trapped air , and then the second port is shut . the strongly hydrophobic nature of the microporous layer on the anode current collector ensures that fluids do not enter the porous diffusion media . once the exhaust port is shut , the oxidant solution is fed into the oxidant port at a flow rate of 0 . 2 ml per minute via a second 1 / 16 ″ outer diameter ptfe tube . the addition of the oxidant completes the electrochemical circuit , and a voltage of approximately one volt is observed across the electronic load . the external electronic load is then programmed to allow current to pass , and the system produces electricity . a schematic of this configuration is shown in fig9 . in a second preferred embodiment , the cell is constructed in a similar fashion , however the pvdf porting plate is replaced with a hastelloy - c porting plate that now serves as the anode current collector , and the porous gas diffusion media and its channel and gas ports are eliminated . a layer of 1 mg / cm 2 platinum black is applied to the porting plate in the active region above the microfluidic channel to serve as the anode catalyst . also , an additional port for fuel is added to the porting plate after the electrolyte port but upstream of the active region along the length of the fluid channel . the copper wire that previously connected the porous anode current collector to the external load now connects the hastelloy porting plate to the external load . in place of gaseous hydrogen fuel , an aqueous solution of one molar formic acid is prepared and loaded into a third glass syringe in a syringe pump . the electrolyte is pumped at a rate of 0 . 6 ml per minute , while the fuel and oxidant streams are both pumped at a rate of 0 . 2 ml per minute . the external load is operated in the same manner as in example 1 . it is recognized that modifications and variations of the invention disclosed herein will be apparent to those of ordinary skill in the art and it is intended that all such modifications and variations be included in the claims . e . r . choban , p . waszczuk , p . j . a . kenis , t . tzedakis , s . k . yoon , and c . kane , “ microfluid device and synthetic methods ,” u . s . pat . no . 7 , 273 , 541 , 2007 . j . l . cohen , d . j . volpe , and h . d . abruna , “ duel electrolyte membraneless microchannel fuel cells ,” u . s . patent application 2006 / 0 228 622 , 2006 . l . cohen , d . j . volpe , d . a . westly , a . pechenik , and h . d . abruna , “ planar membraneless microchannel fuel cell ,” u . s . pat . no . 7 , 435 , 503 , 2008 . s . hasegawa and h . watanabe , “ cell and power generation method ,” u . s . patent application 2007 / 0 181 418 , 2007 . p . j . a . kenis , a . wieckowski , and e . r . choban , “ membraneless electrochemical cell and microfluidic device without ph constraint ,” u . s . pat . no . 7 , 635 , 530 , 2009 . e . kjeang , d . sinton , and n . djilali , “ a microfluidic fuel cell with flow — through porous electrodes ,” u . s . patent application 2009 / 0 092 882 , 2009 . l . j . markoski , j . moore , and j . lyding , “ electrochemical cells comprising laminar flow induced dynamic conducting interfaces , electronic devices comprising such cells , and methods employing same ,” u . s . pat . no . 6 , 713 , 206 , 2004 . l . j . markoski , d . natarajan , and a . primak , “ fuel cell apparatus and method of fabrication ,” u . s . patent application 2005 / 0 202 305 , 2005 . l . j . markoski , d . natarajan , and a . primak , “ electrochemical cells ,” u . s . patent application 2006 / 0 088 744 , 2006 . l . j . markoski , d . natarajan , and a . primak , “ system for flexible in situ control of water in fuel cells ,” u . s . patent application 2007 / 0 190 393 , 2007 . l . j . markoski , p . j . a . kenis , and e . r . choban , “ fuel cells comprising laminar flow induced dynamic conducting interace , electronic devices comprising such cells , and methods employing same ,” u . s . pat . no . 7 , 252 , 898 , 2007 . l . j . markoski , d . natarajan , and a . primak , “ microfluidic fuel cells ,” u . s . patent application 2008 / 0 248 343 , 2008 . l . j . markoski , a . primak , and d . natarajan , “ microfluidic fuel cell electrode system ,” u . s . patent application 2009 / 0 035 644 , 2009 . l . j . ohlsen and j . c . mallari , “ fuel cell systems having internal multistream laminar flow ,” u . s . patent application 2004 / 0 058 217 , 2004 . w . e . pelton and p . j . a . kenis , “ system of distributed electrochemical cells integrated with microelectronic structures ,” u . s . patent application 2007 / 0 020 496 , 2007 . w . e . pelton and p . j . a . kenis , “ distributed electrochemical cells integrated with microelectronic structures ,” u . s . patent application 2007 / 0 026 266 , 2007 . k . s . salloum and j . posner , “ counter - flow membraneless fuel cell ,” u . s . patent application 2010 / 0 112 391 , 2010 . s . bratsch , “ standard electrode potentials and temperature coefficients in water at 298 . 15 k ,” j . phys chem ref data , vol . 18 , no . 1 , pp . 1 - 21 . 1989 . d . r . lide , crc handbook of chemistry and physics , 75 th ed . cleveland , ohio : crc press , 2004 . r . savinell and s . fritts , “ theoretical performance of a hydrogen bromine rechargeable spe fuel - cell ,” j power sources , vol . 22 , no . 3 - 4 , pp . 423 - 440 , 1988 .