Based on the jet stability theory, under the conditions of gas rotation, fluid compressibility and supercavitation, this paper gives the mathematical model describing the thermal instability of supercavitating liquid jet surrounded by a coaxial rotary gas, and the corresponding numerical method for solving the mathematical model is proposed and verified by the data in reference. Then, this paper analyzes the effects of gas–liquid temperature differences and temperature gradients on jet instability, and studies the thermal stability of supercavitating jet. The results show that the maximum disturbance growth rate, the dominant frequency and the maximum disturbance wave numbers increase linearly with the increase of gas–liquid temperature differences. The existence of temperature gradient inside the jet makes the effects of temperature differences on jet instability more obvious. The temperature gradient will inhibit the effect of supercavitation on jet instability, while gas–liquid temperature difference will promote the effect of supercavitation on jet instability.
{"title":"On the thermal instability of supercavitating liquid jet surrounded by coaxial rotary gas","authors":"M. Lü, Z. Ning","doi":"10.1093/jom/ufab024","DOIUrl":"https://doi.org/10.1093/jom/ufab024","url":null,"abstract":"\u0000 Based on the jet stability theory, under the conditions of gas rotation, fluid compressibility and supercavitation, this paper gives the mathematical model describing the thermal instability of supercavitating liquid jet surrounded by a coaxial rotary gas, and the corresponding numerical method for solving the mathematical model is proposed and verified by the data in reference. Then, this paper analyzes the effects of gas–liquid temperature differences and temperature gradients on jet instability, and studies the thermal stability of supercavitating jet. The results show that the maximum disturbance growth rate, the dominant frequency and the maximum disturbance wave numbers increase linearly with the increase of gas–liquid temperature differences. The existence of temperature gradient inside the jet makes the effects of temperature differences on jet instability more obvious. The temperature gradient will inhibit the effect of supercavitation on jet instability, while gas–liquid temperature difference will promote the effect of supercavitation on jet instability.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48601010","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The vector form intrinsic finite element (VFIFE) method is a solution technique for nonlinear structural problems, which describes a continuous body using a set of particles instead of a mathematical function. Thus, a dynamic particle equation can be established by Newton's law of motion, and a viscous or kinetic damping can be introduced to obtain the steady state of the structure. This paper focuses mainly on the development of a stability condition regarding the explicit central difference method used in VFIFE to guarantee the system's convergence. The process is established and evaluated in combination with a dynamic relaxation method with kinetic damping and discrete control theory. Four numerical examples of structure nonlinear problems are used to verify the accuracy, stability and efficiency of the method.
{"title":"Development of an adaptive explicit algorithm for static simulation using the vector form intrinsic finite element method","authors":"Mien-Li Wang, C. Chuang, J. Lee","doi":"10.1093/jom/ufab022","DOIUrl":"https://doi.org/10.1093/jom/ufab022","url":null,"abstract":"\u0000 The vector form intrinsic finite element (VFIFE) method is a solution technique for nonlinear structural problems, which describes a continuous body using a set of particles instead of a mathematical function. Thus, a dynamic particle equation can be established by Newton's law of motion, and a viscous or kinetic damping can be introduced to obtain the steady state of the structure. This paper focuses mainly on the development of a stability condition regarding the explicit central difference method used in VFIFE to guarantee the system's convergence. The process is established and evaluated in combination with a dynamic relaxation method with kinetic damping and discrete control theory. Four numerical examples of structure nonlinear problems are used to verify the accuracy, stability and efficiency of the method.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47944540","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Haiyin Cao, Yu Huang, Y. Rong, Hao Wu, Minghui Guo
In this study, the influence of inlet pocket size on the static performance of non-Newtonian lubricated hole-entry hybrid journal bearings is theoretically analyzed. The oil film of the bearing is discretized into a nonuniform mesh containing the geometric characteristics of the oil inlet pocket, and the inlet pocket is treated as a micro-oil recess. The Reynolds equation is solved by the finite element method based on Galerkin's techniques, and a new solution strategy to solve the recess/pocket pressure is proposed. The power-law model is used to introduce the non-Newtonian effect. The results show that the static performance characteristics of this type of bearing are greatly affected by the pocket size at both zero speed and high speed.
{"title":"Analysis of the effect of oil inlet size on the static performance characteristics of non-Newtonian lubricated hole-entry hybrid journal bearings","authors":"Haiyin Cao, Yu Huang, Y. Rong, Hao Wu, Minghui Guo","doi":"10.1093/jom/ufab020","DOIUrl":"https://doi.org/10.1093/jom/ufab020","url":null,"abstract":"\u0000 In this study, the influence of inlet pocket size on the static performance of non-Newtonian lubricated hole-entry hybrid journal bearings is theoretically analyzed. The oil film of the bearing is discretized into a nonuniform mesh containing the geometric characteristics of the oil inlet pocket, and the inlet pocket is treated as a micro-oil recess. The Reynolds equation is solved by the finite element method based on Galerkin's techniques, and a new solution strategy to solve the recess/pocket pressure is proposed. The power-law model is used to introduce the non-Newtonian effect. The results show that the static performance characteristics of this type of bearing are greatly affected by the pocket size at both zero speed and high speed.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46552802","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, a method is proposed for extracting fracture parameters in isotropic material cracking via a stable generalized/extended finite element method. The numerical results of the stress intensity factors and scaled condition number of the system matrix are presented and compared with different enrichment schemes or those reported in related references. The good agreement and convergence of the results obtained by the developed method with those obtained by other solutions or enrichment schemes proves the applicability of the proposed approach and confirms its capability of efficiently extracting fracture parameters in isotropic materials.
{"title":"Fracture analysis for materials by a stable generalized/extended finite element method","authors":"H. Jia, Y. Zhao, Y. Nie, S. Q. Li","doi":"10.1093/jom/ufab019","DOIUrl":"https://doi.org/10.1093/jom/ufab019","url":null,"abstract":"\u0000 In this paper, a method is proposed for extracting fracture parameters in isotropic material cracking via a stable generalized/extended finite element method. The numerical results of the stress intensity factors and scaled condition number of the system matrix are presented and compared with different enrichment schemes or those reported in related references. The good agreement and convergence of the results obtained by the developed method with those obtained by other solutions or enrichment schemes proves the applicability of the proposed approach and confirms its capability of efficiently extracting fracture parameters in isotropic materials.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48652301","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This study proposes the optimization of a low-level assembly code to reconstruct the flux for a splitting flux Harten–Lax–van Leer (SHLL) scheme on high-end graphic processing units. The proposed solver is implemented using the weighted essentially non-oscillatory reconstruction method to simulate compressible gas flows that are derived using an unsteady Euler equation. Instructions in the low-level assembly code, i.e. parallel thread execution and instruction set architecture in compute unified device architecture (CUDA), are used to optimize the CUDA kernel for the flux reconstruction method. The flux reconstruction method is a fifth-order one that is used to process the high-resolution intercell flux for achieving a highly localized scheme, such as the high-order implementation of SHLL scheme. Many benchmarking test cases including shock-tube and four-shock problems are demonstrated and compared. The results show that the reconstruction method is computationally very intensive and can achieve excellent performance up to 5183 GFLOP/s, ∼66% of peak performance of NVIDIA V100, using the low-level CUDA assembly code. The computational efficiency is twice the value as compared with the previous studies. The CUDA assembly code reduces 26.7% calculation and increases 37.5% bandwidth. The results show that the optimal kernel reaches up to 990 GB/s for the bandwidth. The overall efficiency of bandwidth and computation performance achieves 127% of the predicted performance based on the HBM2-memory roofline model estimated by Empirical Roofline Tool.
{"title":"Implementation of a parallel high-order WENO-type Euler equation solver using a CUDA PTX paradigm","authors":"F. Kuo, J. S. Wu","doi":"10.1093/jom/ufab016","DOIUrl":"https://doi.org/10.1093/jom/ufab016","url":null,"abstract":"\u0000 This study proposes the optimization of a low-level assembly code to reconstruct the flux for a splitting flux Harten–Lax–van Leer (SHLL) scheme on high-end graphic processing units. The proposed solver is implemented using the weighted essentially non-oscillatory reconstruction method to simulate compressible gas flows that are derived using an unsteady Euler equation. Instructions in the low-level assembly code, i.e. parallel thread execution and instruction set architecture in compute unified device architecture (CUDA), are used to optimize the CUDA kernel for the flux reconstruction method. The flux reconstruction method is a fifth-order one that is used to process the high-resolution intercell flux for achieving a highly localized scheme, such as the high-order implementation of SHLL scheme. Many benchmarking test cases including shock-tube and four-shock problems are demonstrated and compared. The results show that the reconstruction method is computationally very intensive and can achieve excellent performance up to 5183 GFLOP/s, ∼66% of peak performance of NVIDIA V100, using the low-level CUDA assembly code. The computational efficiency is twice the value as compared with the previous studies. The CUDA assembly code reduces 26.7% calculation and increases 37.5% bandwidth. The results show that the optimal kernel reaches up to 990 GB/s for the bandwidth. The overall efficiency of bandwidth and computation performance achieves 127% of the predicted performance based on the HBM2-memory roofline model estimated by Empirical Roofline Tool.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49655167","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Corrigendum to: Nonlinear free vibration of size-dependent microbeams with nonlinear elasticity under various boundary conditions","authors":"F. Lin, J. Peng, S. Xue, L. Yang, J. Yang","doi":"10.1093/jom/ufab018","DOIUrl":"https://doi.org/10.1093/jom/ufab018","url":null,"abstract":"","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41338449","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this research, a surface-residual-based nonlocal stress was introduced into nonlocal damage theory to describe the long-range actions among microstructures that were excluded in the definition of Cauchy stress. By using the surface-residual-based nonlocal stress tensor, a thermodynamically consistent nonlocal integral damage model was established to simulate the strain localization behavior for elastic-brittle damage problems. In this model, both the strain and the damage were taken as nonlocal variables in the free energy function, and the integral-type damage constitutive relationships and the evolution equation were derived via thermodynamic laws in order to ensure the self-consistency within the thermodynamic framework. Based on the nonlocal damage formulations using a real nonlocal stress concept, we simulated the strain localization phenomenon in an elastic bar subjected to uniaxial tension. The results showed clear localizing and softening features of strain in the damage zone, and the boundary effects arising from the nonlocal surface residual were illuminated. Furthermore, the strain localization behaviors for different internal characteristic lengths were simulated, through which we found that the characteristic length was comparable to the size of the strain localization zone.
{"title":"A thermodynamically nonlocal damage model using a surface-residual-based nonlocal stress","authors":"Liyang Huang, Yin Yao","doi":"10.1093/jom/ufab017","DOIUrl":"https://doi.org/10.1093/jom/ufab017","url":null,"abstract":"\u0000 In this research, a surface-residual-based nonlocal stress was introduced into nonlocal damage theory to describe the long-range actions among microstructures that were excluded in the definition of Cauchy stress. By using the surface-residual-based nonlocal stress tensor, a thermodynamically consistent nonlocal integral damage model was established to simulate the strain localization behavior for elastic-brittle damage problems. In this model, both the strain and the damage were taken as nonlocal variables in the free energy function, and the integral-type damage constitutive relationships and the evolution equation were derived via thermodynamic laws in order to ensure the self-consistency within the thermodynamic framework. Based on the nonlocal damage formulations using a real nonlocal stress concept, we simulated the strain localization phenomenon in an elastic bar subjected to uniaxial tension. The results showed clear localizing and softening features of strain in the damage zone, and the boundary effects arising from the nonlocal surface residual were illuminated. Furthermore, the strain localization behaviors for different internal characteristic lengths were simulated, through which we found that the characteristic length was comparable to the size of the strain localization zone.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41414758","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this paper, the available formulae for the curvature of plane curve are reviewed not only for the time-like but also for the space-like parameter curve. Two ways to describe the curve are proposed. One is the straight way to obtain the Frenet formula according to the given curve of parameter form. The other is that we can construct the curve by solving the state equation of Frenet formula subject to the initial position, the initial tangent, normal and binormal vectors, and the given radius of curvature and torsion constant. The remainder theorem of the matrix and the Cayley–Hamilton theorem are both employed to solve the Frenet equation. We review the available formulae of the radius of curvature and examine their equivalence. Through the Frenet formula, the relation among different expressions for the radius of curvature formulae can be linked. Therefore, we can integrate the formulae in the engineering mathematics, calculus, mechanics of materials and dynamics. Besides, biproduct of two new and simpler formulae and the available four formulae in the textbook of the radius of curvature yield the same radius of curvature for the plane curve. Linkage of centrifugal force and radius of curvature is also addressed. A demonstrative example of the cycloid is given. Finally, we use the two new formulae to obtain the radius of curvature for four curves, namely a circle. The equivalence is also proved. Animation for 2D and 3D curves is also provided by using the Mathematica software to demonstrate the validity of the present approach.
{"title":"Construction of a curve by using the state equation of Frenet formula","authors":"J. Chen, J. W. Lee, S. Kao, Y. Chou","doi":"10.1093/jom/ufab014","DOIUrl":"https://doi.org/10.1093/jom/ufab014","url":null,"abstract":"\u0000 In this paper, the available formulae for the curvature of plane curve are reviewed not only for the time-like but also for the space-like parameter curve. Two ways to describe the curve are proposed. One is the straight way to obtain the Frenet formula according to the given curve of parameter form. The other is that we can construct the curve by solving the state equation of Frenet formula subject to the initial position, the initial tangent, normal and binormal vectors, and the given radius of curvature and torsion constant. The remainder theorem of the matrix and the Cayley–Hamilton theorem are both employed to solve the Frenet equation. We review the available formulae of the radius of curvature and examine their equivalence. Through the Frenet formula, the relation among different expressions for the radius of curvature formulae can be linked. Therefore, we can integrate the formulae in the engineering mathematics, calculus, mechanics of materials and dynamics. Besides, biproduct of two new and simpler formulae and the available four formulae in the textbook of the radius of curvature yield the same radius of curvature for the plane curve. Linkage of centrifugal force and radius of curvature is also addressed. A demonstrative example of the cycloid is given. Finally, we use the two new formulae to obtain the radius of curvature for four curves, namely a circle. The equivalence is also proved. Animation for 2D and 3D curves is also provided by using the Mathematica software to demonstrate the validity of the present approach.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-06-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44845196","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this study, the hyper-viscoelastic behavior of adhesive films was characterized. A constitutive model was developed by combining the Mooney–Rivlin hyperelastic model and a viscoelastic model expressed in terms of the Prony series to describe the constitutive behavior of the adhesive films. The material parameters of the developed constitutive model were determined through single-step stress relaxation tests conducted for 30 min at four strain levels: 100%, 200%, 300% and 400%. Based on the reduced gradient method, the optimized material parameters were then evaluated by curve fitting the experimental data. To validate the proposed constitutive model, we performed the tensile tests at different strain rates from 5 × 10−4 to 5 × 10−1 s−1 and the multistep stress relaxation tests on the adhesive films. The model predictions and experimental data were in good agreement. Thus, the proposed hyper-viscoelastic constitutive model with parameters determined through single-step stress relaxation tests is effective in characterizing the mechanical behavior of adhesive films.
{"title":"Characterizing the hyper-viscoelastic behavior of adhesive films","authors":"H. Hsu, J. Tsai","doi":"10.1093/JOM/UFAB013","DOIUrl":"https://doi.org/10.1093/JOM/UFAB013","url":null,"abstract":"\u0000 In this study, the hyper-viscoelastic behavior of adhesive films was characterized. A constitutive model was developed by combining the Mooney–Rivlin hyperelastic model and a viscoelastic model expressed in terms of the Prony series to describe the constitutive behavior of the adhesive films. The material parameters of the developed constitutive model were determined through single-step stress relaxation tests conducted for 30 min at four strain levels: 100%, 200%, 300% and 400%. Based on the reduced gradient method, the optimized material parameters were then evaluated by curve fitting the experimental data. To validate the proposed constitutive model, we performed the tensile tests at different strain rates from 5 × 10−4 to 5 × 10−1 s−1 and the multistep stress relaxation tests on the adhesive films. The model predictions and experimental data were in good agreement. Thus, the proposed hyper-viscoelastic constitutive model with parameters determined through single-step stress relaxation tests is effective in characterizing the mechanical behavior of adhesive films.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-06-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45365906","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Chun-Kuei Chen, Sheng-Qi Chen, Wei‐Mon Yan, Wen‐Ken Li, Ta-Hui Lin
The present study is concerned with the experimental impingement of two consecutive droplets on an inclined solid surface. Attention is mainly paid to the effects of impingement timing with various oblique angles (Φ) of the surface on the impact phenomena, which mainly affect the maximum droplet spreading diameter. The investigation considers four impingement scenarios differentiated by impingement timing, namely Case 1: single-droplet impingement; Case 2 of Δt1: the moment when the leading droplet starts spreading along the oblique surface; Case 3 of Δt2: the moment when the leading droplet reaches its maximum spreading; and Case 4 of Δt3: the moment when the leading droplet starts retracting. It is observed that deformation behavior of two successive droplets impacting on the inclined surface experiences a complex asymmetric morphology evolution due to the enhancement of gravity effect and various conditions of the impingement timing. The merged droplet becomes slender with increasing oblique surface angle in the final steady shape, causing the decrease in the value of front and back contact angles. The impingement timing has a significant influence on the change of the maximum height of the merged droplet. The coalesced droplet spreads to the maximum dimensionless width diameter at Δt = Δt2 and the oblique angle of Φ = 45°, but reaches the maximum dimensionless height for Δt = Δt2 at Φ = 30°. The front contact angles converge to a fixed value eventually for all conditions of impingement timing, and the values become lower with the increasing surface inclination.
{"title":"Experimental study on two consecutive droplets impacting onto an inclined solid surface","authors":"Chun-Kuei Chen, Sheng-Qi Chen, Wei‐Mon Yan, Wen‐Ken Li, Ta-Hui Lin","doi":"10.1093/JOM/UFAB012","DOIUrl":"https://doi.org/10.1093/JOM/UFAB012","url":null,"abstract":"\u0000 The present study is concerned with the experimental impingement of two consecutive droplets on an inclined solid surface. Attention is mainly paid to the effects of impingement timing with various oblique angles (Φ) of the surface on the impact phenomena, which mainly affect the maximum droplet spreading diameter. The investigation considers four impingement scenarios differentiated by impingement timing, namely Case 1: single-droplet impingement; Case 2 of Δt1: the moment when the leading droplet starts spreading along the oblique surface; Case 3 of Δt2: the moment when the leading droplet reaches its maximum spreading; and Case 4 of Δt3: the moment when the leading droplet starts retracting. It is observed that deformation behavior of two successive droplets impacting on the inclined surface experiences a complex asymmetric morphology evolution due to the enhancement of gravity effect and various conditions of the impingement timing. The merged droplet becomes slender with increasing oblique surface angle in the final steady shape, causing the decrease in the value of front and back contact angles. The impingement timing has a significant influence on the change of the maximum height of the merged droplet. The coalesced droplet spreads to the maximum dimensionless width diameter at Δt = Δt2 and the oblique angle of Φ = 45°, but reaches the maximum dimensionless height for Δt = Δt2 at Φ = 30°. The front contact angles converge to a fixed value eventually for all conditions of impingement timing, and the values become lower with the increasing surface inclination.","PeriodicalId":50136,"journal":{"name":"Journal of Mechanics","volume":" ","pages":""},"PeriodicalIF":1.7,"publicationDate":"2021-06-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45766840","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}