SEARCH · Engineering Papers

Results for “mass function”

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.

SourceResearch connectionFrom yearThrough year

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

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Bulk Fermi surfaces of the Dirac type-II semimetallic candidate NiTe 2

Here, we present a study on the Fermi surface of the Dirac type-II semimetallic candidate NiTe 2 via the temperature and angular dependence of the de Haas–van Alphen (dHvA) effect measured in single crystals grown through Te flux. In contrast to its isostructural compounds like PtSe 2 , band-structure calculations predict NiTe 2 to display a tilted Dirac node very close to its Fermi level that is located along the Γ to A high-symmetry direction within its first Brillouin zone. The angular dependence of the dHvA frequencies is found to be in agreement with the first-principles calculations when the electronic bands are slightly shifted with respect to the Fermi level (ε F ), and therefore provide support for the existence of a Dirac type-II node in NiTe 2 . Nevertheless, we observed mild disagreements between experimental observations and density functional theory (DFT) calculations as, for example, nearly isotropic and light experimental effective masses. This indicates that the dispersion of the bands is not well captured by DFT. Despite the coexistence of Dirac-like fermions with topologically trivial carriers, samples of the highest quality display an anomalous and large either linear or sublinear magnetoresistivity. Finally, this suggests that Lorentz invariance breaking Dirac-like quasiparticles dominate the carrier transport in this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Electromagnetic proton–neutron mass difference

In this paper, we discuss the Cottingham formula and evaluate the proton–neutron electromagnetic mass difference exploiting the state-of-the-art phenomenological input.We decompose individual contributions to the mass splitting into Born, inelastic and subtraction terms. We evaluate the subtraction-function contribution connecting the input based on experimental data with the operator product expansion matched to QCD which allows us to avoid model dependence and to reduce errors of this contribution. We evaluate inelastic and Born terms accounting for modern low-$Q^2$ data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Modeling Electrochemical and Rheological Characteristics of Suspension-Based Electrodes for Redox Flow Cells

Flowable suspension-based electrodes (FSEs) have gained attention in recent years, as the integration of solid materials into electrochemical flow cells can offer improved performance and flexible operation. However, under conditions that engender favorable electrochemical properties (e.g., high particle loading, high conductivity, high surface area), FSEs can exhibit non-Newtonian characteristics that impose large pumping losses and flow-dependent transport rates. These multifaceted trade-offs motivate the use of models to broadly explore scaling relationships and better understand design rules for electrochemical devices. To this end, we present a one-dimensional model, integrating porous electrode theory with FSE rheology as well as flow-dependent electron and mass transport under pressure-driven flow. We study FSE behavior as a function of material properties and operating conditions, identifying key dimensionless groups that describe the underlying physical processes. We assess flow cell performance by quantifying electrode polarization and relative pumping losses, establishing generalized property-performance relationships for FSEs. Importantly, we expound relevant operating regimes—based on a subset of dimensionless groups—that inform practical operating envelopes, ultimately helping to guide FSE and cell engineering for electrochemical systems.

25 ENERGY STORAGE↗

DOE OSTI2025

Mapping Spatiotemporal Solvent Velocity from Measured Concentration Gradients in a Polarized Electrolyte

The electric-field induced motion of neutral species impedes the efficacy of electrochemical devices. By combining operando X-ray transmission measurements with continuum mechanics, we have developed a methodology for determining the velocity of neutral solvent molecules under an applied field. The X-ray transmission experiments were used to determine ion concentration profiles as a function of space and time in a polymer electrolyte. The unsteady state solvent mass balance equation was solved numerically with experimental concentration profiles to map spatiotemporal solvent velocities. We compare our experimentally derived results with predictions made with concentrated solution theory. We use the cation transference number as the only adjustable parameter to match experimental measurements of both concentration and solvent velocity. Our approach may be used to determine solvent velocity with any operando technique used to measure time-dependent ion concentration profiles.

Abdo, Emily E↗

DOE OSTI2026

Mathematical Modeling of the Potential and Time Dependence of Ir Dissolution from Hydrous Ir Oxide Oxygen Evolution Catalysts

One of the main degradation mechanisms of hydrous iridium oxide acidic oxygen evolution reaction (OER) catalysts is dissolution and loss into the acidic membrane. While degradation models have been proposed, there is a gap in understanding the potential and time dependence of the iridium dissolution reaction and its mechanistic underpinnings. In this work, Ir dissolution rates measured as a function of time and potential via time-resolved inductively-coupled plasma mass spectrometry (ICP-MS) in aqueous acidic electrolyte are used to establish a mathematical model for Ir dissolution. The mathematical model is generated using proposed formation and dissolution reactions for Ir species. Through comparison with the ICP-MS data and existing information on the potential-dependent Ir phase, we find that the potential and time-dependence of dissolution can be modeled as dissolution of an oxide phase, here represented as IrO 2 , with potential dependent kinetics and formation of a passivating species, a process with a rate-limiting step that is not potential dependent. This understanding of the potential dependence of dissolution and passivation kinetics using aqueous electrolyte half-cell measurements can be used to predict the degradation of Ir oxide in operating energy conversion devices relying on the OER, such as proton-exchange membrane water electrolyzers.

Kariuki, Nancy N. [Argonne National Laboratory (AN↗

DOE OSTI2021

Visualizing bleb mass dynamics in single cells using quantitative phase microscopy

Understanding biological responses to directed energy (DE) is critical to ensure the safety of personnel within the Department of Defense. At the Air Force Research Laboratory, we have developed or adapted advanced optical imaging systems that quantify biophysical responses to DE. One notable cellular response to DE exposure is the formation of blebs, or semi-spherical protrusions of the plasma membrane in living cells. In this work, we demonstrate the capacity of quantitative phase imaging (QPI) to both visualize and quantify the formation of membrane blebs following DE exposure. QPI is an interferometric imaging tool that uses optical path length as a label-free contrast mechanism and is sensitive to the non-aqueous mass density, or dry mass, of living cells. Blebs from both CHO-K1 and U937 cells were generated after exposure to a series of 600 ns, 21.2 kV/cm electric pulses. These blebs were visualized in real time, and their dry mass relative to the rest of the cell body was quantified as a function of time. It is our hope that this system will lead to an improved understanding of both DE-induced and apoptotic blebbing.

Steelman, Zachary A.↗

DOE OSTI2023

IDAES-PSE 2.2.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. New features in the 2.2 release include: DiagnosticsToolbox DiagnosticsToolbox (idaes.core.util.model_diagnostics. DiagnosticsToolbox is a new self-guiding toolbox for assisting with diagnosing modeling issues. The toolbox contains a wide range of tools for identifying common modeling issues and conveying these to the user, and more features will be added in the coming months. To try the new toolbox, use the following steps: from idaes.core.util import DiagnosticsToolbox Build an instance of your model Ensure your models has zero degrees of freedom (as best you can) dt = DiagnosticsToolbox(model) dt.report_structural_issues() Follow the prompts Deprecations DegeneracyHunter is being deprecated in favor of the new DiagnosticsToolbox. The features and capabilities of DegeneracyHunter will be moved over to the DiagnositcsToolbox in the coming months. Removal of usage of Pyomo's expr.current due to deprecation upstream. Removed support for radial basis functions in ALAMOpy as these did not work in the current implementation. Those wishing to use radial basis functions should use Pysmo instead. New Capabilities MSContactor model for unit operations involving mass transfer between 2 or more streams, such as membrane, solvent extraction and solids leaching operations. Other Improvements Improved formulation for solubility products with additional parameters for better scaling of terms Better backward compatibility, verification and robustness testing of core models and APIs Added support for Separator models without energy balances

DiagnosticsToolbox↗

DOE OSTI2022

The photon content of the proton in the CT18 global analysis

Recently, two photon PDF sets based on implementations of the LUX ansatz into the CT18 global analysis were released. In CT18lux, the photon PDF is calculated directly using the LUX master formula for all scales, \mu μ . In an alternative realization, CT18qed, the photon PDF is initialized at the starting scale, \mu_0 μ 0 , using the LUX formulation and evolved to higher scales \mu(>\mu_0) μ ( > μ 0 ) with a combined QED+QCD kernel at \mathcal{O}(\alpha),~\mathcal{O}(\alpha\alpha_s) 𝒪 ( α ) , 𝒪 ( α α s ) and \mathcal{O}(\alpha^2) 𝒪 ( α 2 ) . In the small- x x region, the photon PDF uncertainty is mainly induced by the quark and gluon PDFs, through the perturbative DIS structure functions. In comparison, the large- x x photon uncertainty comes from various low-energy, nonperturbative contributions, including variations of the inelastic structure functions in the resonance and continuum regions, higher-twist and target-mass corrections, and elastic electromagnetic form factors of the proton. We take the production of doubly-charged Higgs pairs, (H^{++}H^{--}) ( H + + H − − ) , as an example of scenarios beyond the Standard Model to illustrate the phenomenological implications of these photon PDFs at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

DOE OSTI2024

Nuclear Dependence of Charged Current Inclusive $\bar{\nu}_\mu -$A Cross Section Measurements at 6 GeV at MINERvA

The main advantage of the MINERvA detector is that it uses differently sized nuclear targets, in the same neutrino beam, allowing for a direct study of nuclear medium effects, with reduced flux and detector uncertainties. This thesis presents the first high statistics charged current antineutrino-nucleus inclusive two-dimensional cross section measurements using an antineutrino beam with $\langle E_{\bar\nu}\rangle\sim$ 6 GeV, on iron, lead, carbon and hydrocarbon targets, as a function of bjorken $x$ and four-momentum transfer squared $Q^2$, invariant mass $W$ and $Q^2$ and muon longitudinal and transverse momenta, $p_t$ and $p_z$. These nuclei of different sizes were utilized to measure the nuclear dependence of the cross section. The antineutrino interactions that were utilized for the analysis presented in this thesis, were having an antimuon in the final state with energy 2-20GeV and with an antimuon scattering angle less than $17^o$. This analysis represents one of the largest antineutrino datasets analyzed to date in the energy region of $\langle E_{\bar\nu} \rangle \sim$ 6 GeV. This direct evidence for the nuclear medium effects, can help benchmarking models, which can be used by future neutrino experiments like DUNE, by extrapolating these results to Argon target.

Gaur, Prameet Kumar↗

DOE OSTI2025

Truncating 2D Framework Materials Down to a Single Pore: Synthetic Approaches and Opportunities

Here, in this Accounts article, we summarize our recent work on truncating conjugated two-dimensional framework materials down to a single pore, or a single macrocycle. Conjugated 2D architectures have emerged as one of the most synthetically adaptable motifs for coupling semiconductivity and porosity in metal–organic frameworks (MOFs) and covalent organic frameworks (COFs). However, despite their prevalence, 2D architectures have several limitations. In particular, the strong interlayer π–π stacking can limit both processability and the accessibility of internal active sites. We have found that simple macrocycles preserve key aspects of 2D framework structure and function, including porosity and out-of-plane electrical conductivity, while providing improved processability, surface tunability, and mass transport properties. In this article, we first describe our synthetic approach and general design considerations. Specifically, we show how ditopic analogues of the tritopic ligands commonly found in the synthesis of 2D MOFs and COFs can be used to achieve a diverse library of conjugated macrocycles that resemble fragments of semiconducting frameworks in both form and function. The length of the peripheral side chains, the size of the aromatic core, and the solubility of intermediates are all key variables in favoring selective macrocycle formation over undesired linear polymers and oligomers. Next, we highlight the unique advantages that macrocycles provide, including improved processability, atomically precise surface tunability, and greater active site accessibility. In particular, the identity of the peripheral side chains dramatically impacts both solubility and colloidal stability as well as crystal size and morphology. We further show how the solution processability and nanoscale dimensions of macrocycles can simplify electronic device fabrication and improve electrochemical performance. Finally, we end with a forward-looking discussion on how macrocycles offer a unique bridge between conjugated molecules and extended frameworks, enabling new application areas and fundamental science.

charge transport↗

DOE OSTI2024

An implicit solution for Asay foil trajectories generated by separable, sustained-production ejecta source models

We present a simple implicit solution for the time-dependent trajectory of a thin Asay foil ejecta diagnostic for the general case where the impinging ejecta cloud is generated by a source function characterized by an arbitrary (sustained) time dependence and a time-independent (stationary) particle velocity distribution. In the limit that the source function time dependence becomes a delta function, this solution—which is amenable to rapid numerical calculations of arbitrary accuracy—exactly recovers a previously published solution for the special case of instantaneous ejecta production. We also derive simple expressions for the free-surface arrival (catch-up) time as well as the true ejecta areal mass accumulation on the accelerating foil and place bounds on the level of error incurred when applying instant-production mass solutions to a sustained-production trajectory. We demonstrate these solutions with example calculations for hypothetical source functions spanning a wide range of ejecta production durations, velocity distributions, and temporal behaviors. These calculations demonstrate how the foil trajectory is often insensitive to the temporal dependence of the source function, instead being dominated by the velocity distribution. We quantify this insensitivity using a “compatibility score” metric. Under certain conditions, one may capitalize upon this insensitivity to obtain a good approximation of the second integral of the velocity distribution from the observed foil trajectory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Loop-by-loop differential equations for dual (elliptic) Feynman integrals

We present a loop-by-loop method for computing the differential equations of Feynman integrals using the recently developed dual form formalism. We give explicit prescriptions for the loop-by-loop fibration of multi-loop dual forms. Then, we test our formalism on a simple, but non-trivial, example: the two-loop three-mass elliptic sunrise family of integrals. We obtain an ε-form differential equation within the correct function space in a sequence of relatively simple algebraic steps. In particular, none of these steps relies on the analysis of q-series. Then, we discuss interesting properties satisfied by our dual basis as well as its simple relation to the known ε-form basis of Feynman integrands. The underlying K3-geometry of the three-loop four-mass sunrise integral is also discussed. Finally, we speculate on how to construct a “good” loop-by-loop basis at three-loop.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Precision measurement of the electron energy-loss function in tritium and deuterium gas for the KATRIN experiment

Abstract The KATRIN experiment is designed for a direct and model-independent determination of the effective electron anti-neutrino mass via a high-precision measurement of the tritium $$\upbeta $$ β -decay endpoint region with a sensitivity on $$m_\nu $$ m ν of 0.2 $$\hbox {eV}/\hbox {c}^2$$ eV / c 2 (90% CL). For this purpose, the $$\upbeta $$ β -electrons from a high-luminosity windowless gaseous tritium source traversing an electrostatic retarding spectrometer are counted to obtain an integral spectrum around the endpoint energy of 18.6 keV. A dominant systematic effect of the response of the experimental setup is the energy loss of $$\upbeta $$ β -electrons from elastic and inelastic scattering off tritium molecules within the source. We determined the energy-loss function in-situ with a pulsed angular-selective and monoenergetic photoelectron source at various tritium-source densities. The data was recorded in integral and differential modes; the latter was achieved by using a novel time-of-flight technique. We developed a semi-empirical parametrization for the energy-loss function for the scattering of 18.6-keV electrons from hydrogen isotopologs. This model was fit to measurement data with a 95% $$\hbox {T}_2$$ T 2 gas mixture at 30 K, as used in the first KATRIN neutrino-mass analyses, as well as a $$\hbox {D}_2$$ D 2 gas mixture of 96% purity used in KATRIN commissioning runs. The achieved precision on the energy-loss function has abated the corresponding uncertainty of $$\sigma (m_\nu ^2)< {{10}^{-2}}{\hbox {eV}^{2}}$$ σ ( m ν 2 ) < 10 - 2 eV 2 [1] in the KATRIN neutrino-mass measurement to a subdominant level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Development of New g-function Data for Simulating a Novel Shallow Bore Ground Heat Exchanger

Geothermal heat pump (GHP), which is also referred as ground source heat pump, is the most energy-efficient technology for space heating and cooling. However, the application of GHP is hindered by its high initial cost, of which approximately 30% is for the ground heat exchanger (GHE). In recent years, researchers have developed different types of shallow bore ground heat exchangers (SBGHE) intending to reduce the cost of drilling. The ability to predict the thermal response of an SBGHE is critical for sizing SBGHE. While g-functions have been commonly used for predicting the thermal response of conventional vertical bore ground heat exchangers (VBGHE), they cannot be directly used for predicting the performance of SBGHE because they did not account for the impacts of the seasonal variation of the soil temperature along with the depth of an SBGHE. In addition, an SBGHE has a larger thermal mass within the borehole than the VBGHE due to the larger borehole diameter of SBGHE. This study develops new g-function data pairs for a new design of SBGHE, which is named Underground Thermal Battery (UTB). Impacts of the seasonal variation of soil temperature on the thermal response of UTB were accounted for by superposing a time-dependent soil temperature onto the g-functions calculated with a numerical model that assumes constant undisturbed soil temperature. The TOUGH program was used to predict the thermal response of several configurations involving multiple UTBs. The results indicate that the proposed methodology is appropriate to generate g-functions for the UTB, and the g-function value of UTB is much lower than that of VBGHE in the time range of 15 min to 1 year due to the large thermal mass and convection heat transfer within the UTB.

Shi, Liang↗

DOE OSTI2024

Measurements of the $Z$ 0 /$γ$* cross section and transverse single spin asymmetry in 510 GeV $p$ + $p$ collisions

The differential cross section for Z 0 production, measured as a function of the boson’s transverse momentum (p T ), provides important constraints on the evolution of the transverse momentum dependent parton distribution functions (TMDs). The transverse single spin asymmetry (TSSA) of the Z 0 is sensitive to one of the polarized TMDs, the Sivers function, which is predicted to have the opposite sign in p + p → W/Z + X from that which enters in semi-inclusive deep inelastic scattering. In this Letter, the STAR Collaboration reports the first measurement of the Z 0 /γ* differential cross section as a function of its p T in p+p collisions at a center-of-mass energy of 510 GeV, together with the Z 0 /γ* total cross section. We also report the measurement of Z 0 /γ* TSSA in transversely polarized p+p collisions at 510 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Strong-field induced fragmentation and isomerization of toluene probed by ultrafast femtosecond electron diffraction and mass spectrometry

We investigate the fragmentation and isomerization of toluene molecules induced by strong-field ionization with a femtosecond near-infrared laser pulse. Momentum-resolved coincidence time-of-flight ion mass spectrometry is used to determine the relative yield of different ionic products and fragmentation channels as a function of laser intensity. Ultrafast electron diffraction is used to capture the structure of the ions formed on a picosecond time scale by comparing the diffraction signal with theoretical predictions. Through the combination of the two measurements and theory, we are able to determine the main fragmentation channels and to distinguish between ions with identical mass but different structures. In addition, our diffraction measurements show that the independent atom model, which is widely used to analyze electron diffraction patterns, is not a good approximation for diffraction from ions. Here, we show that the diffraction data is in very good agreement with ab initio scattering calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗