Variational formulation of ideal fluid flows according to gauge principle

On the basis of the gauge principle of field theory, a new variational formulation is presented for flows of an ideal fluid. The fluid is defined thermodynamically by mass density and entropy density, and its flow fields are characterized by symmetries of translation and rotation. The rotational transformations are regarded as gauge transformations as well as the translational ones. In addition to the Lagrangians representing the translation symmetry, a structure of rotation symmetry is equipped with a Lagrangian Λ_A including the vorticity and a vector potential bilinearly. Euler's equation of motion is derived from variations according to the action principle. In addition, the equations of continuity and entropy are derived from the variations. Equations of conserved currents are deduced as the Noether theorem in the space of Lagrangian coordinate a. Without Λ_A, the action principle results in the Clebsch solution with vanishing helicity. The Lagrangian Λ_A yields non-vanishing vorticity and provides a source term of non-vanishing helicity. The vorticity equation is derived as an equation of the gauge field, and the Λ_A characterizes topology of the field. The present formulation is comprehensive and provides a consistent basis for a unique transformation between the Lagrangian a space and the Eulerian x space. In contrast, with translation symmetry alone, there is an arbitrariness in the transformation between these spaces.     

Numerical simulation of the noise generated by a low Mach number, low Reynolds number jet

In this paper we study the sound field produced by a turbulent round jet with a Mach number of 0.6 based on the centerline velocity and the ambient speed of sound c_∞. The turbulent flow field is found by solving the fully compressible Navier–Stokes equations with help of high-order compact finite difference schemes. It is shown that the simulated flow field is in good agreement with experiments. The corresponding sound field has been obtained with help of the Lighthill equation using two different formulations for the Lighthill stress tensor T_ij. In the first formulation of T_ij the fluctuating density is taken into account. In the second formulation the density is assumed to be constant. As an additional check we have also performed an acoustic calculation using a formulation in which a homogeneous wave equation is solved. The boundary conditions for this homogeneous wave equation are obtained from the numerical simulation of the Navier–Stokes equation. The results obtained with both formulations of the Lighthill stress tensor are nearly identical. This implies that an incompressible formulation of the conservations laws could be used to predict jet noise at low Mach numbers.     

A study of electromechanical switching in ferroelectric single crystals

This article documents both modeling and experimental studies developed to investigate the switching behavior of ferroelectric single crystals. The theoretical model makes a priori ansatz that switching follows the evolution of a particular domain pattern. The choice of this configuration is dictated by the requirement that domains remain compatible during evolution, giving rise to a low-energy path for the overall switching. The construction of this pattern is achieved using multirank laminates. It offers an advantage of specifying different types of domain wall movements, leading to a distinction for the switching types. A loading experiment is performed on a barium titanate (BaTiO₃) single crystal with a constant compressive stress and a cyclic electric field. Both 180 and 90 coercive fields are measured as input parameters required for the theoretical framework. The simulation results show good agreement with the observed strains measured by the present and other available experiments. It is found that depolarization has a non-trivial influence on attainable actuation strains.     

Experimental study of liquid sloshing dynamics in a barge carrying tank

An experimental work has been carried out to study the phenomena of sloshing of liquid in partially filled tanks mounted on a barge exposed to regular beam waves. Three liquid fill levels with liquid depth, h_s to length of tank, l ratio (h_s/l) of 0.163, 0.325 and 0.488, are studied. The time histories of sloshing oscillation are measured along the length of container at predefined locations. The nonlinear behaviour of sloshing oscillation is observed for the regular wave excitation. The spectra of the sloshing oscillation and their qualitative assessment are reported. Attempts are made to evaluate the harmonics present in the sloshing oscillation and compare with the results of earlier studies. The effects of wave excitation frequency and wave height on the sloshing oscillation as well as on the response of the barge are studied.     

On the flow in a uniformly porous pipe

This article considers fully laminar flow of an incompressible viscous fluid in a uniformly porous pipe with suction and injection. An exact solution of the Navier–Stokes equations is given. The velocity filed can be expressed in a series form in terms of the modified Bessel function of the first kind of order n. The volume flux across a plane normal to the flow, the vorticity and the stress on the boundary are presented. The flow properties depend on the cross-Reynolds number, Ua/ν, where U is the suction velocity, a is the radius of the pipe and ν is the kinematic viscosity of the fluid. It is found that for large values of the cross-Reynolds number, the flow near the region of the suction shows a boundary layer character. In this region the velocity and the vorticity vary sharply. Outside the boundary layer, the velocity and the vorticity do not show an appreciable change.     

Separation of variables in one problem of motion of the generalized Kowalevski top

In the problem of motion of the Kowalevski top in a double force field the new case of reduction to a Hamiltonian system with two degrees of freedom was pointed out by Kharlamov [Kharlamov, M.P., 2004. Mekh. Tverd. Tela 34, 47–58]. We show that the equations of motion in this case can be separated by the appropriate change of variables, the new variables U,V being hyperelliptic functions of time. The natural phase variables (components of the angular velocity and the direction vectors of the forces with respect to the movable basis) are expressed via U,V explicitly in elementary algebraic functions.     

Linear and nonlinear sloshing in a circular conical tank

Linear and nonlinear fluid sloshing problems in a circular conical tank are studied in a curvilinear coordinate system. The linear sloshing modes are approximated by a series of the solid spheric harmonics. These modes are used to derive a new nonlinear modal theory based on the Moiseyev asymptotics. The theory makes it possible to both classify steady-state waves occurring due to horizontal resonant excitation and visualise nonlinear wave patterns. Secondary (internal) resonances and shallow fluid sloshing (predicted for the semi-apex angles >60) are extensively discussed.     

Autoparametric quasiperiodic excitation

The dynamics of an autonomous conservative three degree of freedom system which exhibits autoparametric quasiperiodic excitation is investigated. The system is a generalization of a classical system known as the “particle in the plane”. The system exhibits a motion, the z=0 mode, whose stability is governed by a linear second order ODE with quasiperiodic coefficients. The behavior of the latter ODE is studied by using three different methods: numerical integration, harmonic balance and perturbation methods.     

Analysis for a screw dislocation accelerating through the shear-wave speed barrier

An analysis is performed for an accelerating screw dislocation through the shear-wave speed barrier. At this instant, the function that determines the interval of the path of the dislocation motion that contributes to the wave front has roots that change from a pair of complex conjugate to a double real, which subsequently splits into two real ones. The analysis is performed at this transition to supersonic that occurs at the double root maximum of the function that defines the interval of the dislocation path that contributes to the field points. It is found that the stress has a log|ξ-ξ^*|/|ξ-ξ^*|^1/2 singularity in the coefficient of the delta function of the forming Mach front, implying that for this phenomenon the Volterra dislocation model has too strong a discontinuity (step-function) in the displacement to be meaningful. A ramp-core displacement dislocation model analysis, which removes the singularity in the stress, is presented. These results can be useful in a multiscale dislocation dynamics modeling with inertia effects.     

Predicting transition in two- and three-dimensional separated flows

This paper is concerned with the numerical prediction of two- and three-dimensional transitional separated flows of turbomachinery interest. The recently proposed single-point transition model based on the use of a laminar kinetic energy transport equation is considered, insofar as it does not require to evaluate any integral parameter, such as boundary-layer thickness, and is thus directly applicable to three-dimensional flows. A well established model, combining a transition-onset correlation with an intermittency transport equation, is also used for comparison. Both models are implemented within a Reynolds-averaged Navier–Stokes solver employing a low-Reynolds-number k–ω turbulence model. The performance of the transition models have been evaluated and tested versus well-documented incompressible flows past a flat plate with semi-circular leading edge, namely: tests T3L2, T3L3, T3L5, and T3LA1 of ERCOFTAC, with different Reynolds numbers and free-stream conditions, the last one being characterized by a non-zero pressure gradient. In all computations, the first model has proven as adequate as or superior to the second one and has been then applied with success to two more complex test cases, for which detailed experimental data are available in the literature, namely: the two- and three-dimensional flows through the T106 linear turbine cascade.     

功能梯度材料动力学问题的POD模型降阶分析

为了快速分析非均质材料结构在复杂载荷作用下的动态响应, 提出一种模型降阶方法, 只需计算结构在简单均质材料情况下的动力学问题, 进而用其计算结果对非均质材料结构进行分析. 首先, 采用结构内部任意一点处的材料参数值作为整个结构的材料参数, 利用有限元分析软件计算该均质材料结构在动态载荷作用下的位移场建立数据库, 该数据库包含计算模型各个节点(自由度为L)在某时间段内L个时刻的位移; 其次, 对数据库中的信息按照时间离散的特定方式组集成瞬像矩阵, 并利用特征正交分解方法对其进行分解, 得到该模型的H个特征正交基底, 选取其中能反应模型主要特征的个(其中\mathrm{~})作为一组最优基底, 通过这组基底建立模型的低阶离散控制方程; 最后, 求解低阶离散微分方程组, 得到功能梯度材料结构在复杂载荷作用下的位移场. 文中分别给出二维和三维算例, 比较了降阶模型和全阶模型计算结果, 验证了该方法的有效性, 并且计算效率能提高1{\mathrm{}}^{1)}2个数量级.      

错列角度对双圆柱涡激振动影响的数值模拟研究

为研究错列角度\alpha对双圆柱涡激振动问题的影响,采用自主研发的基于CIP (constrained interpolation profile)方法的数值模型,对雷诺数\mathit{Re}=100、错列角度\alpha =0^{°}\mathrm{~}90^{°} (间隔15 #x00B0;)的等直径双圆柱涡激振动问题进行数值模拟. 模型在笛卡尔网格系统下建立,采用具有三阶精度的 CIP 方法求解 N-S (Navier--Stokes)方程,采用浸入边界法处理流--固耦合问题,避免了任意拉欧方法下的网格畸变和重叠动网格技术中的大量信息交换问题,保证了模型的计算效率. 重点分析不同错列角度\alpha上下游圆柱的升阻力系数、位移响应、涡脱频率和尾涡模态等. 结果表明:折合速度U_r=2.0\mathrm{~}3.0时,上下游圆柱升阻力随错列角度的增大基本呈单调增大的趋势;U_{\mathrm{r}}=5.0\mathrm{~}8.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“上凸”趋势,在\alpha =15^{°}\mathrm{~}30^{°}取得最大值;U_r=10.0\mathrm{~}13.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“下凹”趋势,在\alpha =30^{°}\mathrm{~}45^{°}取得最小值,且柱体横流向振幅和升力没有明显的对应关系. 最后,结合尾涡模态对以上规律的成因进行分析. 研究结果可为相关海洋工程设计提供参考.      

高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性. 计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数\overline{C_d}、升力系数均方根{C\mathrm{\text{'}}}_l及斯特劳哈尔数\mathit{St}维持在一定范围内波动. 串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0, 2.5, 3.0, 3.5和4.0这五中间距比进行计算. 串列双圆柱绕流中,当\mathit{Re}=2.2\times 10^4时,在内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当时,下游圆柱的阻力系数为负数. 而当\mathit{Re}=3.0\times 10^6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数. 串列双方柱绕流在\mathit{Re}=1.6\times 10^4和\mathit{Re}=1.0\times 10^6两种工况下的临界间距比分别处于和区间内,且当时,两个雷诺数下的下游方柱阻力系数均为负数.      

应力诱发界面迁移下晶内孔洞的演化

随着微电子技术的迅猛发展, 集成电路中内连导线的失效问题引起广泛关注. 内连导线内部孔洞萌生、长大、漂移和失稳变形成狭长裂纹, 从而导致电路的开路失效. 这是内连导线失效的常见形式. 而界面迁移是导致微结构形态演化的主要机制之一. 本文基于界面迁移下微结构演化的经典理论和弱解描述, 建立了应力诱发界面迁移下微结构演化的有限单元法, 并验证了算法的可靠性. 对铜内连导线中晶内孔洞的演化进行了数值模拟, 详细分析了应力、线宽及形态比对晶内孔洞演化的影响. 研究结果表明, 椭圆形晶内孔洞存在生长和收缩两种演化分叉趋势. 通过大量数值分析得到了晶内孔洞演化的临界应力{\stackrel{\text{?}}{\sigma }}_c 、临界线宽{\stackrel{\text{?}}{h}}_c 和临界形态比{\beta }_c . 当\stackrel{\text{?}}{h}\mathrm{?}{\stackrel{\text{?}}{h}}_c, \stackrel{\text{?}}{\sigma }或\stackrel{\text{?}}{h} 时, 晶内孔洞会沿长轴长大; 反之, 晶内孔洞会收缩甚至愈合. 此外, 应力\beta 越大、线宽\stackrel{\text{?}}{\sigma } 越小或形态比\stackrel{\text{?}}{h} 越大, 晶内孔洞越易发生长大, 且孔洞面积增大速度越快; \beta 越小、\stackrel{\text{?}}{\sigma } 越大或\stackrel{\text{?}}{h} 越小, 晶内孔洞越易发生收缩, 且孔洞面积减小速度越快.      

A robust near-wall elliptic-relaxation eddy-viscosity turbulence model for CFD

An eddy-viscosity model based on Durbin’s elliptic relaxation concept is proposed, which solves a transport equation for the velocity scales ratio instead of , thus making the model more robust and less sensitive to grid nonuniformities. Computations of flows and heat transfer in a plane channel, behind a step and in a round impinging jet show all satisfactory results.     

A gradient theory of small-deformation isotropic plasticity that accounts for the Burgers vector and for dissipation due to plastic spin

This study develops a gradient theory of small-deformation viscoplasticity based on: a system of microforces consistent with its peculiar balance; a mechanical version of the second law that includes, via the microforces, work performed during viscoplastic flow; a constitutive theory that accounts for the Burgers vector through a free energy dependent on , with H^p the plastic part of the elastic-plastic decomposition of the displacement gradient. The microforce balance and the constitutive equations, restricted by the second law, are shown to be together equivalent to a nonlocal flow rule in the form of a coupled pair of second-order partial differential equations. The first of these is an equation for the plastic strain-rate in which the stress T plays a basic role; the second, which is independent of T, is an equation for the plastic spin. A consequence of this second equation is that the plastic spin vanishes identically when the free energy is independent of, but not generally otherwise. A formal discussion based on experience with other gradient theories suggests that sufficiently far from boundaries solutions should not differ appreciably from classical solutions, but close to microscopically hard boundaries, boundary layers characterized by a large Burgers vector and large plastic spin should form.Because of the nonlocal nature of the flow rule, the classical macroscopic boundary conditions need be supplemented by nonstandard boundary conditions associated with viscoplastic flow. As an aid to solution, a variational formulation of the flow rule is derived.Finally, we sketch a generalization of the theory that allows for isotropic hardening resulting from dissipative constitutive dependences on .     

Collision operator including large-angle scattering for moderately coupled plasma

A new approach has been developed to treat the large-angle as well as the small-angle binary collisions in high temperature and high density plasmas when the test particle distribution function f_α is even function about the test velocity and the relations of the mass and the velocity between the test particles and the field particles are satisfied with m_αm_β (such as electron–ion collision or Lorentz-gas model) and . With the approach, the Boltzmann collision operator is derived to be suitable for the plasma considered as weakly coupled (Coulomb logarithm ) and moderately coupled , i.e., for the electron–ion coupling constant Γ_ei<0.1. The modified collision operator has a direct and practical connection to the Rosenbluth potentials, the new reduced electron–ion collision operator differs from the original Fokker–Planck operator for Coulomb collisions by terms of order . Moreover, some calculations of relaxation rate and transport properties are given for new reduced electron–ion collision operator that shows corrections.     

A large deformation theory for rate-dependent elastic–plastic materials with combined isotropic and kinematic hardening

We have developed a large deformation viscoplasticity theory with combined isotropic and kinematic hardening based on the dual decompositions F=F^eF^p [Kröner, E., 1960. Allgemeine kontinuumstheorie der versetzungen und eigenspannungen. Archive for Rational Mechanics and Analysis 4, 273–334] and [Lion, A., 2000. Constitutive modelling in finite thermoviscoplasticity: a physical approach based on nonlinear rheological models. International Journal of Plasticity 16, 469–494]. The elastic distortion F^e contributes to a standard elastic free-energy ψ^(e), while , the energetic part of F^p, contributes to a defect energy ψ^(p) – these two additive contributions to the total free energy in turn lead to the standard Cauchy stress and a back-stress. Since F^e=FF^p-1 and , the evolution of the Cauchy stress and the back-stress in a deformation-driven problem is governed by evolution equations for F^p and – the two flow rules of the theory.We have also developed a simple, stable, semi-implicit time-integration procedure for the constitutive theory for implementation in displacement-based finite element programs. The procedure that we develop is “simple” in the sense that it only involves the solution of one non-linear equation, rather than a system of non-linear equations. We show that our time-integration procedure is stable for relatively large time steps, is first-order accurate, and is objective.