Evolution Of Temperature Profiles Research Articles

On learning latent dynamics of the AUG plasma state

In this work, we demonstrate the utility of state representation learning applied to modeling the time evolution of electron density and temperature profiles at ASDEX-Upgrade (AUG). The proposed model is a deep neural network, which learns to map the high dimensional profile observations to a lower dimensional state. The mapped states, alongside the original profile's corresponding machine parameters, are used to learn a forward model to propagate the state in time. We show that this approach is able to predict AUG discharges using only a selected set of machine parameters. The state is then further conditioned to encode information about the confinement regime, which yields a simple baseline linear classifier, while still retaining the information needed to predict the evolution of profiles. We, then, discuss the potential use cases and limitations of state representation learning algorithms applied to fusion devices.

Core transport modelling of the DTT full power scenario using different fuelling strategies

A theory-based integrated modelling work of plasma response to deuterium fuelling in the new Divertor Tokamak Test facility (DTT) is performed, using the 1.5D transport code JETTO with the quasi-linear anomalous transport model QuaLiKiz for the core region. The full power DTT scenario E1 is investigated. It is characterised by 28.8 MW of Electron Cyclotron Resonance Heating, 10 MW of Neutral Beam Injection and 6 MW of Ion Cyclotron Resonance Heating to the plasma. Plasma density and temperature profile evolution is calculated up to the separatrix using two different fuelling methods, gas puffing and pellet injection, and two different seeding gases, argon and neon. To sustain the desired pedestal density level with gas puffing a big amount of neutral flux at the separatrix is needed. The feasibility limits of the pumping system are exceeded, regardless of the type of impurity introduced, thus making the use of pellets mandatory. The simulations performed with pellet injection as fuelling method predict that the pedestal density is well sustained with realistic parameters foreseen for the DTT pellet injector. Strong dependence of the core density on the electron cyclotron resonance (ECR) power deposition profile is found. Trapped Electron Modes dominance, low outward flux and strongly hollow density in the inner core region are foreseen with central peaked ECR power deposition profile. Ion Temperature Gradient modes dominance, inward flux and robust density sustainment on the whole radial interval are predicted for spread ECR power deposition, though with central density close to the ECR cut-off limit and with peaked impurity densities. An intermediate deposition extension is found to sustain the whole density profile and to obtain flatter core densities, as previously predicted for the reference full power DTT scenario by fixed pedestal simulations. The ECR deposition is negligibly modified by refraction changes both during a single pellet cycle and after several pellet cycles, indicating full compatibility between the ECR system and the pellet injection system.

Model-based linear–quadratic–integral controller for simultaneous regulation of the current profile and normalized beta in NSTX-U

Achieving advanced scenarios that are characterized by steady-state operation, stable plasma confinement, and high-performance plasmas is one of the primary objectives of the National Spherical Tokamak eXperiment-Upgrade (NSTX-U). Active control of the plasma may be necessary to achieve these conditions. In particular, control algorithms that can simultaneously optimize the shapes and values of plasma profiles and scalars, respectively, may play a critical role in robustly achieving and sustaining these advanced scenarios. In this work, a model-based optimal control algorithm is developed for feedback control of the current profile in NSTX-U. The linear, finite-dimensional, control model is derived by discretizing and linearizing the magnetic diffusion equation in combination with empirical correlations for electron density, electron temperature, and noninductive current drives. The linear, time-variant model is then used for designing a linear–quadratic–integral (LQI) controller that is capable of regulating both the safety factor and the normalized beta around desired targets. The controller determines the neutral beam injection powers and the overall plasma current that are needed to achieve the desired current profile and normalized beta. The proposed controller is tested in higher-fidelity nonlinear simulations that employ 1D models for the evolutions of both current and temperature profiles using the Control Oriented Transport SIMulator (COTSIM). The closed-loop simulations show the effectiveness of the controller at shaping the safety factor in NSTX-U while achieving the desired normalized beta.

Electrolyte ohmic heating during electrodeposition: The role of coupled kinetic-transport phenomena that lead to morphological and microstructural changes

Understanding the morphological changes to the electrodes during electrochemical processes is important to many fields. These morphological changes often lead to dendrite formation and significant heat generation that can be detrimental to the process, especially the performance and lifetime of batteries. In situ infrared thermometry and microscopy are used to study the evolution of temperature profiles during the microscale electrodeposition of copper from copper sulfate. Under certain conditions, the growth front of the deposit receded and the fastest growing dendrites would suddenly stop growing. Analysis of physicochemical, temperature, and current changes show that many dynamic phenomena are correlated and associated with localized pH changes that change the chemistry of electrodeposition and potential drops occurring across the electrodes that alter the transport around dendrites.

Development of Thomson Scattering Measurement System for Long Duration Discharges on the QUEST Spherical Tokamak

The Thomson scattering control system has been modified to measure the time evolution of the electron density and temperature profiles during long duration discharges in the QUEST spherical tokamak. The system consists of a signal generator and a control circuit. The former accepts a QUEST main trigger and provides multiple triggers, each of which starts a short-term (e.g. 15 s) measurement. The latter provides triggers to synchronize the oscilloscope and laser oscillation during a short-term measurement. The system was used in 1000 s long duration discharges in QUEST, and the temporal evolutions of density and temperature profiles were obtained successfully. It was found the profiles are stationary after about 300 s.

The new sea ice thermodynamics code for the INM RAS Earth System model: The design and comparison of one- and zero-dimensional approaches with the observational data

Abstract This work is devoted to the comparison of one- (1-D) and zero-dimensional (0-D) models of sea ice thermodynamics. 1-D thermodynamics solvers imply the solution of the diffusion equation with penetrating radiation in the moving domain (moving boundary problem), while 0-D implementations neglect the heat capacity of ice and penetrating radiation, that leads to a linear temperature profile by the construction. So far, some climate models use 0-D implementation of thermodynamics in a block of sea ice, while other models use 1-D implementation with a small number of nodes. In this work we present our fully implicit 0-D and 1-D version of snow–ice thermodynamics, which relieves the model of Courant-type time step constraints. Implicitness is achieved by the relaxation method for nonlinear heat capacity and thermal conductivity coefficients with simultaneous search for the temperature that satisfies the boundary conditions. We introduce the experiment with analytical forcing to compare evolution of temperature profiles. This model experiment shows that 0-D model underestimates ice thickness and overestimates the amplitude of temperature profile variation compared to 1-D solution. Despite the large differences in thickness and internal temperatures, the surface temperatures are very close. The models were also validated on the data of the SHEBA field experiment. According to comparison with observation data, the 0-D model predicts ice thickness and spatiotemporal temperature distribution much worse compared to 1-D model, but the amplitude of error in surface temperature is slightly greater for 0-D code. It can be concluded that 0-D model of sea ice thermodynamics is a good tool for predicting surface temperature, but it gives a large error in thickness and temperature distribution.

Unsteady Pressure-Driven Electrokinetic Slip Flow and Heat Transfer of Power-Law Fluid through a Microannulus

To guarantee the transporting efficiency of microdevices associated with fluid transportation, mixing, or separation and to promote the heat transfer performance of heat exchangers in microelectronics, the hydrodynamic behaviors at unsteady and steady states, as well as the thermal characteristics at the steady state in a pressure-driven electrokinetic slip flow of power-law fluid in a microannulus are studied. To present a more reliable prediction, the slip phenomenon at walls and nonlinear rheology of liquid are incorporated. The modified Cauchy momentum equation applicable to all time scales and energy equations, are analytically solved in the limiting case of a Newtonian fluid and numerically solved for power-law fluids. The transient velocity profile, time evolution of flow rate, temperature profile, and heat transfer rate are computed at different flow behavior indices, electrokinetic width, slip lengths, and Brinkman numbers, thereby, the coupling effect of nonlinear rheology, slip hydrodynamics, and annular geometry on flow and thermal behaviors is explored. The unsteady flow takes a longer time to achieve the steady state for shear thinning fluids or greater slip lengths. The flow behavior index and slip length play a significant role in the flow rate and heat transfer performance. The relevant discussion can serve as a theoretical guide for the operation and thermal management of annular geometry-related flow actuation systems.

Economic optimization of enhanced geothermal systems using a fully analytical temperature profile

ABSTRACT In enhanced geothermal systems (EGS), the efficiency of the energy recovery process is highly dependent on the produced fluid’s temperature. In this work, a new analytical scheme is proposed for predicting the spatiotemporal evolution of temperature profile based on the three-dimensional energy conservation equation solved on the vertical plane passing through the wells. The solution is given in an error-function form which suits well for feasibility studies. An economic study is then conducted on the cost of operation based on the analytically described temperature and pressure distribution. The cumulative gain of operation, including the cost of pumping and reduced gain due to thermal breakthrough, is optimized with respect to the well spacing and the water injection rate. Optimum well spacing – injection rate combinations have been obtained for varied project life times. The maximum cumulative gain calculated by this method can be used for capex analyses and project lifetime predictions.

Flexocaloric effect in superelastic materials

We present a combined theoretical–experimental study of flexocaloric effects in superelastic materials exhibiting structural transitions. We study a Ginzburg–Landau model combined with a vibrational model for a beam near a ferroelastic transition loaded with a three-point bending setup. We also perform experiments on a Cu–Al–Ni single crystal undergoing a martensitic transition. We measure bent beam profiles, vertical force vs vertical deflection during a slow isothermal process, time evolution of the bending and unbending amplitudes, and the evolution of temperature profiles. We also compute the evolution of heat source and heat sink profiles. Finally, we study the location of acoustic emission events during the bending/unbending experiment. Our observations are consistent with the model predictions and allow us to identify the main physical parameters relevant for flexocaloric applications.

Simulations of Evolution of Electron Temperature Profiles During Electron Cyclotron Heating of High-Temperature Plasma in Gas-Dynamic Magnetic Trap

Simulations of Evolution of Electron Temperature Profiles During Electron Cyclotron Heating of High-Temperature Plasma in Gas-Dynamic Magnetic Trap

Transport simulations of pre-thermal quench shattered pellet injection in ITER: code verification and assessment of key trends

The material assimilation of the shattered pellet injection (SPI) in the pre-thermal quench phase in ITER is studied numerically by means of one-dimensional (1D) transport simulations. Such a simplified 1D approach is useful to perform extensive and systematic studies of key engineering parameters to optimize the ITER Disruption Mitigation System. The simulation results are compared with two-dimensional (2D) axisymmetric simulations by the nonlinear magnetohydrodynamic code JOREK for 5 neon/95 hydrogen SPI in the 15 MA hydrogen L-mode discharge to clarify the characteristics of SPI assimilation that can be analyzed within the range of the 1D model. Reasonable agreement between the 1D SPI simulation by the INDEX code and the 2D simulations by the JOREK code is found for total ablation rates, the radiation power, and the density and temperature profile evolution. The key process that was studied with the transport code is the onset of the radiative cold front that destabilizes the plasma current profiles. The injection parameters for the neon mixed hydrogen SPI are widely scanned to identify the timescale for the radiative cold front to form. Depending on the relative velocity of the cold front and the SPI fragments, the plasma cooling process can differ significantly. SPI with high injection velocities and large shard sizes results even in an inside-out plasma cooling that leads to hollow temperature profiles in the simulations reported here.

A high-power long lifetime beam dump for the Thomson scattering diagnostic system in the XuanLong-50 experiment.

The Energy iNNovation's XuanLong-50 is a spherical torus experiment with up to 10s plasma operation duration. A 3 J/50Hz pulsed laser is used in the Thomson scattering diagnostic system that is developed to measure the time evolutions of plasma electron temperature and density profiles. The expected laser pulse number is about 7.5 × 106/year with a power load of 150W. To meet at least 1-year lifetime requirement, a Chevron type beam dump with polished molybdenum plates is designed and fabricated, which absorbs the laser beam energy in a 3D structure to reduce the laser fluence deposited on the material surface. To prevent the backscattered stray light from interfering with the Thomson scattering measurements, a 7.5m beam path with folding mirrors is set between the beam dump and the plasma scattering volumes. Details of the beam dump design procedure including the laser beam profile control, multi-pulse laser damage threshold, heat dissipation, Zemax modeling, folding mirror selection, and beam path enclosure are presented together with the testing results.

The Heavy-Ion Beam Diagnostic of the ISTTOK Tokamak-Highlights and Recent Developments.

The unique arrangement of the heavy-ion beam diagnostic in ISTTOK enables one to measure the evolution of temperature, density and pressure-like profiles in normal and AC discharges. The fast chopping beam technique provided the possibility to reduce the noise on the measurements of the plasma pressure-like profile and for the precise control of the plasma column position in real time. The consequent improvements in S/N levels allowed the observation of the effects of runaway beam magnetic energy conversion into plasma local heating. In addition, it made it possible to follow the evolution of the quiescent plasma maintained during AC transitions when the plasma current is null. The use of a new operation mode in the cylindrical energy analyzer provided an improved resolution up to five times in determining the fluctuations of the plasma potential as compared to the normal operation mode. Such analyzer is extremely compact (250 mm × 250 mm × 120 mm) and provides a unique geometry in order to cover the whole plasma diameter. The detector configuration choice gives the possibility for the simultaneous measurements of plasma poloidal magnetic field, plasma pressure-like and plasma potential profiles together with their fluctuations.

Time evolution of temperature profiles retrieved from 13 years of infrared atmospheric sounding interferometer (IASI) data using an artificial neural network

Abstract. The three infrared atmospheric sounding interferometers (IASIs), launched in 2006, 2012, and 2018, are key instruments to weather forecasting, and most meteorological centres assimilate IASI nadir radiance data into atmospheric models to feed their forecasts. The European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) recently released a reprocessed homogeneous radiance record for the whole IASI observation period, from which 13 years (2008–2020) of temperature profiles can be obtained. In this work, atmospheric temperatures at different altitudes are retrieved from IASI radiances measured in the carbon dioxide absorption bands (654–800 and 2250–2400 cm−1) by selecting the channels that are the most sensitive to the temperature at different altitudes. We rely on an artificial neural network (ANN) to retrieve atmospheric temperatures from a selected set of IASI radiances. We trained the ANN with IASI radiances as input and the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis version 5 (ERA5) as output. The retrieved temperatures were validated with ERA5, with in situ radiosonde temperatures from the Analysed RadioSoundings Archive (ARSA) network and with EUMETSAT temperatures retrieved from IASI radiances using a different method. Between 750 and 7 hPa, where IASI is most sensitive to temperature, a good agreement is observed between the three datasets: the differences between IASI on one hand and ERA5, ARSA, or EUMETSAT on the other hand are usually less than 0.5 K at these altitudes. At 2 hPa, as the IASI sensitivity decreases, we found differences up to 2 K between IASI and the three validation datasets. We then computed atmospheric temperature linear trends from atmospheric temperatures between 750 and 2 hPa. We found that in the past 13 years, there is a general warming trend of the troposphere that is more important at the poles and at mid-latitudes (0.5 K/decade at mid-latitudes, 1 K/decade at the North Pole). The stratosphere is globally cooling on average, except at the South Pole as a result of the ozone layer recovery and a sudden stratospheric warming in 2019. The cooling is most pronounced in the equatorial upper stratosphere (−1 K/decade). This work shows that ANN can be a powerful and simple tool to retrieve IASI temperatures at different altitudes in the upper troposphere and in the stratosphere, allowing us to construct a homogeneous and consistent temperature data record adapted to trend analysis.

Unsteady ballistic heat transport in two-dimensional harmonic graphene lattice

We study the evolution of initial temperature profiles in a two-dimensional isolated harmonic graphene lattice. Two heat transfer problems are solved analytically and numerically. In the first problem, the evolution of a spatially sinusoidal initial temperature profile is considered. This profile is usually generated in real experiments based on the transient thermal grating technique. It is shown that at short times the amplitude of the profile decreases by an order magnitude and then it performs small decaying oscillations. A closed-form solution, describing the decay of the amplitude at short times is derived. It shows that due to symmetry of the lattice, the anisotropy of the ballistic heat transfer is negligible at short times, while at large times it is significant. In the second problem, a uniform spatial distribution of the initial temperature in a circle is specified. The profile is the simplest model of graphene heating by an ultrashort localized laser pulse. The corresponding solution has the symmetry of the lattice and many local maxima. Additionally, we show that each atom has two distinct temperatures corresponding to motions in zigzag and armchair directions. Presented results may serve for proper statement and interpretation of laboratory experiments and molecular dynamics simulations of unsteady heat transfer in graphene.

Thermographic in-situ investigation of precise electrochemical machining

Modelling of ECM is a powerful tool to improve the cost- and time intensive tool-development process. While different simulation models for DC-ECM have been published and validated, the virtual prediction of the process variant with oscillating cathode is still hindered by its highly transient nature. In order to develop and later validate a model of PECM, high-resolution information on the temperature profile and gas evolution during one oscillation is required. Temperature and gas concentration both increase during the electric pulses and constantly decrease due to flushing. This paper presents a fixture that allows the thermographic analyses of the PECM process. Consequently, the conditions within the electrolyte channel during PECM are captured which allows a deeper understanding of the process.

Experimental and numerical investigation of the response of GLT beams exposed to fire

This paper presents a combined experimental and numerical investigation of the behavior of glued laminated timber beams when exposed to fire. The influence on the time variation of charring rates based on the evolution of temperature profiles is examined for different fire scenarios and durations as well as different beam’s cross-section sizes. Predictions of charring depths provided by numerical simulations of heat transfer and simplified charring rate models are compared. In the absence of a mass transport representation, a Bayesian inference is introduced to identify the temperature-dependent material parameters for the conventional heat transfer model. A similar approach is adopted when adjusting the selected parameters of the charring rate models to account for variations in actual fire scenarios, which potentially depart from standard fire exposure. When compared to experimental results, both approaches confirmed their predictive capabilities, particularly in the stage of initial design. Since presented in the framework of Bayesian statistics, they open the door to fully stochastic analysis with an emphasis on the load bearing capacity of the studied beams.

Numerical study on the flow characteristics and smoke temperature evolution under double fires condition with a metro train in tunnel

Metro system as the most effective transportation method in urban area, greatly relieve the traffic pressure. The fire safety problem of metro system cannot be ignored as it is easy to develop into a large-area fire due to the relatively closed condition of the metro system. Double fire source is an important process in the development to large-area fire, however, there is still no relevant researches studying the flow characteristics and smoke movement under double fire sources with a metro train in the tunnel. This paper conducted a series of the full-scale numerical simulation to investigate the evolution of longitudinal temperature profile and maximum temperature rise under condition of double fire source with metro train in the tunnel. The velocity and temperature profile along the tunnel ceiling was measured and analyzed. It is found that (1) the longitudinal temperature distribution of tunnel fire smoke is still exponentially attenuated under different fire source separation distance with a metro train in the tunnel, and the decay coefficient decreases with the increase of fire source separation distance; (2) the maximum temperature rise of fire smoke under the tunnel ceiling decreases with the increase of the distance between the fire sources. Correlations were developed to describe the longitudinal temperature profile and maximum temperature rise under condition of double fire sources with a metro train in the tunnel through the analysis of the effect of additional air entrainment area, smoke interaction as well as the blockage effect of the metro train. The new findings in this paper can provide basic understanding of smoke movement under condition of double fire sources with a metro train in the tunnel, the proposed correlations can provide support for the safety operation of the metro system.

CFD Comparative Study Between Different Forms of Solar Greenhouses and Orientation Effect

Only scarce studies that were adopted have considered two properties, the structure safety and energy, where the aero-dynamic and energetic phenomena were taken into account simultaneously in the agricultural greenhouses area. In fact, in this numerical study, the response of the greenhouse has been investigated in outside climate conditions, by considering the orientation relatively to the wind direction velocity and solar trajectory. A resolution of the physical problem combined between the thermal and dynamical fluid flow equations have been based on the Ansys Fluent software. The results showed that the difference between inside and outside air temperature of greenhouse has been strongly affected by the reorientation of the tunnel greenhouse structure, or by the design of the tunnel structure that was adopted in the dome and chapel shape. Moreover, the safety properties of greenhouse structure linked to the drag stress can be developed when based on the interaction fluid-structure analysis. In this view, a temperature profile evolution versus different heights inside greenhouse was highlighted. As well as like continues of our previous study of the drag evolution over tunnel design body proved by the results found in the literature will be compared with chapel and dome designs.

W-band system-on-chip electron cyclotron emission imaging system on DIII-D.

Monolithic, millimeter-wave "system-on-chip" (SoC) technology has been employed in heterodyne receiver integrated circuit radiometers in a newly developed Electron Cyclotron Emission Imaging (ECEI) system on the DIII-D tokamak for 2D electron temperature profile and fluctuation evolution diagnostics. A prototype module operating in the E-band (72 GHz-80 GHz) was first employed in a 2 × 10 element array that demonstrated significant improvements over the previous quasi-optical Schottky diode mixer arrays during the 2018 operational campaign of the DIII-D tokamak. For compatibility with International Thermonuclear Experimental Reactor relevant scenarios on DIII-D, the SoC ECEI system was upgraded with 20 horn-waveguide receiver modules. Each individual module contains a University of California Davis designed W-band (75 GHz-110 GHz) receiver die that integrates a broadband low noise amplifier, a double balanced down-converting mixer, and a ×4 multiplier on the local oscillator (LO) chain. A ×2 multiplier and two IF amplifiers are packaged and selected to further boost the signal strength and downconvert the signal frequency. The upgraded W-band array exhibits >30 dB additional gain and 20× improvement in noise temperature compared with the previous Schottky diode radio frequency mixer input systems; an internal 8 times multiplier chain is used to bring down the LO frequency below 12 GHz, thereby obviating the need for a large aperture for quasi-optical LO coupling and replacing it with coaxial connectors. Horn-waveguide shielding housing avoids out-of-band noise interference on each individual module. The upgraded ECEI system plays an important role for absolute electron temperature evolution and fluctuation measurements for edge and core region transport physics studies.