Newtonian Kinetic Theory and the Ergodic-Nonergodic Transition

In a recent work we have discussed how kinetic theory, the statistics of classical particles obeying Newtonian dynamics, can be formulated as a field theory. The field theory can be organized to produce a self-consistent perturbation theory expansion in an effective interaction potential. In the present work we use this development for investigating ergodic-nonergodic (ENE) transitions in dense fluids. The theory is developed in terms of a core problem spanned by the variables ρ, the number density, and B, a response density. We set up the perturbation theory expansion for studying the self-consistent model which gives rise to a ENE transition. Our main result is that the low-frequency dynamics near the ENE transition is the same for Smoluchowski and Newtonian dynamics. This is true despite the fact that term by term in a density expansion the results for the two dynamics are fundamentally different.     

A laboratory model of a spring thermal bar

A laboratory setup for thermal bar phenomenon simulation is developed and mechanisms behind the thermal bar formation and evolution are experimentally studied. The results of the model experiments are compared with the data from natural observations, and the principal feasibility of laboratory simulation of thermal bar observed in natural water basins is demonstrated.     

Translational motion of domain walls in a strong magnetic field circularly polarized in the basal plane of a uniaxial ferromagnet

A substantially nonlinear theory is developed for the translational motion of domain walls (DWs) in ferromagnets with large easy-axis anisotropy under the influence of a strong magnetic field circularly polarized in the basal plane of the ferromagnet. This theory is a generalization of the well-known theories of DW drift that are limited to an approximation quadratic in the field magnitude. The analytical results are confirmed by computer simulation performed on the basis of the Landau-Lifshitz equations.     

Four-Body Bound-States Systems in Relativistic Scalar Quantum Field Theory: Variational Basis-State Approach

The variational method within the Hamiltonian formalism of quantum field theory has been used in order to investigate the effect of virtual pairs for four-body scalar systems consisting of two particles and two antiparticles of the same mass. The scalar particles and antiparticles interact via a massive or massless mediating scalar field. The ground state energy solutions of Fock-space variational trial states (\(|N{\bar {N}}N{\bar {N}}{ \rangle }+|N{\bar {N}}N{\bar {N}}N{\bar {N}\rangle }\)) of the relativistic wave equations have been studied. We have compared these results with the previous work of four-body system (variational trial states of the form \(|N{ \bar {N}}N{\bar {N}}{\rangle }\)) and we have shown that the inclusion of virtual pairs has a noticeable effect at low coupling and at high coupling the energy of the system is changed by an important amount. In other words, the calculations show that the inclusion of virtual pairs augments the binding energy of the four-body system by a substantial amount at strong coupling. This study can pave the way for some new ideas in order to investigate the effect of virtual pairs, for example, for a bound-states quark-antiquark or tetraquark systems in future.     

Evaluation of a Thermosyphon Heat Pipe Operation and Application in a Waste Heat Recovery System

Efficient and economical utilization of industrial waste heat would result in reduced energy use and thereby contribute to reduction of greenhouse gas emissions to the atmosphere. Two-phase thermosyphon technology has demonstrated the potential capability for waste heat recovery, but it has not been yet utilized in large-scale industrial applications. As a part of an industrial project, various types of thermosyphon heat pipes have been designed and tested for extraction of waste heat and process control in aluminum industry. This article presents the heat and mass transfer model, developed to provide a fast and accurate simulation tool for industrial application of thermosyphon heat pipe technology for waste heat utilization. The mathematical model considers the energy, momentum, and mass transfer equations, in their one-dimensional form, to predict output parameters of the thermosyphon and enable parametric and sensitivity analysis. The mathematical model structure is set up in a way that the least numerical cost and time is spent while the model accuracy is kept at acceptable level for the defined application. To provide experimental data for validation of the simulation model, the proposed thermosyphon was tested experimentally using a test set-up instrumented for this purpose. The simulation results are found to be in good agreement with the experimental data. The developed model and code are viable to be used as a simple and fast tool for modeling, design, and optimization of the thermosyphon as an element in a heat recovery module.     

SGS Reaction rate modelling for MILD combustion based on machine-learning combustion mode classification: Development and a priori study

A neural network (NN) aided model is proposed for the filtered reaction rate in moderate or intense low-oxygen dilution (MILD) combustion. The framework of the present model is based on the partially stirred reactor (PaSR) approach, and the fraction of the reactive structure appearing in the PaSR is predicted using different NN’s, to consider both premixed and non-premixed conditions while allowing the use of imbalanced training data between premixed and non-premixed combustion direct numerical simulation (DNS) data. The key ingredient in the present model is the use of local combustion mode prediction performed by using another NN, which is developed in a previous study. The trained model was then assessed by using two unknown combustion DNS cases, which yields much higher dilution level (more intense MILD condition) and higher Karlovitz number than the DNS cases used as training data. The model performance assessment has been carried out by means of the Pearson’s correlation coefficient and mean squared error. For both the present model and zeroth-order approximated reaction rate, the correlation coefficient with the target values shows relatively high values, suggesting that the trend of predicted field, by the present model and zeroth-order approximation, is well correlated with the actual reaction rate field. This suggests that the use of PaSR equation is promising if the fraction of the reactive structure is appropriately predicted, which is the objective in the present study. On the other hand, substantially lower mean squared error is observed for a range of filter sizes for the present model than that for the zeroth-order approximation. This suggests that the present filtered reaction rate model can account for the SGS contribution reasonably well.     

Mixed element FEM level set method for numerical simulation of immiscible fluids

A new realization of a finite element level set method for simulation of immiscible fluid flows is introduced and validated on numerical benchmarks. The new method involves a mixed discretization of the dependent variables, discretizing the flow variables with non-conforming Rannacher–Turek finite elements while using a simple first order conforming discretization of the level set field. A three step segregated solution approach is employed, first a discrete projection method is used to decouple and compute the velocity and pressure separately, after which the level set field can be computed independently.The developed method is tested and validated on a static bubble test case and on a numerical rising bubble test case for which a very accurate benchmark solution has been established. The new approach is also compared against two commercial simulation codes, Ansys Fluent and Comsol Multiphysics, which shows that the developed method is a magnitude or more accurate and at the same time significantly faster than state of the art commercial codes.     

The effects of self-fields on the electron trajectory and gain in a two-stream electromagnetically pumped free-electronlaser with axial guiding field

A theory of a two-stream free-electron laser in a combined electromagnetic wiggler (EMW) is developed, in which we use an axial-guide magnetic field and take into account the effects of the self-fields. The electron trajectories and the small signal gain are derived. The stability of the trajectories, the characteristics of the linear-gain, and the normalised maximum gain are studied numerically. The results show that there are nine stable groups of orbits in the presence of self-fields instead of seven groups reported in the absence of the self-field. It is also shown that the normalised gains of four groups of the orbits are decreasing and those for the rest of them are increasing with growing \barΩ₀. Furthermore, it is found that the two-stream laser with self-field enhances the maximum gain in comparison with the single stream case.     

A novel coupled level set and volume of fluid method for sharp interface capturing on 3D tetrahedral grids

We present a new three-dimensional hybrid level set (LS) and volume of fluid (VOF) method for free surface flow simulations on tetrahedral grids. At each time step, we evolve both the level set function and the volume fraction. The level set function is evolved by solving the level set advection equation using a second-order characteristic based finite volume method. The volume fraction advection is performed using a bounded compressive normalized variable diagram (NVD) based scheme. The interface is reconstructed based on both the level set and the volume fraction information. The novelty of the method lies in that we use an analytic method for finding the intercepts on tetrahedral grids, which makes interface reconstruction efficient and conserves volume of fluid exactly. Furthermore, the advection of volume fraction makes use of the NVD concept and switches between different high resolution differencing schemes to yield a bounded scalar field, and to preserve both smoothness and sharp definition of the interface. The method is coupled to a well validated finite volume based Navier–Stokes incompressible flow solver. The code validation shows that our method can be employed to resolve complex interface changes efficiently and accurately. In addition, the centroid and intercept data available as a by-product of the proposed interface reconstruction scheme can be used directly in near-interface sub-grid models in large eddy simulation.     

The Universe and the Quantum Computer

It is first pointed out that there is a common mathematical model for the universe and the quantum computer. The former is called the histories approach to quantum mechanics and the latter is called measurement-based quantum computation. Although a rigorous concrete model for the universe has not been completed, a quantum measure and integration theory has been developed which may be useful for future progress. In this work we show that the quantum integral is the unique functional satisfying certain basic physical and mathematical principles. Since the set of paths (or trajectories) for a quantum computer is finite, this theory is easier to treat and more developed. We observe that the sum of the quantum measures of the paths is unity and the total interference vanishes. Thus, constructive interference is always balanced by an equal amount of destructive interference. As an example we consider a simplified two-slit experiment.     

Simulation of a critical state without critical slowing down by a renormalization group multi-grid method

We have implemented and tested a method to eliminate critical slowing down from Monte Carlo (and possibly other) simulations of very large systems, even for a “critical” state, where the correlation length of fluctuations is the size of the sample. Static correlation functions on all length scales (down to microscopic!) may thus be obtained with an amount of work asymptotically growing with size only as in conventional simulations of nearly uncorrelated states. In the simulation we use a finite but large set of approximated renormalized coupling constants, which describe very closely the coarse grained variable interactions of the simulated model. As a test bed, the 2D Ising model has been used. The method has the advantage that critical states of other types of systems can be simulated along the same line.     

Numerical simulation of flow in a screw-type blood pump

This study presents a numerical investigation of the flow field in a screw pump designed to circulate biological fluid such as blood. A simplified channel flow model is used to allow analysis using a Cartesian set of coordinates. Finite analytic (FA) numerical simulation of the flow field inside the channel was performed to study the influence of Reynolds number and pressure gradient on velocity distribution and shear stresses across the channel cross-section. Simulation results were used to predict flow rates, circulatory flow and the shear stresses, which are known to be related to the level of red blood cell damage (hemolysis) caused by the pump. The study shows that high shear levels are confined to small regions within the channel cross-section, but the circulatory nature of the flow causes an increased percentage of blood elements to pass through the high shear regions, and increases the likelihood of cell damage.     

Near-field propagation of terahertz pulses from a large-aperture antenna

The near-field propagation behavior of terahertz (THz) pulses generated by a planar large-aperture photoconducting THz transmitter has been characterized. A simulation model based on Huygens-Fresnel diffraction theory has been developed that permits accurate prediction of the spatiotemporal profiles of the THz beam everywhere and gives excellent agreement with experimental measurements. Two key conclusions emerge from this research, namely, the realization that for practical laboratory setups one is always working in the near-field regime and that the proper temporal shape of the THz field at the antenna is one that rises rapidly but decays slowly.     

Dynamical implications of gluonic excitations in meson-meson systems

We study meson-meson interactions using an extended q ² [`(q)]<sup>2</sup> (g)\bar q^2 (g) basis that allows calculating coupling of an ordinary meson-meson system to a hybrid-hybrid one. We use a potential model matrix in this extended basis which at quark level is known to provide a good fit to numerical simulations of a q <sup>2</sup> [`(q)]²\bar q^2 system in pure gluonic theory for static quarks in a selection of geometries. We use a combination of resonating group method formalism and Born approximation to include the quark motion using wave functions of a q[`(q)]q\bar q potential within a cluster. This potential is taken to be quadratic for ground states and has an additional smeared $\frac{1} {r}$\frac{1} {r} (Gaussian) for the matrix elements between hybrid mesons. For the parameters of this potential, we use values chosen to 1) minimize the error resulting from our use of a quadratic potential and 2) best fit the lattice data for differences of Σ _g and Π _u configurations of the gluonic field between a quark and an antiquark. At the quark (static) level, including the gluonic excitations, was noted to partially replace the need for introducing many-body terms in a multiquark potential. We study how successful such a replacement is at the (dynamical) hadronic level of relevance to actual hard experiments. Thus we study the effects of both gluonic excitations and many-body terms on mesonic transition amplitudes and the energy shifts resulting from the second-order perturbation theory (i.e. from the respective hadron loops). The study suggests introducing both energy and orbital excitations in wave functions of scalar mesons that are modelled as meson-meson molecules or are supposed to have a meson-meson component in their wave functions.     

非对称电极表面微观形貌对交流电渗流速的影响

经典交流电渗理论是利用电场进行非机械式微流体驱动的基础.传统理论交流电渗理论以双电层理论为基础,通过耦合电场方程以及流场方程得到微电极表面交流电渗流速表达式,通常与实验流速相差较大. 以电极表面微观形貌对交流电渗流速的影响为研究目标,定义微电极表面粗糙度为微观形貌特征参数,建立了等效双电层模型,并对传统交流电渗流速公式进行了修正.理论并仿真分析了表面粗糙度对于交流电渗流速的影响,利用非对称电极对交流电渗微流体驱动进行了实验研究,并进行对比分析.结果表明,理论分析与实验结果具有较好的一致性. 关键词： 交流电渗 电极表面粗糙度 等效双电层     

新型二次电子倍增阴极的蒙特卡罗模拟研究

提出了一种新型二次电子倍增阴极强流二极管,并对其进行了动力学理论简化模型和蒙特卡罗数值模拟的对比验证研究。首先,基于设计结构原型,根据二次电子发射特性进行合理简化,建立了动力学模型,获得了电子速度、位移以及渡越时间的解析结果,并结合Vaughan的二次电子产额模型,确定了该新型二次电子倍增阴极强流二极管的理论工作区间;其次,理论分析了施加径向电场的重要意义,并给出了二次电子运动特征参数(最大位移、渡越时间、碰撞能量等)的理论预估结果;最后,对该新型二次电子倍增阴极强流二极管进行了蒙特卡罗模拟研究,获得了电子的运动轨迹、碰撞能量以及二次电子倍增工作区间等物理图像,并将蒙特卡罗数值模拟结果与理论结果进行了比对,两者吻合程度较好,对可能的误差来源进行了分析讨论。理论和模拟结果表明:新型二次电子倍增阴极强流二极管概念可行,工作区间内通过调整施加电场与磁场幅值,可有效达到电子运动状态可控的目标。另外,理论粗估了二次电子倍增饱和条件下的阴极发射电流密度,结果表明:发射电流密度可达kA/cm2水平,具备强流发射特性;增加外加径向场强幅值可有效提升发射电流密度。最后,对该新型二次电子倍增阴极设计步骤和依据进行了讨论。     

声表面横波压力传感器的电极质量负载效应分析方法

针对声表面横波(Surface Transverse Wave,STW)压力传感器,提出了一种考虑叉指电极质量负载的有限元和微扰理论结合分析方法。首先构建包含叉指电极的有限元模型,通过计算电极表面的面电荷得到器件的等效输入导纳,由导纳谐振点处提取STW传播声场,然后根据微扰理论,由STW传播声场以及压力引起的传感器结构的应力应变,得到包含电极质量负载影响的STW传感器的灵敏度。以一种基于ST-90°X石英材料的点压式STW压力传感器为例进行了理论计算,电极厚度为0.17 μm和0.21 μm时,计算灵敏度分别为-54.9 kHz/bar和-21.7 kHz/bar,说明电极质量负载对于传感器有较大影响。制作电极厚度为0.17 μm的传感器进行测试,得到-57 kHz/bar的压力灵敏度,实验结果和理论结果吻合,证明计算方法的有效性。      

Effect of a two-phase wedge-sliding model on the ingredient drift of a stable mixed fluid and its computing method

A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system,and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon.In the model,the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids.For solving above two indexes of the model,a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress,respectively.As an example,the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.