Purpose: To improve numerical simulation of the non-contact tonometry test by using Arbitrary Eulerian-Lagrangian deforming mesh in the coupling between computational fluid dynamics model of an air jet and finite element model of the human eye. Methods: Computational fluid dynamics model simulated impingement of the air puff and consisted of 25920 wedge6 elements and employed Spallart-Allmaras model to simulate capture turbulence of the air jet. The time span of the jet wais 30 ms and maximum Reynolds number

Giant plasmonic surface enhancement has been observed in gold coated micron sized inverse pyramids entrapping a gold nanoparticle. The amplification of both surface enhanced Raman and photoluminescence signals was found to be dependent on the diameter of trapped gold nanoparticle and around 50-fold enhancement was detected for 250nm diameter sample relatively to the 50nm one. Finite differential time domain simulations, performed to determine the near-field distribution in the structure, showed that when the nanoparticle protrudes into the hotspot zone of the void, coupling of electromagnetic field occurs and the plasmon-related near-field enhancement is concentrated into the close vicinity of the nanoparticle, mainly into the close gaps around the tangential points of the curved sphere and the flat pyramid surface. This results in a more than 15 times increase of the near-field intensity, compared to the empty void.

In this paper, the experimental protocol which we disclose is designed to identify the values for both the constant in the Kozeny/Carman model, which relates to the linear component of permeability, and the variable kinetic coefficient in the newly minted Q- modified Ergun model, which relates to the non-linear components of permeability, without involving any new theoretical development. Moreover, kinetic contributions to measured pressure gradient, which are not accounted for in some currently accepted empirical fluid flow equations, such as Poiseuille’s for flow in empty conduits and Kozeny/Carman for flow in packed conduits, but which nevertheless contribute to measured pressure drop and thus hamper the identification of the value of the constant relative to the laminar component, are captured and lumped together into a single variable kinetic parameter-the kinetic coefficient.

This paper proposed an explicit and simple representation of velocity fluctuation and the Reynolds stress tensor in terms of the mean velocity field. The proposed turbulence equations are closed. The proposed formulations reveal that the mean vorticity is the key source of producing turbulence. It is found that there are no velocity fluctuation and turbulence if there were no vorticity. As a natural consequence, the laminar- turbulence transition condition was obtained in a rational way.

This paper attempts to clarify an long-standing issue about the number of unknowns in the Reynolds-Averaged Navier-Stokes equations (RANS). This study shows that all perspectives regarding the numbers of unknowns in the RANS stem from the misinterpretation of the Reynolds stress tensor. The current literature consider that the Reynolds stress tensor has six unknown components; however, this study shows that the Reynolds stress tensor actually has only three unknown components, namely the three components of fluctuation velocity. This understanding might shed a light to understand the well-known closure problem of turbulence.

This work examines atomic and molecular signatures in laser-induced plasma in standard ambient temperature and pressure environments, including background contributions to the spectra that depend on the laser pulse-width. Investigations include solids, gases, and nano-particles. Abel inversions of measured line-of-sight data reveal insight into the radial plasma distribution. For nominal 6 nanosecond laser pulses and for pulse-energies in the range of 100 to 800 mJ, expansion dynamics and turbulence due to shock phenomena are elucidated to address local equilibrium details that are frequently assumed in spatially averaged emission spectroscopy. Chemical equilibrium computations reveal temperature dependence of selected plasma species. Specific interests include atomic hydrogen (H) and cyanide (CN). Atomic hydrogen spectra indicate axisymmetric shell structures and isentropic expansion of the plasma kernel over and above the usual shockwave. The recombination radiation of CN emanates within the first 100 nanoseconds for laser-induced breakdown in a 1:1 CO2:N2 gas mixture when using nanosecond laser pulses to create the micro-plasma. The micro-plasma is generated using 1064 nm, 150 mJ, 6 ns Q-switched Nd:YAG laser radiation. Measurements of the optical emission spectra utilize a 0.64 m Czerny-Turner type spectrometer and an intensified charge-coupled device.

This paper is to summarize the involvement of the stress energy tensor in the study of fluid mechanics. In the first part we will see the implication that carries the stress energy tensor in the framework of general relativity. In the second part, we will study the stress energy tensor under the mechanics of perfect fluids, allowing us to lead third party in the case of Newtonian fluids, and in the last part we will see that it is possible to define space-time as a no-Newtonian fluids.

This work communicates laser-plasma experiments in a gaseous mixture of hydrogen and nitrogen. Time-resolved spectroscopy measures the first four Balmer series hydrogen lines together with selected neutral and ionized nitrogen lines. Optical breakdown plasma is generated in a 1:1 hydrogen:nitrogen mixture at ambient temperature and 0.27-atm pressure. Time-resolved spectroscopy records emitted radiation with spatial resolution along the slit height for the H_α, H_β, H_γ, and H_δ lines. For 13 selected time delays from 0.25 µs to 3.25 µs and 0.025 µ gate-widths, micro-plasma diagnostics is evaluated. Of interest are the peak-separation and width of H_δ and width of H_γ for electron densities in the range of 0.1 to 1.0×10¹⁷ cm^-3, and comparisons with H_β and H_α diagnostics. Integral inversions interrogate spatial distributions of the plasma expansion. Applications include laboratory and stellar astrophysics plasma diagnosis.

We investigate the flow structure and dynamics of moderate-Rayleigh-number ($Ra$) thermal convection in a two-dimensional inclined porous layer. Direct numerical simulations (DNS) confirm the emergence of $\mathit{O}(1)$ aspect-ratio large-scale convective rolls, with one `natural' roll rotating in the counterclockwise direction and one `antinatural' roll rotating in the clockwise direction. As the inclination angle $\phi$ is increased, the background mean shear flow intensifies the natural-roll motion, while suppressing the antinatural-roll motion. Moreover, our DNS reveal---for the first time in single-species porous medium convection---the existence of \emph{spatially-localized} convective states at large $\phi$, which we suggest are enabled by subcritical instability of the base state at sufficiently large inclination angles. To better understand the physics of inclined porous medium convection at different $\phi$, we numerically compute steady convective solutions using Newton iteration and then perform secondary stability analysis of these nonlinear states using Floquet theory. Our analysis indicates that the inclination of the porous layer stabilizes the boundary layers of the natural roll, but intensifies the boundary-layer instability of the antinatural roll. These results facilitate physical understanding of the large-scale cellular flows observed in the DNS at different values of $\phi$.

This article attempts to formulate a mathematical model for a potential explanation regarding the unavoidable impact of a rigid body's peculiar shape on the seamless flow over it. The solid body completely immersed in a Newtonian fluid and respectively has a relative open circuit flow on it will typically experience various observable phenomena like flow separation, flow transition, down-wash, stalling at the higher angle of attack, stalling velocity and how cambered airfoil can typically generate lift at a zero incidence angle. This article respectively represents an understanding of the laminar flow over a rigid body's external surface with due respect to its distinctive shape and size. This working paper formulates a more realistic and simplified mathematical model for open circuit laminar flow over a body, based on the historical data of aerodynamics and theoretical mechanics. This is intended to properly estimate forces on the continuous surface of the body in a laminar flow, to properly explain, understand and predict mentioned phenomena. Most of all the mechanism of streamline formation and its deformation with due regards to flow, shape and size of the body in an open-circuit laminar are formulated mathematically to enhance better design theory which can reduce experimentation while designing a streamlined body.

Recently a new class of experiments in the field of surface gravity waves has emerged. These experiments serve as demonstration of an analogy to a broad variety of phenomena in optics and quantum mechanics. In particular, experiments involving Airy water-wave packets were carried out. The Airy wave packets have attracted tremendous attention in optics and quantum mechanics owing to their unique properties, spanning from an ability to propagate along parabolic trajectories without spreading, and to accumulating a phase that scales with the cubic power of time. Non-dispersive Cosine-Gauss wave packets and self-similar Hermite-Gauss wave packets, also well known in the field of optics and quantum mechanics, were recently studied using surface gravity waves as well. These wave packets demonstrated self-healing properties in water wave pulses as well, preserving their width despite being dispersive. Finally, this new approach also allows to observe diffractive focusing from a temporal slit with finite width.

The article presents the results of tests of minor head losses through PVC and PP elbows for a flow of water and mixtures of water and sand with grain sizes of up to 0.5 mm and concentrations of 5.6 g∙L-1, 10.84 g∙L-1, and 15.73 g∙L-1. The tests were carried out at variable flow velocities for three elbow diameters of 63, 75, and 90 mm. The flow rate, pressure difference in the tested cross-sections, and temperature of the fluids were measured and automatically recorded. The results of the measurements were used to develop mathematical models for determining the minor head loss coefficient as a function of elbow diameter, sand concentration in the liquid, and Reynolds number. The mathematical model was developed by cross validation. It was shown that when the concentration of sand in the liquid was increased by 1.0 g∙L-1, the coefficient of minor head loss through the elbows increased, in the Reynolds number range of 4.6∙104 − 2.1∙105, by 0.3−0.01% for PP63, 0.6−0.03 % for PP75, 1.1−0.06 % for PP90, 0.8−0.01 % for PVC63, 0.8−0.02 % for PVC75, and 0.9−0.04 % for PVC90. An increase in Re from 5∙104 to 2∙106 for elbows with diameters of 63, 75 and 90 mm caused a 7.3 %, 6.8 %, and 6.0 % decrease in the minor head loss coefficient, respectively.

The motivation of this work is to investigate experimentally the influence of non-thermal plasma (NTP) application on the reaction kinetics of atmospheric pressure steam gasification of charcoal using a thermostatically controlled drop tube reactor. A gliding-arc generator provides about 1 kW electrical power NTP. For comparison thermal gasification is investigated under comparable flow and specific energy input conditions providing additional heat to the steam. Optical temperature measurement 20 cm flow down of the NTP zone is utilized to characterize the specific enthalpy of the reactive flow. The composition of produced syngas is measured by a gas analyzer and used for the calculation of gas flow rates. The results show a NTP-enhancement on the production of individual syngas components (H2, CO, CH4), especially on hydrogen production by around 39%. The syngas-based carbon conversion and hydrogen release are calculated from the carbon and hydrogen balance between the correspondent content in syngas and in the feedstock. The NTP promoted the carbon conversion and hydrogen release by 25% and 31%, respectively. The first-order reaction kinetics are determined by data-fitting in an Arrhenius diagram. The plasma enhanced the reaction rate coefficients by 27%. Based on experimental results and other literature, possible plasma-induced reactions are proposed.

This work aimed to characterize a DBD plasma equipment through optical and electrical measurements, seeking to obtain a greater knowledge of the plasma production process and how it behaves through the adopted parameters, such as frequency and voltage applied between electrodes, at a fixed distance of 1.7 mm. In order to measure them, three different characterization techniques were applied. The first method was the Lissajous figures, a technique quite effective for a complete electrical characterization of DBD equipment. The second technique used was the Optical Emission Spectroscopy, a tool used for the diagnosis of plasma, being it possible to identify the excited species produced in filamentary and diffuse discharge in the plasma. And finally, the triple Langmuir probe technique was used to obtain the electron temperature and electron density. Based on this study, it was possible to identify the equipment efficiency in different regimes. The electron temperature measurement for both systems analyzed were 27.96 eV and 20.69 eV to the filamentary and diffuse regimes, respectively. The density of electrons number to these regimes were 1.09 × 10²¹ m⁻³ and 1.56 × 10²¹ m⁻³.

Deformation of the surface of diamagnetic liquid by magnetic field is called the “Moses Effect”. Physics and applications of the direct and inverse Moses effects are reviewed. Experimental techniques enabling visualization of the effects are surveyed. Impact of magnetic field on micro- and macroscopic properties of liquids is addressed. Influence of the surface tension on the shape of the near-surface dip formed in a diamagnetic liquid by magnetic field is reported. Floating of diamagnetic bodies driven by the Moses effect is treated. The effect of the “magnetic memory of water” in its relation to the Moses Effect is discussed. The dynamics of self-healing of near-surface dips due to the Moses Effect is considered.

This paper presents a hydrodynamic numerical model for a pitching wave energy converter (WEC). The model uses potential wave theory and is based on Cummins' equation, with nonlinear hydrostatic restoring stiffness and excitation forces based on instantaneous body position and water surface elevation. The numerical model can include non-linear forces, like quadratic drag, power-take-off and other forces that may account for unknown viscous effects observed in experiments. The paper discusses the applicability and limitations of the code and presents the cases where assumptions and simplifications can be made. The goal is to conclude on the simplest, yet accurate, version of the model by evaluating its accuracy using experimental data. The case study for validation is Floating Power Plant's (FPP's) WEC [1]. In the full-scale commercial project, FPP's device consists of a semisubmersible platform hosting a wind turbine (5-8MW) and 4 WECs, each one connected to the platform by a rotation shaft. Due to the configuration of the platform, strong interactions occur between the WECs and the structure, as they are very closely spaced. In order to validate a numerical model able to simulate these hydrodynamic interactions, wave basin experiments with a similar but simplified setup were performed.

Deformation flows are flows incorporating shear, strain and rotational components. These flows are ubiquitous in the geophysical flows, such as the ocean and atmosphere. They appear near almost any salience, such as isolated coherent structures (vortices and jets), various fixed obstacles (submerged obstacles, continental boundaries). Fluid structures subject to such deformation flows may exhibit drastic changes in motion. In this review paper, we focus on the motion of a small number of coherent vortices embedded in deformation flows. Problems involving isolated one and two vortices are addressed. When considering a single-vortex problem, the main focus is on the evolution of the vortex boundary and its influence on the passive scalar motion. Two vortex problems are addressed with the use of point vortex models, and the resulting stirring patterns of neighbouring scalars are studied by a combination of numerical and analytical methods from the dynamical system theory. Many dynamical effects are reviewed with emphasis on the emergence of chaotic motion of the vortex phase trajectories and the scalars in their immediate vicinity.

Among the most important variables in the design of falling film microreactors (FFMRs) is the liquid film thickness as well as the gas/liquid interfacial area, which dictate the mass and heat transfer rates. In a previous work conducted in our lab the characteristics of a free-falling Newtonian liquid film have been studied and appropriate correlations have been proposed. In this work the geometrical characteristics of a non-Newtonian shear thinning liquid, flowing in an inclined open microchannel, have been experimentally investigated and design correlations that can predict with reasonable accuracy the features of a FFMR have been proposed. The test section used was an open μ-channel with square cross section (Wo=1200 μm) made of brass which can be set to various inclination angles. The liquid film characteristics were measured by a non-intrusive technique that is based on the features of a μ-PIV system. Relevant CFD simulations revealed that the volume average dynamic viscosity over the flow domain is practically the same as the corresponding asymptotic viscosity value, which can thus be used in proposed the design equations. Finally, a generalized algorithm for the design of FFMRs, containing non-Newtonian shear thinning liquids, is suggested.

Homogeneous and heterogeneous chemical reactions in partially ionized magneto-nano-liquid are investigated theoretically using finite element method (FEM). The effects of ions and electrons collisions on the transport of heat and mass are analyzed for both the cases of heterogeneous and homogeneous chemical reactions. The simultaneous effects of dispersion of nanosized particles in partially ionized nano-liquid in the presence of magnetic field are also investigated. Through numerical experiments, it is noted that the temperature of partially ionized nano-liquid increases when electrons collision rate and ion collisions are increased. The transport rate of reacting species decreases when heterogeneous and homogeneous chemical reactions strengths are increased. It is also observed that the effect of electron collisions on the flow in y-direction is the same to that of ion collisions on the flow in y-direction. Homogeneous and heterogeneous chemical reactions have similar effects on concentration of chemically reacting species in qualitative sense. However, in quantitative sense, homogeneous chemical reaction has more significant effect on the concentration reacting species as compared to heterogeneous chemical reaction.

This is a review of efforts to discover action principles for Hydrodynamics and Thermodynamics. The work began with a determined attempt to create sources for Einstein's equations that describe continuous matter - as opposed to point particles. This requires a relativistic action principle. This report is mainly about (non relativistic) Adiabatic Thermodynamics, non equilibrium thermodynamics based on the method of Gibbs. The main relations of classical thermodynamics can be derived from an action principle, in the classical global form and in the local form advocated by Callen. The entropy is defined as the negative derivative of the free energy (free energy density) with respect to the volume (density). This is integrated with an action principle for hydrodynamics that was known to Lagrange in 1760. It is strongly limited to irrotational velocity fields, nevertheless it provides the needed Lagrangian for a large part of thermodynamics. This is the subject of the first 3 sections. It includes a brief review of a new theory of mixtures. The most difficult part of the program has been the discovery of Conservative Hydrodynamics with the required four independent degrees of freedom, with general flows. We give an account of the genesis and the structure of this theory, then a tour of recent applications.

The Landau damping effect was observed in collisionless plasma, as a microscopic resonant mechanism between electromagnetic radiation and the collective modes. In this paper we demonstrate the occurrence of the Landau damping at macroscopic scale in the interaction between water waves and anharmonic lattice of magnetic buoys. By coupling the Navier-Stokes equations for incompressible fluid with the nonlinear dynamics of an anharmonic magnetic lattice we obtain a resonant transfer of momentum and energy between the two systems. The velocity of the flow is obtained in the Stokes approximation with Basset type of drag force. The dynamics of the buoys is calculated in the surfactant approximation for a specific frequency, then we use Fourier analysis to obtain the general time variable interaction. After involving an integral Dirichlet transform we obtain the time dependent expression of the drag force, the interaction waves-lattice with a new term in the form of a Caputo fractional derivative. We compare the results of the model with experiments performed in a wave tank with free floating magnetic buoys under the action of small amplitude gravitational waves. This configuration can be applied in studies for the attenuation with resonant damping of rogue waves, storms or tsunamis.

Filament wound composite pipes are frequently used in the field were transmission of high pressured chemical fluids, disposal of industrial wastes, oil and natural gas transmission takes place. In oil and gas industry, the pipelines transporting heavy crude oil are subjected to variable pressure waves causing fluctuating stress levels in the pipes. Computational Fluid Dynamics Analysis was performed using Ansys 15.0 Fluent software to study the effects of these pressure waves on some specified joints in the pipes. Depending on the type of heavy crude oil being used, the flow behavior indicated a considerable degree of stress levels in certain connecting joints, causing the joints to become weak over a prolonged period of use. In this research comparison of various pipe joints was done by using different material and the output result of the stress levels of the pipe joints were checked so that the life of the pipe joints can be optimized by the change of material.

Plasma being the fourth and most abundant form of matter extensively exists in the universe in the inter-galactic regions. It provides an electrically neutral medium of unbound negative and positive charged particles, which has been produced by subjecting air and various other gaseous mixtures to strengthen the electromagnetic field and by heating compressed air or inert gasses for creating negative and positive charged particles known as ions. Nowadays, many researchers are paying attention to the formation of artificial Plasma and its potential benefits for mankind. The literature is sparsely populated with the applications of Plasma. This paper presents specific methods of generation and applications of Plasma, which benefits humankind in various fields, such as in electrical, mechanical, chemical and medical fields. These applications include hydrogen production from alcohol, copper bonding, semiconductor processing, surface treatment, Plasma polymerization, coating, Plasma display panels, antenna beam forming, nanotechnology, Plasma Torch, Plasma pencils, low-current non-thermal Plasmatron, treatment of prostate cancer, Plasma source ion implantation, cutting by Plasma, Plasma etching, pollution control, neutralization of liquid radioactive waste, etc. Resultantly, worth of Plasma technology in the medical industry is increasing exponentially that is closing the gap between its benefits and cost of equipment used for generating and controlling it.

This paper attempts to clarify an issue regarding the lasting unsolved problem of turbulence, namely the closure problem. This study shows that all perspectives regarding the numbers of un- known quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor. The current literatures have a consensus that the Reynolds stress tensor has six unknowns; however, this study shows that the Reynolds stress tensor actually has only three ones, namely the three components of fluctuation velocity. With this new understanding, the closed turbulence equations for incompressible flows are proposed.

Wetting of deformable substrates has gained significant interest over the past decade due to its extensive applications and uses. This interest has developed due to technological advances which are able to capture interfacial behavior taking place when a liquid droplet is placed on a deformable substrate. Researchers have developed different theories to explain processes taking place in the process of wetting of deformable/soft substrates. For the scope of this review, we will consider the fluid to be Newtonian, partially wetting, and surface forces are incorporated with the help of disjoining/conjoining pressure acting in the vicinity of the apparent, three-phase contact line. The following subjects are briefly reviewed: (i) Equilibrium of droplets on soft substrates. It is shown that properties of the disjoining/conjoining pressure isotherm and properties of the deformable substrate determine both the shape of the liquid droplet and deformation of the substrate; (ii) Equilibrium conditions of droplets on deformable substrates. It is shown that for a droplet to be at equilibrium on a deformable substrate under consideration, Jacobi’s sufficient condition is satisfied; (iii) Hysteresis of contact angle of sessile droplets on deformable substrates. It is shown that as the elasticity of the deformable substrate is increased, both advancing and receding contact angles are reduced.

In this paper, an analytical investigation of the combined effect of through flow and magnetic field on the convective instability in an electrically conducting fluid layer, bounded in a Hele-Shaw cell is presented within the context of linear stability theory. The Galarkin method is utilized to solve the eigenvalue problem. The outcome of the important parameters on the stability of the system is examined analytically as well as graphically. It is observed that the through flow and magnetic field have both stabilizing effects, while the Hele-Shaw number has destabilizing effect on the stability of system. It is also found that the oscillatory mode of convection possible only when the magnetic Prandtl number takes the values less than unity.

The dynamics of a spherical particle in an asymmetric serpentine is studied by finite element method (FEM) simulations in a physically unconstrained system. The two-way coupled time dependent solutions illustrate the path of the particle along a curve where a secondary flow (Dean flow) has developed. The simulated conditions were adjusted to match those of an experiment for which particles were focused under inertial focusing conditions. The obtained rotational modes allowed to infer the influence of the local flow around the particle. We propose a new approach to find the decoupled secondary flow contribution employing a quasi-Stokes flow.

Widely used in hydraulics, the Colebrook equation for flow friction relates implicitly to the input parameters; the Reynolds number, Re and the relative roughness of inner pipe surface, ε/D with the output unknown parameter; the flow friction factor, λ; λ=f(λ, Re, ε/D). In this paper, a few explicit approximations to the Colebrook equation; λ≈f(Re, ε/D), are generated using the ability of artificial intelligence to make inner patterns to connect input and output parameters in explicit way not knowing their nature or the physical law that connects them, but only knowing raw numbers, {Re, ε/D}→{λ}. The fact that the used genetic programming tool does not know the structure of the Colebrook equation which is based on computationally expensive logarithmic law, is used to obtain better structure of the approximations which is less demanding for calculation but also enough accurate. All generated approximations are with low computational cost because they contain a limited number of logarithmic forms used although for normalization of input parameters or for acceleration, but they are also sufficiently accurate. The relative error regarding the friction factor λ, in best case is up to 0.13% with only two logarithmic forms used. As the second logarithm can be accurately approximated by the Padé approximation, practically the same error is obtained also using only one logarithm.

The paper is considering of the basic differential equations of hydrodynamics: the equation of continuity and motion. On the simplest example, it is shown that the equation of continuity in a system with the equation of motion leads to contradictions and erroneous results of modeling. A more correct form of the continuity equation is described. It is shown that the equations of motion can be written in the form of complete differentials. Three possible integral forms of the equations of motion are presented. As a conclusion, the existence and smoothness of the solution of the Navier-Stokes equations are considered.

Due to failure of the continuum hypothesis for higher Knudsen numbers, rarefied gases and microflows of gases are particularly difficult to model. Macroscopic transport equations compete with particle methods, such as DSMC to find accurate solutions in the rarefied gas regime. Due to growing interest in micro flow applications, such as micro fuel cells, it is important to model and understand evaporation in this flow regime. Here, evaporation boundary conditions for the R13 equations, which are macroscopic transport equations with applicability in the rarefied gas regime, are derived. The new equations utilize Onsager relations, linear relations between thermodynamic fluxes and forces, with constant coefficients, that need to be determined. For this, the boundary conditions are fitted to DSMC data and compared to other R13 boundary conditions from kinetic theory and Navier-Stokes-Fourier (NSF) solutions for two one-dimensional steady-state problems. Overall, the suggested fittings of the new phenomenological boundary conditions show better agreement to DSMC than the alternative kinetic theory evaporation boundary conditions for R13. Furthermore, the new evaporation boundary conditions for R13 are implemented in a code for the numerical solution of complex, two-dimensional geometries and compared to NSF solutions. Different flow patterns between R13 and NSF for higher Knudsen numbers are observed.

This work examines atomic and molecular signatures in laser-induced plasma in standard ambient temperature and pressure environments, including background contributions to the spectra that depend on the laser pulse-width. Investigations include solids, gases, and nano-particles. Abel inversions of measured line-of-sight data reveal insight into the radial plasma distribution. For nominal 6 nanosecond laser pulses and for pulse-energies in the range of 100 to 800 milli-Joules, expansion dynamics and turbulence due to shock phenomena are elucidated to address local equilibrium details that are frequently assumed in spatially averaged emission spectroscopy. Chemical equilibrium computations reveal temperature dependence of selected plasma species. Specific interests include atomic hydrogen (H) and cyanide (CN). The atomic H spectra, collected following optical breakdown in ultra-high-pure hydrogen and 9:1 mixtures of ultra-pure hydrogen and nitrogen gases, indicate spherical shell structures and isentropic expansion of the plasma kernel over and above the usual shockwave. The recombination radiation of CN emanates within the first 100 nanoseconds for laser-induced breakdown in a 1:1 CO₂:N₂ gas mixture when using nanosecond laser pulses to create the micro-plasma. The micro-plasma is generated using 1064 nm, 150 mJ, 6 ns Q-switched Nd:YAG laser radiation. Measurements of the optical emission spectra utilize a 0.64 m Czerny-Turner type spectrometer and an intensified charge-coupled device.

The present study provides molecular insight into the mechanisms underlying energy dissipation and lubrication of a smooth contact lubricated by an ionic liquid. We have performed normal and lateral force measurements with a surface forces apparatus and by colloidal probe atomic force microscopy on following model systems: 1-ethyl-3-methyl imidazolium bis-(trifluoro-methylsulfonyl)imide, in dry state and in equilibrium with ambient (humid) air; the surface was either bare mica or functionalized with a polymer brush. The velocity-dependence of the friction force reveals two different regimes of lubrication, boundary-film lubrication, with distinct characteristics for each model system, and fluid-film lubrication above a transition velocity V*. The underlying mechanisms of energy dissipation are evaluated with molecular models for stress-activated slip and flow, respectively. The stress-activated slip assumes that two boundary layers (composed of ions/water strongly adsorbed to the surface) slide pass each; the bond dynamics and the strength of the interaction at the slip plane dictate the change in friction -decreasing, increasing or remaining constant- with velocity in the boundary-film lubrication regime. Above a transition velocity V*, friction monotonically increases with velocity in the three model systems. Here, layers of ions slide past each (“flow”) under a shear stress and friction depends on a shear-activation volume that is significantly affected by confinement. The proposed friction model provides a molecular perspective of the lubrication of smooth contacts by ionic liquids and allows identifying the physical parameters that control friction.

Based on author's previous work [Sun, B. The Reynolds Navier-Stokes Turbulence Equations of Incompressible Flow Are Closed Rather Than Unclosed. Preprints 2018, 2018060461 (doi: 10.20944/preprints201806.0461.v1)], this paper proposed an explicit representation of velocity fluctuation and formulated the Reynolds stress tensor in terms of the mean velocity field. The proposed closed Reynolds Navier-Stokes turbulence formulations reveal that the mean vorticity is the key source of producing turbulence.

This paper shown that turbulence closure problem is not an issue at all. All mistakes in the literature regarding the numbers of unknown quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor, i.e., all literatures have stated that the symmetric Reynolds stress tensor has six unknowns; however, it actually has only three unknowns, i.e., the three components of fluctuation velocity. We shown the integral-differential equations of the Reynolds mean and fluctuation equations have exactly eight equations, which equal to the numbers of quantities in total, namely, three components of mean velocity, three components of fluctuation velocity, one mean pressure and one fluctuation pressure. That is why we claim in this paper, that the Reynolds Navier-Stokes turbulence equations of incompressible flow are closed rather than unclosed. This study may help to solve the puzzle that has eluded scientists and mathematicians for centuries.

Pneumatic nebulizers (as variations based on the Collison nebulizer) have been widely used for producing fine aerosol droplets from a liquid material. The basic working principle of those nebulizers has been qualitatively described as utilization of the negative pressure associated with an expanding gas jet to syphon liquid into the jet stream, then to blow and shear into liquid sheets, filaments, and eventually droplets. Detailed quantitative analysis based on fluid mechanics theory is desirable, to gain in-depth understanding of the liquid aspiration mechanism among other aspects of the Collison nebulizer behaviour. The purpose of present work is to investigate the nature of negative pressure distribution associated with compressible gas jet flow in the Collison nebulizer by a computational fluid dynamics (CFD) analysis, using an OpenFOAM® compressible flow solver. The value of the negative pressure associated with a gas jet flow is examined by varying geometric parameters of the jet expansion channel adjacent to the outlet of jet orifice. Such an analysis can provide valuable insights into fundamental mechanisms for liquid aspiration, helpful for designing improved pneumatic atomizer in the Aerosol Jet® direct-write system for micro-feature, high-aspect-ratio material deposition

Anomalous waves and rogue events are closely associated with irregularities and unexpected events occurring at various levels of physics, such as in optics, in oceans and in the atmosphere. Mathematical modeling of rogue waves is a highly actual field of research, which has evolved over the last decades into a specialized part of mathematical physics. The applications of the mathematical models for rogue events is directly relevant to technology development for prediction of rogue ocean waves, and for signal processing in quantum units. In this survey, a comprehensive perspective of the most recent developments in methods for representing rogue waves is given, along with discussion of the devised and forms and solutions. The standard nonlinear Schrödinger equation, the Hirota equation, the MMT equation and further to other models are discussed, and their properties highlighted. This survey shows that the most recent advancement in modeling rogue waves give models which can be used to establish methods for prediction of rogue waves at open seas, which is important for the safety and activity of marine vessels and installations. The study further puts emphasis on the difference between the methods, and how the resulting models form a basis for representing rogue waves in various forms, solitary or with a wave-background. This review has also a pedagogic component directed towards students and interested non-experts, and forms a complete survey of the most conventional and emerging methods published until recently

The rising of an oil drop in a non-Newtonian viscous solution is studied experimentally. In this case, the shape of the ascending drop is strongly affected by the non-Newtonian properties of the surrounding liquid. We found that the so-called velocity discontinuity phenomena is observed for drops larger than a certain critical size. Beyond the critical velocity, the formation of a long tail is observed, from which small droplets are continuously emitted. We determined that the fragmentation of the tail results mainly from the effect of capillary effects. We explore the idea of using this configuration as a new encapsulation technique, where the size and frequency of droplets can be well predicted.

In this study, an experimentally validated computational model was developed to investigate the hydrodynamics in a rotor-stator vortex RVR agglomeration reactor having a rotating disc at the centre with two shrouded outer plates. A numerical simulation was performed using a simplified form of the reactor geometry to compute the 3D flow field in batch mode operations. Thereafter, the model was validated using data from a 2D Particle Image Velocimetry (PIV) flow analysis performed during the design of the reactor. Using different operating speeds—70, 90, 110 and 130 rpm, the flow fields were computed numerically followed by a comprehensive data analysis. The simulation results showed separated boundary layers on the rotating disc and the stator. The flow field within the reactor is characterized by a rotational plane circular forced vortex flow in which the streamlines are concentric circles with a rotational vortex. Overall, the results of the numerical simulation demonstrate a fairly good agreement between the CFD model and the experimental data as well as the available theoretical predictions. The swirl ratio β was found to be approximately 0.4044, 0.4038, 0.4044 and 0.4043 for operating speeds of N=70, 90, 110 and 130 rpm respectively. In terms of the spatial distribution, the turbulence intensity and kinetic energy are concentrated on the outer region of the reactor while the axial velocity showed a decreasing intensity towards the shroud. However, a comparison of the CFD and experimental predictions of the tangential velocity and the vorticity amplitude profiles shows that these parameters were under-predicted by the experimental analysis which could be attributed to some of the experimental limitations rather than the robustness of the CFD model or numerical code.

We describe results from a Monte-Carlo simulation of Bell-CHSH type correlations in hydrodynamic walkers. We study feasibility of a real life walker test with relevant hydrodynamic parametric ranges. We observe the generic formation of pairs of walkers strongly anti-correlated both in position and momentum. With this source of entangled walkers, we model the insertion of 2 pins in the bath as a notion of measure, akin to the polarizers of photonic Bell tests. This insertion of pins, either static or dynamic, introduces 2 weak field signals. Each field has the physical form of a standing wave Bessel hat, representing the non-local (field mediated) influences of the measure on the walkers. With this representation of the measure, we develop protocol for a Bell game with actual hydrodynamic walkers. We model both static and dynamic insertion of pins in the walker bath. Static pins give us numerical S > 2, as a permissible Bell violation for a non-local (field based) effect. Dynamic insertion of the pins, however, leads to causal space separation of the two arms. We observe the again expected S ≤ 2. We argue for the hydrodynamic implementation and observation of these effects as a walker visualization of Bell inequalities.

This paper presents the results of aerodynamical performance research focused on maintaining the thermal comfort and increasing the energy efficiency of a typical social housing unit located in a high-density urban area. Bioclimatic design strategies are used to develop a sustainable and economic technology in existing housing clusters in Mexico City. A full-scale prototype, built on campus facilities, was used to study the flow conditions around the design. All scaled prototypes implement similar criterion. Furthermore, a scaled prototype is evaluated within a low speed wind tunnel installation. Additionally, numerical simulations were performed at transient state based on previous physical measurements and historical local climatic conditions to determine preferable modifications.

A new formulation for a proposed solution to the 3D Navier-Stokes Equations in cylindrical co-ordinates coupled to the continuity and level set convection equation is presented in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow and a connection between the level set function and composite velocity vector for the additive solution is shown. For the case of a vertical tube configuration with small inclination angle, results are obtained for the level set function defining the interface in both the radial and azimuthal directions. It is found that the curvature dependent part of the problem alone induces sinusoidal azimuthal interfacial waves wheras when the curvature is small oscillating radial interfacial waves occur. The implications of two extremes indicate the importance of looking at the full equations including curvature.

The feasibility of using full ion kinetics, instead of gyrokinetics, in simulating low-frequency Ion-Temperature-Gradient (ITG) instabilities in tokamaks has recently been demonstrated by Sturdevant et al. [Physics of Plasmas 24, 081207 (2017)]. In that work, a variational integrator was developed to integrate the full orbits of ions in toroidal geometry, which proved to be accurate in capturing both the short-time scale cyclotron motion and long time scale drift motion. The present work extends that work in three aspects. First, we implement a relatively simple full orbit integrator, the Boris integrator, and demonstrate that the accuracy of this integrator is also sufficient for simulation of ITG instabilities. Second, the equilibrium magnetic configuration is extended to general toroidal configuration specified numerically, enabling simulation of realistic equilibria reconstructed from tokamak experiments. Third, we extend that work to the nonlinear regime and investigate the nonlinear saturation of ITG instabilities. To verify the new numerical implementation of the orbit integrator and magnetic configuration, the linear electrostatic ITG frequency and growth rate are compared with those given in Sturdevant's work and good agreement is found.

We have recently solved the tumbling-snake model for concentrated polymer solutions and entangled melts in the presence of both steady-state and transient shear and uniaxial elongational flows, supplemented by a variable link tension coefficient. Here, we provide the transient and stationary solutions of the tumbling-snake model under biaxial elongation both analytically, for small and large elongation rates, and via Brownian dynamics simulations, for the case of planar elongational flow over a wide range of rates, times, and the model parameters. We show that both the steady-state and transient first planar viscosity predictions are similar to their uniaxial counterparts, in accord with recent experimental data. The second planar viscosity seems to behave in all aspects similarly to the shear viscosity, if shear rate is replaced by elongation rate.

A four-parameter kinematic model for the position of a fluid parcel in a time-varying ellipse is introduced. For any ellipse advected by an arbitrary linear two-dimensional flow, the rates of change of the ellipse parameters are uniquely determined by the four parameters of the velocity gradient matrix, and vice-versa. This result, termed ellipse/flow equivalence, provides a stronger version of the well-known result that a linear velocity field maps an ellipse into another ellipse. Moreover, ellipse/flow equivalence is shown to be a manifestation of Stokes' theorem. This is done by deriving a matrix-valued relationship, called the geometric Stokes' theorem, that involves a spatial integral over the velocity gradient tensor, thus accounting for the two strain terms in addition to the divergence and vorticity. General expressions for various physical properties of an elliptical ring of fluid are also derived. The ellipse kinetic energy is found to be composed of three portions, associated respectively with the circulation, the rate of change of the moment of inertia, and the variance of parcel angular velocity around the ellipse. A particular innovation is the use of four matrices, termed the IJKL basis, that greatly facilitate the required calculations.

In this paper Brinkman type fluid over an infinite plate between side walls is being investigated. The flow is generated by oscillating shear stress of the bottom plate and the solutions are obtained by using Fourier integral transformation. The obtained results are presented in steady and transient states for both sin and cos shear stresses. The general solutions are reduced to some special cases corresponding, namely to the Brinkman type fluid over an infinite plate and flow of a Newtonian viscous fluid. Graphical illustrations are carried out to have in depth analysis of the involved physical parameters

This paper investigates the tank drainage problem of an isothermal, unsteady, incompressible electrically conducting Power law fluid. Analytic solution have been obtained from governing continuity and momentum equations subject to appropriate boundary conditions by using Perturbation method. The Power law fluid model solution without MHD is retrieved from this proposed model on substitution . Declaration on behalf of velocity profile, volume flux, average velocity, connection of time with respect to length of the tank and requirement of time for whole drainage of fluid are acquired. Special effects of numerous emerging parameter’s on velocity profile vz and depth of the fluid in the tank are graphically presented. Keywords: Tank drainage, Power law MHD fluid, Analytical solution.

In our previous works we have proposed design equations that can predict with reasonable accuracy the transition point from homogeneous to heterogeneous regime as well as the gas holdup and the mean Sauter diameter at the homogeneous regime. The validity of the proposed correlations was checked with data obtained using different geometrical configurations and several Newtonian and non-Newtonian liquids as well as the addition of surfactants. However, in all the experiments the gas phase was atmospheric air. This work investigates the effect of gas phase properties by conducting experiments employing various gases (i.e., air, CO₂, He) that cover a wide range of physical property values. Experiments revealed that only the use of low-density gas (He) has a measurable effect on bubble column performance. More precisely, when the low-density gas (He) is employed, the transition point shifts to higher gas flow rates and the gas holdup decreases, a fact attributed to the lower momentum force exerted by the gas. In view of the new data, the proposed correlations have been slightly modified to include the effect of gas phase properties and it is found that they can predict the aforementioned quantities with an accuracy of ±15%. It has been also proved that CFD simulations are an accurate means for assessing the flow characteristics inside a bubble column.

In this work the efficiency of a new μ-mixer design is investigated. As in this type of devices the Reynolds number is low, mixing is diffusion dominated and it can be enhanced by creating secondary flows. In this study we propose the introduction of helical inserts into a straight tube to create swirling flow. The influence of the insert’s geometrical parameters (pitch and length of the propeller blades) and of the Reynolds number on the mixing efficiency and on the pressure drop are numerically investigated. The mixing efficiency of the device is assessed by calculating a number, i.e. the Index of Mixing Efficiency that quantifies the uniformity of concentration at the outlet of the device. The influence of the design parameters on the mixing efficiency is assessed by performing a series of “computational” experiments, in which the values of the parameter are selected using DOE methodology. Finally using the numerical data, appropriate design equations are formulated, which, for given values ​​of the design parameters, can estimate with reasonable accuracy both the mixing efficiency and the pressure drop of the proposed mixing device.

The Gurney flap (GF) is a miniature lift-enhancement device and is usually mounted at the trailing edge of an airfoil. The GF and has been successfully applied to isolated airfoils, multi-element airfoils, and aircraft wings as well as helicopter rotor due to its attractive features of simplicity, cost-effective and separation control. The GF also has aroused the attention of researchers in the turbomachinery industry. However, limited studies are currently available on the application of a GF to an axial fan.Hence, in this paper, we conduct a wind tunnel and computational fluid dynamic (CFD) investigation on an axial fan profiled with a NACA 65-(12)10 airfoil to evaluate the effect of the GF on the performance of the fan. We also present the detailed flow features of the fan with and without the GF after validating the simulation results with the experimental results. The experimental results show that as the GF is installed higher on the fan blade, it can produce a higher total pressure rise accompanied with a greater loss of efficiency. The installation of the GF also enlarges the work capacity of the fan. Detailed flow field analysis, including the surface pressure distribution, vorticity distribution at the trailing edge and streamline distribution of the fan, is carried out to understand the mechanisms of the effect of the GF on the performance of the fan.

In this study, we design a new refractive index sensor based on a metal–insulator–metal waveguide coupled with a notched ring and vertical rectangular resonators. We use the finite element method to study the propagation characteristics of the sensor. According to the calculation results, the transmission spectrum exhibits a typical Fano resonance shape. The phenomenon of Fano resonance is caused by the coupling between the wide spectrum and narrow spectrum. In the design, the wide spectrum signal is generated by the vertical rectangular resonator, while the narrow spectrum signal is generated by the notched ring resonator. In addition, we varied the structural parameters of the resonators and filled the structure with media of different refractive indices to study the sensing properties. The maximum achieved sensitivity of the sensor reached 1071.4 nm/RIU. The results reveal potential applications of the coupled system in the field of sensors.

The water cushion depth of stilling basin with shallow-water cushion is a key factor that affects the flow regime of hydraulic jump in the basin. However, the specific depth at which the water cushion is considered as “shallow” has not be stated clearly for now, and only conceptual description is provided. This paper attempts to specify the best water cushion depth based on the flow regime of hydraulic jump and underflow speed; namely, in case of critical hydraulic jump in the basin, the best water cushion depth is located where the minimum distance to the bottom plate of the stilling basin is 1/5~1/4 of the water cushion depth. The theoretical analysis indicates, at different inclinations of discharge chute (θ) and depth ratios of inlet (m), instead of monotonic change, the Froude number (Fr) at inlet of the stilling basin with shallow-water cushion firstly reduces and then increases as the flow velocity at discharge chute inlet (V) increases; the parameters of inflection point (critical flow velocity and critical Fr) increase as the inclinations of discharge chute (θ) and depth ratios of inlet (m) increase. Such regularity is the theoretical basis for selecting representative study cases. The reliability of the large eddy simulation calculation results are verified by a model test; in the paper, 30 cases including five different Froude numbers and six shallow-water cushion depths are selected, for calculating the hydraulic factors such as flow profile, flow regime and flow velocity in the stilling basin with shallow-water cushion; and the varying pattern between the best depth of stilling basin with shallow-water cushion (depth-to-length ratio) and the inflow Froude number is obtained which indicates that the best depth of stilling basin with shallow-water cushion varies little as the change of the Froude number before reaching the critical Froude number; however, the best depth-to-length ratio of stilling basin with shallow-water cushion increases as the Froude number increases after the critical Froude number is reached. The study results in this paper are of reference significance to design and calculation of the stilling basin with shallow-water cushion.

In this work, we focus on investigating electrothermal flow in a rotating electric field (ROT-ETF), with primary attention paid to the horizontal traveling-wave electrothermal (TWET) vortex induced at the center of the electric field. The frequency-dependent flow profiles in the microdevice are analyzed using different heat transfer models. Accordingly, we address in particular the importance of electrode cooling in ROT-ETF as metal electrodes of high thermal conductivity while substrate material of low heat dissipation capability are employed to develop such microfluidic chips. Under this circumstance, cooling of electrode array due to external natural convection on millimeter-scale electrode pads for external wire connection occurs and makes the internal temperature maxima shift from the electrode plane to a bit of distance right above the cross-shaped interelectrode gaps, giving rise to reversal of flow rotation from a typical repulsion-type to attraction-type induction vortex, which is in good accordance with our experimental observations of co-field TWET streaming at frequencies on the order of reciprocal charge relaxation time of the bulk fluid. These results point out a way to make a correct interpretation of out-of-phase electrothermal streaming behavior, which holds great potential for handing high-conductivity analytes in modern microfluidic systems.

The Delaware River has made a marked recovery in the half-century since the adoption of the Delaware River Basin Commission (DRBC) Compact in 1961 and passage of the Federal Clean Water Act amendments during the 1970s. During the 1960s, the DRBC set a 3.5 mg/l dissolved oxygen criteria for the river based on an economic analysis that concluded a waste load abatement program designed to meet fishable water quality goals would generate significant recreation and environmental benefits. Scientists with the Delaware Estuary Program have recently called for raising the 1960s DO criteria along the Delaware River from 3.5 mg/l to 5.0 mg/l to protect anadromous American shad and Atlantic sturgeon and address the prospect of rising temperatures, sea levels, and salinity in the estuary. This research concludes through a marginal abatement cost (MAC) analysis that it would be cost effective to raise DO levels to meet a more stringent standard by prioritizing agricultural conservation and wastewater treatment investments in the Delaware River watershed to reduce 90% of the pollutant load 13.6 million kg/year of nitrogen (30 million lb/year) for $160 million at 35% of the $449 million annual cost. The annual least cost to reduce nitrogen loads and raise dissolved oxygen levels to meet more stringent water quality standards in the Delaware River totals $45 million for atmospheric NOX reduction, $130 million for wastewater treatment, $132 million for agriculture conservation, and $141 million for urban stormwater retrofitting. This 21st century least cost analysis estimates that $50 million/year is needed to reduce pollutant loads in the Delaware River to raise dissolved oxygen levels to 4.0 mg/l, $150 million/year is needed to reach 4.5 mg/l, and $449 million/year is needed to reach 5.0 mg/l.

Starting with a hypothetical geostrophic zonal current in an unbounded ocean, the investigation points out the response of this simple system to a thermal forcing, applied to the free surface and consistent with the maintenance of the geostrophic balance. The main result is the formation of a meridional component of the current, according to the Sverdrup relation, such that the full velocity vector rotates clockwise for heating and anticlockwise for cooling to adjust eventually in the initial zonal direction for large depths.

Nanopores with irregularities are promising tools for distinguishing nano-size objects by their shape, but the forces on the object that critically influence its axial and rotational movement are unclear. The physics of the situation was described using the Poisson-Nernst-Planck and Navier-Stokes equations. With uniformly charged object the axial Coulomb and dielectric pressure (which opposes it and is surprisingly important), control the object's axial movement and rotation. Even without external pressure the hydrodynamic pressure is significant (negative at its upper and positive at its lower surface), but its total value is almost zero. If the object is charged only on the upper surface the axial upper Coulomb pressure is near zero close to the center, but negative near its end (the pressure is zero at the lower surface). The total axial dielectric pressure, which is largely dominated by the pressure at the upper surface, is positive along the length of the object becoming pronounced near its end. The axial hydrodynamic pressure is negative and significant at the upper surface (zero at the lower surface), diminishes in value near the object's end, critically influencing its axial movement, which becomes much faster. At its end the axial dielectric pressure prevails, and controls its rotation.

Thin film transistors (TFTs) using In-Ga-Zn Oxide (IGZO) as active layer and the gate insulator was treated with NH3 plasma and N2O plasma, respectively, which is fabricated on flexible PI substrate in this work. The performance of IGZO TFTs with different plasma species and treatment time are investigated and compared. The experiment results show that the plasma species and treatment time play an important role in the threshold voltage, field-effect mobility, Ion/Ioff ratio, sub-threshold swing (SS) and bias stress stability of the devices. The TFT with a 10 seconds NH3 plasma treatment shows the best performance; specifically, threshold voltage of 0.34 V, field-effect mobility of 15.97 cm2/Vs, Ion/Ioff ratio of 6.33×107, and sub-threshold swing of 0.36 V/dec. The proposed flexible IGZO-TFTs in this paper can be used as driving devices in the next-generation flexible displays.

In this work, we proposed a novel design of microfluidic mixer utilizing the amplified Marangoni chaotic advection induced by AC continuous electrowetting of a metal droplet situated in electrolyte solution, due to the linear and quadratic voltage-dependence of flow velocity at small or large voltages, respectively. Unlike previous researchers exploiting the unidirectional surface stress with DC bias at droplet/medium interface for pumping of electrolyte where the resulting flow rate is linearly proportional to the field intensity, dominance of another kind of dipolar flow pattern caused by local Marangoni stress at the drop surface in sufficiently intense AC electric field is demonstrated by both theoretical analysis and experimental observation, which exhibits a quadratic growth trend as a function of the applied voltage. The dipolar shear stress merely appears at larger voltages and greatly enhances the mixing performance by inducing chaotic advection between the neighboring laminar flow. The mixer design developed herein, on the basis of amplified Marangoni chaotic advection around a liquid metal droplet at larger AC voltages, has great potential for chemical reaction and MEMS actuator applications, because of generating high-throughput and excellent mixing performance at the same time.

Today, Colebrook’s equation is mostly accepted as an informal standard for modeling of turbulent flow in hydraulically smooth and rough pipes including transient zone in between. The empirical Colebrook’s equation relates the unknown flow friction factor (λ) with the known Reynolds number (R) and the known relative roughness of inner pipe surface (ε/D). It is implicit in unknown friction factor (λ). Implicit Colebrook’s equation cannot be rearranged to derive friction factor (λ) directly and therefore it can be solved only iteratively [λ=f(λ, R, ε/D)] or using its explicit approximations [λ≈f(R, ε/D)]. Of course, approximations carry in certain error compared with the iterative solution where the highest level of accuracy can be reached after enough number of iterations. The explicit approximations give a relatively good prediction of the friction factor (λ) and can reproduce accurately Colebrook’s equation and its Moody’s plot. Usually, more complex models of approximations are more accurate and vice versa. In this paper, numerical values of parameters in various existing approximations are changed (optimized) using genetic algorithms to reduce maximal relative error. After this improvement computational burden stays unchanged while accuracy of approximations increases in some of the cases very significantly.

In the past few years, micromodels have become a useful tool for visualizing flow phenomena in porous media with pore structures, e.g., the multifluid dynamics in soils or rocks with fractures in natural geomaterials. Micromodels fabricated using glass or silicon substrates incur high material cost; in particular, the microfabrication-facility cost for making a glass or silicon-based micromold is usually high. This may be an obstacle for researchers investigating the two-phase-flow behavior of porous media. A rigid thermoplastic material is a preferable polymer material for microfluidic models because of its high resistance to infiltration and deformation. In this study, cyclic olefin copolymer (COC) was selected as the substrate for the micromodel because of its excellent chemical, optical, and mechanical properties. A delicate micromodel with a complex pore geometry that represents a two-dimensional (2D) cross-section profile of a fractured rock in a natural oil or groundwater reservoir was developed for two-phase-flow experiments. Using an optical visualization system, we visualized the flow behavior in the micromodel during the processes of imbibition and drainage. The results show that the flow resistance in the main channel (fracture) with a large radius was higher than that in the surrounding area with small pore channels when the injection or extraction rates were low. When we increased the flow rates, the extraction efficiency of the water and oil in the mainstream channel (fracture) did not increase monotonically because of the complex two-phase-flow dynamics. These findings provide a new mechanism of residual trapping in porous media.

In this paper, a two-phase model of air shock wave induced by rock-fall in closed goaf was proposed. The model was made up of the uniform motion phase (velocity was close to 0 m•s-1) and the acceleration movement phase. The uniform motion phase was determined by experience and the acceleration movement phase was derived by the theoretical analysis. After this, a series of experiments were performed to verify the two-phase model and obtained the law of the uniform motion phase. By observing, the acceleration movement phase was taking a larger portion and the uniform motion phase was smaller when height of rock-fall was higher. By comparison, experimental results of different falling heights showed good agreements with theoretical analysis, which verifies the effectiveness of the two-phase model. Finally, the model was tested with computational fluid dynamics (CFD) numerical simulation by three groups of different falling height. The two-phase model could provide a reference and basis for estimating the air shock waves' velocity and design the protective measures.

NACA 0015 and NACA 4415 aerofoil are most common four digits and broadly used aerodynamic shape. Both of the shapes are extensively used for various kind of applications including turbine blade, aircraft wing and so on. NACA 0015 is symmetrical and NACA 4415 is unsymmetrical in shape. Consequently, they have big one-of-a-kind in aerodynamic traits at the side of widespread differences of their utility and performance. Both of them undergo the same fluid principle while applied in any fluid medium giving dissimilar outcomes in aerodynamics behavior. On this work, experimental and numerical investigation of each NACA 0015 and NACA 4415 is done to decide their performance. For this purpose, aerofoil section is tested for a prevalence range attack of angle (AOA). The study addresses the performance of NACA 0015 and NACA 4415 and evaluates the dynamics of flow separation, lift, drag, pressure and velocity contour and so on. This additionally enables to layout new optimistic aerofoil, which is critical to enhance the efficiency and performance of an aircraft in terms of lift enhancement and drag reduction.