Thermal food processing optimization: Algorithms and software

The algorithms and graphic user interface software package ?OPT-PROx? are developed to meet food engineering needs related to canned food thermal processing simulation and optimization. The adaptive random search algorithm and its modification coupled with penalty function?s approach, and the finite difference methods with cubic spline approximation are utilized by ?OPT-PROx? package (http://tomakechoice. com/optprox/index.html). The diversity of thermal food processing optimization problems with different objectives and required constraints are solvable by developed software. The geometries supported by the ?OPT-PROx? are the following: (1) cylinder, (2) rectangle, (3) sphere. The mean square error minimization principle is utilized in order to estimate the heat transfer coefficient of food to be heated under optimal condition. The developed user friendly dialogue and used numerical procedures makes the ?OPT-PROx? software useful to food scientists in research and education, as well as to engineers involved in optimization of thermal food processing.

On the flexural vibration of cylinders under axial loads:Numerical and experimental study

The flexural vibration of a homogeneous isotropic linearly elastic cylinder of any aspect ratio is analysed in this paper. Natural frequencies of a cylinder under uniformly distributed axial loads acting on its bases are calculated numerically by the Ritz method with terms of power series in the coordinate directions as approximating functions. The effect of axial loads on the flexural vibration cannot be described by applying infinitesimal strain theory, therefore, geometrically nonlinear strain–displacement relations with second-order terms are considered here. The natural frequencies of free–free, clamped–clamped, and sliding–sliding cylinders subjected to axial loads are calculated using the proposed three-dimensional Ritz approach and are compared with those obtained with the finite element method and the Bernoulli–Euler theory. Different experiments with cylinders axially compressed by a hydraulic press are carried out and the experimental results for the lowest flexural frequency are compared with the numerical results. An approach based on the Ritz formulation is proposed for the flexural vibration of a cylinder between the platens of the press with constraints varying with the intensity of the compression. The results show that for low compressions the cylinder behaves similarly to a sliding–sliding cylinder, whereas for high compressions the cylinder vibrates as a clamped–clamped one.

A molecular dynamics study of swift heavy ion irradiation of amorphous silica: the role of thermal effects

Irradiation with swift heavy ions (SHI), roughly defined as those having atomic masses larger than 15 and energies exceeding 1 MeV/amu, may lead to significant modification of the irradiated material in a nanometric region around the (straight) ion trajectory (latent tracks). In the case of amorphous silica, SHI irradiation originates nano-tracks of higher density than the virgin material (densification). As a result, the refractive index is increased with respect to that of the surroundings. Moreover, track overlapping leads to continuous amorphous layers that present a significant contrast with respect to the pristine substrate. We have recently demonstrated that SHI irradiation produces a large number of point defects, easily detectable by a number of experimental techniques (work presented in the parallel conference ICDIM). The mechanisms of energy transfer from SHI to the target material have their origin in the high electronic excitation induced in the solid. A number of phenomenological approaches have been employed to describe these mechanisms: coulomb explosion, thermal spike, non-radiative exciton decay, bond weakening. However, a detailed microscopic description is missing due to the difficulty of modeling the time evolution of the electronic excitation. In this work we have employed molecular dynamics (MD) calculations to determine whether the irradiation effects are related to the thermal phenomena described by MD (in the ps domain) or to electronic phenomena (sub-ps domain), e.g., exciton localization. We have carried out simulations of up to 100 ps with large boxes (30x30x8 nm3) using a home-modified version of MDCASK that allows us to define a central hot cylinder (ion track) from which heat flows to the surrounding cold bath (unirradiated sample). We observed that once the cylinder has cooled down, the Si and O coordination numbers are 4 and 2, respectively, as in virgin silica. On the other hand, the density of the (cold) cylinder increases with respect to that of silica and, furthermore, the silica network ring size decreases. Both effects are in agreement with the observed densification. In conclusion, purely thermal effects do not explain the generation of point defects upon irradiation, but they do account for the silica densification.

Application of postprocessing techniques methods for noise reduction in cylindrical near field antenna measurements

Three different methods to reduce the noise power in the far-field pattern of an antenna when it is measured in a cylindrical near field system are presented and compared. The first one is based on a modal filtering while the other two are based on spatial filtering, either on an antenna plane or either on a cylinder of smaller radius. Simulated and measured results will be presented.

Caracterización de la impedancia de superficie para el cálculo del acoplo mutuo en un cilindro circular conductor cubierto por una dieléctrico multicapa

A novel formulation for the surface impedance characterization is introduced for the canonical problem of surface fields on a perfect electric conductor (PEC) circular cylinder with a dielectric coating due to a electric current source using the Uniform Theory of Diffraction (UTD) with an Impedance Boundary Condition (IBC). The approach is based on a TE/TM assumption of the surface fields from the original problem. Where this surface impedance fails, an optimization is performed to minimize the error in the SD Green's function between the original problem and the equivalent one with the IBC. This new approach requires small changes in the available UTD based solution with IBC to include the geodesic ray angle and length dependence in the surface impedance formulas. This asymptotic method, accurate for large separations between source and observer points, in combination with spectral domain (SD) Green's functions for multidielectric coatings leads to a new hybrid SD-UTD with IBC to calculate mutual coupling among microstrip patches on a multilayer dielectric-coated PEC circular cylinder. Results are compared with the eigenfunction solution in SD, where a very good agreement is met.

Influence of the composition of natural gas-hydrogen blends (HCNG) on performance an CO2 emissions of a spark ignition engine

Addition of hydrogen to natural gas could be a short-term alternative to nowadays fossil fuels as the emissions of greenhouse gases may be reduced. The aim of this study is to evaluate the performance and emissions of a park ignition engine fuelled with pure natural gas, pure hydrogen and different blends of hydrogen and natural gas (HCNG). The increase of the hydrogen fraction leads to variations in the cylinder pressure and CO2 emissions. In this work, a combustion model based on thermodynamic equations is used considering separated zones for the burned and unburned gases. The results show that the maximum cylinder pressure gets higher as the fraction of hydrogen in the blend increases. The presence of hydrogen in the blend leads to a drecrease in the CO2 emissions. Due to hydrogen properties, leaner fuel-air mixtures can be used along with the appropiate spark timing, leading to an engine emissions improvement without a performance worsening.

Monte Carlo simulation on a system of hard cylinders at a very high packing fraction

Using the Monte Carlo method the behavior of a system of true hard cylinders has been studied. Values of the length-to-breadth ratio L/D and packing fraction η have been chosen similar to those of real nematic liquid crystals. Results include radial distribution function g(r), structure factor S(k), and orientational order parameter M. These results lead to the conclusion that the hard cylinder model may be a useful reference for real mesomorphic phases.

Matrix-free time-stepping methods for the solution of TriGlobal instability problems

La inmensa mayoría de los flujos de relevancia ingenieril permanecen sin estudiar en el marco de la teoría de estabilidad global. Esto es debido a dos razones fundamentalmente, las dificultades asociadas con el análisis de los flujos turbulentos y los inmensos recursos computacionales requeridos para obtener la solución del problema de autovalores asociado al análisis de inestabilidad de flujos tridimensionales, también conocido como problema TriGlobal. En esta tesis se aborda el problema asociado con la tridimensionalidad. Se ha desarrollado una metodología general para obtener soluciones de problemas de análisis modal de las inestabilidades lineales globales mediante el acoplamiento de métodos de evolución temporal, desarrollados en este trabajo, con códigos de mecánica de fluidos computacional de segundo orden, utilizados de forma general en la industria. Esta metodología consiste en la resolución del problema de autovalores asociado al análisis de inestabilidad mediante métodos de proyección en subespacios de Krylov, con la particularidad de que dichos subespacios son generados por medio de la integración temporal de un vector inicial usando cualquier código de mecánica de fluidos computacional. Se han elegido tres problemas desafiantes en función de la exigencia de recursos computacionales necesarios y de la complejidad física para la demostración de la presente metodología: (i) el flujo en el interior de una cavidad tridimensional impulsada por una de sus tapas, (ii) el flujo alrededor de un cilindro equipado con aletas helicoidales a lo largo su envergadura y (iii) el flujo a través de una cavidad abierta tridimensinal en ausencia de homogeneidades espaciales. Para la validación de la tecnología se ha obtenido la solución del problema TriGlobal asociado al flujo en la cavidad tridimensional, utilizando el método de evolución temporal desarrollado acoplado con los operadores numéricos de flujo incompresible del código CFD OpenFOAM (código libre). Los resultados obtenidos coinciden plentamente con la literatura. La aplicación de esta metodología al estudio de inestabilidades globales de flujos abiertos tridimensionales ha proporcionado por primera vez, información sobre la transición tridimensional de estos flujos. Además, la metodología ha sido adaptada para resolver problemas adjuntos TriGlobales, permitiendo el control de flujo basado en modificaciones de las inestabilidades globales. Finalmente, se ha demostrado que la cantidad moderada de los recursos computacionales requeridos para la solución del problema de valor propio TriGlobal usando este método numérico, junto a su versatilidad al poder acoplarse a cualquier código aerodinámico, permite la realización de análisis de inestabilidad global y control de flujos complejos de relevancia industrial. Abstract Most flows of engineering relevance still remain unexplored in a global instability theory context for two reasons. First, because of the difficulties associated with the analysis of turbulent flows and, second, for the formidable computational resources required for the solution of the eigenvalue problem associated with the instability analysis of three-dimensional base flows, also known as TriGlobal problem. In this thesis, the problem associated with the three-dimensionality is addressed by means of the development of a general approach to the solution of large-scale global linear instability analysis by coupling a time-stepping approach with second order aerodynamic codes employed in industry. Three challenging flows in the terms of required computational resources and physical complexity have been chosen for demonstration of the present methodology; (i) the flow inside a wall-bounded three-dimensional lid-driven cavity, (ii) the flow past a cylinder fitted with helical strakes and (iii) the flow over a inhomogeneous three-dimensional open cavity. Results in excellent agreement with the literature have been obtained for the three-dimensional lid-driven cavity by using this methodology coupled with the incompressible solver of the open-source toolbox OpenFOAM®, which has served as validation. Moreover, significant physical insight of the instability of three-dimensional open flows has been gained through the application of the present time-stepping methodology to the other two cases. In addition, modifications to the present approach have been proposed in order to perform adjoint instability analysis of three-dimensional base flows and flow control; validation and TriGlobal examples are presented. Finally, it has been demonstrated that the moderate amount of computational resources required for the solution of the TriGlobal eigenvalue problem using this method enables the performance of instability analysis and control of flows of industrial relevance.

CO2 emissions from a spark ignition engine operatin on natural-hydrogen blends (HCNG)

The addition of hydrogen to natural gas could be a short-term alternative to today’s fossil fuels, as greenhouse gas emissions may be reduced. The aim of this study is to evaluate the emissions and performance of a spark ignition engine fuelled by pure natural gas, pure hydrogen, and different blends of hydrogen and natural gas (HCNG). Increasing the hydrogen fraction leads to variations in cylinder pressure and CO2 emissions. In this study, a combustion model based on thermodynamic equations is used, considering separate zones for burned and unburned gases. The results show that the maximum cylinder pressure rises as the fraction of hydrogen in the blend increases. The presence of hydrogen in the blend leads to a decrease in CO2 emissions. Due to the properties of hydrogen, leaner fuel–air mixtures can be used along with the appropriate spark timing, leading to an improvement in engine emissions with no loss of performance.

Magnetic self-field effects on bare-tether current collection

It has been recently suggested that the magnetic field created by the current in a bare tether could sensibly reduce its electron collection capability in the magnetised ionosphere, a region of closed magnetic surfaces disconnecting the cylinder from infinity. In this paper, the ohmic voltage drop along the tether is taken into account in considering self-field effects. Separate analyses are carried out for the thrust and power generation and drag modes of operation, which are affected in different ways. In the power generation and drag modes, bias decreases as current increases along the tether, starting at the anodic, positively-biased end (upper end in the usual, eastward-flying spacecraft); in the thrust mode of operation, bias increases as current increases along the tether, starting at the lower end. When the ohmic voltage drop is considered, self-field effects are shown to be weak, in all cases, for tape tethers, and for circular cross-section tethers just conductive in a thin outer layer. Self-field effects might become important, in the drag case only, for tethers with fully conductive cross sections that are unrealistically heavy.

An experimental and theoretical study of unsteady flow (gust) effects on structures

The gust wind tunnel at IDR, Universidad Politécnica de Madrid (UPM), has been enhanced and the impact of the modification has been characterized. Several flow quality configurations have been tested. The problems in measuring gusty winds with Pitot tubes have been considered. Experimental results have been obtained and compared with theoretically calculated results (based on potential flow theory). A theoretical correction term has been proposed for unsteady flow measurements obtained with Pitot tubes. The effect of unsteady flow on structures and laying bodies on the ground has been also considered. A theoretical model has been proposed for a semi-circular cylinder and experimental tests have been performed to study the unsteady flow effects, which can help in clarifying the phenomenon.

Measurements and Prediction of Pedestrian Walking Loads on an Aluminium Catwalk

During the last two decades the topic of human induced vibration has attracted a lot of attention among civil engineering practitioners and academics alike. Usually this type of problem may be encountered in pedestrian footbridges or floors of paperless offices. Slender designs are becoming increasingly popular, and as a consequence, the importance of paying attention to vibration serviceability also increases. This paper resumes the results obtained from measurements taken at different points of an aluminium catwalk which is 6 m in length by 0.6 m in width. Measurements were carried out when subjecting the structure to different actions:1)Static test: a steel cylinder of 35 kg was placed in the middle of the catwalk; 2)Dynamic test: this test consists of exciting the structure with singles impulses; 3)Dynamic test: people walking on the catwalk. Identification of the mechanical properties of the structure is an achievement of the paper. Indirect methods were used to estimate properties including the support stiffness, the beam bending stiffness, the mass of the structure (using Rayleigh method and iterative matrix method), the natural frequency (using the time domain and frequency domain analysis) and the damping ratio (by calculating the logarithmic decrement). Experimental results and numerical predictions for the response of an aluminium catwalk subjected to walking loads have been compared. The damping of this light weight structure depends on the amplitude of vibration which complicates the tuning of a structural model. In the light of the results obtained it seems that the used walking load model is not appropriate as the predicted transient vibration values (TTVs) are much higher than the measured ones.

Dynamic & Cyclic behaviour of ballast in the long term as determined in Cedex's track box

The 6 cylinder servo-hydraulic loading system of CEDEX's track box (250 kN, 50 Hz) has been recently implemented with a new piezoelectric loading system (±20 kN, 300 Hz) allowing the incorporation of low amplitude high frequency dynamic load time histories to the high amplitude low frequency quasi-static load time histories used so far in the CEDEX's track box to assess the inelastic long term behavior of ballast under mixed traffic in conventional and high- speed lines. This presentation will discuss the results obtained in the first long-duration test performed at CEDEX's track box using simultaneously both loading systems, to simulate the pass-by of 6000 freight vehicles (1M of 225 kN axle loads) travelling at a speed of 120 km/h over a line with vertical irregularities corresponding to a medium quality lin3e level. The superstructure of the track tested at full scale consisted of E 60 rails, stiff rail pads (mayor que 450 kN/mm), B90.2 sleepers with USP 0.10 N/mm and a 0.35 m thick ballast layer of ADIF first class. A shear wave velocity of 250 m/s can be assumed for the different layers of the track sub-base. The ballast long-term settlements will be compared with those obtained in a previous long-duration quasi- static test performed in the same track, for the RIVAS [EU co-funded] project, in which no dynamic loads where considered. Also, the results provided by a high diameter cyclic triaxial cell with ballast tested in full size will be commented. Finally, the progress made at CEDEX's Geotechnical Laboratory to reproduce numerically the long term behavior of ballast will be discussed.

Drifting plasma collection by a positive biased tether wire in LEO-like plasma conditions: current measurement and plasma diagnostic

BETs is a three-year project financed by the Space Program of the European Commission, aimed at developing an efficient deorbit system that could be carried on board any future satellite launched into Low Earth Orbit (LEO). The operational system involves a conductive tape-tether left bare to establish anodic contact with the ambient plasma as a giant Langmuir probe. As a part of this project, we are carrying out both numerical and experimental approaches to estimate the collected current by the positive part of the tether. This paper deals with experimental measurements performed in the IONospheric Atmosphere Simulator (JONAS) plasma chamber of the Onera-Space Environment Department. The JONAS facility is a 9- m3 vacuum chamber equipped with a plasma source providing drifting plasma simulating LEO conditions in terms of density and temperature. A thin metallic cylinder, simulating the tether, is set inside the chamber and polarized up to 1000 V. The Earth's magnetic field is neutralized inside the chamber. In a first time, tether collected current versus tether polarization is measured for different plasma source energies and densities. In complement, several types of Langmuir probes are used at the same location to allow the extraction of both ion densities and electron parameters by computer modeling (classical Langmuir probe characteristics are not accurate enough in the present situation). These two measurements permit estimation of the discrepancies between the theoretical collection laws, orbital motion limited law in particular, and the experimental data in LEO-like conditions without magnetic fields. In a second time, the spatial variations and the time evolutions of the plasma properties around the tether are investigated. Spherical and emissive Langmuir probes are also used for a more extensive characterization of the plasma in space and time dependent analysis. Results show the ion depletion because of the wake effect and the accumulation of- ions upstream of the tether. In some regimes (at large positive potential), oscillations are observed on the tether collected current and on Langmuir probe collected current in specific sites.

Extracting energy from Vortex-Induced Vibrations: A parametric study

Here, Vortex-Induced Vibrations (VIVs) of a circular cylinder are analyzed as a potential source for energy harvesting. To this end, VIV is described by a one-degree-of-freedom model where fluid forces are introduced from experimental data from forced vibration tests. The influence of some influencing parameters, like the mass ratio m∗ or the mechanical damping ζ in the energy conversion factor is investigated. The analysis reveals that: (i) the maximum efficiency ηM is principally influenced by the mass-damping parameter m∗ζ and there is an optimum value of m∗ζ where ηM presents a maximum; (ii) the range of reduced velocities with significant efficiency is mainly governed by m∗, and (iii) it seems that encouraging high efficiency values can be achieved for high Reynolds numbers.