Freezing is an important operation in biotherapeutics industry. However, water crystallization in solution, containing electrolytes, sugars and proteins, is difficult to control and usually leads to substantial spatial solute heterogeneity. Herein, we address the influence of the geometry of freezing direction (axial or radial) on the heterogeneity of the frozen matrix, in terms of local concentration of solutes and thermal history. Solutions of hemoglobin were frozen radially and axially using small-scale and pilot-scale freezing systems. Concentration of hemoglobin, sucrose and pH values were measured by ice-core sampling and temperature profiles were measured at several locations. The results showed that natural convection is the major source for the cryoconcentration heterogeneity of solutes over the geometry of the container. A significant improvement in this spatial heterogeneity was observed when the freezing geometry was nonconvective, i.e., the freezing front progression was unidirectional from bottom to top. Using this geometry, less than 10% variation in solutes concentration was obtained throughout the frozen solutions. This result was reproducible, even when the volume was increased by two orders of magnitude (from 30 mL to 3 L). The temperature profiles obtained for the nonconvective freezing geometry were predicted using a relatively simple computational fluid dynamics model. The reproducible solutes distribution, predictable temperature profiles, and scalability demonstrate that the bottom to top freezing geometry enables an extended control over the freezing process. This geometry has therefore shown the potential to contribute to a better understanding and control of the risks inherent to frozen storage. (c) 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1212-1221, 2013