Reduced order models for industrial thermo-fluid problems

A new method is presented to generate reduced order models (ROMs) in Fluid Dynamics problems of industrial interest. The method is based on the expansion of the flow variables in a Proper Orthogonal Decomposition (POD) basis, calculated from a limited number of snapshots, which are obtained via Computational Fluid Dynamics (CFD). Then, the POD-mode amplitudes are calculated as minimizers of a properly defined overall residual of the equations and boundary conditions. The method includes various ingredients that are new in this field. The residual can be calculated using only a limited number of points in the flow field, which can be scattered either all over the whole computational domain or over a smaller projection window. The resulting ROM is both computationally efficient(reconstructed flow fields require, in cases that do not present shock waves, less than 1 % of the time needed to compute a full CFD solution) and flexible(the projection window can avoid regions of large localized CFD errors).Also, for problems related with aerodynamics, POD modes are obtained from a set of snapshots calculated by a CFD method based on the compressible Navier Stokes equations and a turbulence model (which further more includes some unphysical stabilizing terms that are included for purely numerical reasons), but projection onto the POD manifold is made using the inviscid Euler equations, which makes the method independent of the CFD scheme. In addition, shock waves are treated specifically in the POD description, to avoid the need of using a too large number of snapshots. Various definitions of the residual are also discussed, along with the number and distribution of snapshots, the number of retained modes, and the effect of CFD errors. The method is checked and discussed on several test problems that describe (i) heat transfer in the recirculation region downstream of a backwards facing step, (ii) the flow past a two-dimensional airfoil in both the subsonic and transonic regimes, and (iii) the flow past a three-dimensional horizontal tail plane. The method is both efficient and numerically robust in the sense that the computational effort is quite small compared to CFD and results are both reasonably accurate and largely insensitive to the definition of the residual, to CFD errors, and to the CFD method itself, which may contain artificial stabilizing terms. Thus, the method is amenable for practical engineering applications. Resumen Se presenta un nuevo método para generar modelos de orden reducido (ROMs) aplicado a problemas fluidodinámicos de interés industrial. El nuevo método se basa en la expansión de las variables fluidas en una base POD, calculada a partir de un cierto número de snapshots, los cuales se han obtenido gracias a simulaciones numéricas (CFD). A continuación, las amplitudes de los modos POD se calculan minimizando un residual global adecuadamente definido que combina las ecuaciones y las condiciones de contorno. El método incluye varios ingredientes que son nuevos en este campo de estudio. El residual puede calcularse utilizando únicamente un número limitado de puntos del campo fluido. Estos puntos puede encontrarse dispersos a lo largo del dominio computacional completo o sobre una ventana de proyección. El modelo ROM obtenido es tanto computacionalmente eficiente (en aquellos casos que no presentan ondas de choque reconstruir los campos fluidos requiere menos del 1% del tiempo necesario para calcular una solución CFD) como flexible (la ventana de proyección puede escogerse de forma que evite contener regiones con errores en la solución CFD localizados y grandes). Además, en problemas aerodinámicos, los modos POD se obtienen de un conjunto de snapshots calculados utilizando un código CFD basado en la versión compresible de las ecuaciones de Navier Stokes y un modelo de turbulencia (el cual puede incluir algunos términos estabilizadores sin sentido físico que se añaden por razones puramente numéricas), aunque la proyección en la variedad POD se hace utilizando las ecuaciones de Euler, lo que hace al método independiente del esquema utilizado en el código CFD. Además, las ondas de choque se tratan específicamente en la descripción POD para evitar la necesidad de utilizar un número demasiado grande de snapshots. Varias definiciones del residual se discuten, así como el número y distribución de los snapshots,el número de modos retenidos y el efecto de los errores debidos al CFD. El método se comprueba y discute para varios problemas de evaluación que describen (i) la transferencia de calor en la región de recirculación aguas abajo de un escalón, (ii) el flujo alrededor de un perfil bidimensional en regímenes subsónico y transónico y (iii) el flujo alrededor de un estabilizador horizontal tridimensional. El método es tanto eficiente como numéricamente robusto en el sentido de que el esfuerzo computacional es muy pequeño comparado con el requerido por el CFD y los resultados son razonablemente precisos y muy insensibles a la definición del residual, los errores debidos al CFD y al método CFD en sí mismo, el cual puede contener términos estabilizadores artificiales. Por lo tanto, el método puede utilizarse en aplicaciones prácticas de ingeniería.

Thermo-Mechanical Treatment Effects on Stress Relaxation and Hydrogen Embrittlement of Cold-Drawn Eutectoid Steels

The effects of the temperature and stretching levels used in the stress-relieving treatment of cold-drawn eutectoid steel wires are evaluated with the aim of improving the stress relaxation behavior and the resistance to hydrogen embrittlement. Five industrial treatments are studied, combining three temperatures (330, 400, and 460 °C) and three stretching levels (38, 50 and 64% of the rupture load). The change of the residual stress produced by the treatments is taken into consideration to account for the results. Surface residual stresses allow us to explain the time to failure in standard hydrogen embrittlement tests

Energy-Entropy-Momentum integration of discrete thermo-visco-elastic dynamics.

A novel time integration scheme is presented for the numerical solution of the dynamics of discrete systems consisting of point masses and thermo-visco-elastic springs. Even considering fully coupled constitutive laws for the elements, the obtained solutions strictly preserve the two laws of thermo dynamics and the symmetries of the continuum evolution equations. Moreover, the unconditional control over the energy and the entropy growth have the effect of stabilizing the numerical solution, allowing the use of larger time steps than those suitable for comparable implicit algorithms. Proofs for these claims are provided in the article as well as numerical examples that illustrate the performance of the method.

Simulating the remote handling of the blanket segments in DEMO fusion reactor with thermo-mechanical meshfree multibody dynamics

The present study shows a first approach to the simulation of the remote handling oper- ation which takes into account the thermal and flexible behavior of the blanket segments and its implications on the remote handling equipment, in order to validate and improve its design.

Induced thermo-mechanical stress in CPV receivers with cycled high intensity light

CPV receivers are made of materials with very different lineal expansion coefficients. Strong variations in DNI due to the passage of clouds can cause sudden temperature changes that creates mechanical stress. For common solder and metal filled polymers the plastic limit could be reached causing substantial fatigue. The best forecast of receiver reliability is therefore achieved by applying an intermittent light source with nominal irradiance level and a number of cycles equal to the expected cloud passages for a given site. The UPM has developed specialized equipment, dubbed the LYSS (Light cYcling Stressing Source), for carrying out such experiments. The small thermal capacity of receivers allows simulating more than 25000 cycles per week. The number of deep transients expected for Madrid in 30 years operation, based on available data, is about 45000. We are currently using the system to cycle a ?Ge/Ag Epoxy/aluminum? receiver, which shows no degradation after 20000 cycles. The equipment can cast up to 200 and 70 W/cm2 on 0.1 and 1 cm2 cells, respectively.

Analysis of thermo-mechanical stress effects on CPV receiver under high intensity light cycles

After construction of the LYSS (Light cYcling Stressing Source) in early 2014, several CPV receivers, with and without secondary optical element (SOE) have been aged under fast transient illumination cycling,. The test plan for Madrid consisted of 50000 cycles. Receivers with poor heat spreaders showed low reliability but those with thicker metal layers passed the test well. The operation of LYSS along 8 months, after more than 250000 cycles, did not show any significant failure, except lamp reposition every 120 hours, in average. The equipment seems valid for unveiling weak receiver designs with respect to intensive illumination, in steady and transient modes.

A temperature-based thermodynamically consistent integration scheme for discrete thermo-elastodynamics

The formulation of thermodynamically consistent (TC) time integration methods was introduced by a general procedure based on the GENERIC form of the evolution equations for thermo-mechanical problems. The use of the entropy was reported to be the best choice for the thermodynamical variable to easily provide TC integrators. Also the employment of the internal energy was proved to not involve excessive complications. However, attempts towards the use of the temperature in the design of GENERIC-based TC schemes have so far been unfruitful. This paper complements the said procedure to attain TC integrators by presenting a TC scheme based on the temperature as thermodynamical state variable. As a result, the problems which arise due to the use of the entropy are overcome, mainly the definition of boundary conditions. What is more, the newly proposed method exhibits the general enhanced numerical stability and robustness properties of the entropy formulation.

A Fluctuation-Dissipation Model for Electrical Noise

This paper shows that today’s modelling of electrical noise as coming from noisy resistances is a non sense one contradicting their nature as systems bearing an electrical noise. We present a new model for electrical noise that including Johnson and Nyquist work also agrees with the Quantum Mechanical description of noisy systems done by Callen and Welton, where electrical energy fluctuates and is dissipated with time. By the two currents the Admittance function links in frequency domain with their common voltage, this new model shows the connection Cause-Effect that exists between Fluctuation and Dissipation of energy in time domain. In spite of its radical departure from today’s belief on electrical noise in resistors, this Complex model for electrical noise is obtained from Nyquist result by basic concepts of Circuit Theory and Thermo- dynamics that also apply to capacitors and inductors.

The role of spatial and temporal radiation deposition in inertial fusion chambers: the case of HiPER¿

The first wall armour for the reactor chamber of HiPER will have to face short energy pulses of 5 to 20 MJ mostly in the form of x-rays and charged particles at a repetition rate of 5–10 Hz. Armour material and chamber dimensions have to be chosen to avoid/minimize damage to the chamber, ensuring the proper functioning of the facility during its planned lifetime. The maximum energy fluence that the armour can withstand without risk of failure, is determined by temporal and spatial deposition of the radiation energy inside the material. In this paper, simulations on the thermal effect of the radiation–armour interaction are carried out with an increasing definition of the temporal and spatial deposition of energy to prove their influence on the final results. These calculations will lead us to present the first values of the thermo-mechanical behaviour of the tungsten armour designed for the HiPER project under a shock ignition target of 48 MJ. The results will show that only the crossing of the plasticity limit in the first few micrometres might be a threat after thousands of shots for the survivability of the armour.

Silica final lenses under HiPER laser fusion reactor operation conditions

We have studied the thermo-mechanical response and atomistic degradation of final lenses in HiPER project. Final silica lenses are squares of 75 × 75 cm2 with a thickness of 5 cm. There are two scenarios where lenses are located at 8 m from the centre: •HiPER 4a, bunches of 100 shots (maximum 5 DT shots <48 MJ at ≈0.1 Hz). No blanket in chamber geometry. •HiPER 4b, continuous mode with shots ≈50 MJ at 10 Hz to generate 0.5 GW. Liquid metal blanket in chamber design.

Influence of cement properties in the reaction rate and mechanical behavior of concrete with high fl y ash content

The use of fly ash (FA) as an admixture to concrete is broadly extended for two main reasons: the reduction of costs that supposes the substitution of cement and the micro structural changes motivated by the mineral admixture. Regarding this second point, there is a consensus that considers that the ash generates a more compact concrete and a reduction in the size of the pore. However, the measure in which this contributes to the pozzolanic activity or as filler is not well defined. There is also no justification to the influence of the physical parameters, fineness of the grain and free water, in its behavior. This work studies the use of FA as a partial substitute of the cement in concretes of different workability (dry and wet) and the influence in the reactivity of the ash. The concrete of dry consistency which serves as reference uses a cement dose of 250 Kg/m 3 and the concrete of fluid consistency utilized a dose of cement of 350 Kg/m 3 . Two trademark of Portland Cement Type 1 were used. The first reached the resistant class for its fineness of grain and the second one for its composition. Moreover, three doses of FA have been used, and the water/binder ratio was constant in all the mixtures. We have studied the mechanical properties and the micro-structure of the concretes by means of compressive strength tests, mercury intrusion porosimetry (MIP) and thermal analysis (TA). The results of compressive strength tests allow us to observe that concrete mixtures with cements of the same classification and similar dosage of binder do not present the same mechanical behavior. These results show that the effective water/binder ratio has a major role in the development of the mechanical properties of concrete. The study of different dosages using TA, thermo-gravimetry and differential thermal analysis, revealed that the portlandite content is not restrictive in any of the dosages studied. Again, this proves that the rheology of the material influences the reaction rate and content of hydrated cement products. We conclude that the available free water is determinant in the efficiency of pozzolanic reaction. It is so that in accordance to the availability of free water, the ashes can react as an active admixture or simply change the porous distribution. The MIP shows concretes that do not exhibit significant changes in their mechanical behavior, but have suffered significant variation in their porous structure

Lifetime of silica final lenses subject to HiPER irradiation conditions

The goal of the European laser fusion project, is to build an engineering facility for repetitive laser operation (HiPER 4a) and later a fusion reactor (HiPER 4b). A key aspect for laser fusion energy is the final optics. At the moment, it is based on silica transmission lenses located 8 m away from the chamber center. Lens lifetime depends on the irradiation conditions. We have used a 48 MJ shock ignition target for calculations. We have studied the thermo-mechanical effects of ions and X-rays on the lenses. Ions lead to lens melting and must therefore be mitigated. On the other hand, X-rays (~1% of the energy) does not produce either a significant temperature rise or detrimental stresses. Finally, we calculated the neutron flux and gamma dose rate on the lenses. Next, based on a simple model we studied the formation of color centers in the sample, which lead to optical absorption. Calculations show that simultaneous neutron and gamma irradiation does not significantly increase the optical absorption during the expected lifetime of the HiPER 4a facility. Under severe conditions (HiPER 4b), operation above 800 K or lens refreshing by thermal annealing treatments seem to assure adequate behavior.

A three-dimensional dynamical model of Amundsenisen icefield, Svalbard.

Amundsenisen is an ice field, 80 km2 in area, located in Southern Spitsbergen, Svalbard. Radio-echo sounding measurements at 20 MHz show high intensity returns from a nearly flat basal reflector at four zones, all of them with ice thickness larger than 500m. These reflections suggest possible subglacial lakes. To determine whether basal liquid water is compatible with current pressure and temperature conditions, we aim at applying a thermo mechanical model with a free boundary at the bed defined as solution of a Stefan problem for the interface ice-subglaciallake. The complexity of the problem suggests the use of a bi-dimensional model, but this requires that well-defined flowlines across the zones with suspected subglacial lakes are available. We define these flow lines from the solution of a three-dimensional dynamical model, and this is the main goal of the present contribution. We apply a three-dimensional full-Stokes model of glacier dynamics to Amundsenisen icefield. We are mostly interested in the plateau zone of the icefield, so we introduce artificial vertical boundaries at the heads of the main outlet glaciers draining Amundsenisen. At these boundaries we set velocity boundary conditions. Velocities near the centres of the heads of the outlets are known from experimental measurements. The velocities at depth are calculated according to a SIA velocity-depth profile, and those at the rest of the transverse section are computed following Nye’s (1952) model. We select as southeastern boundary of the model domain an ice divide, where we set boundary conditions of zero horizontal velocities and zero vertical shear stresses. The upper boundary is a traction-free boundary. For the basal boundary conditions, on the zones of suspected subglacial lakes we set free-slip boundary conditions, while for the rest of the basal boundary we use a friction law linking the sliding velocity to the basal shear stress,in such a way that, contrary to the shallow ice approximation, the basal shear stress is not equal to the basal driving stress but rather part of the solution.

A technical note on application of internally finned tubes in solar parabolic trough absorber pipes

The heterogeneous incoming heat flux in solar parabolic trough absorber tubes generates huge temperature difference in each pipe section. Helical internal fins can reduce this effect, homogenising the temperature profile and reducing thermal stress with the drawback of increasing pressure drop. Another effect is the decreasing of the outer surface temperature and thermal losses, improving the thermal efficiency of the collector. The application of internal finned tubes for the design of parabolic trough collectors is analysed with computational fluid dynamics tools. Our numerical approach has been qualified with the computational estimation of reported experimental data regarding phenomena involved in finned tube applications and solar irradiation of parabolic trough collector. The application of finned tubes to the design of parabolic trough collectors must take into account issues as the pressure losses, thermal losses and thermo-mechanical stress, and thermal fatigue. Our analysis shows an improvement potential in parabolic trough solar plants efficiency by the application of internal finned tubes.

Analysis of internal helically finned tubes for parabolic trough design by CFD tools

This paper has analysed the effect of the utilization of internal finned tubes for the design of parabolic trough collectors with computational fluid dynamics tools. Our numerical approach has been qualified with the computational estimation of reported experimental data regarding phenomena involved in finned tube applications and solar irradiation of parabolic trough collector. The application of finned tubes to the design of parabolic trough collectors must take into account features as the pressure losses, thermal losses and thermo-mechanical stress and thermal fatigue. Our analysis shows an improvement potential in parabolic trough solar plants efficiency by the application of internal finned tubes.