Processing math: 100%
Research article Special Issues

Multiscale simulation of a polymer melt flow between two coaxial cylinders under nonisothermal conditions

  • We successfully extend a multiscale simulation (MSS) method to nonisothermal well-entangled polymer melt flows between two coaxial cylinders. In the multiscale simulation, the macroscopic flow system is connected to a number of microscopic systems through the velocity gradient tensor, stress tensor and temperature. At the macroscopic level, in addition to the momentum balance equation, we consider the energy balance equation, where heat generation plays an important role not only in the temperature distribution but also in the flow profile. At the microscopic level, a dual slip-link model is employed for well-entangled polymers. To incorporate the temperature effect into the microscopic systems, we used the time-temperature superposition rule for the slip-link model, in which the temperature dependence of the parameters is not known; on the other hand, the way to take into account the temperature effect in the macroscopic equations has been well established. We find that the extended multiscale simulation method is quite effective in revealing the relation between nonisothermal polymeric flows for both steady and transient cases and the microscopic states of polymer chains expressed by primitive paths and slip-links. It is also found that the temperature-dependent reptation-time-based Weissenberg number is a suitable measure for understanding the extent of the polymer chain deformation in the range of the shear rate used in this study.

    Citation: Yuji Hamada, Takeshi Sato, Takashi Taniguchi. Multiscale simulation of a polymer melt flow between two coaxial cylinders under nonisothermal conditions[J]. Mathematics in Engineering, 2021, 3(6): 1-22. doi: 10.3934/mine.2021042

    Related Papers:

    [1] Abdulwahab Ali Almazroi . Survival prediction among heart patients using machine learning techniques. Mathematical Biosciences and Engineering, 2022, 19(1): 134-145. doi: 10.3934/mbe.2022007
    [2] Weidong Gao, Yibin Xu, Shengshu Li, Yujun Fu, Dongyang Zheng, Yingjia She . Obstructive sleep apnea syndrome detection based on ballistocardiogram via machine learning approach. Mathematical Biosciences and Engineering, 2019, 16(5): 5672-5686. doi: 10.3934/mbe.2019282
    [3] Lili Jiang, Sirong Chen, Yuanhui Wu, Da Zhou, Lihua Duan . Prediction of coronary heart disease in gout patients using machine learning models. Mathematical Biosciences and Engineering, 2023, 20(3): 4574-4591. doi: 10.3934/mbe.2023212
    [4] Jingchi Jiang, Xuehui Yu, Yi Lin, Yi Guan . PercolationDF: A percolation-based medical diagnosis framework. Mathematical Biosciences and Engineering, 2022, 19(6): 5832-5849. doi: 10.3934/mbe.2022273
    [5] Abigail Ferreira, Rui Lapa, Nuno Vale . A retrospective study comparing creatinine clearance estimation using different equations on a population-based cohort. Mathematical Biosciences and Engineering, 2021, 18(5): 5680-5691. doi: 10.3934/mbe.2021287
    [6] Atsushi Kawaguchi . Network-based diagnostic probability estimation from resting-state functional magnetic resonance imaging. Mathematical Biosciences and Engineering, 2023, 20(10): 17702-17725. doi: 10.3934/mbe.2023787
    [7] Michele L. Joyner, Cammey Cole Manning, Whitney Forbes, Valerie Bobola, William Frazier . Modeling Ertapenem: the impact of body mass index on distribution of the antibiotic in the body. Mathematical Biosciences and Engineering, 2019, 16(2): 713-726. doi: 10.3934/mbe.2019034
    [8] Wenzhu Song, Lixia Qiu, Jianbo Qing, Wenqiang Zhi, Zhijian Zha, Xueli Hu, Zhiqi Qin, Hao Gong, Yafeng Li . Using Bayesian network model with MMHC algorithm to detect risk factors for stroke. Mathematical Biosciences and Engineering, 2022, 19(12): 13660-13674. doi: 10.3934/mbe.2022637
    [9] Anastasia-Maria Leventi-Peetz, Kai Weber . Probabilistic machine learning for breast cancer classification. Mathematical Biosciences and Engineering, 2023, 20(1): 624-655. doi: 10.3934/mbe.2023029
    [10] Tingxi Wen, Hanxiao Wu, Yu Du, Chuanbo Huang . Faster R-CNN with improved anchor box for cell recognition. Mathematical Biosciences and Engineering, 2020, 17(6): 7772-7786. doi: 10.3934/mbe.2020395
  • We successfully extend a multiscale simulation (MSS) method to nonisothermal well-entangled polymer melt flows between two coaxial cylinders. In the multiscale simulation, the macroscopic flow system is connected to a number of microscopic systems through the velocity gradient tensor, stress tensor and temperature. At the macroscopic level, in addition to the momentum balance equation, we consider the energy balance equation, where heat generation plays an important role not only in the temperature distribution but also in the flow profile. At the microscopic level, a dual slip-link model is employed for well-entangled polymers. To incorporate the temperature effect into the microscopic systems, we used the time-temperature superposition rule for the slip-link model, in which the temperature dependence of the parameters is not known; on the other hand, the way to take into account the temperature effect in the macroscopic equations has been well established. We find that the extended multiscale simulation method is quite effective in revealing the relation between nonisothermal polymeric flows for both steady and transient cases and the microscopic states of polymer chains expressed by primitive paths and slip-links. It is also found that the temperature-dependent reptation-time-based Weissenberg number is a suitable measure for understanding the extent of the polymer chain deformation in the range of the shear rate used in this study.


    Obstructive sleep apnea (OSA) is featured by repeated upper airway narrowing or obstruction, causing a substantial reduction or cessation of airflow during sleep [1,2]. These disruptions in breathing lead to intermittent hypercapnia and hypoxemia, causing metabolic dysregulation, endothelial dysfunction, and systemic inflammatory responses [2,3]. These pathophysiologic changes have contributed to various major comorbidities [2,3,4,5]. OSA is a global health concern and economic burden, with estimates of the incidence of moderate or severe OSA reaching 6%-17% in the general adult population and up to 49% in the elderly population [4,6], with most of the cases being undiagnosed [7].

    No single physical examination finding can be used to diagnose or exclusively account for OSA. The gold standard diagnostic test for OSA relies on in-lab all-night polysomnography (PSG) [7,8]. However, this method is costly, and patients require continuous sleep monitoring by medical staff [2]; thus, limited availability of PSG may cause delays in diagnosing OSA and an increased disease burden [7,9].

    Several body measurements predispose patients to OSA. Body weight, body mass index (BMI), central body fat distribution, and large neck and waist circumferences are possible risk factors [2,6,10]. In addition, increased age and males are also at risk for OSA [2,6]. Despite the epidemiological reports demonstrating that age, sex, and BMI are associated with OSA severity, an integrated predicting model that uses only these three factors to assess OSA risk has not yet been developed.

    Machine learning approaches have expanded the ability to design efficient predictive and diagnostic tools and have gained popularity in many areas of medical research and clinical applications [11,12]. Logistic regression (LR) is one of the computational methods for probing multivariate data to detect the independent effects of a variable after adjusting for the contributions of the others [13]. An artificial neural network (ANN) is a supervised learning algorithm with advantages in analyzing nonlinear data, enhancing data interpretation, and designing effective predictive and diagnostic models to assist different clinical tasks [11,14].

    The present study investigated the clinical predictive validity of age, sex, and BMI for sleep apnea and applied LR and ANN algorithms to develop a risk-predicting model using these three parameters as predictors.

    This retrospective observational study used clinical data from the medical records of participants who had undergone PSG at the Sleep Center of Shuang Ho Hospital, Taipei Medical University, Taiwan, between January 2015 and January 2020.

    The inclusion criteria for this study were: (1) participants aged≥18 years and (2) who had undergone a full-night PSG. Participants were excluded if (1) they had received any OSA treatments or been admitted to the hospital for any invasive surgery for OSA, (2) they had undergone PSG for continuous positive airway pressure (CPAP) titration, (3) the total recording time was < 4 h [8], and (4) the total sleep time was < 3 h [15].

    Before the PSG examination, a baseline screening questionnaire was used to assess each participant's basic information, medication, and surgical history. The participants were measured for height, weight, and BMI (kg/m2) at the time of registration.

    All the participants underwent a full-night PSG recorded by qualified technicians. All the sleep signals were obtained using an Embla N7000 Recording System (Medcare Flaga, Reykjavik, Iceland) and Sandman Elite sleep diagnostic software (Natus Neurology Inc., Middleton, WI, USA).

    Experienced sleep technologists scored all sleep apnea and hypopnea events following the American Academy of Sleep Medicine 2012 updated guideline [8,16]. Each participant's apnea-hypopnea index (AHI) was calculated according to the total number of apnea and hypopnea events divided by the total sleeping time. Hypopnea was defined as a 30-89% reduction in nasal airflow for at least 10s, with arousal or at least a decrease of 3% oxygen saturation. Apnea was defined as a more than 90% reduction in oral airflow for at least 10s with or without arousal or decreased saturation.

    The severity of OSA was classified as none (AHI < 5), mild (AHI = 5-15), moderate (AHI = 15-30), and severe (AHI ≥ 30) [2,8]. To develop the predictive models for the risk of moderate-to-severe OSA, we divided the participants into none-and-mild (AHI < 15) and moderate-to-severe (AHI ≥ 15) OSA groups for further analysis.

    Taipei Medical University-Joint Institutional Review Board (TMU-JIRB: N201911007) approved the protocol of this study. Waivers of informed consent were approved by the TMU-JIRB (TMU-JIRB: N201911007) for this retrospective study involving the secondary analysis of existing data.

    All analyses were performed using STATISTICA software (TIBCO Software, Tulsa, OK, USA). Continuous data are presented as mean and standard deviation, whereas categorical variables are expressed as frequencies and percentages. We used Fisher's exact test to determine associations between two categorical variables and Pearson's chi-squared (X2) test to determine the statistical significance between the means of two continuous variables in the participants with none-to-mild and moderate-to-severe OSA. The one-way analysis of variance (ANOVA) was used to determine the statistical significance between the means of ≥ 3 independent groups. A two-tailed p-value less than 0.05 was considered statistically significant.

    All the LR and ANN models were established and modeled using STATISTICA software. In the models, age and BMI were set as continuous variables, and sex was designated as a categorical variable. For our ANNs, the employed architecture was feed-forward, combined with a back-propagation algorithm with multilayer perceptron (MLP). A bias term was used for each neural layer in the MLP. We set the number of neurons in the hidden layer from 1 to 50. We initially partitioned the eligible PSG data into a training-testing set and an independent validation set, containing ~80% and ~20% of the data, respectively. The data from the training-testing set were used to determine the generalizability of the LR and ANN analyses using five-fold cross-validation (Figure 1A). All participants were stratified by age, sex, BMI, and OSA severity to maintain consistent proportions of demographic data across the training-testing and validation sets. We reported the mean AUCs for the training and testing sets of the LR and ANN models through five-fold cross-validation. The model with the best AUC in the five testing sets was identified to have the best predictive performance and used as the final model. The final models were then applied to the independent validation set, and we calculated the accuracy, sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and AUCs of the validation set derived from the confusion matrix [17] (Figure 1A). The classification threshold was set to 0.5 for all models.

    Figure 1.  Study flowchart and ANN model. (A) Flowchart showing the study cohort selection and division. (B) Schematic representation of the structure of the artificial neural network (ANN) model. AHI, apnea-hypopnea index; ANN, artificial neural network; BMI, body mass index; OSA, obstructive sleep apnea.

    During the study period, 11, 100 participants underwent PSG in the Sleep Center. Of these, 1, 665 participants were excluded from the baseline analyses for the following reasons: they had been admitted to the hospital for invasive surgery or CPAP titration (n = 247), the total recording time < 4h (n = 1, 386), the total sleep time < 3h (n = 32), missing BMI information (n = 7) or missing AHI information (n = 6) (Figure 1A). Overall, 9, 422 participants (2918 women and 6504 men, mean age: 49.2 ± 14.0 years) who had undergone PSG were enrolled in this study. Of these, 5, 920 (62.8%) had an AHI of ≥15, indicating moderate-to-severe OSA, and 3, 502 (37.2%) had an AHI < 15, indicating none-to-mild OSA.

    A higher proportion of men have moderate-to-severe OSA than women (Table 1). Compared with participants with none-to-mild OSA, participants with moderate-to-severe OSA were older and had a higher BMI. In addition, participants with moderate-to-severe OSA had a shorter total sleep time, lower sleep efficiency, shorter rapid eye movement duration, and lower oxygen saturation during the PSG recording than those with none-to-mild OSA (Table 1).

    Table 1.  Demographic characteristics of the cohort.
    Variables Moderate-to-severe OSA
    (AHI ≥ 15)
    None-to-mild OSA
    (AHI < 15)
    p-value
    n 5, 920 3, 502
    Age 50.8 ± 13.5 46.6 ± 14.4 < 0.0001
    Female, n (%) 1266 (21.4) 1652 (47.2) < 0.0001
    Height (cm) 167.2 ± 8.5 164.4 ± 8.9 < 0.0001
    Weight (kg) 79.7±16.1 66.1 ± 13.1 < 0.0001
    BMI (kg/m2) 28.4 ± 4.8 24.3 ± 3.8 < 0.0001
    Polysomnography
    Total recording time (min) 366.6±11.9 367.8±13.4 < 0.0001
    Total sleep time (min) 275.8±60.2 281.4±61.1 < 0.0001
    Sleep efficiency (%) 75.2±16.4 76.5±16.5 0.0002
    Non-REM sleep (min) 240.6±51.7 242.2±51.3 0.133
    REM sleep (min) 35.2±22.5 39.1±23.7 < 0.0001
    AHI 44.7±23.8 6.4±4.5 < 0.0001
    SO2 mean (%) 94.5±2.6 96.6±1.3 < 0.0001
    Data are expressed as mean ± standard deviation. The p-value represents the comparison between the AHI ≥ 15 and AHI < 15 groups.

     | Show Table
    DownLoad: CSV

    A multivariate analysis was applied to adjust for age, sex, and BMI to determine the independent predictors of moderate-to-severe OSA. Table 2 presents the variables' crude and adjusted odds ratio (OR) values in the LR model. After adjustment for the associated variables, the regression analysis revealed that the female sex was associated with a 0.28-fold reduced risk of moderate-to-severe OSA (95% CI: 0.25-0.32; p < 0.0001), and every unit increase in age or BMI was associated with a 1.04-fold (95%CI: 1.04-1.05; p < 0.0001) or 1.27-fold (95%CI: 1.26-1.29; p < 0.0001) increase in the likelihood of moderate-to-severe OSA, respectively.

    Table 2.  Crude and adjusted odds ratios from the logistic regression analyses identifying associations between selected variables and OSA.
    Unadjusted Model Adjusted Model
    Variables OR 95% CI OR 95% CI
    Age 1.02*† 1.02-1.03 1.04*† 1.04-1.05
    Female 0.30* 0.28-0.33 0.28* 0.25-0.32
    BMI 1.27*† 1.25-1.28 1.27*† 1.26-1.29
    The odds ratio (OR) and 95% confidence intervals (CI) represent the association between the variables and the presence of moderate-to-severe OSA (AHI ≥ 15). *p < 0.0001; †Odds ratio of per unit changes.

     | Show Table
    DownLoad: CSV

    As mentioned, age, sex, and BMI were used as input parameters for developing the LR and ANN models for predicting moderate-to-severe OSA. As shown in Figure 1A, the 7328 participants in the training-testing group were randomly partitioned into training and testing sets for five-fold cross-validation. Table 3 showed that age, sex, BMI, and OSA severity distribution remained identical across the five testing sets and the independent validation sets.

    Table 3.  Description of the input and output variables within the testing datasets for cross-validation and the independent validation set of the models.
    Variable Testing set 1 Testing set 2 Testing set 3 Testing set 4 Testing set 5 Independent validation set p-value
    n 1466 1465 1466 1465 1466 2094
    AHI≥15, n (%) 920 (62.8) 919 (62.7) 920 (62.8) 919 (62.7) 920 (62.8) 1322 (63.1) 1.000
    Age 49.0 ± 0.4 49.0 ± 0.4 49.9 ± 0.4 49.1 ± 0.4 49.0 ± 0.4 49.3 ± 0.3 0.503
    Female, n (%) 483 (33.0) 438 (29.9) 472 (32.2) 445 (30.4) 431 (29.4) 649 (31.0) 0.266
    BMI (kg/m2) 26.9 ± 4.9 26.9 ± 4.9 26.9 ± 5.0 26.9 ± 5.0 26.8 ± 4.8 26.9 ± 4.9 0.973
    The p-value represents the multiple comparisons across the five testing and the independent validation sets.

     | Show Table
    DownLoad: CSV

    Following adequate training, the five testing sets that achieved the best prediction performance after five-fold cross-validation were ANN models with 3, 14, 25, 27, and 29 hidden neurons, respectively. The mean AUCs of the training sets were 0.808 ± 0.004 (Figure 2A) and 0.808 ± 0.015 for the testing sets (Figure 2B). Among the generated ANN models, the model containing 25 hidden neurons obtained the highest testing AUC and was selected as the final model. After being applied to the independent validation set, the ANN-based model achieved the accuracy, sensitivity, specificity, PPV, and NPV of 76.4%, 87.7%, 56.9%, 77.7%, 73.0%, with the AUC of 0.807, respectively. Results of the LR model revealed a mean AUC of 0.802 ± 0.004 for the training sets (Figure 2C) and 0.801 ± 0.015 for the testing sets (Figure 2D). After being validated to the independent validation set, the accuracy, sensitivity, specificity, PPV, and NPV were 76.3%, 87.5%, 57.0%, 77.7%, and 72. 7% with an AUC of 0.806 for the LR-based model. Figure 3 shows the ROC curve and AUC of the final ANN and LR models applied to the independent validation set.

    Figure 2.  AUC-defined performance of the predicting models. The ROC curves and AUCs of the (A) training, (B) testing sets of the logistic regression (LR) model, and (C) training, (D) testing sets of the artificial neural network (ANN) model for predicting moderate-to-severe OSA in participants who have undergone polysomnography by using age, sex, and BMI as the predictors. The AUC values represent the mean ± standard deviation of the AUCs of five training and testing sets.
    Figure 3.  The ROC curves and AUC value of the final LR and ANN models applied to the independent validation set (n = 2094). The pink line represents the ROC curve of the artificial neural network (ANN) model, and the blue line represents the ROC curve of the logistic regression (LR) model.

    The present study validated age, sex, and BMI as parameters for use in LR- and ANN models to estimate the risk of OSA in participants who had undergone PSG. The established LR- and ANN-based model integrating the three parameters achieved a notable prediction performance with an AUC of 0.806 and 0.807 for predicting moderate-to-severe OSA in the validation set. The relatively high AUC values indicated that the models could distinguish participants between groups of interest, namely, participants with none-to-mild and those with moderate-to-severe OSA. Patients with OSA are prone to various complications and require interventions [2,7]; our models may contribute to further understanding these three variables' fundamental role in assessing OSA risk.

    OSA has been recognized as a global health and economic burden, and early identification and diagnosis of OSA is a crucial issue in preventive medicine [4,6,7]. Since standard diagnostic tests for OSA rely on expensive in-laboratory PSG, which can be challenging to schedule in high-demand situations, identifying crucial factors and developing a simple and reliable tool to estimate the risk of OSA is clinically relevant [7,9].

    Anatomical and physiological differences in the oropharynx, larynx and upper respiratory tract, changes in respiratory physiology associated with different sleep structures, and differences in the distribution of adipose tissue have been suggested as reasons for sex differences in developing OSA [18]. Deterioration of the collapsibility of the pharyngeal airway and the sensitivity of the respiratory control system may contribute to the pathogenesis of OSA with aging [2,18]. Hormonal changes in menopausal women have been linked to an increased risk of OSA and a decreased sex difference in OSA with age [19,20]. The prevalence of OSA has been reported to increase 1.4 to 3.2-fold with each 10-year increment in age, and men had a 1.7 to 3.0-fold higher odds of developing OSA than women [2,6]. These findings support the rationale for using age and sex as parameters in a generic prediction model to predict OSA risk.

    Obesity is a well-recognized risk factor for developing OSA, and a higher BMI is associated with greater severity of OSA for both sexes [19,21]. As already mentioned, anthropometric profiles, such as body weight, BMI, central body fat distribution, and neck and waist circumference, are predisposed to the risk of OSA [2,6,10]. In patients with OSA, occlusion of the nasal and oro-pharynx may occur during sleep when the tongue and epiglottis move backward and lean together. Individuals with obesity have excess soft tissue around the collapsible pharynx, making the narrow airway more prone to collapse [2,3,21]. In line with the findings of other studies, we observed that the participants with moderate-to-severe OSA had a higher BMI. Although other anthropometric profiles have been associated with OSA, neck and waist circumference were considered surrogates for BMI [4,21,22]. The study by Huxley et al. has also suggested that "measures of central obesity were superior to BMI as discriminators of cardiovascular risk factors, although the differences were small and unlikely to be of clinical relevance" [23]. Furthermore, the same study demonstrated that "combining BMI with any measure of central obesity did not improve the discriminatory capability of the individual measures" [23]. In this study, we proposed BMI as a simple, well-known, and readily available predictor of OSA for the prediction models.

    Despite the epidemiological reports indicating that age, sex, and BMI are associated with OSA, an integrated predicting model using only these three attributes to assess OSA risk has not yet been developed. Advances in machine learning techniques have improved the ability of researchers and clinical applications to analyze massive amounts of data and develop appropriate diagnostic and predictive methods [24,25,26,27]. Several studies have applied machine learning algorithms to predict OSA risk. However, most of these studies have incorporated more complex parameters, such as sleep questionnaires [28,29,30], at-home PSG [15], cephalometric images [31], neck and waist circumference measurements [28,29,32,33], and facial-cervical measurements [28,32], into the predictive models.

    In this context, our study used a large hospital-based sample of adults, allowing accurate comparison between none-to-mild and moderate-to-severe cases. The generated LR- and ANN-boosted models exhibited relatively high AUCs with adequate accuracy, sensitivity, PPV, and NPV, revealing the ability of these models to distinguish between the groups of interest, namely individuals with moderate-to-severe and none-to-mild OSA. Furthermore, our models used simple and easy-to-obtain parameters (age, sex, and BMI), which may provide a direction for developing predictive models under conditions of limited and constrained information.

    Both LR and ANN are widely used models in biomedicine. This study's LR and ANN models achieved similar predictive validity in predicting OSA in patients undergoing PSG, which may be related to our data structure. ANN has the advantage of analyzing complex, multidimensional and nonlinear variables [34,35,36,37]. When the predictor-outcome relationships in the dataset are relatively simple and linear, LR using linear combinations of variables may achieve performance that approximates that of ANN [36,38]. This consideration was contextualized in the conclusion by Dreiseitl and Ohno-Machado [38] that "there is no single algorithm that performs better than all other algorithms on any given data set and application area." and also summarized by Sargent [34] that "both methods should continue to be used and explored in a complementary manner."

    The present study has several limitations. First, in our cohort, the incidence of moderate-to-severe OSA was 62.8%, which is higher than the reported prevalence of OSA in the general population [4,19,20,39]. In clinical practice, patients referred for a PSG exam are usually expected to have OSA or another sleep disorder; therefore, the participant population in this study cannot be representative of the general population with a very low prevalence of OSA. Second, the data set used in this study was limited to a single center and included only Asian ethnic groups. Therefore, the trained model may not be generalizable beyond the institution or other ethnic groups [28,29]. Third, this retrospective observational study did not exclude undiagnosed medical conditions associated with OSA, such as neurological, cardiovascular, and pulmonary disorders, which may affect the model's predictive performance. Our model did not contain detailed anthropometric imaging or measurements, which may limit the model's ability to distinguish the disease-specific causes of OSA. Finally, the relatively low specificity for the models (57.0% and 56.9%) suggests a higher rate of false positives, meaning that a healthy person is more likely to receive an incorrect diagnosis of OSA using the proposed tool. As Jonas et al. [30] concluded, "There is uncertainty about the accuracy or clinical utility of all potential screening tools." a PSG to verify OSA severity and confirm the diagnosis is warranted for patients with clinically suspected sleep apnea.

    Our study validated the diagnostic relevance of age, sex, and BMI for OSA. The generated LR- and ANN-based model effectively predicted OSA risk using these three simple parameters. Our results, therefore, may improve insights into OSA-related essential risk factors and provide knowledge to facilitate individualized risk stratification and formulate diagnostic decisions for sleep apnea.

    This work was funded by the National Science and Technology Council (NSTC 111-2314-B-038-132-MY3) to CC Chung.

    The authors declare there is no conflict of interest.



    [1] Y. Masubuchi, Simulating the flow of entangled polymers, Annu. Rev. Chem. Biomol. Eng., 5 (2014), 11-33.
    [2] T. Sato, A review on transport phenomena of polymeric liquids, J. Soc. Rheol. Jpn., 48 (2020), 1-14.
    [3] R. B. Bird, R. C. Armstrong, O. Hassager, Dynamics of polymeric liquids, Volume 1 fluid mechanics, 3 Eds., John Wiley & Sons, Inc, 1987.
    [4] M. M. Cross, Rheology of non-Newtonian fluids: A new flow equation for pseudoplastic systems, J. Colloid. Sci., 20 (1965), 417-437.
    [5] J. D. Ferry, Viscoelastic properties of polymers, 3 Eds., John Wiley & Sons, Inc, 1980.
    [6] H. H. Chiang, C. A. Hieber, K. K. Wang, A unified simulation of the filling and postfilling stages in injection molding. Part Ⅰ: Formulation, Polym. Eng. Sci., 31 (1991), 116-124.
    [7] R. E. Otmani, M. Zinet, M. Boutaous, H. Benhadid, Numerical simulation and thermal analysis of the filling stage in the injection molding process: Role of the mold-polymer interface, J. Appl. Polym. Sci., 121 (2011), 1579-1592.
    [8] A. Khalifeh, J. R. Clermont, Numerical simulations of non-isothermal three-dimensional flows in an extruder by a finite-volume method, J. Non-Newton. Fluid Mech., 126 (2005), 7-22.
    [9] P. S. B. Zdanski, M. Vaz Jr, Three-dimensional polymer melt flow in sudden expansions: Nonisothermal flow topology, Int. J. Heat Mass Tran., 52 (2009), 3585-3594.
    [10] A. K. Doufas, A. J. McHugh, C. Miller, Simulation of melt spinning including flow-induced crystallization Part Ⅰ. Model development and predictions, J. Non-Newton. Fluid Mech., 92 (2000), 27-66.
    [11] G. W. M. Peters, F. P. T. Baaijens, Modeling of non-isothermal viscoelastic flows, J. Non-Newton. Fluid Mech., 68 (1997), 205-224.
    [12] K. Kunisch, X. Marduel, Optimal control of non-isothermal viscoelastic fluid flow, J. Non-Newton. Fluid Mech., 88 (2000), 261-301.
    [13] A. Wachs, J. R. Clermont, Non-isothermal viscoelastic flow computations in an axisymmetric contraction at high Weissenberg numbers by a finite volume method, J. Non-Newton. Fluid Mech., 95 (2000), 147-184.
    [14] M. Laso, H. C. Öttinger, Calculation of viscoelastic flow using molecular models: the CONNFFESSIT approach, J. Non-Newton. Fluid Mech., 47 (1993), 1-20.
    [15] R. Keunings, Micro-macro methods for the multiscale simulation of viscoelastic flow using molecular models of kinetic theory, Rheology Reviews, 2004 (2004), 67-98.
    [16] M. A. Hulsen, A. P. G. van Heel, B. H. A. A. van den Brule, Simulation of viscoelastic flows using Brownian configuration fields, J. Non-Newton. Fluid Mech., 70 (1997), 79-101.
    [17] P. Halin, G. Lielens, R. Keunings, V. Legat, The Lagrangian particle method for macroscopic and micro-macro viscoelastic flow computations, J. Non-Newton. Fluid Mech., 79 (1998), 387-403.
    [18] W. E, B. Engquist, The heterogeneous multiscale methods, Commun. Math. Sci., 1 (2003), 87-132.
    [19] S. Yasuda, R. Yamamoto, A model for hybrid simulations of molecular dynamics and computational fluid dynamics, Phys. Fluids, 20 (2008), 113101.
    [20] T. Murashima, T. Taniguchi, Multiscale Lagrangian fluid dynamics simulation for polymeric fluid, J. Polym. Sci. B Polym. Phys., 48 (2010), 886-893.
    [21] T. Murashima, T. Taniguchi, Multiscale simulation of history-dependent flow in entangled polymer melts, Europhys. Lett., 96 (2011), 18002.
    [22] T. Murashima, T. Taniguchi, Flow-history-dependent behavior of entangled polymer melt flow analyzed by multiscale simulation, J. Phys. Soc. Jpn., 81 (2012), 1-7.
    [23] T. Sato, K. Harada, T. Taniguchi, Multiscale simulations of flows of a well-entangled polymer melt in a contraction-expansion channel, Macromolecules, 52 (2019), 547-564.
    [24] T. Sato, T. Taniguchi, Multiscale simulations for entangled polymer melt spinning process, J. NonNewton. Fluid Mech., 241 (2017), 34-42.
    [25] M. Doi, J. Takimoto, Molecular modeling of entanglement, Phil. Trans. R. Soc. Lond. A, 361 (2003), 641-652.
    [26] K. Kremer, G. S. Grest, Dynamics of entangled linear polymer melts: A molecular-dynamics simulation, J. Chem. Phys., 92 (1990), 5057-5086.
    [27] S. Yasuda, R. Yamamoto, Synchronized molecular-dynamics simulation via macroscopic heat and momentum transfer: An application to polymer lubrication, Phys. Rev. X, 4 (2014), 041011.
    [28] S. Yasuda, R. Yamamoto, Synchronized molecular-dynamics simulation for the thermal lubrication of a polymeric liquid between parallel plates, Comput. Fluids, 124 (2016), 185-189.
    [29] S. Yasuda, Synchronized molecular-dynamics simulation of the thermal lubrication of an entangled polymeric liquid, Polymers, 11 (2019), 131.
    [30] T. Sato, T. Taniguchi, Rheology and entanglement structure of well-entangled polymer melts: A slip-link simulation study, Macromolecules, 52 (2019), 3951-3964.
    [31] F. Habla, A. Woitalka, S. Neuner, O. Hinrichsen, Development of a methodology for numerical simulation of non-isothermal viscoelastic fluid flows with application to axisymmetric 4:1 contraction flows, Chem. Eng. J., 207-208 (2012), 772-784.
    [32] H. Pol, S. Banik, B. L. Azad, S. Thete, P. Doshi, A. Lele, Nonisothermal analysis of extrusion film casting process using molecular constitutive equations, Rheol. Acta, 53 (2014), 85-101.
    [33] M. Zatloukal, Measurements and modeling of temperature-strain rate dependent uniaxial and planar extensional viscosities for branched LDPE polymer melt, Polymer, 104 (2016), 258-267.
    [34] X. Zhuang, J. Ouyang, W. M. Li, Y. G. Li, Three-dimensional simulations of non-isothermal transient flow and flow-induced stresses during the viscoelastic fluid filling process, Int. J. Heat Mass Tran., 104 (2017), 374-391.
    [35] P. Gao, X. Wang, J. Ouyang, Numerical investigation of non-isothermal viscoelastic filling process by a coupled finite element and discontinuous Galerkin method, Int. J. Heat Mass Tran., 140 (2019), 227-242.
    [36] M. Doi, S. Edwards, The theory of polymer dynamics, Oxford University Press, 1986.
    [37] G. Marrucci, Dynamics of entanglements: A nonlinear model consistent with the Cox-Merz rule, J. Non-Newton. Fluid Mech., 62 (1996), 279-289.
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5597) PDF downloads(690) Cited by(2)

Figures and Tables

Figures(12)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog