GURTIN-TYPE VARIATIONAL PRINCIPLES FOR DYNAMICS OF A NON-LOCAL THERMAL EQUILIBRIUM SATURATED POROUS MEDIUM

Based on the porous media theory and by taking into account the efects of the pore fuid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fuid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fuid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small de- formation of the solid phase, small velocity of the fuid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fuid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles, especially Hu-Washizu type variational principles, for the initial boundary value problems of dy- namic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fuid-saturated elastic porous media, which have been discussed previously.

NONLINEAR J-E CHARACTERISTICS IN THE ELECTRIC-THERMAL EQUILIBRIUM STATE FOR HIGH DENSITY POLYETHYLENE CONDUCTIVE COMPOSITES

The nonlinear J-E characteristics under self-heating equilibrium for conductive composites based on high density polyethylene were studied. The results show that there are identical conduction mechanisms under self-heating equilibrium for the composites with various initial resistivities determined by filler content or ambient temperature. The nonlinear conduction behavior was involved in the limited microstructure transformations of the conducting network induced by electrical field applied and the corresponding self-heating effect. A reversible thermal fuse (RTF) model was suggested to interpret the physical origin of the nonlinear J-E characteristics.

Sensitivity Analysis of Thermal Equilibrium Parameters of MIKE 11 Model:A Case Study of Wuxikou Reservoir in Jiangxi Province of China

Sensitivity analysis of thermal equilibrium parameters in the reservoir module of MIKE 11 model was conducted for the Wuxikou Reservoir in Jiangxi Province of China in order to apply the module to the environmental impact assessment to accurately predict water temperature of reservoirs.Results showed that radiation parameter A and evaporation-first parameter were much more sensitive than other parameters.The values of the radiation parameter A ranged from 0.10 to 0.34.The values of evaporation-first parameter varied from 0 to 10.The sensitivity of solar absorption parameters was less than that of evaporation parameter,of which light attenuation values ranged from 0.5 to 0.7,and this parameter would not impact model results if it was more than 2.Constants in Beer's law ranged from 0.2 to 0.7.Radiation parameter B was not more sensitive than evaporation parameter and its reasonable range was higher than 0.48.The fitting curves showed consistent changing tendency for these parameters within the reasonable ranges.Additionally,all the thermal equilibrium parameters had much more important effects on surface water temperature than deep water temperature.Moreover,if no observed data could be obtained,the local empirical value would be used to input to the MIKE 11 model to simulate the changes in the discharged outflow-water temperature qualitatively.

An Extension of One-Period Nash Equilibrium Model in Non-Life Insurance Markets

This paper deals with an extension of the one-period model in non-life insurance markets (cf. [1]) by using a transition probability matrix depending on some economic factors. We introduce a multi-period model and in each period the solvency constraints will be updated. Moreover, the model has the inactive state including some uninsured population. Similar results on the existence of premium equilibrium and sensitivity analysis for this model are presented and illustrated by numerical results.

A non-equilibrium molecular dynamics study of the thermal conductivity of uranium dioxide

The thermal conductivity of uranium dioxide in the temperature range of 300–2400 K was estimated by non-equilibrium molecular dynamics calculation using Fourier law.The Kawamura function was adopted as the interatomic potential function.The calculated thermal conductivities are found to be strongly dependent on the temperature of the simulation cube.The thermal conductivity simulation results are compared with the experiment results and agreed well with the experimental results when the temperature is above 600 K.The thermal conductivities scale effects are found to be existed in uranium dioxide nanometer thin film.The approximate mean free paths of phonons in different temperatures can be examined.

An Improved Three-Dimensional Non-Equilibrium Mixing Pool Model

This paper presents an improved three-dimensional non-equilibrium mixing pool model.It is a simplified form of the original model and is more practical for applications.The simulation re-sults show that the industrial scale distillation tray columns can be described closely by the improvedmodel.The effects of model parameters,such as the number of mixing pools,the point efficiencyand flow pattern,on separation are analyzed quantitatively.

Three-Temperature Model Applied to Thermochemical Non-Equilibrium Reentry Flows in 2D—Seven Species

In this work, a study involving the fully coupled Euler and Navier-Stokes reactive equations is performed. These equations, in conservative and finite volume contexts, employing structured spatial discretization, on a condition of thermochemical non-equilibrium, are analyzed. High-order studies are accomplished using the Spectral method of Streett, Zang, and Hussaini. The high enthalpy hypersonic flows around a circumference, around a reentry capsule, along a blunt body, and along a double ellipse in two-dimensions are simulated. The Van Leer, Liou and Steffen Jr., and Steger and Warming flux vector splitting algorithms are applied to execute the numerical experiments. Three temperatures, which are the translational-rotational temperature, the vibrational temperature, and the electron temperature, are used to accomplish the numerical comparisons. Excellent results were obtained with minimum errors inferior to 6.0%. The key contribution of this work is the correct implementation of a three temperature model coupled with the implementation of three algorithms to perform the numerical simulations, as well the description of energy exchange mechanisms to perform more realistic simulations.

Thermal quantum discord in Heisenberg models with Dzyaloshinski-Moriya interaction

We study the quantum discord of the bipartite Heisenberg model with the Dzyaloshinski-Moriya （DM） interaction in thermal equilibrium state and discuss the effect of the DM interaction on the quantum discord. The quantum entanglement of the system is also discussed and compared with quantum discord. Our results show that the quantum discord may reveal more properties of the system than quantum entanglement and the DM interaction may play an important role in the Heisenberg model.

Non-equilibrium ignition criterion for magnetized deuterium–tritium fuel

In this paper, non-equilibrium ignition conditions for magnetized cylindrical deuterium–tritium plasma in the presence of an axial magnetic field have been investigated. It is expected that temperature imbalance between ions and electrons as well as the axial magnetic field will relax the threshold of ignition conditions.Therefore, ignition conditions for this model are derived numerically involving the energy balance equation at the stagnation point. It has been derived using parametric space including electron and ion temperature(T_e, T_i), areal density(q R), and seed magnetic field-dependent free parameters of B/q, mB, and BR. For B/ρ < 10~6 G cm^3 g^(-1),mB < 4 × 10~4 G cm g^(-1), and BR <3 × 10~5 G cm, the minimum fuel areal density exceeds between ρR >0.002 g cm^(-2), ρR> 0.25 g cm^(-2), and ρR > 0.02 g cm^(-2),respectively. The practical equilibrium conditions also addressed which is in good agreement with the corresponding one-temperature magnetized mode proposed in previous studies. Moreover, it has been shown that the typical criterion of BR ≥(6.13–4.64) × 10~5 G cm would be expectable. It is also confirmed that the minimum product of areal density times fuel temperature in equilibrium model is located in the range of T = 6–8 keV for all these free parameters, depending on the magnitude of the magnetic field. This is the entry point for the non-equilibrium model consistent with equilibrium model.

The Global Stability Analysis of a Mathematical Cellular Model of Hepatitis C Virus Infection with Non-Cytolytic Process

The aim of this work is to analyse the global dynamics of an extended mathematical model of Hepatitis C virus (HCV) infection in vivo with cellular proliferation, spontaneous cure and hepatocyte homeostasis. We firstly prove the existence of local and global solutions of the model and establish some properties of this solution as positivity and asymptotic behaviour. Secondly we show, by the construction of appropriate Lyapunov functions, that the uninfected equilibrium and the unique infected equilibrium of the mathematical model of HCV are globally asymptotically stable respectively when the threshold number and when .

Mathematical Model and Non-Pharmaceutical Control of the Coronavirus 2019 Disease in Madagascar

For Madagascar, with the uncertainty over vaccines against the novel coronavirus 2019 and its variants, non-pharmaceutical approach is widely used. Our objective is to propose a mathematical control model which will serve as a tool to help decision-makers in the strategy to be implemented to better face the pandemic. By separating asymptomatic cases which are often not reported and symptomatic who are hospitalized after tests;we develop a mathematical model of the propagation of covid-19 in Madagascar, by integrating control strategies. We study the stability of the model by expressing the basic reproduction number using the next-generation matrix. Simulation with different parameters shows the effects of non-pharmaceutical measures on the speed of the disease spread. By integrating a control parameter linked to compliance with barrier measures in the virus propagation equation, we were able to show the impacts of the implementation of social distancing measures on the basic reproduction number. The strict application of social distancing measures and total confinement is unfavorable for economic situation even if they allow the contamination to be reduced quickly. Without any restrictions, the disease spreads at high speed and the peak is reached fairly quickly. In this condition, hospitals are overwhelmed and the death rate increases rapidly. With 50% respect for non-pharmaceutical strategies such as rapid detection and isolation of positive cases and barrier gestures;the basic reproduction number R0 can go down from 3 to 1.7. The pressures on the economic and social situation are rather viable. It is the most suitable for the Malagasy health system. The results proposed are a way to control the spread of the disease and limit its devastation in a country like Madagascar.

焊接热循环下Nb原子非平衡晶界偏聚行为的动力学建模

采用Gleeble-3500型热模拟机,对3种不同Nb含量(0.180%、0.085%、0.025%)的HSLA钢在不同冷却速率下的焊接热循环过程进行模拟,借助电子显微探针(EPMA)技术对Nb原子的非平衡晶界偏聚行为进行表征,得到了不同Nb含量和冷却速率下晶界处Nb原子偏聚量,基于非平衡晶界偏聚动力学理论,计算得到空位形成能为1.403 eV,同时确定了Nb含量分别为0.180%和0.085%钢中Nb原子非平衡偏聚的临界时间常数分别为7.59×10^(3)和6.0×10^(3)。进一步将Nb原子-空位复合体的扩散系数作为变量,计算并修正得到了其与等效时间的函数关系,最终建立了焊接热循环条件下Nb原子非平衡偏聚的修正动力学方程。

主动式气膜冷却对高超声速飞行器等离子体鞘套的影响

以钝锥模型为研究对象,采用热化学非平衡计算方法并结合可压缩N-S方程、SST k-ω湍流模型,对钝锥模型肩部台阶喷流和头部逆向喷流流场进行了数值模拟。在考虑不同喷流压强、喷口位置和数量等因素的基础上,分析了主动式气膜冷却对等离子体鞘套的影响。结果表明:高超声速飞行器采用主动式气膜冷却技术时,喷口的数量、位置及喷流压强对等离子体密度均具有显著的影响。肩部的切向喷流可有效抑制模型壁面附近的等离子体密度,进而可能对高频电磁波的传输和目标雷达散射截面(RCS)产生影响。头部逆向喷流可显著改变等离子体的分布情况,不同的逆向喷流参数配置会导致明显的差异。

Thermal Analysis of Vehicular Twin-Tube Hydraulic Gas-Precharged Shock Absorbers

In this study of temperature rising in vehicular twin-tube hydraulic gas-precharged shock absorbers,thermodynamic analyses were conducted via simulations.Equations on heat conduction,heat convection as well as radiation were derived by applying certain laws governing heat transfer;an equivalent thermal resistance network model of a shock absorber undergoing heat transfer was established innovatively;moreover,the shock absorber’s thermodynamic model of control volume system was built by using the first law of thermodynamics;and finally,time required for shock absorber to reach thermal equilibrium and corresponding value of steady temperature were calculated by programming.In this way,a lower thermal equilibrium temperature will be achieved,hence help to improve reliability of shock absorbers in work by offering low ambient temperature,by reducing amplitudes and frequencies of external incentives exerted on them and by increasing flow rate of ambient air passing around them.

Two-Dimensional Mathematical Model of Tidal Current and Sediment for Oujiang Estuary and Wenzhou Bay

A 2-D mathematical model of tidal current and sediment has been developed for the Oujiang Estuary and the Wenzhou Bay. This model accomodates complicated features including multiple islands, existence of turbidity, and significant differ-ence in size distribution of bed material. The governing equations for non-uniform suspended load and bed load transport are presented in a boundary-fitted orthogonal curvilinear coordinate system. The numerical solution procedures along with their initial conditions, boundary conditions, and movable boundary technique are presented. Strategies for computation of the critical condition of deposition or erosion, sediment transport capacity, non-uniform bed load discharge, etc. are suggested. The model verification computation shows that, the tidal levels computed from the model are in good agreement with the field data at the 18 tidal gauge stations. The computed velocities and flow directions also agree well with the values measured along the totally 52 synchronously observed verticals distributed over 8 cross sections. The coraputed tidal water throughputs through the Huangda'ao cross section are close to the measured data. And the computed values of bed deformation from Yangfushan to the estuary outfall and in the outer-sea area are in good agreement with the data observed from 1986 to 1992. The changes of tidal volumes through the estuary, velocities in different channels and the bed form due to the influence of the reclamation project on the Wenzhou shoal are predicted by means of this model.

Physics-based analysis and simulation model of electromagnetic interference induced soft logic upset in CMOS inverter

The instantaneous reversible soft logic upset induced by the electromagnetic interference（EMI） severely affects the performances and reliabilities of complementary metal–oxide–semiconductor（CMOS） inverters. This kind of soft logic upset is investigated in theory and simulation. Physics-based analysis is performed, and the result shows that the upset is caused by the non-equilibrium carrier accumulation in channels, which can ultimately lead to an abnormal turn-on of specific metal–oxide–semiconductor field-effect transistor（MOSFET） in CMOS inverter. Then a soft logic upset simulation model is introduced. Using this model, analysis of upset characteristic reveals an increasing susceptibility under higher injection powers, which accords well with experimental results, and the influences of EMI frequency and device size are studied respectively using the same model. The research indicates that in a range from L waveband to C waveband, lower interference frequency and smaller device size are more likely to be affected by the soft logic upset.

Direct modeling for computational fluid dynamics

All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes （NS） equations. The current computational fluid dynamics （CFD） focuses on the numerical solution of partial differential equations （PDEs）, and its aim is to get the accurate solution of these governing equations. Under such a CFD practice, it is hard to develop a unified scheme that covers flow physics from kinetic to hydrodynamic scales continuously because there is no such governing equation which could make a smooth transition from the Boltzmann to the NS modeling. The study of fluid dynamics needs to go beyond the traditional numer- ical partial differential equations. The emerging engineering applications, such as air-vehicle design for near-space flight and flow and heat transfer in micro-devices, do require fur- ther expansion of the concept of gas dynamics to a larger domain of physical reality, rather than the traditional dis- tinguishable governing equations. At the current stage, the non-equilibrium flow physics has not yet been well explored or clearly understood due to the lack of appropriate tools. Unfortunately, under the current numerical PDE approach, it is hard to develop such a meaningful tool due to the absence of valid PDEs. In order to construct multiscale and multiphysics simulation methods similar to the modeling process of con- structing the Boltzmann or the NS governing equations, the development of a numerical algorithm should be based on the first principle of physical modeling. In this paper, instead of following the traditional numerical PDE path, we introduce direct modeling as a principle for CFD algorithm develop- ment. Since all computations are conducted in a discretized space with limited cell resolution, the flow physics to be mod- eled has to be done in the mesh size and time step scales. Here, the CFD is more or less a direct construction of dis- crete numerical evolution equations, where the mesh size and time step will play dynamic roles in the modeling process. With the variation of the ratio between mesh size and local particle mean free path, the scheme will capture flow physics from the kinetic particle transport and collision to the hydro- dynamic wave propagation. Based on the direct modeling, a continuous dynamics of flow motion will be captured in the unified gas-kinetic scheme. This scheme can be faithfully used to study the unexplored non-equilibrium flow physics in the transition regime.

Modelling of Electrical Conductivity of a Silver Plasma at Low Temperature

During the working of electrical fuses, inside the fuse element the silver ribbon first begins to melt, to vaporize and then a fuse arc appears between the two separated parts of the element. Second, the electrodes are struck and the burn-back phenomenon takes place. Usually, the silver ribbon is enclosed inside a cavity filled with silica sand. During the vaporization of the fuse element, one can consider that the volume is fixed so that the pressure increase appears to reach pressures higher than atmospheric pressure. Thus, in this paper two pressures, 1 atm and 10 atm, are considered. The electrical field inside the plasma can reach high values since the distance between the cathode surface and the anode surface varies with time. That is to say from zero cm to one cm order. So we consider various electrical fields： 102 V/m, 103 V/m, 5×103 V/m, 104 V/m at atmospheric pressure and 105 V/m at a pressure of 10 atm. This study is made in heavy species temperature range from 2,400 K to 10,000 K. To study the plasma created inside the electric fuse, we first need to determine some characteristics in order to justify some hypotheses. That is to say： are the classical approximations of the thermal plasmas physics justified？ In other words： plasma frequency, the ideality of the plasma, the Debye-Hückel approximation and the drift velocity versus thermal velocity. These characteristics and assumptions are discussed and commented on in this paper. Then, an evaluation of non-thermal equilibrium versus considered electrical fields is given. Finally, considering the high mobility of electrons, we evaluate the electrical conductivities.

Analysis and FDTD Modeling of the Influences of Microwave Electromagnetic Waves on Human Biological Systems

The interactions of electromagnetic waves with the human body are complex and depend on several factors related to the characteristics of the incident wave, including its frequency, its intensity, the polarization of the tissue encountered, the geometry of the tissue and its electromagnetic properties. That’s to say, the dielectric permittivity, the conductivity and the type of coupling between the field and the exposed body. A biological system irradiated by an electromagnetic wave is traversed by induced currents of non-negligible density;the water molecules present in the biological tissues exposed to the electromagnetic field will begin to oscillate at the frequency of the incident wave, thus creating internal friction responsible for the heating of the irradiated tissues. This heating will be all the more important as the tissues are rich in water. This article presents the establishment from a mathematical and numerical analysis explaining the phenomena of interaction and consequences between electromagnetic waves and health. Since the total electric field in the biological system is unknown, that is why it can be determined by the Finite Difference Time Domain FDTD method to assess the electromagnetic power distribution in the biological system under study. For this purpose, the detailed on the mechanisms of interaction of microwave electromagnetic waves with the human body have been presented. Mathematical analysis using Maxwell’s equations as well as bio-heat equations is the basis of this study for a consistent result. Therefore, a thermal model of biological tissues based on an electrical analogy has been developed. By the principle of duality, an electrical model in the dielectric form of a multilayered human tissue was used in order to obtain a corresponding thermal model. This thermal model made it possible to evaluate the temperature profile of biological tissues during exposure to electromagnetic waves. The simulation results obtained from computer tools show that the temperature in the biological tissue is a linear function of the duration of exposure to microwave electromagnetic waves.

Radiation mean free path for non-LTE fully ionized medium

１ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｒａｄｉａｔｉｏｎｍｅａｎｆｒｅｐａｔｈ（ＲＭＦＰ）［１，２］ｉｓａｖｅｒｙｉｍｐｏｒｔａｎｔｐｈｙｓｉｃａｌｐａｒａｍｅｔｅｒｉｎｃａｌｃｕｌａｔｉｏｎｓｏｆｒａｄｉａｔｉｏｎｈｙｄｒｏｄｙｎａｍｉｃｓａｎｄｔｒａｎ...