Computational fluid dynamics model of rhythmic motion of charged droplets between parallel electrodes

2017 ◽  
Vol 822 ◽  
pp. 31-53 ◽  
Author(s):  
Rudolf Flittner ◽  
Michal Přibyl

A mathematical model of rhythmic motion of a charged droplet between two parallel electrodes is developed in this study. The work is motivated by recent experimental findings that report oscillatory behaviour of water in oil droplets under a direct current electric field. The model considers the presence of a charged droplet placed in a dielectric medium. The droplet is immediately attracted to the electrode with the opposite polarity. When approaching the electrode, the electric charge is electrochemically reversed within the droplet, which is then repelled to the other electrode. The entire process can periodically repeat. The model is able to track a deformable liquid–liquid interface, the dynamics of the wetting process at the electrodes and the dynamics of electrochemical charge transfer between the droplet and the electrodes. The dependences of the oscillation frequency, charge acquired by the droplet and charging time on several model parameters (surface charge density on electrodes, kinetic parameter of charging, droplet–electrode contact angle, droplet size, liquid permittivity) are examined. Qualitative agreement of the model predictions with available experimental data is obtained, e.g. the oscillation frequency increases with growing electric field strength or droplet size. Our model represents the first successful attempt to predict oscillatory motion of aqueous droplets by a pseudo-three-dimensional two-phase approach. Our model also strongly supports the theory that the oscillatory motion relies on the combination of electrochemical charge injection at the electrodes and electrostatic attraction/repulsion processes.

Energies ◽  
2018 ◽  
Vol 11 (9) ◽  
pp. 2399 ◽  
Author(s):  
Fengbo Yang ◽  
Xinyu Xue ◽  
Chen Cai ◽  
Zhu Sun ◽  
Qingqing Zhou

In recent years, multirotor unmanned aerial vehicles (UAVs) have become more and more important in the field of plant protection in China. Multirotor unmanned plant protection UAVs have been widely used in vast plains, hills, mountains, and other regions, and become an integral part of China’s agricultural mechanization and modernization. The easy takeoff and landing performances of UAVs are urgently required for timely and effective spraying, especially in dispersed plots and hilly mountains. However, the unclearness of wind field distribution leads to more serious droplet drift problems. The drift and distribution of droplets, which depend on airflow distribution characteristics of UAVs and the droplet size of the nozzle, are directly related to the control effect of pesticide and crop growth in different growth periods. This paper proposes an approach to research the influence of the downwash and windward airflow on the motion distribution of droplet group for the SLK-5 six-rotor plant protection UAV. At first, based on the Navier-Stokes (N-S) equation and SST k–ε turbulence model, the three-dimensional wind field numerical model is established for a six-rotor plant protection UAV under 3 kg load condition. Droplet discrete phase is added to N-S equation, the momentum and energy equations are also corrected for continuous phase to establish a two-phase flow model, and a three-dimensional two-phase flow model is finally established for the six-rotor plant protection UAV. By comparing with the experiment, this paper verifies the feasibility and accuracy of a computational fluid dynamics (CFD) method in the calculation of wind field and spraying two-phase flow field. Analyses are carried out through the combination of computational fluid dynamics and radial basis neural network, and this paper, finally, discusses the influence of windward airflow and droplet size on the movement of droplet groups.


2021 ◽  
Author(s):  
Andrea Manzoni ◽  
Aronne Dell'Oca ◽  
Martina Siena ◽  
Alberto Guadagnini

<p>We consider transient three-dimensional (3D) two-phase (oil and water) flows, taking place at the core-scale. In this context, we aim at exploiting the full information content associated with available information of (i) the 3D distribution of oil saturation and (ii) the overall pressure difference across the rock sample, to estimate the set of model parameters. We consider a continuum-scale description of the system behavior upon relying on the widely employed Brooks-Corey model for the characterization of relative permeabilities and on the capillary pressure correlation introduced by Skjaeveland et al. (2000). To provide a transparent way of assessing the results of the inversion, we rely on a synthetic reference scenario. The latter is intended to mimic having at our disposal 3D and section-averaged distributions of (time-dependent) oil saturations of the kind that can be acquired during typical laboratory experiments. These are in turn corrupted by way of a random noise, to address the influence of experimental uncertainties. We focus on diverse scenarios encompassing imbibition and drainage conditions. We employ two population-based optimization algorithms, i.e., (i) the particle swarm optimization (PSO); and (ii) the differential evolution (DE), which enable one to effectively tackle the high-dimensionality parameters space (i.e., 12 dimensions in our setting) we consider. Model calibration results are of satisfactory quality for the majority of the tested scenarios, whereas the DE algorithm is associated with highest effectiveness.</p><p><strong>References</strong></p><p>S.M. Skjaeveland; L.M. Siqveland; A. Kjosavik; W.L. Hammervold Thomas; G.A. Virnovsky (2000). Capillary Pressure Correlation for Mixed-Wet Reservoirs SPE Res Eval & Eng 3 (01): 60–67. https://doi.org/10.2118/60900-PA</p>


2020 ◽  
Vol 2020 ◽  
pp. 1-17
Author(s):  
Jianxiu Qin ◽  
Huiqiang Zhang

Combustion instabilities in a small MMH/NTO liquid rocket engine used for satellite attitude and course control are numerically investigated. A three-dimensional Navier-Stokes code is developed to simulate two-phase spray combustion for cases with five different droplet Sauter Mean Diameters. As the droplet size increases from 30 microns to 80 microns, pressure oscillations are stronger with larger amplitudes. But an increase of the droplet size in the range of 80 microns to 140 microns indicates a reduction in the amplitudes of pressure oscillations. This trend is the same as the Hewitt criterion. The first tangential (1T) mode and the first longitudinal (1L) mode self-excited combustion instabilities are captured in the 60-micron and 80-micron cases. Abrupt spikes occur in the mass fraction of MMH and coincide with abrupt spikes in the mass fraction of NTO at the downstream regions just adjacent to the impinging points. Thus, local combustible high-dense mixtures are formed, which result in quasiconstant volume combustion and abrupt pressure spikes. The propagation and reflection of pressure waves in the chamber stimulate the combustion instability. When the droplet size is too small or too large, it is difficult to form local high-dense premixtures and combustion is stable in the chamber.


1999 ◽  
Vol 384 ◽  
pp. 59-91 ◽  
Author(s):  
S. M. LEE ◽  
I. S. KANG

A three-dimensional analysis is performed to investigate the effects of an electric field on the steady deformation and small-amplitude oscillation of a bubble in dielectric liquid. To deal with a general class of electric fields, an electric field near the bubble is approximately represented by the sum of a uniform field and a linear field. Analytical results have been obtained for steady deformation and modification of oscillation frequency by using the domain perturbation method with the angular momentum operator approach.It has been found that, to the first order, the steady shape of a bubble in an arbitrary electric field can be represented by a linear combination of a finite number of spherical harmonics Yml, where 0[les ]l[les ]4 and [mid ]m[mid ][les ]l. For the oscillation about the deformed steady shape, the overall frequency modification from the value of free oscillation about a spherical shape is obtained by considering two contributions separately: (i) that due to the deformed steady shape (indirect effect), and (ii) that due to the direct effect of an electric field. Both the direct and indirect effects of an electric field split the (2l+1)-fold degenerate frequency of Yml modes, in the case of free oscillation about a spherical shape, into different frequencies that depend on m. However, when the average is taken over the (2l+1) values of m, the frequency splitting due to the indirect effect via the deformed steady shape preserves the average value, while the splitting due to the direct effect of an electric field does not.The oscillation characteristics of a bubble in a uniform electric field under the negligible compressibility assumption are compared with those of a conducting drop in a uniform electric field. For axisymmetric oscillation modes, deforming the steady shape into a prolate spheroid has the same effect of decreasing the oscillation frequency in both the drop and the bubble. However, the electric field has different effects on the oscillation about a spherical shape. The oscillation frequency increases with the increase of electric field in the case of a bubble, while it decreases in the case of a drop. This fundamental difference comes from the fact that the electric field outside the bubble exerts a suppressive surface force while the electric field outside the conducting drop exerts a pulling force on the surface.


2014 ◽  
Vol 602-605 ◽  
pp. 674-679
Author(s):  
Hao Chen ◽  
Li Li Xie

This paper presents the 3D dynamic crack growth simulation of unreinforced baked brick shear wall by using particle discretization scheme finite element method (PDS-FEM), which is efficient and capable of computing bifurcation/branching in cracking. The technology of fast modelling of bricks and cements by applying VB script in AUTOCAD is illustrated briefly. The shear wall including mortar joints is modelled in detail. The model parameters are calibrated by using standard static tests. Since the computation cost is high in structural level fracture analysis, parallel computation technology is employed. Finally, with two-phase failure criterion of mortar under multi-dimension stress state, the performance of low and high loading speed is compared. The numerical results verify the availability of dynamic fracture analysis of masonry structure by using PDS-FEM.


2018 ◽  
Vol 8 (12) ◽  
pp. 2642 ◽  
Author(s):  
Wei Xu ◽  
Cancan Zhang ◽  
Xinyuan Ji ◽  
Hongyan Xing

The measurement of the atmospheric electric field is of great significance for the study of thunderstorm cloud charge models. Traditional electric field meters can only measure the vertical component of the atmospheric electric field, and thus it is difficult to invert the structure of the thunderstorm cloud. A three-dimensional atmospheric electric field meter was developed to simultaneously measure the horizontal and vertical components of the atmospheric electric field in this paper. The effective measurement linearly relates the measured electric field to the induced voltage, and the nonlinear equations of the three-dimensional atmospheric electric field and the thunderstorm cloud-charging model parameters were derived. The particle swarm optimization algorithm (PSO) and the three-dimensional atmospheric electric field were used to invert the thunderstorm clouds. Experimental observations of the three-dimensional electric field in a cloud during a thunderstorm were analyzed. Combined with the typical charged structure model, parameters such as the charge and relative distance of the thunderstorm cloud were determined. The results showed that the value of the inversion fitness function reached 0.7288, and the charge structure was even. The measurement of the three-dimensional atmospheric electric field provides a new means of observation for the study of atmospheric electricity.


1997 ◽  
Vol 473 ◽  
Author(s):  
Heng-Chih Lin ◽  
Edwin C. Kan ◽  
Toshiaki Yamanaka ◽  
Simon J. Fang ◽  
Kwame N. Eason ◽  
...  

ABSTRACTFor future CMOS GSI technology, Si/SiO2 interface micro-roughness becomes a non-negligible problem. Interface roughness causes fluctuations of the surface normal electric field, which, in turn, change the gate oxide Fowler-Nordheim tunneling behavior. In this research, we used a simple two-spheres model and a three-dimensional Laplace solver to simulate the electric field and the tunneling current in the oxide region. Our results show that both quantities are strong functions of roughness spatial wavelength, associated amplitude, and oxide thickness. We found that RMS roughness itself cannot fully characterize surface roughness and that roughness has a larger effect for thicker oxide in terms of surface electric field and tunneling behavior.


1986 ◽  
Vol 51 (5) ◽  
pp. 1001-1015 ◽  
Author(s):  
Ivan Fořt ◽  
Vladimír Rogalewicz ◽  
Miroslav Richter

The study describes simulation of the motion of bubbles in gas, dispersed by a mechanical impeller in a turbulent low-viscosity liquid flow. The model employs the Monte Carlo method and it is based both on the knowledge of the mean velocity field of mixed liquid (mean motion) and of the spatial distribution of turbulence intensity ( fluctuating motion) in the investigated system - a cylindrical tank with radial baffles at the wall and with a standard (Rushton) turbine impeller in the vessel axis. Motion of the liquid is then superimposed with that of the bubbles in a still environment (ascending motion). The computation of the simulation includes determination of the spatial distribution of the gas holds-up (volumetric concentrations) in the agitated charge as well as of the total gas hold-up system depending on the impeller size and its frequency of revolutions, on the volumetric gas flow rate and the physical properties of gas and liquid. As model parameters, both liquid velocity field and normal gas bubbles distribution characteristics are considered, assuming that the bubbles in the system do not coalesce.


Sign in / Sign up

Export Citation Format

Share Document