A numerical model based on fundamental engineering principles was developed and validated to establish a relationship between process parameters and the final properties of woodbased composite boards. The model simulates the mat formation, then compresses the reconstituted mat to its final thickness in a virtual press. The number of interacting variables during the hot-compression process is prohibitively large to assess a wide variety of data by experimental means. Therefore, the main advantage of the model based approach that the effect of the hot-compression parameters on the final properties of wood-based composite boards can be monitored without extensive experimentation. The mat formation part of the model is based on the Monte Carlo simulation technique to reproduce the spatial structure of the mat. The dimensions and the density of each flake are considered as random variables in the model, which follow certain probability density distributions. The parameters of these distributions are derived from data collected on industrial flakes by using an image analysis technique. The model can simulate the structure of a threelayer oriented strandboard (OSB) mat as well as the structure of random fiber networks. A grid is superimposed on the simulated mat and the number of flakes, the thickness, and the density of the mat at each grid point are computed. Additionally, the model predicts the change in several void volume fractions within the mat and the contact area between the flakes during consolidation. The void volume fractions are directly related to the physical properties of the mat, such as thermal conductivity, diffusivity, and permeability, and the contact area is an indicator of the effectively bonded area within the mat. The heat and mass transfer part of the model predicts the change of air content, moisture content, and temperature at designated mesh points in the cross section of the mat during the hotcompression. The water content is subdivided into vapor and bound water components. The free water component is not considered in the model due to the low (typically 6-7 %) initial moisture content of the flakes. The gas phase (air and vapor) moves by bulk flow and diffusion, while the bound water only moves by diffusion across the mat. The heat flow occurs by conduction and convection. The spatial derivatives of the resulting coupled partial differential equations are discretized by finite differences. The resulting ordinary differential equation in time is solved by a differential-algebraic system solver (DDASSL). The internal environment within the mat can be predicted among different initial and boundary conditions by this part of the hot-compression model. In the next phase of the research, the viscoelastic (time, temperature, and moisture dependent) response of the flakes was modeled using the time-temperature-moisture superposition principle of polymers. A master curve was created from data available in the literature, which describes the changing relaxation modulus of the flakes as a function of moisture and temperature at different locations in the mat. Then the flake mat was compressed in a virtual press. The stress-strain response is highly nonlinear due to the cellular structure of the mat. Hooke's Law was modified with a nonlinear strain function to account for the behavior of the flake mat in transverse compression. This part of the model gives insight into the vertical density profile formation through the thickness of the mat. Laboratory boards were produced to validate the model. A split-plot experimental design, with three different initial mat moisture contents (5, 8.5, 12 %), three final densities (609, 641, 673 kg êm3 or 38, 40, 42 lb ê ft3), two press platen temperatures (150, 200 °C), and three different press closing times (40, 60, 80 s) was applied to investigate the effect of production parameters on the internal mat conditions and the formation of the vertical density profile. The temperature and gas pressure at six locations in the mat, and the resultant density profiles of the laboratory boards, were measured. Adequate agreement was found between the model predicted and the experimentally measured temperature, pressure, and vertical density profiles. The complete model uses pressing parameters (press platen temperature, press schedule) and mat properties (flake dimensions and orientation, density distribution, initial moisture content and temperature) to predict the resulting internal conditions and vertical density profile formation within the compressed board. The density profile is related to all the relevant mechanical properties (bending strength, modulus of elasticity, internal bond strength) of the final board. The model can assist in the optimization of the parameters for hot-pressing woodbased composites and improve the performance of the final panel. / Ph. D.

The biogeochemical processes responsible for production and emission of ammonia from stored liquid dairy manure are governed by environmental factors (e.g. manure temperature, moisture) and manure characteristics (e.g. total ammoniacal nitrogen concentration, pH). These environmental factors and manure characteristics vary spatially as a result of spatially heterogeneous physical, chemical, and biological properties of manure. Existing process-based models used for estimating ammonia emission consider stored manure as a homogeneous system and do not consider these spatial variations leading to inaccurate estimations. In this study, a one-dimensional compartmental biogeochemical model was developed to (i) estimate spatial variation of temperature and substrate concentration (ii) estimate spatial variations and rates of biogeochemical processes, and (iii) estimate production and emission of ammonia from stored scraped liquid dairy manure. A one-dimension compartmentalized modeling approach was used whereby manure storage is partitioned into several sections in vertical domain assuming that the conditions are spatially uniform within the horizontal domain. Spatial variation of temperature and substrate concentration were estimated using established principles of heat and mass transfer. Pertinent biogeochemical processes were assigned to each compartment to estimate the production and emission of ammonia. Model performance was conducted using experimental data obtained from National Air Emissions Monitoring Study conducted by the United States Environmental Protection Agency. A sensitivity analysis was performed and air temperature, manure pH, wind speed, and manure total ammoniacal nitrogen concentration were identified as the most sensitive model inputs. The model was used to estimate ammonia emission from a liquid dairy manure storage of a dairy farm located in Rockingham and Franklin counties in Virginia. Ammonia emission was estimated under different management and weather scenarios: two different manure storage periods from November to April and May to October using historical weather data of the two counties. Results suggest greater ammonia emissions and manure nitrogen loss for the manure storage period in warm season from May to October compared to the storage period in cold season from November to April. / Ph. D.

<p>Co-axial dripping is one of the many ways to make drops with a core-shell structure for encapsulated materials. However, in systems where the capsule components are not density matched or surfactants are not used, the shell will eventually thin and break if not solidified in time. If the shell fails before solidifying, the core will leak out and result in a non-functional capsule. This study assumes that all capsules will fail once the core has reached 80% eccentricity, meaning a shell region has thinned to 20% of its original thickness (~70 µm). In reality, rupture of the shell depends more on stochastic defects and disturbances, but locally decreasing the shell thickness will increase the probability of capsule rupture. With this assumption, the survival time of a core-shell drop is inversely proportional to the relative velocity of the inner drop, where the greater this relative velocity, the faster the shell phase will thin. Stoke's law is generally used to approximate the speed of a sphere in a fluid. However, this study demonstrates that Stoke's law is insufficient for predicting the inner drop's motion for a compound drop. This is due to internal flows that develop within all fluid drops because of shear forces on the drop’s external face during freefall. For core-shell drops, prior studies report how the inner drop velocity can change in magnitude and direction as a function of its eccentricity, meaning its position within the outer drop. Since previous studies did not analyze this core-shell drop relationship with a 50 vol% core and a high viscosity shell, a model was built in COMSOL Multiphysics to understand how the claims from literature would apply to a previous encapsulation study (Betancourt, 2021). The model was also put through a series of validation tests that confirmed the model’s ability to accurately represent the speed and direction of inner drop motion. The final model configuration was then used to identify the transition point between buoyancy-driven and internal flow-driven failure modes observed during the production of core-shell drops in a previous encapsulation study for phase change materials (Betancourt, 2021). The model results showed how the estimated inner drop velocity was significantly reduced once accounting for the internal flows within the shell phase of a compound drop. While this study does help characterize the motion of an inner drop and could be used to find a material system with a favorable velocity profile, it is still recommended to use an in-air curing system to produce concentric capsules. Achieving a concentric capsule would still require this co-axial dripping setup to be modified significantly. </p> <p>Betancourt-Jimenez, D., Wells, B., Youngblood, J. P., & Martinez, C. J. (2021). Encapsulation of biobased fatty acid amides for phase change material applications. <em>Journal of Renewable and Sustainable Energy</em>, <em>13</em>(6), 064101. https://doi.org/10.1063/5.0072105</p>

<p>  </p> <p>Heat pipes and vapor chambers are passive two-phase heat transport devices that are used for thermal management in electronics. The passive operation of a heat pipe is facilitated by capillary wicking of the working fluid through a porous wick, and thus is subject to an operational limit in terms of the maximum pressure head that the wick can provide. This operational limit, often termed as the capillary limit, is the maximum heat input at which the pressure drop in the wick is balanced by the maximum capillary pressure head; operating a heat pipe or a vapor chamber above the capillary limit at steady-state leads to dryout. It thus becomes important to predict the performance of heat pipes and vapor chambers and explore the parametric design space to provide guidelines for minimized thermal resistance while satisfying this capillary limit. An increasingly critical aspect is to predict the transient thermal response of vapor chambers. Moreover, heat pipes and vapor chambers are extensively being used in electronic systems where the power input is dictated by the end-user activity and is expected to even exceed the capillary limit for brief time intervals. Thus, it is imperative to understand the behavior of heat pipes and vapor chambers when operated at steady and transient heat loads above the capillary limit as dryout occurs. However, review of the literature on heat pipe performance characterization reveals that the regime of dryout operation has been virtually unexplored, and thus this thesis aims to fill this critical gap in understanding.</p> <p>The design for minimized thermal resistance of a vapor chamber or a heat pipe is guided by the relative contribution of thermal resistance due to conduction across the evaporator wick and the saturation temperature gradient in the vapor core. In the limit of very thin form factors, the contribution from the vapor core thermal resistance dominates the overall thermal resistance of the vapor chamber; recent work has focused on working fluid selection to minimize overall thermal resistance in this limit. However, the wick thermal resistance becomes increasingly significant as its thickness increases to support higher heat inputs while avoiding the capillary limit. A thermal resistance network model is thus utilized to investigate the importance of simultaneously considering the contributions of the wick and vapor core thermal resistances. A generalized approach is proposed for vapor chamber design which allows <em>simultaneous</em> selection of the working fluid and wick that provides minimum overall thermal resistance for a given geometry and operating condition. While the thermal resistance network model provides a convenient method for exploring the design space, it cannot be used to predict 3-D temperature fields in the vapor chamber. Moreover, such thermal resistance network models cannot predict transient performance and temperature evolution for a vapor chamber. Therefore, an easy-to-use approach is proposed for mapping of vapor chamber transport to the heat diffusion equation using a set of appropriately defined effective anisotropic thermophysical properties, thus allowing simulation of vapor chamber as a sold conduction block. This effective anisotropic properties approach is validated against a time-stepping analytical model and is shown to have good match for both spatial and transient temperature predictions.</p> <p>Moving the focus from steady-state and transient operation of vapor chambers, a comprehensive characterization of heat pipe operation above capillary limit is performed. Different user needs and device workloads can lead to highly transient heat loads which could exceed the notional capillary limit for brief time intervals. Experiments are performed to characterize the transient thermal response of a heat pipe subjected to heat input pulses of varying duration that exceed the capillary limit. Transient dryout events due to a wick pressure drop exceeding the maximum available capillary pressure can be detected from an analysis of the measured temperature signatures. It is discovered that under such transient heating conditions, a heat pipe can sustain heat loads higher than the steady-state capillary limit for brief periods of time without experiencing dryout. If the heating pulse is sufficiently long as to induce transient dryout, the heat pipe may experience an elevated steady-state temperature even after the heat load is reduced back to a level lower than the capillary limit. The steady-state heat load must then be reduced to a level much below the capillary limit to fully recover the original thermal resistance of the heat pipe. The recovery process of heat pipes is further investigated, and a mechanism is proposed for the thermal hysteresis observed in heat pipe performance after dryout. A model for <em>steady-state</em> heat pipe transport is developed based on the proposed mechanism to predict the parametric trends of thermal resistance following recovery from dryout-induced thermal hysteresis, and the model is mechanistically validated against experiments. The experimental characterization of the recovery process demonstrates the existence of a maximum hysteresis curve, which serves as the worst-case scenario for thermal hysteresis in heat pipe after dryout. Based on the learnings from the experimental characterization, a new procedure is introduced to experimentally characterize the steady-state dryout performance of a heat pipe.</p> <p>To recover the heat pipe performance under steady-state, it has been shown that the heat input needs to be lowered down or <em>throttled</em> significantly below the capillary limit. However, due to the highly transient nature of power dissipation from electronic devices, it becomes imperative to characterize heat pipe recovery from dryout under transient operations. Hence, power-throttling assisted recovery of heat pipe from dryout has been characterized under transient conditions. A minimum throttling time interval, defined as time-to-rewet, is identified to eliminate dryout induced thermal hysteresis using power throttling. Dependence of time-to-rewet on throttling power is explored, and guidelines are presented to advise the throttling need and choice of throttling power under transient conditions. </p> <p>The experimental characterization of heat pipe operation at pulse loads above the capillary limit and power throttling following the pulse load helped define the dryout and recovery performance of a heat pipe. Next, a physics-based model is developed to predict the heat pipe <em>transient</em> thermal response under dryout-inducing pulse load and power throttling assisted recovery. This novel model considers wick as a partially saturated media with spatially and temporally varying liquid saturation, and accounts for the effect of wick partial saturation in heat pipe transport. The model prediction are validated against experiments with commercial heat pipe samples, and it is shown that the model can accurately predict dryout and recovery characteristics, namely time-to-dryout, time-to-rewet, and dryout-induced thermal hysteresis, for heat pipes with a range of wick types, heat pipe lengths and pulse loads above the capillary limit. </p> <p>The work discussed in this thesis opens certain questions that are expected to guide further research in this area. First, the thermal hysteresis mechanism proposed could be further validated with direct visualization of the liquid in a vapor chamber. To achieve this, X-ray microscopy is proposed as a viable option for the imaging <em>in situ</em> wetting dynamics in a vapor chamber. Second, the model developed to predict the dryout and recovery characteristics of the heat pipe can be used to design heat pipe with improved performance under pulse loads and power throttling. Third, novel wick designs can be explored that utilize the understanding developed of governing mechanisms for recovery from dryout, and can eliminate thermal hysteresis at powers closer to capillary limit. Fourth, the modeling approach can be extended to predict dryout and recovery trends in vapor chamber since the heat transfer pathways in a vapor chamber are different than those of a heat pipe. Fifth, and lastly it was observed several times during experiments that some of the heat pipe samples would exhibit complete dryout (sudden catastrophic rise in temperature and thermal resistance at the point of dryout) whereas other samples would exhibit partial dryout (noticeable but small increase in thermal resistance at dryout) at operating powers just above the capillary limit. Exploring and explaining the cause of complete dryout, in particular, would be an extremely valuable contribution to the heat pipe research. </p> <p>The work discussed in this thesis has led to the comprehensive development of a functional and mechanistic understanding of heat pipe operation above the notional capillary limit. The experimental procedures developed in this work are utilized to characterize a heat pipe performance under dryout and recovery. The models based on the mechanistic understanding developed from experimental characterization of dryout and recovery operation of a heat pipe have been experimentally validated and are useful for predicting heat pipe performance under dryout-inducing pulse loads and power-throttling.   </p>

<p>  </p> <p>Cavitation is one of the most common causes of failures in axial piston machines. Due to the detrimental effects that cavitation has on unit performance, it is of important consideration both in the design of new units and in defining the operational limits of existing market products. The work in this thesis aimed to contribute to the current knowledge in both areas, with a focus on design considerations with respect to cavitation scalability, and on operating conditions by measuring cavitation severity under separate and combined inciting parameters. Though the application of unit scaling is common in industry for the design of pump families, there have been no comprehensive attempts to quantify whether cavitation in fluid power units may be adequately accounted for in published scaling laws. In this thesis, the scalability of cavitation phenomena was examined through a CFD scaling study performed using a modified version of the Full Cavitation Model.  Results indicate that linear scaling is consistent in maintaining volumetric efficiency performance within 1% across scaled units up to eight times larger or smaller than the baseline. However, the gas and vapor volume distributions vary significantly between scaled units, due largely to the linear non-scalability of fluid inertia and turbulent factors. Physical exchange between phases within a working fluid was shown to be time-dependent, such that the scaled-down unit exhibits bubble collapse rates up to 30% and 150% greater than the baseline and scaled-up units, respectfully. Considering these effects, the presented work demonstrates a potential for increased cavitation damage area when downscaling a unit and reduced risk in upscaling, despite the scaling law being a reliable indicator for volumetric efficiency. </p> <p>To define a more complete study of cavitation under a variety of operating conditions and inciting parameters, this a new experimental procedure and testing circuit was proposed with focus on repeatability by controlled pressure drops and preliminary quantification of inlet fluid quality. By measuring cavitation conditions under pressure starvation, incomplete filling, and combinations thereof, the direct effect of different inception methods on unit performance was shown to be readily identifiable. Through visualization of the inlet flow, reduction in inlet pressure levels was correlated to fluid cloudiness levels and bubble size, with transparency loss at 0.0 bar_g and transition from bubbly to plug flow at -0.4 bar_g. Incomplete filling-induced cavitation was also shown to be detectable by inlet flow conditions, with a distinct change in bubble coalescence rate when operating under shaft speeds greater than or equal to fill speed for a given inlet pressure. </p>

<p>The present study investigates the capability of computational fluid dynamics (CFD) extensively to predict hydrodynamics and heat transfer characteristics of FC-72 flow boiling in a 2.5-mm ´ 5.0-mm rectangular channel and experimentally explores system instabilities: <em>density wave oscillation</em> (DWO), <em>pressure drop oscillation</em> (PDO) and <em>parallel channel instability</em> (PCI) in a micro-channel heat sink containing 38 parallel channels having a hydraulic diameter of 316-μm. </p> <p>The computational method performs transient analysis to model the entire flow field and bubble behavior for subcooled flow boiling in a rectangular channel heated on two opposite walls at high heat flux conditions of about 40% – 80% of <em>critical heat flux</em> (CHF).  The 3D CFD solver is constructed in ANSYS Fluent in which the <em>volume of fluid</em> (VOF) model is combined with a <em>shear stress transport</em> (SST) <em>k</em>-<em>ω</em> turbulent model, a surface tension model, and interfacial phase change model, along with a model for effects of shear-lift and bubble collision dispersion to overcome a fundamental weakness in modeling multiphase flows.  Detailed information about bubble distribution in the vicinity of the heated surface, thermal conduction inside the heating wall, local heat fluxes passing through the solid-fluid interface, and velocity and temperature profiles, which are not easily observed or measured by experiments, is carefully evaluated.  The simulation results are compared to experimental data to validate the solver’s ability to predict the flow configuration with single/double-side heating.  The added momentum by shear-lift is shown to govern primarily the dynamic behavior of tiny bubbles stuck on the heated bottom wall and therefore has a more significant impact on both heat transfer and heated wall temperature.  By including bubble collision dispersion force, coalescence of densely packed bubbles in the bulk region is significantly inhibited, with more giant bubbles even incurring additional breakup into smaller bubbles and culminating in far less vapor accumulation along the top wall.  Including these momentums is shown to yield better agreement with local interfacial behavior along the channel, overall flow pattern, and heat transfer parameters (wall temperature and heat transfer coefficient) observed and measured in experiments.  The computational approach is also shown to be highly effective at predicting local phenomena (velocity and temperature profiles) not easily determined through experiments.  Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms, including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with the experiment.  Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations.  Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method.  Predicted wall temperature is reasonably uniform in the middle of the heated length but increases in the entrance region due to sensible heat transfer in the subcooled liquid and decreases toward the exit, primarily because of flow acceleration resulting from increased void fraction.  When it comes to analyzing heat transfer mechanisms at extremely high heat flux via CFD, predicted flow pattern, bubble behavior, and heat transfer parameters (such as wall temperature excursion and thermal energy concentration) clearly represent phenomena of premature CHF, which take place slightly earlier than actual operating conditions.  But, despite these slight differences, the present computational work does demonstrate the ability to effectively predict the severe degradation in heat transfer performance commonly encountered at heat fluxes nearing CHF.  </p> <p>Much of the published literature addressing flow instabilities in thermal management systems employing micro-channel modules are focused on instability characteristics of the module alone, and far fewer studies have aimed at understanding the relationship between these characteristics and compressive volume in the flow loop external to the module.  From a practical point of view, developers of micro-channel thermal management systems for many modern applications are in pursuit of practical remedies that would significantly mitigate instabilities and their impact on cooling performance.  Experiments are executed using FC-72 as a working fluid with a wide range of mass velocities and a reasonably constant inlet subcooling of ~15°C.  The flow instabilities are reflected in pressure fluctuations detected mainly in the heat sink’s upstream plenum.  Both inlet pressure and pressure drop signals are analyzed in pursuit of amplitude and frequency characteristics for different mass velocities and over a range of heat fluxes.  The current experimental study also examines the effects of compressible volume location in a closed pump-driven flow loop designed to deliver FC-72 to a micro-channel test module having 38 channels with 315-μm hydraulic diameter.  Three accumulator locations are investigated: upstream of the test module, downstream of the test module, and between the condenser and pump.  Both high-frequency temporal parameter data and high-speed video records are analyzed for ranges of mass velocity and heat flux, with inlet subcooling held constant at ~15°C.  PDO is shown to dominate when the accumulator is situated upstream, whereas PCI is dominant for the other two locations.  Appreciable confinement of bubbles in individual channels is shown to promote rapid axial bubble growth.  The study shows significant variations in the amount of vapor generated and dominant flow patterns among channels, a clear manifestation of PCI, especially for low mass velocities and high heat fluxes.  It is also shown effects of the heat sink’s instabilities are felt in other components of the flow loop.  The parametric trends for PCI are investigated with the aid of three different types of stability maps which show different abilities at demarcating stable and unstable operations.  PDO shows severe pressure oscillations across the micro-channel heat sink, with rapid bubble growth and confinement, elongated bubble expansion in both directions, flow stagnation, and flow reversal (including vapor backflow to the inlet plenum) constituting the principal sequence of events characterizing the instability.  Spectral analysis of pressure signals is performed using Fast Fourier Transform, which shows PDO extending the inlet pressure fluctuations with the same dominant frequency to other upstream flow loop components, with higher amplitudes closer to the pump exit.  From a practical system operation point of view, throttling the flow upstream of the heat sink eliminates PDO but renders PCI dominant, and placing the accumulator in the liquid flow segment of the loop between the condenser and pump ensures the most stable operation.</p>

<p dir="ltr">As aircraft engine manufacturers continue to embark on their pursuit of higher-efficiency, lower-emissions gas turbines, a prevailing theme in the industry has been the increase of the engine bypass ratio. As the optimization space for engine bypass ratios trends towards smaller and smaller engine core sizes, the feasibility of centrifugal compressors as the final stage in an axial-centrifugal compressor becomes apparent due to their performance advantages at smaller scales.</p><p dir="ltr">This study performed an investigation into the aerodynamics of a stationary diffusion system intended for use with a final stage aeroengine centrifugal compressor using experimental and numerical techniques. Experimental work was performed at the Purdue Compressor Research Lab at Purdue University’s Maurice J. Zucrow Laboratories. Data were collected from several iterations of rapidly prototyped, additively manufactured diffuser and deswirl parts with internal instrumentation features. Furthermore, computational work on the stage was conducted using the Ansys Turbosystem.</p><p dir="ltr">The goal of this research is to identify trends in stationary diffusion system designs and the geometric features that cause them. Furthermore, the ability of steady computational fluid dynamics methods to predict these changes was evaluated using two turbulence models to produce a simulation of the compressor flow field. When used in conjunction with one another, the efficient use of computational methods to identify an optimal design and rapid prototyping to validate it leads to the determination of the best diffusion system design at a lower cost and time requirement than what is otherwise currently possible.</p><p dir="ltr">The different geometries which were tested identified the negative effects of spanwise vane contouring on the diffuser performance and the ability of endwall divergence to augment the pressure recovery performance of a design at the expense of increased losses. A full understanding of the effect of each design parameter is enabled by iterative inclusion or exclusion of certain design parameters. Furthermore, the use of computational fluid dynamics showed that the BSLEARSM turbulence model performs reasonably well in predicting the build-to-build performance trends. However, neither the BSLEARSM nor the SST turbulence model were able to accurately identify performance trends for the deswirl. For this reason, additional work is warranted to identify an optimal set of parameters to characterize the high axial and meridional turning present in this component.</p>

The thermal efficiency of gas turbines increases with the temperature of the gas entering its turbine component. To enable high inlet temperatures, even those that far exceed the melting point of the turbine materials, the turbine must be cooled. One way is by internal cooling, where cooler air passes through U-ducts embedded inside turbine vanes and blades. Since the flow and heat transfer in these ducts are highly complicated, computational fluid dynamics (CFD) based on RANS have been used extensively to explore and assess design concepts. However, RANS have been found to be unreliable – giving accurate results for some designs but not for others. In this study, large-eddy simulations (LES) were performed for a U-duct with a trapezoidal cross section to assess four widely used RANS turbulence models: realizable k-ε (k-ε), shear-stress transport (SST), Reynolds stress model with linear pressure strain (RSM-LPS), and the seven-equation stress-omega full Reynolds stress model (RSM).<div><br></div><div>When examining the capability of steady RANS, two versions of the U-duct were examined, one with a staggered array of pin fins and one without pin fins. Results obtained for the heat-transfer coefficient (HTC) were compared with experimental measurements. The maximum relative error in the predicted “averaged” HTC was found to be 50% for k-ε and RSM-LPS, 20% for SST, and 30% for RSM-τω when there are no pin fins and 25% for k-ε, 12% for the SST and RSM-τω when there are pin fins. When there are no pin fins, all RANS models predicted a large separated flow region downstream of the turn, which the experiment does show to exist. Thus, all models predicted local distributions poorly. When there were pin fins, they behaved like guide vanes in turning the flow and confined the separation around the turn. For this configuration, all RANS models predicted reasonably well.<br></div><div><br></div><div>To understand why RANS cannot predict the HTC in the U-duct after the turn when there are no pin fins, LES were performed. To ensure that the LES is benchmark quality, verification and validation were performed via LES of a straight duct with square cross section where data from experiments and direct numerical simulation (DNS) are available. To ensure correct inflow boundary condition is provided for the U-duct, a concurrent LES is performed of a straight duct with the same trapezoidal cross section and flow conditions as the U-duct. Results obtained for the U-duct show RANS models to be inadequate in predicting the separation due to their inability to predict the unsteady separation about the tip of the turn. To investigate the limitations of the RANS models, LES results were generated for the turbulent kinetic energy, Reynolds-stresses, pressure-strain rate, turbulent diffusion, pressure diffusion, turbulent transport, and velocity-temperature correlations with focus on understanding their behavior induced by the turn region of the U-duct. As expected, the Boussinesq assumption was found to be incorrect, which led to incorrect predictions of Reynolds stresses. For RSM-τω, the modeling of the pressure-strain rate was found to match LES data well, but huge error was found on modeling the turbulent diffusion. This huge error indicates that the two terms in the turbulent diffusion – pressure diffusion and turbulent transport – should be modeled separately. Since the turbulent transport was found to be ignorable, the focus should be on modeling the pressure diffusion. On the velocity-temperature correlations, the existing eddy-diffusivity model was found to be over simplified if there is unsteady separation with shedding. The generated LES data could be used to provide the guidance for a better model.<br></div>

Este trabalho propõe um modelo de transferência de calor e massa para prever a formação de gelo em um container preenchido com legumes congelados. O problema físico é modelado como um meio poroso composto pelo próprio produto e o ar em seu entorno. O regime de convecção natural é assumido dentro do container, o qual promove o transporte de massa. Como uma primeira validação, o modelo é simulado considerando diferentes temperaturas de ar externo, causadas por flutuações da vizinhança. Resultados para quatro ciclos de temperaturas foram comparados, variando separadamente a temperatura média do ar, amplitude e frequência de oscilação. De modo geral, é observado que a temperatura do produto se comporta assim como era esperado e este resultado é diretamente associado à formação de gelo dentro do container. A formação de gelo cresce com uma maior amplitude de oscilação, porém decresce com um aumento na frequência e na temperatura média. Os parâmetros do modelo foram obtidos para dois diferentes produtos: fatias de cenouras congeladas e vagens congeladas, ambos em meio ao ar. As definições de parâmetros são oriundas de revisão bibliográfica, medições experimentais e simulações numéricas. Os parâmetros encontrados para a caracterização desses meios porosos foram similares para ambos os produtos, mesmo eles possuindo diferentes geometrias. A validação experimental foi feita para as fatias de cenoura considerando dois ciclos de temperatura O modelo numérico é capaz de prever o campo de velocidades do ar, as temperaturas do produto e a formação de gelo local. Os resultados foram validados em relação a um grupo independente de resultados numéricos, tal comparação apresentou uma boa concordância. A circulação de ar encontrada é, de fato, devido à convecção natural. O comportamento da temperatura dos produtos simulados concorda com os valores medidos e os valores de temperaturas diferem por menos de 12%. Com respeito à formação de gelo, o modelo é capaz de prevê-la corretamente nas regiões mais suscetíveis a este fenômeno. Porém, a quantidade de gelo formado prevista pelo modelo (1,56 g/semana) é menor do que a experimental (4,67 g/semana), apesar de serem de mesma ordem de magnitude. O efeito de cada parâmetro no modelo é estudado visando detectar maneiras de aprimorar o modelo. Foi encontrado que os parâmetros mais importantes para a formação de gelo total são a difusividade de massa efetiva e o coeficiente de transferência de calor convectivo dentro do container. Ajustando estes parâmetros duas vezes foi possível encontrar resultados melhores com respeito à formação de gelo (3,09 g/semana). / A model of heat and mass transfer is proposed in order to predict frost formation into a closed container filled with frozen vegetables. The physical problem is modeled as a macroporous media composed by the product itself and the surrounding air. Natural convection air flow is assumed into the container, who promotes water mass transport. As a first validation, the model is simulated for several exterior air temperatures, under environmental fluctuations (boundary conditions). Results of four temperature cycles were compared, varying average air temperature, amplitude and frequency of oscillation, one by one. As a general result, it is observed that the product temperature behavior is as expected, and it is directly associated with frost formation into the container. Frost formation increases with large amplitude of oscillation, but decreases with higher frequencies and higher mean temperatures. Model parameters were obtained for two assembling: frozen slices of carrots and air, and frozen extra thin green beans and air. Parameter definition and evaluation combines literature review, measurements and numerical simulation. In general, parameters which characterize these porous media were similar for both products, even though they display different geometries. The experimental validation is performed for carrot slices with two temperature cycles The numerical model is able to predict air velocity field, air and product temperatures, and local frost formation. Results are validated in respect to a set of independent experimental results that shown a good agreement. Air flow circulation is as expected due to natural convection. Product temperature simulated behavior agrees with measurements, and temperature values differ by less than 12%. Respect to frost formation predictions, the model predicts correctly the most susceptible regions to frost formation. However, the quantity of frost formed predicted by the model (1.56 g/ week)is lower than the experimental one (4.67g/week), despite being of the same order of magnitude. The effect of each parameter in the model is study in order to detect how to improve the model. The most important parameters affecting total frost formation are effective mass diffusivity and convective heat coefficient into the storage container. Adjusting these parameters to twice, better results in terms of frost formation could be obtained (3.09 g/ week).

Periciliary liquid (PCL) is a critical component of the respiratory system for maintaining mucus clearance. As PCL homeostasis is affected by evaporation and mechanical forces, which are in turn affected by various breathing conditions, lung morphology and ventilation distribution, the complex process of PCL depth regulation in vivo is not fully understood. We propose an integrative approach to couple a thermo-fluid computational fluid dynamics (CFD) model with an epithelial cell model to study the dynamics of PCL depth using subject-specific human airway models based on multi-detector row computed-tomography (MDCT) volumetric lung images. The thermo-fluid CFD model solves three-dimensional (3D) incompressible Navier-Stokes and transport equations for temperature and water vapor concentration with a realistic energy flux based boundary condition imposed at airway wall. A corresponding one-dimensional (1D) thermo-fluid CFD model is also developed to provide necessary information to the 3D model. Both 1D and 3D models are validated with experimental measurements, and the temperature and humidity distributions in the airways are investigated. Correlations for the dimensionless parameters of Nusselt number and Sherwood number are proposed for characterizing heat and mass transfer in the airways. As one of the key applications of the thermo-fluid CFD model, the water loss rates in the both 1D and 3D airway models are studied. It is found that the secondary flows formed at the bifurcations elevate the regional heat and mass transfer during inspiration and hence the water loss rate, which can only be observed in the 3D models. Among the three human airway models studied in both 1D and 3D, little inter-subject variability is observed for the distributions of temperature and humidity. However, the inter-subject variability could be dramatic for the distribution of water loss rate, as it is greatly affected by airway diameter and regional ventilation. A method is proposed to construct an ion-channel conductance model for both normal and cystic fibrosis (CF) epithelial cells, which couples an existing fluid secretion model with an existing nucleotide and nucleoside metabolism model (collectively named epithelial cell model). The epithelial cell models for both normal and CF are capable of predicting PCL depth based on mechanical stresses and evaporation, and are validated with a wide range of experimental data. With these two models separately validated and tested, the integrated model of the thermo-fluid CFD model and epithelial cell model is applied to MDCT-based human airway models of three CF subjects and three normal subjects to study and compare PCL depth regulation under regular breathing conditions. It is found that evaporative water loss is the dominant factor in PCL homeostasis. Between three types of mechanical forces, cyclic shear stress is the primary factor that triggers ATP release and increases PCL depth. In addition, it is found that that greater diameters of the airways in the 4th-7th generations in CF subjects decrease evaporative water loss, resulting in similar PCL depth as normal subjects. Under regular breathing conditions, the average PCL depths of normal and CF is around 6 to 7 µm, with mechanical forces play a greater role in regulating CF PCL depth. Comparing to 7.68 µm normal base level (considered as optimum PCL depth), this average PCL depth is about 8 to 21% lower. This might suggest that mechanical forces alone cannot entirely balance evaporative water loss, and other mechanisms might be involved.