A Continuum - Molecular Dynamics Hybrid Method for Micro- and Nano-Fluid Flow

Author(s):  
Ali Kharazmi ◽  
Reza Kamali

A computer program based on a Molecular Dynamics-Continuum hybrid numerical method has been developed in which the Navier-Stokes equations are solved in the continuum region and the atomistic molecular dynamics in molecular region. The prepared algorithm and the computer code are capable of computing flows in micro and nano-scale geometries. The coupling between the continuum equations and the molecular dynamics is constructed through constrained dynamics within an overlap region where both molecular and continuum equations are solved simultaneously. An Overlap region is introduced in two directions to improve the choice of using molecular region in smaller areas. The proposed method is used to simulate steady and start-up Couette flow showing quantitative agreement with results from analytical solutions and full molecular dynamics simulations.

1999 ◽  
Vol 5 (1) ◽  
pp. 17-33 ◽  
Author(s):  
Y. S. Choi ◽  
S. H. Kang

A computer code predicting the flows through the centrifugal compressor with the radial vaneless diffuser was developed and applied to investigate the detailed flowfields, i.e., secondary flows and jet-wake type flow pattern in design and off-design conditions. Various parameters such as slip factors, aerodynamic blockages, entropy generation and two-zone modeling which are widely used in design and performance prediction, were discussed.A control volume method based on a general curvilinear coordinate system was used to solve the time-averaged Navier–Stokes equations and SIMPLER algorithm was used to solve the pressure linked continuity equation. The standardk-εturbulence model was used to obtain the eddy viscosity. Performance of the code was verified using the measured data for the Eckardt impeller.


Author(s):  
V.V. Vyshinsky ◽  
K.T. Zoan

The paper introduces an engineering method for assessing the aerodynamic effect of disturbed atmosphere on an aircraft. As a source of vortex structures, we can consider vortex wind wakes that arise when the atmospheric wind flows around the landscape, large structures, moving or stationary aircraft-carrying platforms, vortex wakes behind aircraft, etc. In this study, we consider the situation when a light transport aircraft and an aircraft of the MC-21 type get into the vortex wake behind the super-heavy aircraft A-380 when flying along the glide path. A coherent vortex structure behind the A-380 is formed by the grid method within the framework of the boundary value problem for the Reynolds-averaged Navier —Stokes equations. The evolution and stochastics of the far wake are carried out using the author’s computer code written in the MATLAB system, within the framework of discrete vortices with a Rankine core. The assessment of the increment of forces and moments from the effect of the vortex system on the aircraft was carried out using the panel method.


2020 ◽  
Vol 635 ◽  
pp. A184 ◽  
Author(s):  
B. Dias ◽  
J. B. Scoggins ◽  
T. E. Magin

Context. Composition, mass, and trajectory parameters of meteors can be derived by combining observations with the meteor physics equations. The fidelity of these equations, which rely on heuristic coefficients, significantly affects the accuracy of the properties inferred. Aims. Our objective is to present a methodology that can be used to compute the luminosity of meteor entry based on detailed flow simulations in the continuum regime. Methods. The methodology consists in solving the Navier–Stokes equations using state-of-the-art physico-chemical models for hypersonic flows. It includes accurate boundary conditions to simulate the surface evaporation of the molten material and coupled flow-radiation effects. Such detailed simulations allow for the calculation of heat-transfer coefficients and luminous efficiency, which can be incorporated into the meteor physics equations. Finally, we integrate the radiative transfer equation over a line of sight from the ground to the meteor to derive the luminosity magnitude. Results. We use the developed methodology to simulate the Lost City bolide and to derive the luminosity magnitude, obtaining good agreement between numerical results and observations. The computed color index is more prominent than the observations. This is attributed to a lack of refractory elements such as Ca in the modeled flow that might originate from the vaporization of droplets in the trail, a phenomenon currently not included in the model.


Author(s):  
Arash Karimipour ◽  
Masoud Afrand

Forced convection of water–Cu nanofluid in a two-dimensional microchannel is studied numerically. The microchannel wall is divided into three parts. The entry and exit ones are kept insulated while the middle one has more temperature than the inlet fluid. The whole of microchannel is under the influence of a magnetic field with uniform strength of B0. Slip velocity and temperature jump are involved along the microchannel walls for different values of slip coefficient such as B = 0.001, B = 0.01, and B = 0.1 for Re = 10, Re = 50, and Re = 100. Navier–Stokes equations are discretized and numerically solved by a developed computer code in FORTRAN. Results are presented as the velocity, temperature, and Nusselt number profiles. Moreover, the effect of magnetic field on slip velocity and temperature jump is investigated for the first time in the present work. Larger Hartmann number, Reynolds number, and volume fraction correspond to more heat transfer rate; however, the effects of Ha and ϕ are more significant at higher Re.


2018 ◽  
Vol 10 (6) ◽  
pp. 168781401878365 ◽  
Author(s):  
Zhaoyong Mao ◽  
Jingang Bai

The development of underwater vehicles is facing the problem of sustainable energy supply. This study introduces a small water turbine, the Lenz turbine, for energy generation from the ocean currents which will provide energy for the underwater vehicles. Computational fluid dynamics simulations of the effect of geometric parameters, including the blade radius, chord length, and pitch angle, on the performance of the turbine are conducted. The Reynolds-Averaged Navier–Stokes equations are numerically solved with a sliding mesh method. Thirteen sets of tests in total are performed at different values of tip-speed ratios. The tests are divided into three groups to study the effect of the three parameters mentioned above, separately. The obtained power coefficients, coefficient of torque, and the dynamic torque on a blade are then compared in each group of tests. Pressure contours and velocity contours are given to explain the reason how the geometric parameters affect the rotor performance.


1980 ◽  
Vol 58 ◽  
pp. 293-299
Author(s):  
Lawrence D. Cloutman

AbstractThe solar granulation has been simulated by numerical solution of the multidimensional, time-dependent, nonlinear Navier-Stokes equations applied to the solar atmosphere. Granules may be explained as buoyantly rising bubbles created at the level where T = 8000 K, and which have collapsed into vortex rings. The calculation is in quantitative agreement with observations and has a number of implications for solar physics and convection theory.


2001 ◽  
Vol 444 ◽  
pp. 383-407 ◽  
Author(s):  
ERCAN ERTURK ◽  
THOMAS C. CORKE

The leading-edge receptivity to acoustic waves of two-dimensional parabolic bodies was investigated using a spatial solution of the Navier–Stokes equations in vorticity/streamfunction form in parabolic coordinates. The free stream is composed of a uniform flow with a superposed periodic velocity fluctuation of small amplitude. The method follows that of Haddad & Corke (1998) in which the solution for the basic flow and linearized perturbation flow are solved separately. We primarily investigated the effect of frequency and angle of incidence (−180° [les ] α2 [les ] 180°) of the acoustic waves on the leading-edge receptivity. The results at α2 = 0° were found to be in quantitative agreement with those of Haddad & Corke (1998), and substantiated the Strouhal number scaling based on the nose radius. The results with sound waves at angles of incidence agreed qualitatively with the analysis of Hammerton & Kerschen (1996). These included a maximum receptivity at α2 = 90°, and an asymmetric variation in the receptivity with sound incidence angle, with minima at angles which were slightly less than α2 = 0° and α2 = 180°.


Author(s):  
Saira F. Pineda ◽  
Arjan M. Kamp ◽  
D. Legendre ◽  
Armando J. Blanco

Flow constituted by drops appears in a wide range of natural, biological and engineering situations. For example, liquid-liquid two phase flow inside capillaries constitutes a model commonly used to represent fluid flow in a petroleum reservoir. The typical modeling approach considers inertial forces negligible in comparison to viscous forces, allowing the use of Stokes equation to study flow dynamics. Very few numerical simulations have been made considering inertial effects. In this project, the flow of a periodic train of drops in a viscous suspending fluid, due to the influence of a fixed pressure gradient, was studied by numerical simulation considering the full Navier-Stokes equations. A numerical approach based on a Volume of Fluid (VOF) formulation was employed using JADIM software, developed by the Institut de Mécanique des Fluides de Toulouse, France. JADIM solves Navier-Stokes equations using a VOF finite volume method, second order in space and time using structured mesh. This two-fluid approach without reconstruction of the interface allows simulating two-phase flows with complex interface shapes. Densities of the drops equal to those of the suspending fluid and a constant interface tension were assumed. The effect of drop size, viscosity ratio, interfacial forces and system pressure gradient on the flow dynamics was studied. Parameters values were chosen to be representative for some particular viscous oil. The result validation shows an excellent agreement between both numerical results. However, there are relative differences between them due to the increase in flow velocity when drop relative size increase and validity of Stokes approach is questionable. Results show non-symmetric eddies in the continuum phase, in a referential frame fixed to the drop. The shape of eddies is strongly influenced by viscosity radio. Drop mobility decreases with increasing size. Additionally, drop mobility also decreases when the viscosity ratio increases. Extra pressure gradient of the system due to the presence of the drop shows a strong dependency on the size ratio between the drop and the pore. For size ratio lower than 0.5, the extra pressure gradient required to move the continuum phase is small. However, when drop to micro-channel ratio exceeds 0.5, the extra pressure gradient significantly increases when the drop size increases. Also, viscosity ratio affects on the system pressure loss, especially in cases where the viscosity ratio is high. The analysis of the capillary number effect on the dynamics of the two-phase system shows that it does not influence drop mobility for the drop sizes considered.


2021 ◽  
Vol 1 (4(68)) ◽  
pp. 40-45
Author(s):  
M. Biliaiev ◽  
V. Biliaieva ◽  
O. Berlov ◽  
V. Kozachyna

The problem of estimating the level of air pollution in the working areas near the coal pile is considered. The task is to develop a CFD model that allows to predict the level of air dust pollution, taking into account the process of wetting the surface of the coal pile. To model the process of coal dust transfer in the air, a twodimensional mass transfer equation is used, which takes into account coal dust transfer due to convection and diffusion. The Navier-Stokes equations are used to calculate the air flow field near the coal pile. Finite-difference schemes of splitting are used for numerical integration of modeling equations. Computer code is developed on the basis of created CFD model. The developed code can be used to analyze the effectiveness of the coal surface wetting to reduce dust pollution of work areas near coal piles. The results of a computational experiment are presented.


Author(s):  
Juan B. V. Wanderley ◽  
Andre´ Ramiro ◽  
Thiago Reis ◽  
Antonio Carlos Fernandes ◽  
Carlos Levi

The highly viscous flow problem of roll damping of a FPSO is investigated by means of numerical solution of the unsteady two-dimensional Navier-Stokes equations. The finite volume method using non-structured grid is used to solve the integral form of the governing equations. The cross section of the FPSO hull with an initial roll displacement is let free to oscillate in roll in an initially still fluid. The numerical simulation provides a realistic picture of the physics of the phenomenon, capturing the vortex formation around the bilge keel. Numerical results from roll free decay tests are compared with experimental data showing a fairly good qualitative and quantitative agreement of the roll damping.


Sign in / Sign up

Export Citation Format

Share Document