Constraining storm-scale forecasts of deep convective initiation with surface weather observations

Thesis (Ph.D.)--University of Washington, 2016-06

Effect of addition of maltodextrin on drying kinetics and stickiness of sugar and acid-rich foods during convective drying: experiments and modelling

The effect of addition of maltodextrin on drying kinetics of drops containing fructose, glucose, sucrose and citric acid individually and in mixtures was studied experimentally using single drop drying experiments and numerically by solving appropriate mass and heat transfer equations. The numerical predictions agreed with the experimental moisture and temperature histories within 5-6% average relative (absolute) errors and average differences of +/- 1degreesC, respectively. The stickiness of these drops was determined using the glass transition temperature (T-g) of the drops' surface layer as an indicator. The experimental stickiness histories followed the model predictions with reasonable accuracy. A safe drying (non-sticky) regime in a spray drying environment has been proposed, and used to estimate the optimum amount of addition of maltodextrin for successful spray drying of 120 micron diameter droplets of fruit juices. (C) 2003 Elsevier Ltd. All rights reserved.

Theoretical investigation of convective instability in inclined and fluid-saturated three-dimensional fault zones.

The convective instability of pore-fluid flow in inclined and fluid-saturated three-dimensional fault zones has been theoretically investigated in this paper. Due to the consideration of the inclined three-dimensional fault zone with any values of the inclined angle, it is impossible to use the conventional linear stability analysis method for deriving the critical condition (i.e., the critical Rayleigh number) which can be used to investigate the convective instability of the pore-fluid flow in an inclined three-dimensional fault zone system. To overcome this mathematical difficulty, a combination of the variable separation method and the integration elimination method has been used to derive the characteristic equation, which depends on the Rayleigh number and the inclined angle of the inclined three-dimensional fault zone. Using this characteristic equation, the critical Rayleigh number of the system can be numerically found as a function of the inclined angle of the three-dimensional fault zone. For a vertically oriented three-dimensional fault zone system, the critical Rayleigh number of the system can be explicitly derived from the characteristic equation. Comparison of the resulting critical Rayleigh number of the system with that previously derived in a vertically oriented three-dimensional fault zone has demonstrated that the characteristic equation of the Rayleigh number is correct and useful for investigating the convective instability of pore-fluid flow in the inclined three-dimensional fault zone system. The related numerical results from this investigation have indicated that: (1) the convective pore-fluid flow may take place in the inclined three-dimensional fault zone; (2) if the height of the fault zone is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone stabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; (3) if the thickness of the stratum is used as the characteristic length of the system, a decrease in the inclined angle of the inclined fault zone destabilizes the three-dimensional fundamental convective flow in the inclined three-dimensional fault zone system; and that (4) the shape of the inclined three-dimensional fault zone may affect the convective instability of pore-fluid flow in the system. (C) 2004 Published by Elsevier B.V.

Comparisons of diffusive and viscous contributions to transport coefficients of light gases in single-walled carbon nanotubes

We examine here the relative importance of different contributions to transport of light gases in single walled carbon nanotubes, using methane and hydrogen as examples. Transport coefficients at 298 K are determined using molecular dynamics simulation with atomistic models of the nanotube wall, from which the diffusive and viscous contributions are resolved using a recent approach that provides an explicit expression for the latter. We also exploit an exact theory for the transport of Lennard-Jones fluids at low density considering diffuse reflection at the tube wall, thereby permitting the estimation of Maxwell coefficients for the wall reflection. It is found that reflection from the carbon nanotube wall is nearly specular, as a result of which slip flow dominates, and the viscous contribution is small in comparison, even for a tube as large as 8.1 nm in diameter. The reflection coefficient for hydrogen is 3-6 times as large as that for methane in tubes of 1.36 nm diameter, indicating less specular reflection for hydrogen and greater sensitivity to atomic detail of the surface. This reconciles results showing that transport coefficients for hydrogen and methane, obtained in simulation, are comparable in tubes of this size. With increase in adsorbate density, the reflection coefficient increases, suggesting that adsorbate interactions near the wall serve to roughen the local potential energy landscape perceived by fluid molecules.

Assessment of P availability in heavily fertilized soils using the diffusive gradient in thin films (DGT) technique

Phosphorus-availability tests typically provide an indication of quantity of P available (Colwell bicarbonate-extractable P), or of the intensity of supply (0.01 M CaCl2-extractable P). The soil's capacity to buffer P is more difficult to assess, and is generally estimated using a P-adsorption curve. The diffusive gradient in thin films (DGT) approach may provide a simpler means of assessing a soil's ability to maintain soil solution P. Optimal extraction conditions were found to be 24 h exposure of DGT samplers to saturated soil. The DGT approach was evaluated on a range of 24 soils, some of which had high Colwell- (>100 mu g g(-1)) and Bray 1- (>30 mu g g(-1)) extractable P content, but showed a tomato (Lycopersicon esculentum Mill.) yield response to the addition of P fertilizer. The DGT approach provided an excellent separation of soils on which tomato showed a yield response, from those where fertilizer P did not increase dry-matter yield. Phosphorus accumulation was strongly correlated with soil solution P concentration and anion exchange resin-extractable P, but showed poor correlation with Colwell- or Bray 1-extractable P. The DGT P accumulation rate of 3.62 x 10(-7) to 4.79 x 10(-5) mol s(-1) m(-3) for the soils tested was comparable to the uptake rate of roots of tomato plants that were adequately supplied with P (2.25 x 10(-5) mol s(-1) m(-3)).

Experimental and analytical study of the effect of contact angle on liquid convective heat transfer in microchannels

In this paper we examine the effect of contact angle (or surface wettability) on the convective heat transfer coefficient in microchannels. Slip flow, where the fluid velocity at the wall is non-zero, is most likely to occur in microchannels due to its dependence on shear rate or wall shear stress. We show analytically that for a constant pressure drop, the presence of slip increases the Nusselt number. In a microchannel heat exchanger we modified the surface wettability from a contact angle of 20 degrees-120 degrees using thin film coating technology. Apparent slip flow is implied from pressure and flow rate measurements with a departure from classical laminar friction coefficients above a critical shear rate of approximately 10,000 s(-1). The magnitude of this departure is dependant on the contact angle with higher contact angles surfaces exhibiting larger pressure drop decreases. Similarly, the non-dimensional heat flux is found to decrease relative to laminar non-slip theory, and this decrease is also a function of the contact angle. Depending on the contact angle and the wall shear rate, variations in the heat transfer rate exceeding 10% can be expected. Thus the contact angle is an important consideration in the design of micro, and even more so, nano heat exchangers. (c) 2006 Elsevier Ltd. All rights reserved.

CFD modelling of the fast pyrolysis of an in-flight cellulosic particle subjected to convective heat transfer

The pyrolysis of a freely moving cellulosic particle inside a 41.7mgs -1 source continuously fed fluid bed reactor subjected to convective heat transfer is modelled. The Lagrangian approach is adopted for the particle tracking inside the reactor, while the flow of the inert gas is treated with the standard Eulerian method for gases. The model incorporates the thermal degradation of cellulose to char with simultaneous evolution of gases and vapours from discrete cellulosic particles. The reaction kinetics is represented according to the Broido–Shafizadeh scheme. The convective heat transfer to the surface of the particle is solved by two means, namely the Ranz–Marshall correlation and the limit case of infinitely fast external heat transfer rates. The results from both approaches are compared and discussed. The effect of the different heat transfer rates on the discrete phase trajectory is also considered.

Mechanisms of convective drying in thick beds of solids

The literature on heat and mass transfer mechanisms in the convective drying of thick beds of solids has been critically reviewed. Related mathematical models of heat transfer are also considered. Experimental and theoretical studies were made of the temperature distribution within beds, and of drying rates, with various materials undergoing convective drying. The experimental work covered thick beds of hygroscopic and non-hygroscopic materials (glass beads of different diameters, polystyrene pellets, activated alumina and wood powder) at air temperatures of 54°C to 84°C. Tests were carried out in a laboratory drying apparatus comprising a wind tunnel through which the air, of controlled temperature and humidity, was passed over a sample suspended from a balance. Thermocouples were inserted at different depths within the sample bed. The temperature distribution profiles for both hygroscopic and non-hygroscopic beds exhibited a clear difference between the temperatures at the surface and bottom during the constant rate period. An effective method was introduced for predicting the critical moisture content. During the falling rate the profiles showed the existence of a receding evaporation plane; this divided the system into a hotter dry zone in the upper section and a wet zone near the bottom. A graphical procedure was established to predict accurately the position of the receding evaporation front at any time. A new mathematical model, based on the receding evaporation front phenomenon, was proposed to predict temperature distributions throughout a bed during drying. Good agreement was obtained when the model was validated by comparing its predictions with experimental data. The model was also able to predict the duration of each drying stage. In experiments using sample trays of different diameters, the drying rate was found to increase with a decrease in the effective length of the bed surface. During the constant rate period with trays of a small effective length, i.e. less than 0.08 m, an 'inversion' in temperature distribution occurred in the bed; the bottom temperature increased and became greater than that of the surface. Experimental measurements were verified in several ways to ensure this phenomenon was real. Theoretical explanations are given for both the effective length and temperature inversion phenomena.

Transition in convective flows heated internally

The stability of internally heated convective flows in a vertical channel under the influence of a pressure gradient and in the limit of small Prandtl number is examined numerically. In each of the cases studied the basic flow, which can have two inflection points, loses stability at the critical point identified by the corresponding linear analysis to two-dimensional states in a Hopf bifurcation. These marginal points determine the linear stability curve that identifies the minimum Grashof number (based on the strength of the homogeneous heat source), at which the two-dimensional periodic flow can bifurcate. The range of stability of the finite amplitude secondary flow is determined by its (linear) stability against three-dimensional infinitesimal disturbances. By first examining the behavior of the eigenvalues as functions of the Floquet parameters in the streamwise and spanwise directions we show that the secondary flow loses stability also in a Hopf bifurcation as the Grashof number increases, indicating that the tertiary flow is quasi-periodic. Secondly the Eckhaus marginal stability curve, that bounds the domain of stable transverse vortices towards smaller and larger wavenumbers, but does not cause a transition as the Grashof number increases, is also given for the cases studied in this work.

Charge effects on hindrance factors for diffusion and convection of solute in pores I

The transport of a spherical solute through a long circular cylindrical pore filled with an electrolyte solution is studied numerically, in the presence of constant surface charge on the solute and the pore wall. Fluid dynamic analyses were carried out to calculate the flow field around the solute in the pore to evaluate the drag coefficients exerted on the solute. Electrical potentials around the solute in the electrolyte solution were computed based on a mean-field theory to provide the interaction energy between the charged solute and the pore wall. Combining the results of the fluid dynamic and electrostatic analyses, we estimated the rate of the diffusive and convective transport of the solute across the pore. Although the present estimates of the drag coefficients on the solute suggest more than 10% difference from existing studies, depending on the radius ratio of the solute relative to the pore and the radial position of the solute center in the pore, this difference leads to a minor effect on the hindrance factors. It was found that even at rather large ion concentrations, the repulsive electrostatic interaction between the charged solute and the pore wall of like charge could significantly reduce the transport rate of the solute.

Charge effects on the hindered transport of macromolecules across the endothelial surface glycocalyx layer

A fluid mechanical and electrostatic model for the transport of solute molecules across the vascular endothelial surface glycocalyx layer (EGL) was developed to study the charge effect on the diffusive and convective transport of the solutes. The solute was assumed to be a spherical particle with a constant surface charge density, and the EGL was represented as an array of periodically arranged circular cylinders of like charge, with a constant surface charge density. By combining the fluid mechanical analyses for the flow around a solute suspended in an electrolyte solution and the electrostatic analyses for the free energy of the interaction between the solute and cylinders based on a mean field theory, we estimated the transport coefficients of the solute across the EGL. Both of diffusive and convective transports are reduced compared to those for an uncharged system, due to the stronger exclusion of the solute that results from the repulsive electrostatic interaction. The model prediction for the reflection coefficient for serum albumin agreed well with experimental observations if the charge density in the EGL is ranged from approximately -10 to -30 mEq/l.