fegomezl/Brinicle: Modeling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modeled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behavior of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.