Liquid Characteristics Under Melting/Solidification Conditions Using Energy Conserving Dissipative Particle Dynamics

Author(s):  
Erik Johansson ◽  
Toru Yamada ◽  
Jinliang Yuan ◽  
Bengt Sundén ◽  
Yutaka Asako ◽  
...  

In this paper, energy conserving Dissipative Particle Dynamics (DPDe) is used to study liquid characteristics when the walls are kept at a melting temperature. The formulation of the phase change problem is based on the latent heat model available in the literature. It is incorporated into the DPDe method to simulate a one-dimensional solid-liquid moving boundary problem. The solution domain is considered to be a two-dimensional Cartesian box where DPDe particles are randomly distributed. Periodic boundary conditions are applied in the flow direction and solid DPDe particles are placed as additional layers on the top and bottom of the domain. The DPDe result was compared with the available analytical solution and the effects of the DPDe parameters and thermal characteristics are discussed.

Author(s):  
M. S. Zaman ◽  
M. G. Satish

It is crucial to understand how one fluid is displaced by another at different temperature through a capillary, as many industrial and reservoir enhanced recovery methods fall into this category. Dissipative particle dynamics (DPD) method has been successfully applied to model mesoscale behaviors of many processes. In this paper, DPD method with energy conservation has been applied to model non-isothermal fluid displacement in capillary tube. Validation of the in-house computer code written in C# is carried out by modeling isothermal no-slip fluid flow. Simulation of non-isothermal fluid displacement using energy conserving DPD gives insight about the parameters affecting the flow.


Author(s):  
Anuj Chaudhri ◽  
Jennifer R. Lukes

The energy conserving formulation of the Dissipative Particle Dynamics (DPD) mesoscale method for multicomponent systems is analyzed thoroughly. A new framework is established by identifying the dimensionless groups using general scaling factors. When the scaling factors are chosen based on the solvent in a multicomponent system, the reduced system of equations can easily be solved computationally. Simulation results are presented for one dimensional transient and steady-state heat conduction in a random DPD solid, which compare well with existing published and analytical solutions. This model is extended to two dimensions and shows excellent agreement with the analytical solution.


Author(s):  
Toru Yamada ◽  
Anurag Kumar ◽  
Yutaka Asako ◽  
Mohammad Faghri

Dissipative particle dynamics (DPD) with energy conservation was applied to simulate forced convection in parallel-plate channels with boundary conditions of constant wall temperature (CWT) and constant wall heat flux (CHF). DPD is a coarse-grained version of molecular dynamics. An additional governing equation for energy conservation was solved along with conventional DPD equations where inter-particle heat flux accounts for changes in mechanical and internal energies when particles interact with surrounding particles. The solution domain was considered to be two-dimensional with periodic boundary condition in the flow direction. Additional layers of particles on top and bottom of the channel were utilized to apply no-slip velocity and temperature boundary conditions. The governing equations for energy conservation were modified based on periodic fully developed velocity and temperature conditions. The results were shown via velocity and temperature profiles across the channel cross section. The Nusselt numbers were calculated from the temperature gradient at the wall using a second order accurate forward difference approximation. The results agreed well with the exact solutions to within 2.3%.


Author(s):  
Toru Yamada ◽  
Yutaka Asako ◽  
Mohammad Faghri ◽  
Bengt Sundén

The liquid flow in sub-micron channels is simulated using multi-body dissipative particle dynamics (MDPD) to study the effect of the surface tension between liquid and wall surface on the flow in sub-micon scale. The solution domain is considered to be two-dimensional, where DPD particles are randomly distributed. Periodic boundary condition is employed in the flow direction and the solid walls are created by distributing DPD particles in the additional layers on the top and bottom of the domain. The different surface tensions between liquid and wall surface are obtained by changing the interaction parameters between the liquid and wall DPD particles. The ratio of Capillary number (Ca) to Reynolds number (Re) is used to relate the DPD units to the physical units. The results are shown in the form of slip length and the effect of the surface tension on the liquid flow in sub-micron channels is discussed.


Author(s):  
Toru Yamada ◽  
Yutaka Asako ◽  
Mohammad Faghri ◽  
Chungpyo Hong

The effective thermal conductivity of Al2O3/water and CuO/water nanofluids were modeled by numerically solving steady heat flow in one-dimensional channels. This was accomplished by using energy conserving dissipative particle dynamics (DPDe). The effects of the interfacial thermal resistance and the Brownian motion of nanoparticles were incorporated in the model by modifying the conductive interaction parameter in the energy equation. The results were presented in the form of the thermal conductivity of nanofluids as functions of particle volume fraction and temperature, and were compared with the available experimental and analytical results. The present model agreed well with the experimental results for Al2O3/water nanofluid while there were discrepancies between the model and the results for CuO/water nanofluid.


Sign in / Sign up

Export Citation Format

Share Document