Dynamics Method Research Articles

Molecular dynamics study on the thermal properties of DGEBA/DETA/Ag/SWCNT-Ag composite materials.

Traditional conductive adhesives based on epoxy resin system often encounter problems such as high brittleness and low heat resistance. Therefore, it is particularly important to improve the thermal and mechanical properties of the conductive adhesive. In this study, the effects of SWCNT-Ag and SWCNT fillers on the thermal properties of DGEBA/DETA/Ag conductive adhesive system were studied by using molecular dynamics to construct different cross-linking models. The final results show that the addition of SWCNT and SWCNT-Ag can significantly improve the thermal properties of the conductive adhesive. However, the nanosilver particles on the surface of SWCNT-Ag act as a bridge for the connection between SWCNT and Ag in the conductive adhesive. Therefore, SWCNT-Ag has a more positive impact on the thermal properties of DGEBA/DETA/Ag conductive adhesive system. In this paper, the influence of SWCNT-Ag on the thermal properties of traditional DGEBA/DETA/Ag conductive adhesive system was studied by using Materials Studio software. The volume shrinkage, glass transition temperature, thermal expansion coefficient, and thermal conductivity of the material were calculated based on COMPASS force field. The thermal conductivity is calculated by using reverse non-equilibrium molecular dynamics method. Finally, it is found that SWCNT-Ag has a positive effect on the thermal properties of the conductive adhesive system by comparing several groups of calculation data.

Two dynamically discovered compact object candidate binary systems from LAMOST low-resolution survey

ABSTRACT We report two binary systems, LAMOST J035540 + 381550 (hereafter J035540) and LAMOST J035916 + 400732 (hereafter J035916), identified through the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS). Each of these two systems contains an M-type star orbiting with a invisible compact object candidate. Follow-up spectroscopic observations of Palomar 200-in. telescope (P200) enhance radial velocity measurements. We use radial velocities from LAMOST and P200, as well as light curves from Zwicky Transient Facility (ZTF) to constrain orbital parameters. The masses of the visible M-type stars are estimated by fitting the MIST isochrones and spectral energy distributions. The mass functions for the unseen companions are: $0.22\pm 0.01\,{\rm M}_{\odot }$ for J035540 and $0.16\pm 0.01\, {\rm M}_{\odot }$ for J035916. With the orbital and stellar parameters derived above and assuming the orbital inclination is 90° (edge-on), we find that the minimum masses of the invisible companions exceeds that of the visible stars. The single-lined feature and the dynamical evidence suggest the presence of compact objects. J035540’s ZTF light curve, modelled with phoebe, yields a compact object mass of $0.70^{+0.12}_{-0.05}\, {\rm M}_{\odot }$. For J035916, ellipsoidal modulation analysis constrains the light-curve amplitude, yielding a compact object mass range of $0.57-0.90\, {\rm M}_{\odot }$. The mass estimates indicate that both are likely white dwarfs. These findings underscore the efficiency of optical time-domain surveys and dynamical methods in identifying faint, massive white dwarfs, along with other compact objects in binaries.

Exciton Transport in the Nonfullerene Acceptor O-IDTBR from Nonadiabatic Molecular Dynamics.

Theory, computation, and experiment have given strong evidence that charge carriers in organic molecular crystals form partially delocalized quantum objects that diffuse very efficiently via a mechanism termed transient delocalization. It is currently unclear how prevalent this mechanism is for exciton transport. Here we carry out simulation of singlet Frenkel excitons (FE) in a molecular organic semiconductor that belongs to the class of nonfullerene acceptors, O-IDTBR, using the recently introduced FE surface hopping nonadiabatic molecular dynamics method. We find that FE are, on average, localized on a single molecule in the crystal due to sizable reorganization energy and moderate excitonic couplings. Yet, our simulations suggest that the diffusion mechanism is more complex than simple local hopping; in addition to hopping, we observe frequent transient delocalization events where the exciton wave function expands over 10 or more molecules for a short period of time in response to thermal excitations within the excitonic band, followed by de-excitation and contraction onto a single molecule. The transient delocalization events lead to an increase in the diffusion constant by a factor of 3-4, depending on the crystallographic direction as compared to the situation where only local hopping events are considered. Intriguingly, O-IDTBR appears to be a moderately anisotropic 3D "conductor" for excitons but a highly anisotropic 2D conductor for electrons. Taken together with previous simulation results, two trends seem to emerge for molecular organic crystals: excitons tend to be more localized and slower than charge carriers due to higher internal reorganization energy, while exciton transport tends to be more isotropic than charge transport due to the weaker distance dependence of excitonic versus electronic coupling.

Effects of the temperature, strain rate, and loading conditions on the deformation behaviors and mechanical properties of the Ni/Ni3Al superalloy.

Using the molecular dynamics method, we comprehensively studied the effects of temperature, strain rate, and loading conditions on the deformation behaviors and the mechanical properties of the Ni/Ni3Al superalloy. Our investigation revealed that, an increase of the deformation temperature led to a significant improvement of plastic deformation capacity of the system, but the tensile strength and elastic modulus decreased. And the tensile strength and plastic deformation capacity of the system drastically increased with the strain rate. At high deformation temperature and strain rate, the loading conditions had a large effect on the deformation behaviors and the mechanical properties of the system. The difference of the mechanical properties of the system was mainly due to the different deformation mechanism of the system under different deformation temperature, strain rate and loading conditions. Our study offered a theoretical framework for explaining the difference of the mechanical properties for the Ni/Ni3Al superalloy at different service conditions.

Investigation on the Blade Number of Pico-scale Crossflow Turbine for Low Head by Numerical Method

Pico-scale crossflow turbines (CFT) can be an alternative solution to meet electrical energy needs, especially in remote rural areas. CFT is recommended because of its suitability in low head (< 5 m) conditions and fluctuating discharge conditions. One of the parameters that influences the performance of a CFT is the number of blades of the runner. CFT was discovered in 1903 and is still developing; however, the study of the physical phenomena of flow due to the blade number on the energy conversion process has yet to be comprehensively depicted. Therefore, this study aims to analyze the effect of the blade's number of runners on CFT performance using the computational fluid dynamics (CFD) method. The CFD method can visualize the flow field more detail than analytical and experimental. The CFD method is run with a moving mesh feature (transient) and pressure-based solver with a head condition of 3 m. The blades number studied were 16, 18, 22, 24, 26, and 30. Based on the results, the relationship of the CFT efficiency to blade number is described using a second-order multiple regression polynomial, and runner rotation is parabolic. Based on the performance curve, the CFT with 26 blades has the highest performance for low-head conditions.

Investigation on the Blade Number of Pico-scale Crossflow Turbine for Low Head by Numerical Method

Pico-scale crossflow turbines (CFT) can be an alternative solution to meet electrical energy needs, especially in remote rural areas. CFT is recommended because of its suitability in low head (< 5 m) conditions and fluctuating discharge conditions. One of the parameters that influences the performance of a CFT is the number of blades of the runner. CFT was discovered in 1903 and is still developing; however, the study of the physical phenomena of flow due to the blade number on the energy conversion process has yet to be comprehensively depicted. Therefore, this study aims to analyze the effect of the blade's number of runners on CFT performance using the computational fluid dynamics (CFD) method. The CFD method can visualize the flow field more detail than analytical and experimental. The CFD method is run with a moving mesh feature (transient) and pressure-based solver with a head condition of 3 m. The blades number studied were 16, 18, 22, 24, 26, and 30. Based on the results, the relationship of the CFT efficiency to blade number is described using a second-order multiple regression polynomial, and runner rotation is parabolic. Based on the performance curve, the CFT with 26 blades has the highest performance for low-head conditions.

Real-time imaging of intracellular deformation dynamics in vibrated adherent cell cultures.

Mechanical vibration has been shown to regulate cell proliferation and differentiation in vitro and in vivo. However, the mechanism of its cellular mechanotransduction remains unclear. Although the measurement of intracellular deformation dynamics under mechanical vibration could reveal more detailed mechanisms, corroborating experimental evidence is lacking due to technical difficulties. In this study, we aimed to propose a real-time imaging method of intracellular structure deformation dynamics in vibrated adherent cell cultures and investigate whether organelles such as actin filaments connected to a nucleus and the nucleus itself show deformation under horizontal mechanical vibration. The proposed real-time imaging was achieved by conducting vibration isolation and making design improvements to the experimental setup; using a high-speed and high-sensitivity camera with a global shutter; and reducing image blur using a stroboscope technique. Using our system, we successfully produced the first experimental report on the existence of the deformation of organelles connected to a nucleus and the nucleus itself under horizontal mechanical vibration. Furthermore, the intracellular deformation difference between HeLa and MC3T3-E1 cells measured under horizontal mechanical vibration agrees with the prediction of their intracellular structure based on the mechanical vibration theory. These results provide new findings about the cellular mechanotransduction mechanism under mechanical vibration.

Performance analysis of coaxial shear static mixer for hydrogen blending into natural gas

Blending a specific proportion of hydrogen into a natural gas (NG) pipeline network is an effective strategy for achieving long-distance, large-scale, and cost-effective hydrogen transport. Static mixers played a crucial role in this process. In this study, a coaxial shear static mixer with ring-shaped structures instead of turbulence elements is proposed. The mixer exhibited favorable mixing performance and low pressure loss. The flow characteristics of NG and hydrogen within the mixer were investigated with computational fluid dynamics (CFD) methods. The effects of the number of cavities, hydrogen pipe diameter, hydrogen blending ratio and pressure on the mixer performances were investigated, taking into account both mixing uniformity and pressure loss. When the number of cavities in the coaxial shear mixer increases to 48, the mixing uniformity reaches 95.43% (SMX is 95%), and the pressure loss is 67.02 Pa (SMX is 500 Pa).

Comparative study of nudged elastic band and molecular dynamics methods for diffusion kinetics in solid-state electrolytes

Abstract Considerable efforts are being made to transition current lithium-ion and sodium-ion batteries towards the use of solid-state electrolytes. Computational methods, specifically nudged elastic band (NEB) and molecular dynamics (MD) methods, provide powerful tools for the design of solid-state electrolytes. The MD method is usually the choice for studying the materials involving complex multiple diffusion paths or having disordered structures. However, it relies on simulations at temperatures much higher than working temperature. This paper studies the reliability of the MD method using the system of Na diffusion in MgO as a benchmark. We carefully study the convergence behavior of the MD method and demonstrate that total effective simulation time of 12 ns can converge the calculated diffusion barrier to about 0.01 eV. The calculated diffusion barrier is 0.31 eV from both methods. The diffusion coefficients at room temperature are 4.3×10-9 and 2.2×10-9 cm2/s, respectively, from the NEB and MD methods. Our results justify the reliability of the MD method, even though high temperature simulations have to be employed to overcome the limitation on simulation time.

Complex Fluid Models of Mixed Quantum–Classical Dynamics

Several methods in nonadiabatic molecular dynamics are based on Madelung’s hydrodynamic description of nuclear motion, while the electronic component is treated as a finite-dimensional quantum system. In this context, the quantum potential leads to severe computational challenges and one often seeks to neglect its contribution, thereby approximating nuclear motion as classical. The resulting model couples classical hydrodynamics for the nuclei to the quantum motion of the electronic component, leading to the structure of a complex fluid system. This type of mixed quantum–classical fluid models has also appeared in solvation dynamics to describe the coupling between liquid solvents and the quantum solute molecule. While these approaches represent a promising direction, their mathematical structure requires a certain care. In some cases, challenging higher-order gradients make these equations hardly tractable. In other cases, these models are based on phase-space formulations that suffer from well-known consistency issues. Here, we present a new complex fluid system that resolves these difficulties. Unlike common approaches, the current system is obtained by applying the fluid closure at the level of the action principle of the original phase-space model. As a result, the system inherits a Hamiltonian structure and retains energy/momentum balance. After discussing some of its structural properties and dynamical invariants, we illustrate the model in the case of pure-dephasing dynamics. We conclude by presenting some invariant planar models.

Synthesis, Docking, and Machine Learning Studies of Some Novel Quinolinesulfonamides-Triazole Hybrids with Anticancer Activity.

In the presented work, a series of 22 hybrids of 8-quinolinesulfonamide and 1,4-disubstituted triazole with antiproliferative activity were designed and synthesised. The title compounds were designed using molecular modelling techniques. For this purpose, machine-learning, molecular docking, and molecular dynamics methods were used. Calculations of the pharmacokinetic parameters (connected with absorption, distribution, metabolism, excretion, and toxicity) of the hybrids were also performed. The new compounds were synthesised via a copper-catalysed azide-alkyne cycloaddition reaction (CuAAC). 8-N-Methyl-N-{[1-(7-chloroquinolin-4-yl)-1H-1,2,3-triazol-4-yl]methyl}quinolinesulfonamide was identified in in silico studies as a potential strong inhibitor of Rho-associated protein kinase and as a compound that has an appropriate pharmacokinetic profile. The results obtained from in vitro experiments confirm the cytotoxicity of derivative 9b in four selected cancer cell lines and the lack of cytotoxicity of this derivative towards normal cells. The results obtained from silico and in vitro experiments indicate that the introduction of another quinolinyl fragment into the inhibitor molecule may have a significant impact on increasing the level of cytotoxicity toward cancer cells and indicate a further direction for future research in order to find new substances suitable for clinical applications in cancer treatment.

Simulation of ordering and diffusion in CrxMoNbTaVW alloys within the framework of the N-body approach when specifying interatomic interactions

Purpose. Development of interatomic potentials for atomistic simulation of alloys based on refractory metals of the V-Cr-Nb-Mo-Ta-W system and atomistic simulation using these ordering and diffusion potentials in CrxMoNbTaVW alloys. Methods. The development of the interatomic potentials of the V-Cr-Nb-Mo-Ta-W system was carried out within the framework of the 𝑁-body approach; To optimize the parameters of the potentials, the results of calculations within the framework of the electron density functional theory using the VASP software package were used as target values; the simulations of ordering and diffusion was carried out using methods of molecular dynamics and the developed We have previously used the combined method of molecular dynamics and the Monte Carlo method (MD+MK). Results. Within the framework of the 𝑁-body approach, potentials are constructed that complement the V-Nb-Mo-W potential system that we built earlier to the V-Cr-Nb-Mo-Ta-W system. The constructed potentials predict the characteristics of alloys in good agreement with experimental data, CALPHAD data and density functional theory data. Using these potentials by the MD+MC method, CrxMoNbTaVW alloys were modeled, where 𝑥 = 0, 0.5, 1, 2 and 3, in the temperature range from 500 ºC to 2300 ºC. The MD+MC calculations and CALPHAD data agreed in the temperature and concentration regions containing one phase of BCC. At a temperature of 1000ºC, MD + MC calculations show the presence of a single BCC phase, which contradicts CALPHAD data, but is consistent with experimental data. The simulation showed that the atomic structure of the CrMoNbTaVW model alloy is self-sufficient to realize the diffusion of components without artificial introduction of vacancies. The absolute values of the diffusion coefficients of the components and the values of the effective activation energy of diffusion in a solid solution of CrMoNbTaVW were calculated for the first time by the method of molecular dynamics. Analysis of the calculated diffusion characteristics indicates that a diffusion mechanism is implemented in the CrMoNbTaVW alloy, which includes coordinated movement of atoms of various grades. Conclusion. The development of methods for atomistic simulations of alloys based on refractory metals of the V-Cr-Nb-Mo-Ta-W system and the results of atomistic simulation of ordering and diffusion in CrxMoNbTaVW alloys have shown their importance for explaining and predicting heat resistance caused by diffusion processes in multicomponent alloys at high temperatures.

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Simulation and analysis of movement trajectories of fresh tea leaf particles based on CFD‐EDEM coupling

AbstractTo study the air‐suction separation mechanism of fresh tea leaves, this article adopts the coupling simulation method of discrete elements and computational fluid dynamics to conduct numerical simulation and experimental research on the movement of fresh tea leaves adsorbed by negative pressure air flow. In the gas–solid coupling model, three‐dimensional scanning technology is applied to establish a three‐dimensional model of fresh tea leaves with one bud and two leaves. Through the Fluent module in CFD, the negative pressure adsorption gas phase of the simulation process was established, bench tests were established, and parameters of the coupling model were optimized. Finally, the minimum error between the simulated value and the test value was 2.85% and the maximum error was 7.43%, indicating that the coupling model established had certain accuracy and could be theoretically analyzed according to the coupling model established. It is found that when the static pressure values of the adsorption surface are −120 Pa, −340 Pa, and −500 Pa respectively, the initial distance from the adsorption surface is 250–550 mm. It is found that fresh tea leaves with one bud and two leaves all fall to the pipe before adsorption, and then move to the adsorption surface. Under the same static pressure value of the adsorption surface, the initial falling distance of tea particles from the adsorption surface (250–550 mm) has a linear function relationship with the distance between the turning point and the adsorption surface. When the initial distance between falling tea particles and the adsorption surface is the same, the distance between the turning point and the adsorption surface is an exponential function when the initial position is 250–350 mm from the adsorption surface, and the distance between the turning point and the adsorption surface is a linear function when the initial position is 400–550 mm from the adsorption surface.

Time-dependent dynamical energy analysis via convolution quadrature

Dynamical Energy Analysis was introduced in 2009 as a novel method for predicting high-frequency acoustic and vibrational energy distributions in complex engineering structures. In this paper we introduce the first time-dependent Dynamical Energy Analysis method. Time-domain models are important in numerous applications including sound simulation in room acoustics, predicting shock-responses in structural mechanics and modelling electromagnetic scattering from conductors. The first step is to reformulate Dynamical Energy Analysis in the time-domain by means of a convolution integral operator. We are then able to employ the Convolution Quadrature method to provide a link between the previous frequency-domain implementations of Dynamical Energy Analysis and fully time-dependent solutions by means of the Z-transform. By combining a modified multistep Convolution Quadrature approach for the time discretisation, together with Galerkin and Petrov-Galerkin methods for the space and momentum discretisations, respectively, we are able to accurately track the propagation of high-frequency transient signals through phase-space. The implementation here is detailed for finite two-dimensional spatial domains and we demonstrate the versatility of our approach by performing a range of numerical experiments for regular, non-convex and irregular geometries as well as different types of wave source.

Research on Multi-Objective Optimization of Low Pulsation Unloading Damping Groove of Axial Piston Pump

The high- and low-pressure switching of the axial piston pump is realized by the structure of the valve plate, and the buffer-groove structure on the valve plate is very important to reduce the pressure shock and flow fluctuation. In order to optimize the structural parameters of the buffer tank of the piston pump, the influence of the triangular-groove structure on the outlet pressure–flow characteristics was analyzed, and the low-pulsation buffer-groove structure was designed. First, the relationship between the structure of the triangular buffer tank and the output pressure and flow characteristics of the pump was analyzed theoretically. The Computational Fluid Dynamics (CFD) method was used to calculate the pressure and flow characteristics of the whole pump flow field of the triangular-groove structure, and the influence of the buffer-groove structure parameters on the outlet flow pulsation characteristics was studied. The multi-objective optimization algorithm was used to optimize the structure parameters of the buffer tank, and the optimized structure significantly reduced the outlet pressure–flow pulsation of the piston pump. The results show that, after optimization, the width angle of the front and rear triangular groove is 82.3°, the depth angle is 12.7°, the optimized pressure pulsation rate is 0.3%, and the optimized flow pulsation rate is 13.7%.

Study on Grain Boundary Mechanical Behaviors of Polycrystalline γ-TiAl Using Molecular Dynamics Simulations

In this study, the molecular dynamics (MD) method was used to study the tensile deformation of polycrystalline γ-TiAl with complex and random grain orientations. Firstly, the tensile deformation was simulated with different average grain sizes (8.60 nm, 6.18 nm, and 4.50 nm) and strain rates (1 × 108 s−1, 5 × 108 s−1, and 1 × 109 s−1). The results show that the peak stress increases with an increase in tensile strain rate, and the peak stress decreases as the grain size decreases, showing an inverse Hall–Petch effect. Upon observing atomic configuration evolution during tensile deformation, it is found that the grain boundary is seriously distorted, which indicates obvious grain boundary sliding occurring. With a further increase in the loading, some dislocations nucleate at the grain boundaries and propagate towards the interior of the grains along the grain boundaries, which demonstrates that dislocation motion is the primary coordination of the mechanical process of the grain boundaries. The dislocation density near the grain boundaries continues to increase, leading to the generation of micro-cracks and eventually causing material failure. Another interesting phenomenon is that the grains rotate, and the specific rotation angle values of each grain are quantitatively calculated. Grain rotation relaxes the stress concentration near the grain boundaries and plays a toughening role. Consequently, the plastic deformation behaviors of polycrystalline γ-TiAl are achieved through the grain boundary mechanical process, that is, grain boundary sliding and grain rotation.

Impact of sediment uniformity from steep tributary on transport in main Channels: A CFD-DEM study

During flash flood events, a prominent size segregation phenomenon often occurs in the mixture of sediment particles of varying diameters moving with the floodwaters. However, the impact of the inherent inhomogeneity of the sediment on its transport process continues to be unclear. This study evaluates the sorting phenomenon in the transport behavior of flood-impacted, loosely deposited sediment particles originating from steep tributary feeders in mountain channels. To achieve this objective, a computational framework integrating Computational Fluid Dynamics and Discrete Element Method (CFD-DEM) was utilized to directly simulate the detailed paths of all particles. After validating the CFD-DEM model, we investigated the influence of the sediment particles’ coefficient of uniformity (Cu) on sediment transport characteristics. The results suggest that under different Cu values, coarse particles initially accelerate rapidly but later slow down and converge to a consistent speed, whereas fine particles demonstrate a more uniform speed increase and eventually maintain a constant velocity. With rising Cu, the peak degree of sorting between coarse and fine particles first increases and then stabilizes. The temporal variation of the dimensionless sediment transport rate of any sediment component at a given section can be represented using a modified dual-branch Weibull function. The overall transport process can be determined through linear superposition based on mass proportion. Additionally, we found that the heterogeneous characteristics of the sediment provide a mode of overall sediment transport where fine particles roll at the bottom, while coarse particles slide and saltate at the top, leading to reduced frictional dissipation in the movement of coarse particles. This phenomenon is commonly observed in mountain flood channels with slopes less than 60 degrees. These findings bolster the explanation for the poor application of existing sediment transport rate formulas in flash flood channels.

Molecular dynamics study on tribological properties of AlCrFeCoNi HEA at different temperatures

The microscopically wear mechanism of AlCrFeCoNi HEA effected by high temperatures still remains unclear. In this work, the tribological properties of AlCrFeCoNi HEA and tungsten carbide (WC) at temperatures ranging from 250, 400, 550, 700, 850 °C was investigated using molecular dynamics methods. Additionally, the shear strain, contact stress, and energy change at different temperatures were analyzed to explore the atomic wear mechanism. The simulated results indicated that the coefficient of friction and wear rate decrease with the increasing temperature, since the dislocation types decrease and the kinetic energy and total system energy increase. And high temperature can inhibit the synergistic effect of strain accumulation and atomic wear on the surface of AlCrFeCoNi HEAs.

A molecular dynamics study on pore structure: Performance comparison between metal foam and artificial mesh porous surface

With the development of surface engineering, porous surfaces have emerged as a significant research subject in boiling heat transfer. The latter, in turn, plays a crucial role in various industries such as power plants, distillation plants, and microelectronic technology. In this paper, the Molecular Dynamics method is adopted to investigate the wicking dynamics and boiling dynamics of two porous surfaces: foam, which exhibits randomly distributed pores, and mesh, composed of ordered square wires with relatively uniform pore sizes. Three wettability, namely hydrophilic, neutral, and hydrophobic wetting states, are assigned to the two porous surfaces de-coupling the effect of wettability from surface structure. Results reveal that, during the wicking process, the foam surface shows better wetting ability as it absorbs liquid under both hydrophilic and neutral wettability. Comparatively, the mesh surface has the fastest wicking speed under hydrophilic wettability yet it becomes non-wetting under neutral wettability. During the boiling process, the boiling dynamics differ greatly under three wettability. More importantly, the difference in surface structure makes the foam surface possess a better heat transfer whereas the mesh surface causes gentle pressure variation. Our findings provide insights into the design of artificial porous surfaces for certain purpose and their potential application.