Nonlinear oscillations of a charged drop levitated in gravity and electrostatic fields

Analytically, on the basis of asymptotic methods, the problem of the nonlinear oscillations of a charged ideal incompressible electroconductive fluid drop levitated at rest in gravity and homogeneous electrostatic fields is solved in the quadratic approximation in two small parameters: the initial drop shape deformation amplitude and the stationary eccentricity of the equilibrium drop shape in the electrostatic field. The calculations are performed in fractional powers of the nonlinear oscillation amplitude. The nonlinear corrections to the oscillation frequencies are always negative and already present in the second-order approximation due to the stationary deformation of the drop in the external fields rather than nonlinear interaction between the modes. In the case considered, in contrast to the nonlinear oscillations of a free charged drop, the expression for the generator of the nonlinearly oscillating drop shape contains terms proportional to the oscillation amplitude to the power 3/2.     

Nonlinear Oscillations of a Charged Drop with Multimodal Initial Deformation of Its Equilibrium Shape

For axisymmetric oscillations of a charged drop, in the quadratic approximation in the amplitude of an arbitrary initial deformation of its equilibrium spherical shape, an expression describing its surface shape as a function of time is obtained. Regularities in the formation of the spectrum of modes excited in the second order as a result of intermodal interaction are analyzed. It is shown that nonlinear oscillations of the drop surface occur in the neighborhood of an elongated spheroid-like figure, not a sphere as followed from the linear analysis.     

Nonlinear analysis of the finite-amplitude oscillations of a charged liquid layer on a rigid spherical core for the multimodal initial deformation of the free surface

Within the framework of the analytical asymptotic procedure, in the quadratic approximation in the small dimensionless amplitude of the capillary oscillation of a charged electroconductive fluid layer on the surface of a rigid spherical core, the problem of calculation of the free surface shape at an arbitrary moment is solved for a multimodal initial deformation. The conditions under which the internal nonlinear resonant interaction between the liquid layer oscillation modes and the intermodal energy exchange are realized are analyzed for both the degenerate and the secondary combinational resonances.     

Equilibrium shapes of a drop pendant in an electrostatic field

The stationary shapes of a conducting fluid drop in the gap between the plates of a plane capacitor are studied. The drop is held on the upper plate by the surface tension forces. The self-consistent problem of the determination of the drop shape and the charge distribution over its surface is solved. Estimates are obtained for the maximum volume of the stationary drop at the given fluid parameters and electric field strength.     

Further considerations of two-dimensional condensation drop profiles and departure sizes

The problem of the equilibrium shape and departure size of two-dimensional dropwise condensation drops on a vertical surface, presented in an earlier work, is extended to include advancing contact angles to 180°. The equation of the surface of the drop is obtained by minimizing (for a given volume) the total energy of the drop, consisting of surface and gravitational energy, using the techniques of variational calculus. The solution is tractable once the advancing contact angle is known, and is taken as an approximation to the axial meridian profile of a threedimensional drop. The receding contact angle is obtained as part of the solution. The drop size is specified by imposing its vertical length in contact with the wall. A maximum value of this length exists which provides a real solution, and this is taken as the departure size of the drop. It is shown that the general departure shape for an advancing contact angle of 180° includes the cases for all advancing contact angles.     

Annular surface-charged jet in an external electric field

The problem of the shape of an annular charged jet in an external electric field is considered. In the strong-field approximation, the formulation of the problem is reduced to a system of ordinary differential equations. The asymptotics of the solutions at large distances from the point of outflow are investigated. The condition of formation of the jet profile non-increasing along the streamwise coordinate is found.     

The influence of the curvature dependence of the surface tension on the geometry of electrically charged menisci

We evaluate how the curvature dependence of surface tension affects the shape of electrically charged interfaces between a perfectly conducting fluid and its vapour. We consider two cases: i) spherical droplets in equilibrium with their vapour; ii) menisci pending in a capillary tube in presence of a conducting plate at given electric potential drop. Tolman-like dependence of surface tension on curvature becomes important when the “nucleation radius” is comparable with the interface curvature radius. In case i) we prove existence of the equilibrium minimal radius and estimate its dependence on the electric fields and Tolman-like curvature effects. In case ii) the menisci are subject to the gravitational force, surface tension and electrostatic fields. We determine the unknown surface of the menisci to which the potential is assigned using an iterative numerical method and show that Tolman-like corrections imply: 1) a variation of the height (up to 10% in some cases) of the tip of the menisci; 2) a decrease of the maximum electrical potential applicable to the menisci before their break-down amounting to 40V over 800V in the considered cases. We conjecture that these effects could be used in new experiments based on electric measurements to determine the dependence of the equilibrium surface tension on curvature. Received January 19, 1998     

Formation and stability of a charged jet in a strong electric field

The shape of a charged jet is determined in the approximation of a strong electric field. The stability of the jet with respect to both axisymmetric and nonaxisymmetric perturbations of the sinusoidal type is investigated in the linear approximation. The domains of predominance of the axisymmetric and bending modes and the longitudinal partition mode are determined. Experimental data on the longitudinal partition of a polymeric jet into several daughter jets are given.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

Profile and departure size of condensation drops on vertical surfaces

The problem of the equilibrium shape and departure size of two-dimensional dropwise condensation drops on a vertical surface is considered. The equation of the surface of the drop is obtained by minimizing (for a given volume) the total energy of the drop, consisting of surface and gravitational energy, using the techniques of variational calculus. The solution is tractable once the advancing contact angle is known, and is taken as an approximation to the axial meridian profile of a three-dimensional drop. The receding contact angle is obtained as part of the solution. The drop size is specified by imposing its vertical length. Upon increasing this vertical length, a point is reached at which no real solution exists, and this is taken as the departure size of the drop. Comparison with measured departure sizes under various body forces from standard to 100 times earth gravity are good.     

Ferrite-filled cavity resonators

An approximation technique is developed for the electromagnetic resonances and electric fields inside a cavity of arbitrary shape whose walls are perfectly conducting and which is filled with a lossless ferrite. Operator notation is introduced and it is proved that the operator for this problem is self-adjoint. A variational expression is introduced and this functional is minimized by employing the Rayleigh-Ritz technique. The solution is in the form of a matrix eigenvalue equation. The general formulas are specialized to the case of a ferrite-filled spherical cavity resonator and some of the lower-order mode resonances are calculated. The technique is briefly contrasted with other approximation techniques which are found in the literature.     

A non-local implicit gradient-enhanced model for unstable behaviors of pseudoelastic shape memory alloys in tensile loading

In this paper, a gradient-enhanced 3-D phenomenological model for shape memory alloys using the non-local theory is developed based on a 1-D constitutive model. The method utilizes a non-local field variable in its constitutive framework with an implicit gradient formulation in order to achieve results independent of the finite element discretization. An efficient numerical approach to implement the non-local gradient-enhanced model in finite element codes is proposed. The model is used to simulate stress drop at the onset of transformation, and its performance is evaluated using different experimental data. The potential of the presented numerical approach for behavior of shape memory alloys in eliminating mesh-dependent simulations is validated by conducting various localization problems. The numerical results show that the developed model can simulate the observed unstable behaviors such as stress drop and deviation of local strain from global strain during nucleation and propagation of martensitic phase.     

Mass transfer between particles and a flow in the presence of a volume chemical reaction

The exterior problem of the mass transfer between a spherical drop and a linear shear flow in the presence of a first-order volume reaction is solved in the diffusion boundary layer approximation. A simple approximate expression for calculating the average Sherwood number for a drop or solid particle of arbitrary shape is proposed. At large Péclet numbers this expression is applicable to any type of flow over the entire range of variation of the reaction rate constant. The problem of diffusion to a spherical drop in a translational Stokesian flow in the presence of a first-order volume reaction was investigated in [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 109–113, November–December, 1987.     

基于内聚力理论的二维二次界面单元在ABAQUS中的UEL程序实现

依据有限元理论,结合内聚力模型法则,推导出二维二次粘结界面单元在大位移情况下的数值格式,得到用形函数表示的单元位移模式、载荷向量和刚度矩阵,并进行了离散化。基于ABAQUS软件的自定义扩展模块,编制了相应的用户单元子程序UEL,通过数值算例验证了该程序的准确性和有效性。这一成果能为在ABAQUS软件中开展相关数值研究,以及开发其他类型的内聚力界面有限单元提供思路和参考。     

薄板分析的线性基梯度光滑伽辽金无网格法

薄板问题的控制方程为四阶微分方程,因而当采用伽辽金法进行分析时,形函数需要满足C$^{1}$连续性要求,且至少使用二次基函数才能保证方法的收敛性.无网格形函数虽然易于满足C$^{1}$连续性要求,但由于不是多项式,其二阶导数的计算较为复杂耗时,同时也对刚度矩阵的数值积分提出了更高的要求.本文提出了一种薄板分析的线性基梯度光滑伽辽金无网格法,该方法的基础是线性基无网格形函数的光滑梯度.在梯度光滑构造的理论框架内,无网格形函数的二阶光滑梯度可以表示为形函数一阶梯度的线性组合,因而可以提高形函数二阶梯度的计算效率.分析表明,线性基无网格形函数的光滑梯度不仅满足其固有的线性梯度一致性条件,还满足本属于二次基函数对应的额外高阶一致性条件,因此能够恰当地运用到薄板结构的伽辽金分析.此外,插值误差分析也很好地验证了线性基无网格光滑梯度的收敛特性.算例结果进一步表明,线性基梯度光滑伽辽金无网格法的收敛率与传统二次基伽辽金无网格法相当,但精度更高,同时刚度矩阵所需的高斯积分点数明显减少.      

智能抛物面天线多目标最优控制

大型天线必须保持非常精确的形状。由于载荷变化及各种随机干扰,为了保持精确的形状,这类结构宜采用主动控制。针对智能抛物面天线结构,建立了对天线反射面的最佳吻合抛物面的RMS精度和作动器能耗为综合目标的多目标优化模型。优化模型以结构强度和作动器性能作为约束条件,模型转化为二次规划问题进行求解,实现智能抛物面天线的准静态最优控制。算例表明,可以用较少的作动器,实现大型天线结构的精密控制。     

Spreading of a liquid over a plane rigid substrate in an electric field

The influence of an electric field on spreading of a thin conducting liquid layer over a plane rigid substrate is investigated theoretically. The conductivity of the liquid is assumed to be so low that the effect of the magnetic field of the currents generated in the liquid under the action of the electric field can be neglected. The spreading is assumed to be so slow that the quasi-steady approximation can be used to calculate the electric field strength which can be considered to be equal to zero inside the liquid. Equations that describe variations in the layer shape are obtained in the lubrication theory approximation. The general formulation of the problem is considered. The solution of the problem is obtained in parametric form when the effect of the gravity force and the surface tension can be neglected. Variations in the layer thickness along the substrate are so smooth that the charge distribution over its surface can be assumed to be the same as that over the substrate surface in the absence of the liquid.     

Radiation of electromagneticwaves by a spheroidal charged drop oscillating in a uniform electrostatic field

The dispersion relation is derived for a spheroidal charged drop oscillating in a uniform electrostatic field and radiating electromagnetic waves in the first and second orders with respect to the dimensionless oscillation amplitude and the steady-state deformation, respectively. For an individual drop and a model cloud the intensity of the electromagnetic radiation and the frequency bandwidth are estimated as functions of the drop dimensions and charge and the external electric field strength.     

Analysis of physical optics far field inverse scattering for the limited data case using radon theory and polarization information

In the inverse scattering problem for perfectly conducting objects, the reconstruction of the shape and size of a convex body from its cross-sectional areas has been formulated as a Radon problem of shape reconstruction from projections. It is shown here that the Physical Optics Inverse Scattering time-domain and frequency-domain identities form a Radon-Fourier transform pair, and the problem of target reconstruction from incomplete data is common to both. The mathematical aspects of reconstruction from projections are examined using concepts of Radon's theory, and the sparse data problemis analyzed. The limited aperture problem is solved via the Radon transform approach utilizing properties of Ludwig's theorems on support. In order to obtain more accurate cross-sectional areas for reconstruction, polarization utilization is investigated to correct the polarization-independent deficiency in the Physical Optics approximation. The technique is applied to the case of a sphere-capped cylinder and the results show substantial improvement over previous scalar approaches.     

Hydrodynamic interaction between a flat surface and an evaporating drop in its own superheated vapor in the case of small Reynolds and Knudsen numbers

An expression for the force of interaction between a flat surface and an evaporating drop moving along the normal to this surface is obtained in the approximation of the hydrodynamic lubrication theory. The gap between the surface and the drop is small. The effects of the slip, the temperature jump, and the evaporation rate of the drop on the time of variation of this gap are considered under the assumption that the temperature of the flat surface exceeds the boiling temperature of the drop.