Temperature-dependent Material Properties Research Articles

Temperature-dependent dynamics of electrokinetic conservative and reactive transport in porous media: A model-based analysis

Electrokinetic techniques employ direct current electric fields to enhance the transport of amendments in low permeability porous media and have been demonstrated effective for in situ remediation of both organic contaminants and heavy metals. The application of electric potential gradients give rise to coupled chemical, hydraulic and electric fluxes, which are at the basis of the main transport mechanisms: electromigration and electroosmosis. Previous research has highlighted the significant impacts of charge interactions and fluid composition, including temperature-dependent properties such as electrolyte conductivity and density, on these transport phenomena. However, current models of electrokinetic applications often assume isothermal conditions and overlook the production of heat resulting from Joule heating.This study provides a detailed model-based investigation, systematically exploring the effects of temperature on electrokinetic conservative and reactive transport in porous media. By incorporating temperature-dependent material properties and progressively investigating the impact of temperature on each transport mechanism, we analyze the effects of temperature variations in both 1D and 2D systems. The study reveals how temperature dynamically influences the physical, chemical and electrostatic processes controlling electrokinetic transport. A temperature increase results in a higher speed of amendments delivery by both electromigration and electroosmosis and increases the kinetics of degradation reactions. The simulations also reveal a feedback mechanism in which higher aqueous conductivity results in increased Joule heating, leading to a faster temperature rise and, subsequently, to higher electrolyte conductivity. Finally, we estimate the electric energy requirements of the system at varying temperatures and show how these changes impact the rate of contaminant removal.

Deconvolution of heat sources for application in thermoelectric micro four-point probe measurements

Joule heating is a primary phenomenon responsible for increasing temperature in electronic devices, and consequently, decreasing the lifetime and performance of electronic devices. However, this unwanted Joule heating can itself be used as a local source of heat to map the temperature-dependent material properties. Recently, micro four-point probe (M4PP) showed a promising potential for characterizing the temperature coefficient of resistance (TCR) and the Seebeck coefficient using Joule heating resulting from the measurement current. Here, we use M4PP to estimate the micrometer scale, relative temperature profile resulting from a single heat source. We introduce a triplet of four-point voltages measured at the second harmonic frequency, to deconvolute the thermoelectric voltage from the three individual heat sources involved. This paper tests and documents the validity of the proposed scheme in the 1–40 μm range on highly doped single crystal Si at 300 K, supporting predominantly Fourier heat transport at these scales. The method of deconvolution of heat sources reduces the complexity in evaluation of length- and time-dependent measurements, specifically used in the characterization of thermoelectric properties. The proposed method may also facilitate a more profound understanding of heat transport on the mesoscopic scale.

Thermoelastic topology optimization for structures with temperature-dependent material properties

Thermoelastic topology optimization for structures with temperature-dependent material properties

Thermal-Mechanical Effect and Removal Mechanism of Ti-6Al-4V During Laser-Assisted Grinding

The low density and high corrosion resistance of titanium alloy make it a material with various applications in the aerospace industry. However, because of its high specific strength and poor thermal conductivity, there are problems such as high cutting force, poor surface integrity, and high cutting temperature during conventional machining. As an advanced processing method with high efficiency and low damage, laser-assisted machining can improve the machinability of titanium alloy. In this study, a picosecond pulse laser-assisted scratching (PPLAS) method considering both the temperature-dependent material properties and ultrashort pulse laser’s characteristics is first proposed. Then, the effects of laser power, scratching depth, and scratching speed on the distribution of stress and temperature field are investigated by simulation. Next, PPLAS experiments are conducted to verify the correctness of the simulation and reveal the removal behavior at various combinations of laser power and scratching depths. Finally, combined with simulated and experimental results, the removal mechanism under the two machining methods is illustrated. Compared with conventional scratching (CS), the tangential grinding force is reduced by more than 60% and the material removal degree is up to 0.948 during PPLAS, while the material removal is still primarily in the form of plastic removal. Grinding debris in CS takes the form of stacked flakes with a “fish scale” surface, whereas it takes the form of broken serrations in PPLAS. This research can provide important guidance for titanium alloy grinding with high surface quality and low surface damage.

A multiphase modeling for investigating temperature history, flow field and solidification microstructure evolution of FeCoNiCrTi coating by laser cladding

In this work, a multiphase numerical model was established to investigate the non-isothermal flow and microstructure evolution during laser cladding. The model considered the beam-powder interaction, temperature dependent material properties, the effect of buoyancy and surface tension on convection. The accuracy of the model was verified comparing the simulated coating geometric morphology with the experimental results. Simulation indicated that the temperature and flow of the molten pool were interactive and mutually causal. On the one hand, the temperature field affected the tension gradient and the liquid flow direction. On the other hand, the velocity field significantly affected the heat transfer in the molten pool. Specifically, the maximum fluid velocity was 9.58 mm/s and the Peclet number (PeT) value was less than 5 at the initial stage, indicating that conduction dominated in heat transfer. Furthermore, the PeT value was greater than 200 when the maximum velocity reached 344 mm/s. The heat transfer in the molten pool was mainly thermal convection. In addition, the distribution diagram of solidification parameters (G and R) was calculated to quantitatively establish the internal correlation between solidification parameters and microstructure combing with experimental observation, and revealed the columnar-to-equiaxed transition (CET). Eventually, the microhardness results showed that HEA coating could improve the Ti6Al4V surface microhardness. Simultaneously, the CET behavior significantly affected the microhardness of the coating.

Thermoelastic free vibration analysis of functionally graded conical shell based on trigonometric higher-order shear deformation theory

The fundamental frequency of a rotating cantilevered porous functionally graded (FG) twisted conical shell with varying thickness along the longitudinal direction is calculated using a trigonometric higher-order shear deformation theory under thermal loading. Finite element method is employed for this purpose. The shell is discretized using eight-noded isoparametric shell elements with seven degrees of freedom per node. Using a simple power law across the transverse direction, the temperature-dependent material properties of the FG shell are determined. The nonlinear temperature distribution across the thickness direction is calculated using the one-dimensional Fourier heat conduction equation. The dynamic equation of motion is derived using Lagrange's equation. Finally, a parametric investigation of the effects of taper ratio, porosity, pretwist angle, temperature, and rotational speed on the fundamental frequency of the porous FG rotating conical shell is performed. It is also discussed how such characteristics affect mode shapes.

Investigation on temperature limitation and failure mechanism of SiC MOSFETs under avalanche conditions

To determine the latent degradation of SiC MOSFETs under single-pulse avalanche conditions, the stability of static characteristics and device capacitances is thoroughly studied by conducting Unclamped Inductive Switching (UIS) tests. No parametric shift is observed even though more than 90 % of the maximum avalanche energy is dissipated, indicating that devices fail in an abrupt mode. The avalanche failure turns out to be reaching the critical temperature. By introducing some corrections, including the temperature-dependent material properties, the aluminium layer, and the model of its melting process, the established model provides more accurate temperature results. The derived critical temperature is around 1200 K. All the devices suffer a fusion of metallization, which is considered to be the main mechanism of avalanche failure at the rated current region. With a sufficiently high current, the avalanche failure could also be accelerated by the turning-on of parasitic BJT.

Modelling and simulation of surface formation in electrical discharge machining based on thermo-hydraulic coupling

Surface formation in electrical discharge machining (EDM) involves melting, evaporation, re-solidification and crater overlapping. However, almost all existing models concerning surface formation in EDM only consider the thermal process and ignore hydrodynamic phenomena during machining. Herein, a novel three-dimensional thermo-hydraulic coupling model incorporating the arbitrary Lagrangian–Eulerian scheme was developed to address the above problem. The model comprehensively considered the surface deformation caused by melt flow, temperature-dependent material properties, phase changes and the effects of practical forces on melt pools to simulate surface formation on an anodic workpiece at a discharge current of 2 A. Simulation results demonstrate that during single-pulse discharge, molten material flowed radially to the crater periphery. The radial melt flow was driven by evaporation recoil pressure and thermocapillary force. The effect of evaporation recoil pressure on melt flow was gradually counteracted by that of surface tension after discharge ignition, and thermocapillary force became the dominant factor that maintained the existence of bulges around craters. During multiple discharges, the bulges generated by different discharges underwent superposition owing to the flow of molten material and subsequently formed surface peaks. The simulation results also provide the measurements of crater diameter, surface roughness and recast layer thickness, which are consistent with experimental results. The proposed modelling method can provide a useful reference for further investigations on surface integrity in EDM with minimal computing resource requirements.

Thermo-mechanical analysis of functionally graded nickel/titanium alloys under different process parameters of directed energy deposition

The build-up of residual stress and porosity plays an important role in the mechanical properties of as-built parts during directed energy deposition (DED). However, the thermo-mechanical behavior of functionally graded nickel/titanium alloys (FGNTAs) in the DED process remains unclear. In this study, a computationally efficient thermo-mechanical coupled model of DED-printed FGNTAs with different porosities is developed. The temperature-dependent effective material properties of FGNTAs are first determined based on the rule of mixture. Subsequently, the temperature, residual stress, and deformation fields within as-built FGNTAs are predicted using the differential quadrature method. The predicted results show excellent agreement with the predictions and experimental data in the literature. Furthermore, the effect of key parameters on the thermo-mechanical responses of DED-printed FGNTAs is investigated, including energy density, porosity distributions, part dimensions, temperature settings, and structural modes. The results show that the residual stress and vertical deflection increase with the higher energy input, length-to-width ratio, and evener porosity distributions, while decreasing with larger porosity, preheating as well environment temperatures, external loads, and power index. These results can provide a useful guideline for the design and manufacturing of high-quality FGNTAs during DED.

Investigation of statistical energy analysis coupling loss factors of vibro-acoustic systems under thermal environment

Vibro-acoustic systems in aviation engineering could operate in a thermal environment and suffer from high-frequency dynamic excitations. The high-frequency vibro-acoustic problems are commonly solved by the Statistical Energy Analysis method, in which the Coupling Loss Factor is the key parameter for modeling vibro-acoustic systems. This work investigates the CLFs of vibro-acoustic systems under thermal environment. First, based on the modal data considering different thermal effects, the modal CLFs between structural modes and acoustic modes are predicted by the Statistical modal Energy distribution Analysis. Then, the CLFs between the structural subsystem and the acoustic subsystem are calculated. The simulation models employed in this research include the coupled plate-cavity system and the cylindrical shell-cavity system. The influence of thermal effects on the CLFs is discussed. For the investigated cases, in which the boundaries of structural susbsystem are set to simply supported, the results demonstrate that the CLFs decrease with increasing temperature when only the temperature-dependent material properties are considered. However, when only the additional stress stiffness is considered, the CLFs experience a more dramatic decrease compared to when only the temperature-dependent material properties are taken into account. Moreover, when both effects are considered, the influence of the additional stress stiffness has a leading effect on the CLFs. The proposed method can be conveniently applied to various vibro-acoustic systems operating in the thermal environment.

Temperature-Dependent Vibration Characteristics of Porous FG Material Plates Utilizing FSDT

This study focuses on the vibration characteristics of the functionally graded materials (FGMs) porous plate. The plate is to be supported on different boundary constraints with linearly varying thicknesses. Existence of different porosity (void) pattern, within the materials, are taken into consideration using the power (P-FGM) and sigmoidal (S-FGM) gradation laws. The current methodology was developed utilizing the FSDT (first-order shear deformation theory) under thermal environmental conditions. Variations of temperature like as uniform, linearly varying, and nonlinear distribution patterns were examined by including temperature-dependent and independent material properties. The equations of motion including all the effects are derived from Hamilton’s principle and, subsequently solved using the Galerkin’s Vlasov method for various plate boundary conditions. Finally, the analytical outcomes are verified numerically, with the existing works. Furthermore, the study demonstrates that the fundamental frequency of porous FGM tapered plate is very close to the result obtained by the other researchers. Moreover, a detailed examination has been carried out to reveal the effect of various factors such as volume exponent index ([Formula: see text], side-to-thickness ratio ([Formula: see text]/[Formula: see text], and temperature effect ([Formula: see text]. In addition to this, some new benchmark results are obtained for free vibration analysis of tapered plates under a thermal environment.

Hygrothermally-Induced Vibration Analysis of Porous FGM Rectangular Mindlin Plates Resting on Elastic Foundation

The nonlinear vibration response of rectangular plates made of functionally graded porous materials (FGPs) induced by hygrothermal loading is investigated in this article using a numerical approach. The effect of elastic foundation on the vibrations is taken into account according to the Winkler–Pasternak model. Hygroscopic stresses produced due to nonlinear rise in moisture concentration are also considered. The temperature-dependent material properties of plate are computed based on the modified Voigt’s rule of mixture and Touloukian experiments for even and uneven distribution patterns of porosity. Within the framework of the first-order shear deformation plate theory and von-Kármán nonlinearity, Hamilton’s principle is utilized in order to derive the equations of motions. To achieve the temporal evolution of maximum lateral deflection of hygrothermally-induced plates, the generalized differential quadrature (GDQ) and Newmark integration methods are employed. Selected numerical results are presented to study the influences of temperature distribution, porosity volume fraction, moisture concentration, geometrical parameters, elastic foundation parameters and FG index on the geometrically nonlinear vibrations of FG porous plates with various boundary conditions.

MHz free electron laser x-ray diffraction and modeling of pulsed laser heated diamond anvil cell

A new diamond anvil cell experimental approach has been implemented at the European x-ray Free Electron Laser, combining pulsed laser heating with MHz x-ray diffraction. Here, we use this setup to determine liquidus temperatures under extreme conditions, based on the determination of time-resolved crystallization. The focus is on a Fe-Si-O ternary system, relevant for planetary cores. This time-resolved diagnostic is complemented by a finite-element model, reproducing temporal temperature profiles measured experimentally using streaked optical pyrometry. This model calculates the temperature and strain fields by including (i) pressure and temperature dependencies of material properties, and (ii) the heat-induced thermal stress, including feedback effect on material parameter variations. Making our model more realistic, these improvements are critical as they give 7000 K temperature differences compared to previous models. Laser intensities are determined by seeking minimal deviation between measured and modeled temperatures. Combining models and streak optical pyrometry data extends temperature determination below detection limit. The presented approach can be used to infer the liquidus temperature by the appearance of SiO2 diffraction spots. In addition, temperatures obtained by the model agree with crystallization temperatures reported for Fe–Si alloys. Our model reproduces the planetary relevant experimental conditions, providing temperature, pressure, and volume conditions. Those predictions are then used to determine liquidus temperatures at experimental timescales where chemical migration is limited. This synergy of novel time-resolved experiments and finite-element modeling pushes further the interpretation capabilities in diamond anvil cell experiments.

Extension of the zone method for the design of circular and tubular concrete columns subjected to a standard fire

PurposeThe paper aims to deal with the enhancement of a simplified method for the design of concrete columns subject to fire toward applications on circular and tubular cross-sections. The original zone method, developed by Hertz as a plastic design method, has been extended by Achenbach for the use as a nonlinear method. This proposed extended zone method (EZM) is verified by checking the theoretical background and is successfully validated by the recalculation of laboratory tests.Design/methodology/approachThe zone method assumes a reduction of a cross-section by a “damaged” zone. The remaining area is modeled with the constant, temperature-dependent material properties. The equations for the calculation of the damaged zone to model the loss of cross-section resistance or stiffness are derived. The proposed equations are validated by the recalculation of laboratory test and compared to the results of the advanced method (AM).FindingsIt can be shown that the EZM is suitable for the check of the fire resistance of circular concrete columns and leads to a safe and economic design. The method provides a suitable alternative to more sophisticated AM. The further extension toward tubular spun columns is discussed und is the object of the ongoing research.Originality/valuePresented enhancement extends the range of applications of the EZMs toward circular and tubular cross sections, which has previously not been examined.

Modulation of Turbulence Flux Budgets by Varying Fluid Properties in Heated High Prandtl Number Flow

The accurate computational modelling of momentum and heat transfer in turbulent wall-bounded flow is still a challenging task, especially for fluids with high molecular Prandtl numbers. Due to the steep temperature gradients, the strong variation of fluid properties becomes an additional important issue. The present Direct Numerical Simulation study specially addresses this problem considering a representative working fluid in terms of a real heat transfer oil with temperature-dependent material properties. Near the heated wall, the increase in molecular viscosity significantly attenuates the turbulent convective fluxes of both momentum and heat, despite the markedly increased enthalpy fluctuations. As the thermal conductivity decreases markedly less with temperature than the viscosity, the dampened turbulent mixing significantly reduces the Nusselt number, whereas the skin friction coefficient only marginally differs from the constant fluid property case. The attenuation of the radial velocity fluctuation consistently reduces the production terms in the budgets of the turbulent momentum and heat fluxes, while the production of the enthalpy fluctuations is only shifted away from the wall without significant change in peak level. The also attenuated redistribution in the turbulent axial stress budget reduces the pressure-strain based steering of turbulence towards isotropy. Facing a consequently more anisotropic turbulent motion poses a further difficulty to standard turbulence modelling mostly developed for isotropic turbulence.

Nonlinear vibration characteristics of elastically supported FRP cylindrical shells under temperature gradient conditions: Theoretical and experimental studies

The nonlinear free vibrations of elastically supported fiber-reinforced polymer cylindrical shells (FRPCSs) under temperature gradient conditions are investigated using both theoretical and tested techniques. A dynamic model of the FRPCSs with consideration of such complex thermal conditions is firstly established based on the first-order shear deformation theory in conjunction with Hamilton’s principle, the Galerkin method, the artificial spring technique, etc. Numerical results at room temperature and in a uniform thermal environment with classical or arbitrary boundary conditions are utilized to give a rough validation of this model. Subsequently, a thermal-vibration system is set up to measure the nonlinear natural frequencies of two FRPCS specimens subjected to four different temperature gradient conditions in free-free boundary constraints, and iterative calculations are performed to identify the temperature-dependent material properties. Finally, the detailed comparisons of calculated and experimental results provide a solid validation of the proposed model. This study offers a practical model tool to forecast the nonlinear free vibrations of the FRPCSs in a complex thermal environment, which can be readily adjusted and extended to other forms of composite shells. Also, the predicted and measured results can help assess the structural thermal-vibration behaviors when the temperature gradient effect needs to be considered. Communicated by Makoto Ohsaki.

On quantifying uncertainty in lightning strike damage of composite laminates: A hybrid stochastic framework of coupled transient thermal-electrical simulations

Lightning strike damage can severely affect the thermo-mechanical performance of composite laminates. It is essential to quantify the effect of lightning strikes considering the inevitable influence of material and geometric uncertainties for ensuring the operational safety of aircraft. This paper presents an efficient support vector machine (SVM)-based surrogate approach coupled with computationally intensive transient thermal-electrical finite element simulations to quantify the uncertainty in lightning strike damage. The uncertainty in epoxy matrix thermal damage and electrical responses of unprotected carbon/epoxy composite laminates is probabilistically quantified considering the stochasticity in temperature-dependent multi-physical material properties and ply orientations. Further, the SVM models are exploited for variance-based global sensitivity analysis to investigate the input parameters' relative influence on the lightning strike-induced damage behavior. Due to the adoption of a coupled SVM-based simulation approach here, it has become possible to carry out a comprehensive uncertainty quantification leading to complete probabilistic descriptions of the electrical and lightning damage parameters despite the requirement of performing a large number of computationally intensive function evaluations. The results reveal that source-uncertainty of the unprotected laminates significantly influences the epoxy matrix decomposition, electrical current density and electric potential, wherein longitudinal electrical conductivity is most sensitive to stochastic variations followed by other electrical, thermal and geometric parameters.

(Invited) Bubble Trouble: Simulation of Molten Oxide Electrolysis in Reduced Gravity

Molten Oxide Electrolysis (MOE) is a leading contender for processing extra-terrestrial minerals because it produces oxygen and liquid metal without the use of consumables. During electrolysis a gravitational field stratifies gas, electrolyte, and metal from each other by density. Processing on the moon, with ~1/6g, or Mars, with ~1/3g, will reduce the buoyancy force slowing bubble velocity and lowering convection intensity. Low thermal conductivity, (~1Wm-1K-1) and high viscosity, (~1Pa⋅s) of molten regolith contribute to thermal and chemical transport being dominated by convection, even as convection velocity decreases. We explore MOE system dynamics in Earth’s (1g), Mars’ (1/3g), and the Moon’s (1/6g) gravity through simulations that include composition and temperature dependent material properties and two phase flow. We implement the level set method to track the interface between the bubble laden flow near the anode and the bubble free flow far from the anode. By grouping bubble dynamics into the material properties of the fluid near the anode we relax the requirements for a fine mesh to track individual bubble-fluid interfaces, and are able to use a mesh that is finer than that required for a bubbly flow mixture model. Thus, all physics can be captured on one mesh. This is particularly relevant because bubble sizes for proposed systems are 1/10 to 1/100 of the system scale. While other simulations have been presented and experiments performed in micro and zero gravity for gas evolving electrolysis systems, none have investigated the downward facing electrode, heat transfer, or species diffusion. The focus of this work is to address this gap in the literature through simulation. Here, multiphysics simulation is supported by inspection of dimensionless groups and observations from other electrolysis systems. This work demonstrates that combining heat and mass transfer with temperature and composition dependent material properties is critical to designing MOE and other electrolysis systems for earth and beyond.

Thermally nonlinear analysis of propagating cracks under generalized thermal shock

In this paper, the thermally nonlinear response of stationary and moving cracks in isotropic media subjected to a hyperbolic thermal shock in the framework of the eXtended Finite Element Method (XFEM) is examined. For this purpose, the thermally nonlinear fully coupled Lord-Shulman (LS) model is considered as governing equations. In addition, two cases of material properties, temperature-dependent (TD) and temperature-independent (TID), are assumed. The geometry, including a predefined crack, is modeled using the XFEM for spatial discretization. For modeling the dynamic crack growth within the XFEM, the Maximum Tangential Stress (MTS) criterion is modified based on the stress field near the tip of a growing crack to calculate the dynamic crack growth angle. Besides, the nonlinear form of the Newmark integration scheme is implemented for the numerical integration in the time domain. Finally, a method based on one of the Barnet-Lothe tensors for isotropic materials is employed to extract the stress intensity factors (SIFs). In the provided examples, the effects of thermal nonlinearity, the temperature-dependency of material properties, and the LS relaxation time on the crack propagation path are explored.

Multiscale thermomechanical assessment of silicon carbide-based nanocomposites in solar energy harvesting applications

Increasing the operating temperature of concentrating solar power plants beyond 650 °C is crucial for increasing their efficiency and energy storage capability, and thus pushing their market penetration. However, one of the main drawbacks that arise with an increase in temperature is material degradation. Ceramic materials are an interesting alternative to nickel–chromium superalloys due to their high thermal resistance and low expansion capabilities. In this work, we characterize the thermomechanical behavior of a silicon carbide-based ceramic material for receiver tubes up to 750 °C. For this purpose, we develop a multiscale simulation framework, combining studies at different scales of observation (i.e., material, and component) using different simulation tools. At the material scale, we generate virtual microstructures that are subject to thermomechanical tests using the lattice model. The material properties obtained at this scale of observation are subsequently used at component level simulations, which are carried out by means of the finite element method. This allows the exploration of new materials designed in silico and the prediction of their performance in solar energy applications. In this work, we use our methodology to assess the thermoelastic behavior of SiSiC tubes at high temperatures and obtain temperature-dependent material properties. At the material level, we report elastic moduli in the range of 250–500 MPa and coefficients of thermal expansion in the range of 6–12·10−6 1/°C, under operation conditions. Regarding the component, we found that SiSiC tubes with a thickness-to-diameter ratio of 0.15 or larger are safe at operating temperatures up to 750 °C.