RESUMEN
An unstable density stratification between two fluids mixes spontaneously under the effect of gravity, a phenomenon known as Rayleigh-Taylor (RT) turbulence. If the two fluids are immiscible, for example, oil and water, surface tension prevents intermixing at the molecular level. However, turbulence fragments one fluid into the other, generating an emulsion in which the typical droplet size decreases over time as a result of the competition between the rising kinetic energy and the surface energy density. Even though the first phenomenological theory describing this emulsification process was derived many years ago, it has remained elusive to experimental verification, hampering our ability to predict the fate of oil in applications such as deep-water spills. Here, we provide the first experimental and numerical verification of the immiscible RT turbulence theory, unveiling a unique turbulent state that originates at the oil-water interface due to the interaction of multiple capillary waves. We show that a single, non-dimensional, and time-independent parameter controls the range of validity of the theory. Our findings have wide-ranging implications for the understanding of the mixing of immiscible fluids. This includes in particular oil spills, where our work enables the prediction of the oil-water interface dynamics that ultimately determine the rate of oil biodegradation by marine bacteria.
RESUMEN
Drop-fiber interactions are fundamental to the operation of technologies such as atmospheric fog capture, oil filtration, refrigeration, and dehumidification. We demonstrate that by twisting together two fibers, a sliding drop's flow path can be controlled by tuning the ratio between its size and the twist wavelength. We find both experimentally and numerically that twisted fiber systems are able to asymmetrically stabilize drops, both enhancing drop transport speeds and creating a rich array of new flow patterns. We show that the passive flow control generated by twisting fibers allows for woven nets that can be "programmed" with junctions that predetermine drop interactions and can be anticlogging. Furthermore, it is shown that twisted fiber structures are significantly more effective at capturing atmospheric fog compared to straight fibers.
RESUMEN
This study explores the impact of electrostatic interactions and hydrogen bonding on tear film stability, a crucial factor for ocular surface health. While mucosal and meibomian layers have been extensively studied, the role of electrolytes in the aqueous phase remains unclear. Dry eye syndrome, characterized by insufficient tear quantity or quality, is associated with hyperosmolality, making electrolyte composition an important factor that might impact tear stability. Using a model buffer solution on a silica glass dome, we simulated physiologically relevant tear film conditions. Sodium chloride alone induced premature dewetting through salt crystal nucleation. In contrast, trace amounts of solutes with hydroxyl groups (sodium phosphate dibasic, potassium phosphate monobasic, and glucose) exhibited intriguing phenomena: quasi-stable films, solutal Marangoni-driven fluid influx increasing film thickness, and viscous fingering due to Saffman-Taylor instability. These observations are rationalized by the association of salt solutions with increased surface tension and the propensity of hydroxyl-group-containing solutes to engage in significant hydrogen bonding, altering local viscosity. This creates a viscosity contrast between the bulk buffer solution and the film region. Moreover, these solutes shield the glass dome, counteracting sodium chloride crystallization. These insights not only advance our understanding of tear film mechanics but also pave the way for predictive diagnostics in dry eye syndrome, offering a robust platform for personalized medical interventions based on individual tear film composition.
Asunto(s)
Electrólitos , Enlace de Hidrógeno , Lágrimas , Lágrimas/química , Electrólitos/química , Humanos , Viscosidad , Cloruro de Sodio/química , Fosfatos/química , Tensión Superficial , Electricidad Estática , Síndromes de Ojo Seco/metabolismo , Humectabilidad , Compuestos de PotasioRESUMEN
The feeding mechanisms of animals constrain the spectrum of resources that they can exploit profitably. For floral nectar eaters, both corolla depth and nectar properties have marked influence on foraging choices. We report the multiple strategies used by honey bees to efficiently extract nectar at the range of sugar concentrations and corolla depths they face in nature. Honey bees can collect nectar by dipping their hairy tongues or capillary loading when lapping it, or they can attach the tongue to the wall of long corollas and directly suck the nectar along the tongue sides. The honey bee feeding apparatus is unveiled as a multifunctional tool that can switch between lapping and sucking nectar according to the instantaneous ingesting efficiency, which is determined by the interplay of nectar-mouth distance and sugar concentration. These versatile feeding mechanisms allow honey bees to extract nectar efficiently from a wider range of floral resources than previously appreciated and endow them with remarkable adaptability to diverse foraging environments.
Asunto(s)
Boca , Néctar de las Plantas , Abejas , Animales , Lengua , Carbohidratos , AzúcaresRESUMEN
We show that unconstrained asymmetric dissolving solids floating in a fluid can move rectilinearly as a result of attached density currents which occur along their inclined surfaces. Solids in the form of boats composed of centimeter-scale sugar and salt slabs attached to a buoy are observed to move rapidly in water with speeds up to 5 mm/s determined by the inclination angle and orientation of the dissolving surfaces. While symmetric boats drift slowly, asymmetric boats are observed to accelerate rapidly along a line before reaching a terminal velocity when their drag matches the thrust generated by dissolution. By visualizing the flow around the body, we show that the boat velocity is always directed opposite to the horizontal component of the density current. We derive the thrust acting on the body from its measured kinematics and show that the propulsion mechanism is consistent with the unbalanced momentum generated by the attached density current. We obtain an analytical formula for the body speed depending on geometry and material properties and show that it captures the observed trends reasonably. Our analysis shows that the gravity current sets the scale of the body speed consistent with our observations, and we estimate that speeds can grow slowly as the cube root of the length of the inclined dissolving surface. The dynamics of dissolving solids demonstrated here applies equally well to solids undergoing phase change and may enhance the drift of melting icebergs, besides unraveling a primal strategy by which to achieve locomotion in active matter.
RESUMEN
For quantum computing (QC) to emerge as a practically indispensable computational tool, there is a need for quantum protocols with end-to-end practical applications-in this instance, fluid dynamics. We debut here a high-performance quantum simulator which we term QFlowS (Quantum Flow Simulator), designed for fluid flow simulations using QC. Solving nonlinear flows by QC generally proceeds by solving an equivalent infinite dimensional linear system as a result of linear embedding. Thus, we first choose to simulate two well-known flows using QFlowS and demonstrate a previously unseen, full gate-level implementation of a hybrid and high precision Quantum Linear Systems Algorithms (QLSA) for simulating such flows at low Reynolds numbers. The utility of this simulator is demonstrated by extracting error estimates and power law scaling that relates [Formula: see text] (a parameter crucial to Hamiltonian simulations) to the condition number [Formula: see text] of the simulation matrix and allows the prediction of an optimal scaling parameter for accurate eigenvalue estimation. Further, we include two speedup preserving algorithms for a) the functional form or sparse quantum state preparation and b) in situ quantum postprocessing tool for computing nonlinear functions of the velocity field. We choose the viscous dissipation rate as an example, for which the end-to-end complexity is shown to be [Formula: see text], where [Formula: see text] is the size of the linear system of equations, [Formula: see text] is the solution error, and [Formula: see text] is the error in postprocessing. This work suggests a path toward quantum simulation of fluid flows and highlights the special considerations needed at the gate-level implementation of QC.
RESUMEN
A quantitative understanding of the coupled dynamics of flow and particles in aerosol and droplet transmission associated with speech remains elusive. Here, we summarize an effort that integrates insights into flow-particle dynamics induced by the production plosive sounds during speech with skin-integrated electronic systems for monitoring the production of these sounds. In particular, we uncover diffusive and ballistic regimes separated by a threshold particle size and characterize the Lagrangian acceleration and pair dispersion. Lagrangian dynamics of the particles in the diffusive regime exhibit features of isotropic turbulence. These fundamental findings highlight the value in skin-interfaced wireless sensors for continuously measuring critical speech patterns in clinical settings, work environments, and the home, based on unique neck biomechanics associated with the generation of plosive sounds. We introduce a wireless, soft device that captures these motions to enable detection of plosive sounds in multiple languages through a convolutional neural network approach. This work spans fundamental flow-particle physics to soft electronic technology, with implications in monitoring and studying critical speech patterns associated with aerosol and droplet transmissions relevant to the spread of infectious diseases.
Asunto(s)
Electrónica , Habla , Aerosoles , Tamaño de la Partícula , Movimiento (Física)RESUMEN
Most acute coronary syndromes are due to a sudden luminal embolism caused by the rupturing or erosion of atherosclerotic plaques. Prevention and treatment of plaque development have become an effective strategy to reduce mortality and morbidity from coronary heart disease. It is now generally accepted that plaques with thin-cap fibroatheroma (TCFA) are precursors to rupturing and that larger plaques and high-risk plaque features (including low-attenuation plaque, positive remodeling, napkin-ring sign, and spotty calcification) constitute unstable plaque morphologies. However, plaque vulnerability or rupturing is a complex evolutionary process caused by a combination of multiple factors. Using a combination of medicine, engineering mechanics, and computer software, researchers have turned their attention to computational fluid mechanics. The importance of fluid mechanics in pathological states for promoting plaque progression, inducing plaque tendency to vulnerability, or even rupture, as well as the high value of functional evaluation of myocardial ischemia has become a new area of research. This article reviews recent research advances in coronary plaque fluid mechanics, aiming to describe the concept, research implications, current status of clinical studies, and limitations of fluid mechanic's characteristic parameters: wall shear stress (WSS), axial plaque shear (APS), and fractional flow reserve (FFR). Previously, most computational fluid dynamics were obtained using invasive methods, such as intravascular ultrasound (IVUS) or optical coherence tomography (OCT). In recent years, the image quality and spatial resolution of coronary computed tomography angiography (CCTA) have greatly improved, making it possible to compute fluid dynamics by noninvasive methods. In the future, the combination of CCTA-based anatomical stenosis, plaque high-risk features, and fluid mechanics can further improve the prediction of plaque development, vulnerability, and risk of rupturing, as well as enabling noninvasive means to assess the degree of myocardial ischemia, thereby providing an important aid to guide clinical decision-making and optimize treatment.
RESUMEN
Fungal pellets are hierarchical systems that can be found in an ample variety of applications. Modeling transport phenomena in this type of systems is a challenging but necessary task to provide knowledge-based processes that improve the outcome of their biotechnological applications. In this work, an upscaled model for total mass and momentum transport in fungal pellets is implemented and analyzed, using elements of the volume averaging and adjoint homogenization methods departing from the governing equations at the microscale in the intracellular and extracellular phases. The biomass is assumed to be composed of a non-Newtonian fluid and the organelles impervious to momentum transport are modeled as a rigid solid phase. The upscaled equations contain effective-medium coefficients, which are predicted from the solution of adjoint closure problems in a three-dimensional periodic domains representative of the microstructure. The construction of these domains was performed for Laccaria trichodermophora based on observations of actual biological structures. The upscaled model was validated with direct numerical simulations in homogeneous portions of the pellets core. It is shown that no significant differences are observed when the dolipores are open or closed to fluid flow. By comparing the predictions of the average velocity in the extracellular phase resulting from the upscaled model with those from the classical Darcy equation (i.e., assuming that the biomass is a solid phase) the contribution of the intracellular fluid phase was evidenced. This work sets the foundations for further studies dedicated to transport phenomena in this type of systems.
Asunto(s)
Modelos Biológicos , Biomasa , Simulación por ComputadorRESUMEN
PURPOSE: The etiology of ureteral dilation in primary nonrefluxing, nonobstructing megaureters is still not well understood. Impaired ureteral peristalsis has been theorized as one of the contributing factors. However, ureteral peristalsis and its "normal" function is not well defined. In this study, using mathematical modeling techniques, we aim to better understand how ureteral peristalsis works. This is the first model to consider clinically observed, back-and-forth, cyclic wall longitudinal motion during peristalsis. We hypothesize that dysfunctional ureteral peristalsis, caused by insufficient peristaltic amplitudes (e.g., circular muscle dysfunction) and/or lack of ureteral wall longitudinal motion (e.g., longitudinal muscle dysfunction), promotes peristaltic reflux (i.e., retrograde flow of urine during an episode of peristalsis) and may result in urinary stasis, urine accumulation, and consequent dilation. METHODS: Based on lubrication theory in fluid mechanics, we developed a two-dimensional (planar) model of ureteral peristalsis. In doing so, we treated ureteral peristalsis as an infinite train of sinusoidal waves. We then analyzed antegrade and retrograde flows in the ureter under different bladder-kidney differential pressure and peristalsis conditions. RESULTS: There is a minimum peristaltic amplitude required to prevent peristaltic reflux. Ureteral wall longitudinal motion decreases this minimum required amplitude, increasing the nonrefluxing range of peristaltic amplitudes. As an example, for a normal bladder-kidney differential pressure of 5 cmH2 O, ureteral wall longitudinal motion increases nonrefluxing range of peristaltic amplitude by 65%. Additionally, ureteral wall longitudinal motion decreases refluxing volumetric flow rates. For a similar normal bladder pressure example of 5 cmH2 O, refluxing volumetric flow rate decreases by a factor of 18. Finally, elevated bladder pressure, not only increases the required peristaltic amplitude for reflux prevention but it increases maximum refluxing volumetric flow rates. For the case without wall longitudinal motion, as bladder-kidney differential pressure increases from 5 to 40 cmH2 O, minimum required peristaltic amplitude to prevent reflux increases by 40% while the maximum refluxing volumetric flow rate increases by approximately 100%. CONCLUSION: The results presented in this study show how abnormal ureteral peristalsis, caused by the absence of wall longitudinal motion and/or lack of sufficient peristaltic amplitudes, facilitates peristaltic reflux and retrograde flow. We theorize that this retrograde flow can lead to urinary stasis and urine accumulation in the ureters, resulting in ureteral dilation seen on imaging studies and elevated infection risk. Our results also show how chronically elevated bladder pressures are more susceptible to such refluxing conditions that could lead to ureteral dilation.
Asunto(s)
Uréter , Obstrucción Ureteral , Humanos , Peristaltismo/fisiología , Dilatación , Uréter/fisiología , Vejiga UrinariaRESUMEN
BACKGROUND: A percutaneous ventricular assist device (pVAD) is an effective method to treat heart failure, but its complications, mainly hemolysis and thrombus formation, cannot be ignored. Accurate evaluation of hemolysis and thrombus formation in pVAD is essential to guide the development of pVAD and reduce the incidence of complications. METHODS: This study optimized the numerical model to predict hemolysis and thrombus formation in pVAD. The hemolysis model is based on the power law function, and the multi-component thrombus prediction model is improved by introducing the von Willebrand factor. RESULTS: The error between the numerical simulation and the hydraulic performance experiment is within 5%. The numerical results of hemolysis are in good agreement with those of in vitro experiments. Meanwhile, the thrombus location predicted by the numerical model is the same as that found in the in vivo experiment. CONCLUSION: The numerical model suggested in this study may therefore accurately assess the possible hemolytic and thrombotic dangers in pVAD, making it an effective tool to support the development of pVAD.
Asunto(s)
Insuficiencia Cardíaca , Corazón Auxiliar , Trombosis , Humanos , Hemólisis , Corazón Auxiliar/efectos adversos , Insuficiencia Cardíaca/cirugía , Simulación por Computador , Trombosis/etiologíaRESUMEN
We present a mathematical model of an experiment in which cells are cultured within a gel, which in turn floats freely within a liquid nutrient medium. Traction forces exerted by the cells on the gel cause it to contract over time, giving a measure of the strength of these forces. Building upon our previous work (Reoch et al. in J Math Biol 84(5):31, 2022), we exploit the fact that the gels used frequently have a thin geometry to obtain a reduced model for the behaviour of a thin, two-dimensional cell-seeded gel. We find that steady-state solutions of the reduced model require the cell density and volume fraction of polymer in the gel to be spatially uniform, while the gel height may vary spatially. If we further assume that all three of these variables are initially spatially uniform, this continues for all time and the thin film model can be further reduced to solving a single, non-linear ODE for gel height as a function of time. The thin film model is further investigated for both spatially-uniform and varying initial conditions, using a combination of analytical techniques and numerical simulations. We show that a number of qualitatively different behaviours are possible, depending on the composition of the gel (i.e., the chemical potentials) and the strength of the cell traction forces. However, unlike in the earlier one-dimensional model, we do not observe cases where the gel oscillates between swelling and contraction. For the case of initially uniform cell and gel density, our model predicts that the relative change in the gels' height and length are equal, which justifies an assumption previously used in the work of Stevenson et al. (Biophys J 99(1):19-28, 2010). Conversely, however, even for non-uniform initial conditions, we do not observe cases where the length of the gel changes whilst its height remains constant, which have been reported in another model of osmotic swelling by Trinschek et al. (AIMS Mater Sci 3(3):1138-1159, 2016; Phys Rev Lett 119:078003, 2017).
Asunto(s)
Nutrientes , Polímeros , Geles , SemillasRESUMEN
A mathematical model of infection, inflammation and immune response in an idealized bronchial tree is developed. This work is based on a model from the literature that is extended to account for the propagation dynamics of an infection between the airways. The inflammation affects the size of the airways, the air flows distribution in the airway tree, and, consequently, the oxygen transfers to blood. We test different infections outcomes and propagation speed. In the hypotheses of our model, the inflammation can reduce notably and sometimes drastically the oxygen flow to blood. Our model predicts how the air flows and oxygen exchanges reorganize in the tree during an infection. Our results highlight the links between the localization of the infection and the amplitude of the loss of oxygen flow to blood. We show that a compensation phenomena due to the reorganization of the flow exists, but that it remains marginal unless the power produced the ventilation muscles is increased. Our model forms a first step towards a better understanding of the dynamics of bronchial infections.
Asunto(s)
Pulmón , Oxígeno , Humanos , Pulmón/fisiología , Inflamación , Modelos Biológicos , Modelos AnatómicosRESUMEN
We have demonstrated the effectiveness of reinforcement learning (RL) in bluff body flow control problems both in experiments and simulations by automatically discovering active control strategies for drag reduction in turbulent flow. Specifically, we aimed to maximize the power gain efficiency by properly selecting the rotational speed of two small cylinders, located parallel to and downstream of the main cylinder. By properly defining rewards and designing noise reduction techniques, and after an automatic sequence of tens of towing experiments, the RL agent was shown to discover a control strategy that is comparable to the optimal strategy found through lengthy systematically planned control experiments. Subsequently, these results were verified by simulations that enabled us to gain insight into the physical mechanisms of the drag reduction process. While RL has been used effectively previously in idealized computer flow simulation studies, this study demonstrates its effectiveness in experimental fluid mechanics and verifies it by simulations, potentially paving the way for efficient exploration of additional active flow control strategies in other complex fluid mechanics applications.
Asunto(s)
Algoritmos , Simulación por Computador , Ambiente , Aprendizaje/fisiología , Modelos Biológicos , Refuerzo en Psicología , Humanos , Fenómenos FísicosRESUMEN
Crustacean filter feeders capture oil droplets with the use of their ramified appendages. These appendages behave as paddles or sieves, based on the system's Reynolds number. Here, we used high-speed videography, scanning electron microscopy and fluid mechanics to study the capturing mechanisms of crude oil droplets and the filtering appendage's wettability by two species of barnacles (Balanus glandula and Balanus crenatus) and of the freshwater cladoceran Daphnia magna. Our results show that barnacle appendages behave as paddles and capture droplets in their boundary layers at low Reynolds number. At high Reynolds number, droplets are most likely to be captured via direct interception. There is an intermediate range of Reynolds number where droplets can be captured by both mechanisms at the same time. Daphnia magna captures droplets in the boundary layers of the third and fourth pair of thoracic legs with a metachronal motion of the appendages. All studied surfaces were revealed to be highly lipophobic, demonstrating captured oil droplets with high contact angles. We also discuss implications of such capture mechanisms and wettability on potential ingestion of crude oil by filter feeders. These results further our understanding of the capture of crude oil by filter feeders, shedding light on the main entry point of oil in marine food webs.
Asunto(s)
Contaminación por Petróleo , Petróleo , Thoracica , Animales , Daphnia , Cadena AlimentariaRESUMEN
Biological tissues are composed of cells surrounded by the extracellular matrix (ECM). The ECM can be thought of as a fibrous polymer network, acting as a natural scaffolding to provide mechanical support to the cells. Reciprocal mechanical and chemical interactions between the cells and the ECM are crucial in regulating the development of tissues and maintaining their functionality. Hence, to maintain in vivo-like behaviour when cells are cultured in vitro, they are often seeded in a gel, which aims to mimic the ECM. In this paper, we present a mathematical model that incorporates cell-gel interactions together with osmotic pressure to study the mechanical behaviour of biological gels. In particular, we consider an experiment where cells are seeded within a gel, which gradually compacts due to forces exerted on it by the cells. Adopting a one-dimensional Cartesian geometry for simplicity, we use a combination of analytical techniques and numerical simulations to investigate how cell traction forces interact with osmotic effects (which can lead to either gel swelling or contraction depending on the gel's composition). Our results show that a number of qualitatively different behaviours are possible, depending on the composition of the gel (i.e. its chemical potentials) and the strength of the cell traction forces. A novel prediction of our model is that there are cases where the gel oscillates between swelling and contraction; to our knowledge, this behaviour has not been reported in experiments. We also consider how physical parameters like drag and viscosity affect the manner in which the gel evolves.
Asunto(s)
Matriz Extracelular , Modelos Teóricos , Matriz Extracelular/fisiología , Geles/análisis , ViscosidadRESUMEN
The Navier slip condition describes the motion of a liquid relative to a neighboring solid surface, with its characteristic Navier slip length being a constitutive property of the solid-liquid interface. Measurement of this slip length is complicated by its small magnitude, expected to be in the nanometer range based on molecular simulations. Here, we report an experimental technique that interrogates the Navier slip length on individual nanoparticles immersed in liquid with subnanometer precision. Proof-of-principle experiments on individual, citrate-stabilized, gold nanoparticles in water give a constant slip length of 2.7 ± 0.6 nm (95% C.I.), independent of particle size. Achieving this feature of size independence is central to any measurement of this constitutive property, which is facilitated through the use of individual particles of varying radii. This demonstration motivates studies that can now validate the wealth of existing molecular simulation data on slip.
Asunto(s)
Oro , Nanopartículas del Metal , Tamaño de la Partícula , Propiedades de Superficie , AguaRESUMEN
The study of the evolution of the atmosphere requires careful consideration of multicomponent gaseous flows under gravity. The gas dynamics under an external force field is usually associated with an intrinsic multiscale nature due to large particle density variation along the direction of force. A wonderfully diverse set of behaviors of fluids can be observed in different flow regimes. This poses a great challenge for numerical algorithms to accurately and efficiently capture the scale-dependent flow physics. In this paper, a well-balanced unified gas-kinetic scheme (UGKS) for a gas mixture is developed, which can be used for the study of cross-scale multicomponent flows under an external force field. The well-balanced scheme here indicates the capability of a numerical method to evolve a gravitational system under any initial condition to the hydrostatic equilibrium and to keep such a solution. Such a property is crucial for an accurate description of multicomponent gas evolution under an external force field, especially for long-term evolving systems such as galaxy formation. Based on the Boltzmann model equation for gas mixtures, the UGKS leverages the space-time integral solution to construct numerical flux functions and, thus, provides a self-conditioned mechanism to recover typical flow dynamics in various flow regimes. We prove the well-balanced property of the current scheme formally through theoretical analysis and numerical validations. New physical phenomena, including the decoupled transport of different gas components in the transition regime, are presented and studied.
RESUMEN
The human eye is a complex biomechanical structure with a range of biomechanical processes involved in various physiological as well as pathological conditions. Fluid flow inside different domains of the eye is one of the most significant biomechanical processes that tend to perform a wide variety of functions and when combined with other biophysical processes play a crucial role in ocular drug delivery. However, it is quite difficult to comprehend the effect of these processes on drug transport and associated treatment experimentally because of ethical constraints and economic feasibility. Computational modeling on the other hand is an excellent means to understand the associated complexity between these aforementioned processes and drug delivery. A wide range of computational models specific to different types of fluids present in different domains of the eye as well as varying drug delivery modes has been established to understand the fluid flow behavior and drug transport phenomenon in an insilico manner. These computational models have been used as a non-invasive tool to aid ophthalmologists in identifying the challenges associated with a particular drug delivery mode while treating particular eye diseases and to advance the understanding of the biomechanical behavior of the eye. In this regard, the author attempts to summarize the existing computational and mathematical approaches proposed in the last two decades for understanding the fluid mechanics and drug transport associated with different domains of the eye, together with their application to modify the existing treatment processes.
Asunto(s)
Sistemas de Liberación de Medicamentos/métodos , Oftalmopatías/tratamiento farmacológico , Ojo/fisiopatología , Modelos Biológicos , Administración Oftálmica , Disponibilidad Biológica , Fenómenos Biomecánicos , Simulación por Computador , Ojo/metabolismo , Oftalmopatías/fisiopatología , Humanos , Distribución TisularRESUMEN
Mathematical modeling of surface deformation during cryopreservation by vitrification is presented in this study. The specific problem under consideration is of a cryoprotective agent (CPA) solution vitrifying in a vial, following previously obtained cryomacroscopy observations. A multiphysics solution is proposed in this study, combining coupled effects associated with heat transfer, fluid mechanics, and solid mechanics. Consistent with previous investigations, this study demonstrates that surface deformation is the result of material flow, which is the combined outcome of temperature gradients developed during the inward cooling process, the tendency of the material to change its volume with temperature, and the exponential increase in material viscosity with the decreasing temperature. During this process, the behavior of the CPA changes from liquid to a solid-like amorphous material, where the arrested flow in the vitrified state results in mechanical stresses. Results of this study show a good qualitative agreement of surface deformation with previously obtained experimental data, and support prior investigations to explain fracture tendencies propagating from the deformed surface. Results of this study also highlight the effect of heat convection in the CPA at the early stage of cooling.