Extraction of volatile oil from aromatic plants with supercritical carbon dioxide: experiments and modeling

An overview of the studies carried out in our laboratories on supercritical fluid extraction (SFE) of volatile oils from seven aromatic plants: pennyroyal (Mentha pulegium L.), fennel seeds (Foeniculum vulgare Mill.), coriander (Coriandrum sativum L.), savory (Satureja fruticosa Beguinot), winter savory (Satureja montana L.), cotton lavender (Santolina chamaecyparisus) and thyme (Thymus vulgaris), is presented. A flow apparatus with a 1 L extractor and two 0.27 L separators was built to perform studies at temperatures ranging from 298 to 353 K and pressures up to 30.0 MPa. The best compromise between yield and composition compared with hydrodistillation (HD) was achieved selecting the optimum experimental conditions of extraction and fractionation. The major differences between HD and SFE oils is the presence of a small percentage of cuticular waxes and the relative amount of thymoquinone, an oxygenated monoterpene with important biological properties, which is present in the oils from thyme and winter savory. On the other hand, the modeling of our data on supercritical extraction of volatile oil from pennyroyal is discussed using Sovova's models. These models have been applied successfully to the other volatile oil extractions. Furthermore, other experimental studies involving supercritical CO2 carried out in our laboratories are also mentioned.

Mathematical modeling and simulation of heat flow through pultrusion die assembly system for thermoset composites

Pultrusion is an industrial process used to produce glass fibers reinforced polymers profiles. These materials are worldwide used when performing characteristics, such as great electrical and magnetic insulation, high strength to weight ratio, corrosion and weather resistance, long service life and minimal maintenance are required. In this study, we present the results of the modelling and simulation of heat flow through a pultrusion die by means of Finite Element Analysis (FEA). The numerical simulation was calibrated based on temperature profiles computed from thermographic measurements carried out during pultrusion manufacturing process. Obtained results have shown a maximum deviation of 7%, which is considered to be acceptable for this type of analysis, and is below to the 10% value, previously specified as maximum deviation. © 2011, Advanced Engineering Solutions.

Experimental and modeling studies on reverse electrodialysis for sustainable power generation

A potentially renewable and sustainable source of energy is the chemical energy associated with solvation of salts. Mixing of two aqueous streams with different saline concentrations is spontaneous and releases energy. The global theoretically obtainable power from salinity gradient energy due to World’s rivers discharge into the oceans has been estimated to be within the range of 1.4-2.6 TW. Reverse electrodialysis (RED) is one of the emerging, membrane-based, technologies for harvesting the salinity gradient energy. A common RED stack is composed by alternately-arranged cation- and anion-exchange membranes, stacked between two electrodes. The compartments between the membranes are alternately fed with concentrated (e.g., sea water) and dilute (e.g., river water) saline solutions. Migration of the respective counter-ions through the membranes leads to ionic current between the electrodes, where an appropriate redox pair converts the chemical salinity gradient energy into electrical energy. Given the importance of the need for new sources of energy for power generation, the present study aims at better understanding and solving current challenges, associated with the RED stack design, fluid dynamics, ionic mass transfer and long-term RED stack performance with natural saline solutions as feedwaters. Chronopotentiometry was used to determinate diffusion boundary layer (DBL) thickness from diffusion relaxation data and the flow entrance effects on mass transfer were found to avail a power generation increase in RED stacks. Increasing the linear flow velocity also leads to a decrease of DBL thickness but on the cost of a higher pressure drop. Pressure drop inside RED stacks was successfully simulated by the developed mathematical model, in which contribution of several pressure drops, that until now have not been considered, was included. The effect of each pressure drop on the RED stack performance was identified and rationalized and guidelines for planning and/or optimization of RED stacks were derived. The design of new profiled membranes, with a chevron corrugation structure, was proposed using computational fluid dynamics (CFD) modeling. The performance of the suggested corrugation geometry was compared with the already existing ones, as well as with the use of conductive and non-conductive spacers. According to the estimations, use of chevron structures grants the highest net power density values, at the best compromise between the mass transfer coefficient and the pressure drop values. Finally, long-term experiments with natural waters were performed, during which fouling was experienced. For the first time, 2D fluorescence spectroscopy was used to monitor RED stack performance, with a dedicated focus on following fouling on ion-exchange membrane surfaces. To extract relevant information from fluorescence spectra, parallel factor analysis (PARAFAC) was performed. Moreover, the information obtained was then used to predict net power density, stack electric resistance and pressure drop by multivariate statistical models based on projection to latent structures (PLS) modeling. The use in such models of 2D fluorescence data, containing hidden, but extractable by PARAFAC, information about fouling on membrane surfaces, considerably improved the models fitting to the experimental data.

Optimization of filling systems for low pressure by Flow 3D

Dissertação de mestrado integrado em Engenharia Mecânica

Optimization of radial jets mixing in cross-flow of combustion chambers using computational fluid dynamics

Cross-Flow, Radial Jets Mixing, Temperature Homogenization, Optimization, Combustion Chamber, CFD

Modeling the inelastic behavior of heat exchangers accounting for fluid-structure interactions

Magdeburg, Univ., Fak. für Maschinenbau, Diss., 2013

An approximate solution method for boundary layer flow of a power law fluid over a flat plate

The work in this paper deals with the development of momentum and thermal boundary layers when a power law fluid flows over a flat plate. At the plate we impose either constant temperature, constant flux or a Newton cooling condition. The problem is analysed using similarity solutions, integral momentum and energy equations and an approximation technique which is a form of the Heat Balance Integral Method. The fluid properties are assumed to be independent of temperature, hence the momentum equation uncouples from the thermal problem. We first derive the similarity equations for the velocity and present exact solutions for the case where the power law index n = 2. The similarity solutions are used to validate the new approximation method. This new technique is then applied to the thermal boundary layer, where a similarity solution can only be obtained for the case n = 1.

Computational fluid dynamics of the right ventricular outflow tract and of the pulmonary artery: a bench model of flow dynamics.

OBJECTIVES: The reconstruction of the right ventricular outflow tract (RVOT) with valved conduits remains a challenge. The reoperation rate at 5 years can be as high as 25% and depends on age, type of conduit, conduit diameter and principal heart malformation. The aim of this study is to provide a bench model with computer fluid dynamics to analyse the haemodynamics of the RVOT, pulmonary artery, its bifurcation, and left and right pulmonary arteries that in the future may serve as a tool for analysis and prediction of outcome following RVOT reconstruction. METHODS: Pressure, flow and diameter at the RVOT, pulmonary artery, bifurcation of the pulmonary artery, and left and right pulmonary arteries were measured in five normal pigs with a mean weight of 24.6 ± 0.89 kg. Data obtained were used for a 3D computer fluid-dynamics simulation of flow conditions, focusing on the pressure, flow and shear stress profile of the pulmonary trunk to the level of the left and right pulmonary arteries. RESULTS: Three inlet steady flow profiles were obtained at 0.2, 0.29 and 0.36 m/s that correspond to the flow rates of 1.5, 2.0 and 2.5 l/min flow at the RVOT. The flow velocity profile was constant at the RVOT down to the bifurcation and decreased at the left and right pulmonary arteries. In all three inlet velocity profiles, low sheer stress and low-velocity areas were detected along the left wall of the pulmonary artery, at the pulmonary artery bifurcation and at the ostia of both pulmonary arteries. CONCLUSIONS: This computed fluid real-time model provides us with a realistic picture of fluid dynamics in the pulmonary tract area. Deep shear stress areas correspond to a turbulent flow profile that is a predictive factor for the development of vessel wall arteriosclerosis. We believe that this bench model may be a useful tool for further evaluation of RVOT pathology following surgical reconstructions.

Fluid dynamical prediction of changed v1 flow at energies available at the CERN Large Hadron Collider

Substantial collective flow is observed in collisions between lead nuclei at Large Hadron Collider (LHC) as evidenced by the azimuthal correlations in the transverse momentum distributions of the produced particles. Our calculations indicate that the global v1-flow, which at RHIC peaked at negative rapidities (named third flow component or antiflow), now at LHC is going to turn toward forward rapidities (to the same side and direction as the projectile residue). Potentially this can provide a sensitive barometer to estimate the pressure and transport properties of the quark-gluon plasma. Our calculations also take into account the initial state center-of-mass rapidity fluctuations, and demonstrate that these are crucial for v1 simulations. In order to better study the transverse momentum flow dependence we suggest a new "symmetrized" v1S(pt) function, and we also propose a new method to disentangle global v1 flow from the contribution generated by the random fluctuations in the initial state. This will enhance the possibilities of studying the collective Global v1 flow both at the STAR Beam Energy Scan program and at LHC.

Patterns of gravity induced aggregate migration during casting of fluid concretes

In this paper, aggregate migration patterns during fluid concrete castings are studied through experiments, dimensionless approach and numerical modeling. The experimental results obtained on two beams show that gravity induced migration is primarily affecting the coarsest aggregates resulting in a decrease of coarse aggregates volume fraction with the horizontal distance from the pouring point and in a puzzling vertical multi-layer structure. The origin of this multi layer structure is discussed and analyzed with the help of numerical simulations of free surface flow. Our results suggest that it finds its origin in the non Newtonian nature of fresh concrete and that increasing casting rate shall decrease the magnitude of gravity induced particle migration. (C) 2012 Elsevier Ltd. All rights reserved.

A framework for hybrid simulations of two-phase flow in porous media

Rough a global coarse problem. Although these techniques are usually employed for problems in which the fine-scale processes are described by Darcy's law, they can also be applied to pore-scale simulations and used as a mathematical framework for hybrid methods that couples a Darcy and pore scales. In this work, we consider a pore-scale description of fine-scale processes. The Navier-Stokes equations are numerically solved in the pore geometry to compute the velocity field and obtain generalized permeabilities. In the case of two-phase flow, the dynamics of the phase interface is described by the volume of fluid method with the continuum surface force model. The MsFV method is employed to construct an algorithm that couples a Darcy macro-scale description with a pore-scale description at the fine scale. The hybrid simulations results presented are in good agreement with the fine-scale reference solutions. As the reconstruction of the fine-scale details can be done adaptively, the presented method offers a flexible framework for hybrid modeling.

Fluid dynamical prediction of changed v1-flow at energies available at the CERN Large Hadron Collider

Substantial collective flow is observed in collisions between lead nuclei at Large Hadron Collider (LHC) as evidenced by the azimuthal correlations in the transverse momentum distributions of the produced particles. Our calculations indicate that the global v1-flow, which at RHIC peaked at negative rapidities (named third flow component or antiflow), now at LHC is going to turn toward forward rapidities (to the same side and direction as the projectile residue). Potentially this can provide a sensitive barometer to estimate the pressure and transport properties of the quark-gluon plasma. Our calculations also take into account the initial state center-of-mass rapidity fluctuations, and demonstrate that these are crucial for v1 simulations. In order to better study the transverse momentum flow dependence we suggest a new"symmetrized" vS1(pt) function, and we also propose a new method to disentangle global v1 flow from the contribution generated by the random fluctuations in the initial state. This will enhance the possibilities of studying the collective Global v1 flow both at the STAR Beam Energy Scan program and at LHC.

Extension of Murray's law using a non-Newtonian model of blood flow.

BACKGROUND: So far, none of the existing methods on Murray's law deal with the non-Newtonian behavior of blood flow although the non-Newtonian approach for blood flow modelling looks more accurate. MODELING: In the present paper, Murray's law which is applicable to an arterial bifurcation, is generalized to a non-Newtonian blood flow model (power-law model). When the vessel size reaches the capillary limitation, blood can be modeled using a non-Newtonian constitutive equation. It is assumed two different constraints in addition to the pumping power: the volume constraint or the surface constraint (related to the internal surface of the vessel). For a seek of generality, the relationships are given for an arbitrary number of daughter vessels. It is shown that for a cost function including the volume constraint, classical Murray's law remains valid (i.e. SigmaR(c) = cste with c = 3 is verified and is independent of n, the dimensionless index in the viscosity equation; R being the radius of the vessel). On the contrary, for a cost function including the surface constraint, different values of c may be calculated depending on the value of n. RESULTS: We find that c varies for blood from 2.42 to 3 depending on the constraint and the fluid properties. For the Newtonian model, the surface constraint leads to c = 2.5. The cost function (based on the surface constraint) can be related to entropy generation, by dividing it by the temperature. CONCLUSION: It is demonstrated that the entropy generated in all the daughter vessels is greater than the entropy generated in the parent vessel. Furthermore, it is shown that the difference of entropy generation between the parent and daughter vessels is smaller for a non-Newtonian fluid than for a Newtonian fluid.

Fluid dynamical prediction of changed v1-flow at energies available at the CERN Large Hadron Collider

Substantial collective flow is observed in collisions between lead nuclei at Large Hadron Collider (LHC) as evidenced by the azimuthal correlations in the transverse momentum distributions of the produced particles. Our calculations indicate that the global v1-flow, which at RHIC peaked at negative rapidities (named third flow component or antiflow), now at LHC is going to turn toward forward rapidities (to the same side and direction as the projectile residue). Potentially this can provide a sensitive barometer to estimate the pressure and transport properties of the quark-gluon plasma. Our calculations also take into account the initial state center-of-mass rapidity fluctuations, and demonstrate that these are crucial for v1 simulations. In order to better study the transverse momentum flow dependence we suggest a new"symmetrized" vS1(pt) function, and we also propose a new method to disentangle global v1 flow from the contribution generated by the random fluctuations in the initial state. This will enhance the possibilities of studying the collective Global v1 flow both at the STAR Beam Energy Scan program and at LHC.

Pain-evoked alterations on gingival blood flow and gingival crevicular fluid (GCF) neuropeptide SP and collagenase-2 (MMP-8) levels

Previous studies have demonstrated that clinical pulpal pain can induce the expression of pro-inflammatory neuropeptides in the adjacent gingival crevice fluid (GCF). Vasoactive agents such as substance P (SP) are known to contribute to the inflammatory type of pain and are associated with increased blood flow. More recent animal studies have shown that application of capsaicin on alveolar mucosa provokes pain and neurogenic vasodilatation in the adjacent gingiva. Pain-associated inflammatory reactions may initiate expression of several pro- and anti-inflammatory mediators. Collagenase-2 (MMP-8) has been considered to be the major destructive protease, especially in the periodontitis-affected gingival crevice fluid (GCF). MMP-8 originates mostly from neutrophil leukocytes, the first line of defence cells that exist abundantly in GCF, especially in inflammation. With this background, we wished to clarify the spatial extensions and differences between tooth-pain stimulation and capsaicin-induced neurogenic vasodilatation in human gingiva. Experiments were carried out to study whether tooth stimulation and capsaicin stimulation of alveolar mucosa would induce changes in GCF MMP-8 levels and whether tooth stimulation would release neuropeptide SP in GCF. The experiments were carried out on healthy human volunteers. During the experiments, moderate and high intensity painful tooth stimulation was performed by a constant current tooth stimulator. Moderate tooth stimulation activates A-delta fibres, while high stimulation also activates C-fibres. Painful stimulation of the gingiva was achieved by topical application of capsaicin-moistened filter paper on the mucosal surface. Capsaicin is known to activate selectively nociceptive C-fibres of stimulated tissue. Pain-evoked vasoactive changes in gingivomucosal tissues were mapped by laser Doppler imaging (LDI), which is a sophisticated and non-invasive method for studying e.g. spatial and temporal characteristics of pain- and inflammation-evoked blood flow changes in gingivomucosal tissues. Pain-evoked release of MMP-8 in GCF samples was studied by immunofluorometric assay (IFMA) and Western immunoblotting. The SP levels in GCF were analysed by Enzyme immunoassay (EIA). During the experiments, subjective stimulus-evoked pain responses were determined by a visual analogue pain scale. Unilateral stimulation of alveolar mucosa and attached gingiva by capsaicin evoked a distinct neurogenic vasodilatation in the ipsilateral gingiva, which attenuated rapidly at the midline. Capsaicin stimulation of alveolar mucosa provoked clear inflammatory reactions. In contrast to capsaicin stimuli, tooth stimulation produced symmetrical vasodilatations bilaterally in the gingiva. The ipsilateral responses were significantly smaller during tooth stimulation than during capsaicin stimuli. The current finding – that tooth stimulation evokes bilateral vasodilatation while capsaicin stimulation of the gingiva mainly produces unilateral vasodilatation – emphasises the usefulness of LDI in clarifying spatial features of neurogenic vasoactive changes in the intra-oral tissues. Capsaicin stimulation of the alveolar mucosa induced significant elevations in MMP-8 levels and activation in GCF of the adjacent teeth. During the experiments, no marked changes occurred in MMP-8 levels in the GCF of distantly located teeth. Painful stimulation of the upper incisor provoked elevations in GCF MMP-8 and SP levels of the stimulated tooth. The GCF MMP-8 and SP levels of the non-stimulated teeth were not changed. These results suggest that capsaicin-induced inflammatory reactions in gingivomucosal tissues do not cross the midline in the anterior maxilla. The enhanced reaction found during stimulation of alveolar mucosa indicates that alveolar mucosa is more sensitive to chemical irritants than the attached gingiva. Analysis of these data suggests that capsaicin-evoked neurogenic inflammation in the gingiva can trigger the expression and activation of MMP-8 in GCF of the adjacent teeth. In this study, it is concluded that experimental tooth pain at C-fibre intensity can induce local elevations in MMP-8 and SP levels in GCF. Depending on the role of MMP-8 in inflammation, in addition to surrogated tissue destruction, the elevated MMP-8 in GCF may also reflect accelerated local defensive and anti-inflammatory reactions.