Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mesoscale systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Turbulence and Cloudiness in Cumulus-Topped Marine Boundary Layers

Turbulence and cloudiness in cumulus-topped marine boundary layers are studied using data from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site. The analysis includes eight periods of nonprecipitating shallow cumulus clouds, spanning 141 h and encompassing 603 individual clouds. On average, the cumulus had bases at 558 m, were 99 m thick, had a chord length of 500 m, and exhibited an hourly base-layer cloudiness of 12%. Changes in cloud fraction were primarily driven by variations in cloud number not chord length. High-resolution Doppler lidar and cloud radar observations were used to estimate updraft and downdraft mass fluxes in the cumulus base layer under both clear and cloudy regions. Cloudy updrafts contributed only ∼23% of the total updraft mass flux at cloud base, indicating that most upward mass flux originated from clear-air eddies. Updraft strength, rather than updraft fraction, was found to predominantly control both clear and cloudy updraft mass fluxes. Cloud-base cloudiness strongly correlated with cloudy updraft mass flux but showed negligible correlation with clear-air updraft mass flux. Clear-air downdraft mass flux exhibited a strong relationship with the ratio of surface buoyancy to inversion strength. Mesoscale analysis revealed that moist patches had greater cloudiness, updraft mass flux, and vertical velocity variance compared to dry patches. Additionally, mesoscale base-layer cloud fraction was highly and significantly correlated with cloudy updraft mass flux in both dry and moist environments. Results presented herein have implications for cumulus parameterization development along with climatological and model evaluation studies conducted at the ENA site.

cumulus clouds↗
DOE OSTI2026

Lidar Measurements and High-Resolution Mesoscale Modeling of Coastally Trapped Disturbances off the Coast of California

Coastally Trapped disturbances (CTDs) are shifts in wind direction from the pre-dominant direction to equatorward to poleward for a period of time. These CTDs occur during the warm season off the California coast and impact coastal weather conditions and planned offshore wind plants. This study assesses the characteristics of CTD events as observed by lidar and other offshore buoys, then evaluates the ability of modeling systems to capture the correct characteristics, leveraging model output from the High-Resolution Rapid Refresh (HRRR) operational modeling system and the NOW-23 (National Offshore Wind) model dataset. CTDs were analyzed for October 2020 and May through to October of 2021, identifying 18 unique CTD events, confirmed by a nearby National Data Buoy Center (NDBC) buoy. The HRRR model captured most of these events, but the NOW-23 model output contained only 12 events. Composites of the wind, temperature, and pressure perturbations pre-, during, and post-event demonstrated the diminishment in wind speed, particularly for the alongshore component. Although the NOW-23 model captured the alongshore wind component and pressure perturbations well, the cross-shore wind component and temperature perturbations varied substantially. When the turbulent kinetic energy deviation and wind shear was positive across all levels pre-event, the NOW-23 modeling system was less likely to capture the CTD event. In contrast, the events that were captured by the model tended to have negative wind shear aloft pre-event.

coastal meteorology↗
DOE OSTI2022

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Mesoscopic-scale complexity in macroscopically-uniform plastic flow of an Al 0.3 CoCrFeNi high-entropy alloy

The concept of a smooth and homogeneous plastic flow of solids is nowadays constantly challenged by various observations of the self-organization of crystal defects on mesoscopic scales pertaining, e.g., to acoustic emission or the evolution of the local strain field. Such investigations would be of particular interest for High-Entropy Alloys (HEAs) characterized by extremely complex microstructures. However, the complexity of their deformation has been only studied in relation to the jerky flow. The present work explores a mesoscale complexity of the macroscopically-smooth plastic flow of an Al 0.3 CoCrFeNi alloy. Here, an inclusive approach involving various experimental methods and complementary nonlinear analyses allowed to reveal transitions between distinct dynamical regimes, which may be assimilated with noises of different colors, in particular the blue noise that is rarely observed in complex systems, and various kinds of reddened behavior. Such findings provide significant new insights into the micro/macro transition in the deformation behavior of HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE)

Convective clouds play an important role in Earth’s climate system and are a known source of extreme weather. Gaps in our understanding of convective vertical motions, microphysics, and precipitation across a full range of aerosol and meteorological regimes continue to limit our ability to predict the occurrence and intensity of these cloud systems. To improve predictability, the National Science Foundation (NSF) sponsored a large field experiment entitled “Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE).” ESCAPE took place between 30 May and 30 September 2022 in the vicinity of Houston, Texas, because this area frequently experiences isolated deep convection that interacts with the region’s mesoscale circulations and its range of aerosol conditions. ESCAPE focused on collecting observations of isolated deep convection through innovative sampling and developing novel analysis techniques. This included the deployment of two research aircraft, the National Research Council of Canada Convair-580 and the Stratton Park Engineering Company Learjet, which combined conducted 24 research flights from 30 May to 17 June. On the ground, three mobile X-band radars and one mobile Doppler lidar truck equipped with soundings were deployed from 30 May to 28 June. From 1 August to 30 September 2022, a dual-polarization C-band radar was deployed and operated using a novel, multisensor agile adaptive sampling strategy to track the entire life cycle of isolated convective clouds. Analysis of the ESCAPE observations has already yielded preliminary findings on how aerosols and environmental conditions impact the convective life cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗
DOE OSTI2023

Derivation of a Multi-Species Mix Model

Multi-species flow takes place when two or more gases, plasmas or fluids are flowing and interacting with one another through particle collisions, chemical reactions, turbulence, and thermonuclear processes; this list is not exhaustive. It has a wide range of applications in geological, biological and engineering systems. In what follows species refers to different species of atoms or molecules that may or may not be partially or fully ionized. The interactions between these species takes place at an atomic or molecular level. The species are electrically neutral. These are very important points, because the coupling between species that are interacting on a mesoscale level or are immiscible would be very different than what is being considered in this work. In addition, in section 4, we present a heuristic argument to enhance the collisional mean free path to incorporate the effects of turbulent eddies by use of Prandtl’s assumption that the integral scale length of a turbulent boundary layer is proportional to the width of the turbulent layer. This leads to coupling terms that are non-Newtonian, because they are a function of the dynamics of the flow. It is also assumed that the turbulent stress terms are much smaller than the species partial pressures. This may be a good assumption if the species are hot, dense plasmas that exist in ICF targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Multiscale characterization of phase change materials for building thermal energy storage applications

Phase change materials (PCMs) store and release large amounts of thermal energy because of their high latent energy storage capacity. However, long-term cyclic stability, supercooling and performance-scalability are some of the major challenges for their use in building thermal energy storage (TES) applications. Here, in this study, we present a comprehensive multiscale characterization of two commercially available organic PCMs, Puretemp 18 and Puretemp 23. At the microscale, differential scanning calorimetry (DSC) was used to characterize phase change temperature, specific heat, and latent heat. At the mesoscale, a heat flow meter apparatus (HFMA), following the ASTM C1784 standard, was employed to measure the phase change temperature, specific heat, and latent heat properties. A comparative analysis of latent heat as a function of temperature was conducted by integrating the DSC and HFMA results. At the macroscale, the thermal performance and cyclic stability of the TES system was evaluated using Puretemp 23. The TES system consisted of a finned tube heat exchanger with a storage volume of 0.0189 m 3 (5 gal), which represents a compact, real-world TES solution suitable for building energy storage. The results showed consistent thermal stability of the PCM over 200 cycles, and the supercooling temperature remained within 0.2 °C, which was not detected in smaller-scale characterization methods. Additionally, the macroscale testing methodology of the PCM revealed that the TES is able to charge and discharge stored latent energy within 2 h under a temperature differential of 16.67 °C measured between the inlet water temperature and the phase transition temperature of the PCM. The proposed multiscale PCM characterization method provides a systematic basis for comparing important thermal storage properties while also investigating the scalability, reliability and integration challenges in large scale TES applications.

Latent heat↗
DOE OSTI2021

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Phase-field simulation of domain size effect on dielectric and piezoelectric responses in K 0.5 Na 0.5 NbO 3 epitaxial thin films with superdomain structures

Size effects of mesoscale ferroelectric domains on the macroscopic dielectric and piezoelectric responses in domain-engineered bulk piezocrystals have been extensively studied for more than two decades. However, less is known about the domain size effects in ferroelectric epitaxial thin films, especially for films consisting of low-symmetry ferroelectric phases and exhibits hierarchical superdomain structures. Herein, using phase-field simulations, we systemically evaluate the effective out-of-plane dielectric and piezoelectric coefficients, κ$^{*}_{33}$ and d$^{*}_{33}$, as a function of the domain periods for two types of superdomain structures in ferroelectric K 0.5 Na 0.5 NbO 3 epitaxial thin films. In one type of the superdomain structures, we find that more than 70% increase of κ$^{*}_{33}$ and nearly 20% increase of d$^{*}_{33}$ can be achieved by tuning the domain period by a few tens of nanometers. Dissimilar behaviors are found when the domain period varies along different lateral directions, suggesting anisotropic domain size effects in thin films. By analyzing the local dielectric and piezoelectric responses from each domain variants and domain walls, we reveal that the domain size effect is governed by the variation of out-of-plane polarization inside the domains. Moreover, we also demonstrate enhanced domain size effects by modulating the misfit strains and temperature to approach the polymorphic phase boundaries, suggesting tunability of the size effect by doping and strain engineering. Furthermore, our results reveal multimodal domain size dependence of dielectric and piezoelectric responses in low-symmetry ferroelectric epitaxial thin films, implying that domain size engineering can be used to tune macroscopic properties of thin-film ferroelectrics, similar to their bulk counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗
DOE OSTI2021

The interpenetration polymer network in a cement paste–waterborne epoxy system

The formation of the interpenetrating polymer network (IPN) structure within the cement-polymer system has been revealed by experiments from 2 dimensional to 3-dimensional scales. However, the microstructure design and performance prediction of IPN as a function of specific C-S-H/polymer components or ingredients' parameters e.g. water to cement ratios (w/c), polymer to cement ratios (p/c), monomer components, degree of polymerization (DP), etc. is by far not available. Here we developed a mesoscale model for IPN visualization based on the Flory-Huggins interaction theory which was applied to cement science for the first time. All the ingredients in the micro-structure were considered as soluble beads with rational sizes based on their properties obtained by molecular dynamic methods. The interaction parameters of each bead were then determined based on their element ratios and chemical backbones. The model was validated with the waterborne epoxy-cement material (WECM) in which a novel waterborne epoxy resin (WEP) was prepared and impregnated. The verification results on the 2D-3D scale show that the developed model predicts the WECM's IPN structures by rule and line concerning DP, w/c, and p/c in the mixture. The C-S-H beads were progressively scaled in size which alters the C-S-H texture to have completely different dissolution characteristics. Beads of diameter ~6 Å are more unstable and soluble which enable them to form a continuous phase in water or a composite structure with WEP. In contrast, beads with diameters larger than 10 Å have different properties with stronger nucleation effects. The results also suggest the impregnation content of WEP in cement-based material should be limited to 10% vol. to prevent the polymer IPN from decomposing into discrete clusters. The application of Flory-Huggins theory in cement-based composites demonstrates great potential in performance prediction and microstructure design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Observational Study of the Thermodynamics and Morphological Characteristics of a Midlatitude Continental Cold Pool Event

Abstract Convectively generated cold pools are important to the Earth system as they exert strong controls on deep convective-storm initiation, intensity, and life cycle. Despite their importance, efforts to introduce such cold pool controls into weather and climate models lack guidance and/or physical constraints from cold pool observations. This work presents a detailed, purely observational analysis of a cold pool event that took place on 23–24 May 2011 in north-central Oklahoma. The characteristics of the cold pool, and the spatiotemporal evolution of the hydrometeors and dynamics in the proximity of the cold pool, are studied with high-resolution observations. The unprecedented dataset used in this work to study cold pool characteristics includes an enhanced network of surface weather stations, a high-temporal-frequency sounding array, and the NEXRAD and Atmospheric Radiation Measurement (ARM) Southern Great Plains radar networks. The potential use of NEXRAD surveillance scans to estimate height and propagation speed of the leading edge of the cold pool (LECP) is presented in this work. Manual identification and tracking of the LECP from NEXRAD imagery shows a spatial and temporal heterogeneity of the LECP properties. Surprisingly, over its detected life cycle, the LECP speed remains almost constant, even though the strength of the cold pool diminishes in time and its height varies. Radar analysis shows that pulses of graupel and hail within downdrafts in the convective system generating the cold pool appeared to be related to temporary increases in the LECP height.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving the representation of isopycnal mixing in E3SM (Final Report)

The DOE-supported Energy Exascale Earth System Model (E3SM) is a major attempt by the Department of Energy to develop a new Earth System Model using unstructured grids, which allow for the model to concentrate resolution where it is most needed and avoid numerical artifacts experienced by previous Earth System Models (ESMs) using regular grids in which the lateral resolution becomes extremely fine as the resolution shrinks at the poles. This effort required rewriting many of the algorithms previously developed for regular grids. One algorithm that was problematic in the first version of E3SM was the representation of mixing due to turbulent ocean eddies with scales smaller than the model grid. These mesoscale eddies are the primary way by which tracers are stirred along density surfaces in the ocean interior. In particular, the process of isopycnal mixing exchanges fresher, more oxygenated waters from polar regions with saltier, nutrient-rich and oxygen poor waters from tropical regions. Previous work in our group has shown that this mechanism is vitally important for bringing oxygen into poorly ventilated tropical regions (Gnanadesikan et al., 2012, 2013) and plays a significant role in determining the uptake of anthropogenic carbon dioxide (Gnanadesikan et al., 2015). However, in E3SMv1 this process was turned off, as turning it on caused the model to become unstable. Additionally, the rate of isopycnal mixing is directly proportional to a mixing coefficient $A_{Redi}$ whose value varies from less than 400 m 2 /s to 2000 m 2 /s across contemporary climate models. The main thrusts of this proposal therefore were to: 1.) Improve the numerics of isopycnal mixing in E3SM; 2.) Identify ways in which climate models are sensitive to the isopycnal mixing coefficient; and 3.) Develop better representations of the isopycnal mixing coefficient.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Flower‐Type Organized Trade‐Wind Cumulus: A Multi‐Day Lagrangian Large Eddy Simulation Intercomparison Study

Shallow cumulus cloud fields in subtropical marine trade wind environments, particularly over the tropical Atlantic Ocean, show distinct organizational patterns. Among these, Flower‐type clouds are characterized by expansive stratiform cloud patches surrounded by regions of scattered convection. The objectives of this study were (a) to construct a case study of a time period during the EUREC 4 A/ATOMIC field campaign when Flower‐type organization was observed, (b) to evaluate the fidelity of a multi‐model ensemble of large eddy simulations of that case, and (c) to analyze the interaction between cloud and precipitation processes and mesoscale organization in the simulations. The simulations follow a quasi‐Lagrangian trajectory, allowing mesoscale features to develop over time in a domain that follows the boundary‐layer airmass. The results show a broad agreement in simulated thermodynamic properties across different LES codes, with Flower‐type cloud patches appearing within hours of each other. The consensus among models is consistent with observations made during the EUREC 4 A/ATOMIC field campaign on the specific day of interest. The cloud structure reveals three distinct peaks in the joint probability densities of cloud base and cloud top height, with the dominant peak at any given time influenced by the stage of cloud organization. The simulated cloud system evolution reveals consistent occurrence of maxima in liquid water path and rain rate before Flower reaches its maximum length scale. Targeted sensitivity tests reveal a weak relationship between Cloud Droplet Number concentration and the extent/degree/type of organization.

EUREC4A↗
DOE OSTI2023

CEGANN: CRYSTAL EDGE GRAPH ATTENTION NEURAL NETWORK

SF-22-156 Machine learning (ML) models and applications in materials design and discovery typically involve the use of feature representations or descriptors followed by a learning algorithm that maps them to user desired properties of interest. Most popular mathematical formulation-based descriptors are not unique across atomic environments and suffer from transferability issues across different application domains and/or material classes. The CEGANN code provides a unified interface to facilitate material characterization across materials across multiple scales (from atomic to mesoscale) and diverse classes of materials ranging from metals oxides, non-metals, and even hierarchical materials such as zeolites and semi ordered materials such as mesophases. CEGANN implements a Graph Attention Network (GAT) type convolution architecture. The details of network architecture can be found in the paper https://doi.org/10.48550/arXiv.2207.10168. The software comes with pretrained examples and dataset for the classification of the following representative systems: (1) Structure-level representation such as space group (2) Structural dimensionality (e.g., bulk, 2D, clusters etc.) (3) Grain boundary identification (4) Nucleation and growth of a zeolite polymorph (5) Characterization of binary mesophases and their phase transitions (6) Growth of ice. The code is written in python programming language.

CHAN, HENRYT↗