Extracorporeal bioartificial liver devices (BAL) are perhaps among the most promising technologies for the treatment of liver failure, but significant technical challenges remain in order to develop systems with sufficient processing capacity and of manageable size. One key limitation is that during BAL operation, when the device is exposed to plasma from the patient, hepatocytes are prone to accumulate intracellular lipids and exhibit poor liver-specific functions. Based on hepatic intermediary metabolism, we have utilized mathematical programming techniques to optimize the biochemical environment of hepatocyte cultures towards the desired effect of increased albumin and urea synthesis. To investigate the feasible range of optimal hepatic function, we have obtained a Pareto optimal set of solutions corresponding to liver-specific functions of urea and albumin secretion in the metabolic framework using multiobjective optimization. The importance of amino acids in the supplementation and the criticality of the metabolic pathways have been investigated using logic-based programming techniques. Since the metabolite measurements are bound to be patient specific, and hence subject to variability, uncertainty has to be integrated with system analysis to improve the prediction of hepatic function. We have used the concept of two stage stochastic programming to obtain robust solutions by considering extracellular variability. The proposed analysis represents a new systematic approach to analyze behavior of hepatocyte cultures and optimize different operating parameters for an extracorporeal device based on real-time conditions. (c) 2005 Wiley Periodicals, Inc.