European Journal of Mechanics B-fluids最新文献

{"title":"Derivation and numerical resolution of 2D shallow water equations for multi-regime flows of Herschel–Bulkley fluids","authors":"","doi":"10.1016/j.euromechflu.2024.07.010","DOIUrl":"10.1016/j.euromechflu.2024.07.010","url":null,"abstract":"<div><p>This paper presents mathematical modelling and simulation of thin free-surface flows of viscoplastic fluids with a Herschel–Bulkley rheology over complex topographies with basal perturbations. Using the asymptotic expansion method, depth-averaged models (lubrication and shallow water type models) are derived for 3D (three-dimensional) multi-regime flows on non-flat inclined topographies with varying basal slipperiness. Starting from the Navier–Stokes equations, two flow regimes corresponding to different balances between shear and pressure forces are presented. Flow models corresponding to these regimes are calculated as perturbations of the zeroth-order solutions. The classical reference models in the literature are recovered by considering their respective cases on a flat-inclined surface. In the second regime case, a pressure term is non-negligible. Mathematically, it leads to a corrective term to the classical regime equations. Flow solutions of the two regimes are compared; the difference appears in particular in the vicinity of sharp changes of slopes. Nonetheless, both regime models are compared with experiments and are found to be in good agreement. Furthermore, numerical examples are shown to illustrate the robustness of the present shallow water models to simulate viscoplastic flows in 3D and over an inclined topography with local perturbations in basal elevation and basal slipperiness. The derived models are adequate for direct (engineering and geophysical) applications to real-world flow problems presenting Herschel–Bulkley rheology like lava and mud flows.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-08-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142151785","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Evaluation on different volume of fluid methods in unstructured solver under the optimized condition","authors":"","doi":"10.1016/j.euromechflu.2024.07.016","DOIUrl":"10.1016/j.euromechflu.2024.07.016","url":null,"abstract":"<div><p>We compared the accuracy of volume of fluid (VOF) methods in unstructured solvers using the following five different methods: 1 - the algebraically compressive VOF method, 2 – simple coupled VOF method with Level Set (S-CLSVOF) method, 3 - interface-compressing VOF method incorporated with Laplacian filter (VOFL), 4 - isoAdvector method, and 5 - isoAdvector method incorporated with Laplacian filter (isoAdvectorL) by incorporating them into OpenFOAM®, an open-source software. To evaluate these methods under proper conditions, we compared the calculation accuracy using the optimized parameters, which are explored by Bayesian optimization. The test cases for advection accuracy of volume fraction and for imbalance of surface tension force in static multiphase fluid fields were considered. In this study, we found that the compression parameters and maximum Courant number should be adjusted to obtain high accuracy simulation according to the simulation condition in VOF and S-CLSVOF method. In VOFL and isoAdvectorL methods, the spurious current can be extremely reduced, which means that these methods are suitable for slow flow with higher Laplace number conditions.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141951680","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Thermal convection of a liquid metal under an alternating magnetic field","authors":"","doi":"10.1016/j.euromechflu.2024.07.015","DOIUrl":"10.1016/j.euromechflu.2024.07.015","url":null,"abstract":"<div><p>The objective of this work is to measure the heat transfer of a liquid metal in a cylindrical cell under the conjugate effects of a temperature difference and a Lorentz force generated by an alternating current in a coil. The experimental results are compared to recent direct numerical simulations (DNS) (Guillou et al., 2022). 25 experiments are performed for a large range of frequency <span><math><mi>f</mi></math></span>, ac intensity amplitude <span><math><msub><mrow><mi>I</mi></mrow><mrow><mn>0</mn></mrow></msub></math></span> and temperature difference between the top and bottom walls <span><math><mrow><mi>Δ</mi><msub><mrow><mi>T</mi></mrow><mrow><mn>0</mn></mrow></msub></mrow></math></span>: <span><math><mrow><mn>15</mn><mo>≤</mo><mi>f</mi><mo>≤</mo><mn>1000</mn><mspace></mspace><mi>Hz</mi></mrow></math></span>, <span><math><mrow><mn>2</mn><mo>≤</mo><msub><mrow><mi>I</mi></mrow><mrow><mn>0</mn></mrow></msub><mo>≤</mo><mn>67</mn></mrow></math></span> A and <span><math><mrow><mn>6</mn><mo>≤</mo><mi>Δ</mi><msub><mrow><mi>T</mi></mrow><mrow><mn>0</mn></mrow></msub><mo>≤</mo><mn>11</mn></mrow></math></span> K. In these experiments, the Hartmann number <span><math><mrow><mi>H</mi><mi>a</mi></mrow></math></span>, the shielding parameter <span><math><msub><mrow><mi>S</mi></mrow><mrow><mi>ω</mi></mrow></msub></math></span> and Rayleigh number <span><math><mrow><mi>R</mi><mi>a</mi></mrow></math></span> vary in the following range: <span><math><mrow><mn>6</mn><mo>≤</mo><mi>H</mi><mi>a</mi><mo>≤</mo><mn>200</mn></mrow></math></span>, <span><math><mrow><mn>1</mn><mo>≤</mo><msub><mrow><mi>S</mi></mrow><mrow><mi>ω</mi></mrow></msub><mo>≤</mo><mn>70</mn></mrow></math></span>, <span><math><mrow><mn>2</mn><mo>.</mo><mn>3</mn><mo>×</mo><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mn>6</mn></mrow></msup><mo>≤</mo><mi>R</mi><mi>a</mi><mo>≤</mo><mn>4</mn><mo>.</mo><mn>1</mn><mo>×</mo><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mn>6</mn></mrow></msup></mrow></math></span>. The experiments with an ac magnetic field are compared with the Rayleigh–Bénard convection (RBC) experiments under the same thermal conditions. Three rings of thermocouples allow characterizing the fluid temperature distribution during the convection. The heat flux at the bottom and top walls are also measured. We observe a very good agreement between the experimental results and the DNS results. As previously shown by numerical simulations, a master curve of <span><math><mrow><mi>N</mi><mi>u</mi><mo>/</mo><mi>P</mi><msub><mrow><mi>e</mi></mrow><mrow><mi>ω</mi></mrow></msub></mrow></math></span> vs. <span><math><mrow><msub><mrow><mi>Q</mi></mrow><mrow><mi>J</mi></mrow></msub><mo>/</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>c</mi></mrow></msub></mrow></math></span> allows predicting the evolution of the heat transfer under different conditions of temperature difference and Lorentz force. Here <span><math><mrow><mi>N</mi><mi>u</mi></mrow></math></span> and <span><math><mrow><mi>P</mi><m","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141882527","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"High-fidelity simulations of Richtmyer–Meshkov flows triggered by a forward-pentagonal bubble with different Atwood numbers","authors":"","doi":"10.1016/j.euromechflu.2024.07.013","DOIUrl":"10.1016/j.euromechflu.2024.07.013","url":null,"abstract":"<div><p>In fluid dynamics, the Atwood number is a dimensionless parameter that quantifies the density difference between two fluids. It is calculated as <span><math><mrow><mi>A</mi><mi>t</mi><mo>=</mo><mrow><mo>(</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>−</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow></math></span>, where <span><math><msub><mrow><mi>ρ</mi></mrow><mrow><mn>1</mn></mrow></msub></math></span> and <span><math><msub><mrow><mi>ρ</mi></mrow><mrow><mn>2</mn></mrow></msub></math></span> represent the densities of the respective fluids. This research employs high-fidelity numerical simulations to examine the Atwood number impacts on Richtmyer–Meshkov (RM) flows triggered by a shocked forward-pentagonal bubble. Five distinct gases — <span><math><msub><mrow><mtext>SF</mtext></mrow><mrow><mn>6</mn></mrow></msub></math></span>, Kr, Ar, Ne, and He — are considered within the forward-pentagonal bubble, encompassed by <span><math><msub><mrow><mtext>N</mtext></mrow><mrow><mn>2</mn></mrow></msub></math></span> gas. In these simulations, a third-order discontinuous Galerkin approach is applied to solve a two-dimensional set of compressible Navier–Stokes-Fourier (NSF) equations for two-component gas flows. To discretize space, hierarchical modal basis functions based on orthogonal-scaled Legendre polynomials are employed. This approach simplifies the NSF equations into a set of ordinary differential equations over time, which are solved using an explicit third-order SSP Runge–Kutta algorithm. The numerical results highlight the notable impact of the Atwood number on the evolution of RM flows in the shocked forward-pentagonal bubble, a phenomenon not previously reported in the literature. The Atwood number exerts a significant influence on the flow patterns, leading to intricate wave formations, shock focusing, jet generation, and interface distortion. Moreover, a comprehensive analysis of the these impact elucidates the mechanisms driving vorticity formation during the interaction process. Additionally, the study conducts a thorough quantitative examination of the Atwood number impacts on the flow fields based on integral quantities and interface features.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0997754624001110/pdfft?md5=3f2bce669d02ed1d1a574cae47d7a3d3&pid=1-s2.0-S0997754624001110-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141844163","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Recent advances in the analysis of turbulent superstructures","authors":"Jörg Schumacher, Wolfgang Schröder","doi":"10.1016/j.euromechflu.2024.07.014","DOIUrl":"https://doi.org/10.1016/j.euromechflu.2024.07.014","url":null,"abstract":"","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.6,"publicationDate":"2024-07-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141784318","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A study on blast wave diffractions and the dynamics of associated vortices inside different grooves kept at various lateral distances from the shock tube","authors":"","doi":"10.1016/j.euromechflu.2024.07.012","DOIUrl":"10.1016/j.euromechflu.2024.07.012","url":null,"abstract":"<div><p>Diffraction is a fundamental phenomenon that occurs when blast or shock waves pass over sudden discontinuous surfaces. It generates a complex flow field consisting of diffracted waves, expansion waves, slipstream, contact surface, and an unstable shear layer, in addition to emitting acoustic waves. In this study, we investigated the diffraction of a blast wave passing over a series of grooved structures with different aspect ratios and geometrical shapes (rectangular, circular, and triangular) using high-speed shadowgraph images. The blast wave Mach number considered in our investigation is 1.34. The grooves feature leading-edge geometrical variations such as rectangular, circular arc, and wedge shapes positioned at various lateral locations from the exit of the shock tube. The aspect ratios of the rectangular grooves vary from 0.33, 0.5, and 0.67. The circular and triangular grooves have an aspect ratio of 0.33. The trajectories and velocities of the primary vortex are calculated by tracking the location of the vortex in the shadowgraph images. Our observations revealed that a large portion of the incident blast wave is abducted inside the groove as the aspect ratio increases in rectangular grooves, resulting in better attenuation of the blast wave. The grooves, which have circular shapes, produced weaker diffraction, which resulted in delayed and weak primary vortex. The triangular grooves produced the strongest primary vortex and have the highest attenuation characteristics among other grooves. The strength and trajectory of the primary vortex formed over the grooves strongly depend on the aspect ratio and the curvature of the leading edge for a given Mach number. Vortices generated from rectangular and triangular grooves exhibit considerable strength and longevity.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141784317","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Detailed 3D URANS analysis of two-phase flow in an airlift pump","authors":"","doi":"10.1016/j.euromechflu.2024.07.011","DOIUrl":"10.1016/j.euromechflu.2024.07.011","url":null,"abstract":"<div><p>An airlift pump is a vertical tube that utilizes the buoyant effects of a gas to lift a liquid. Unlike a standard mechanical pump, the liquid flow rate through the airlift pump is not directly controlled; rather, it depends on the supplied gas flow rate, the tube length and diameter, and the relative height of the liquid supply free surface (submergence ratio). The present study uses the commercial CFD code ANSYS CFX to model the isothermal, 3D, transient flow in an airlift pump using water and air. The model applies pressure boundary conditions at both ends of the tube and specifies the mass flow rate of air through multiple openings in the side of the tube. The bottom of the tube is an inlet of water only and the outlet is a two-phase flow opening. A time-dependent, homogeneous, VOF two-phase RANS CFD modelling approach is used with the air treated as an ideal gas. This work found that a complete 3D domain was necessary for consistent prediction of the airlift performance and physically realistic two-phase flow structures. Statistical analysis of the two-phase flow structures was applied to characterize airlift pump instability and better understand the physics of the airlift pump.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141691324","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Modeling drag coefficients of spheroidal particles in rarefied flow conditions","authors":"","doi":"10.1016/j.euromechflu.2024.07.008","DOIUrl":"10.1016/j.euromechflu.2024.07.008","url":null,"abstract":"<div><p>Transport of particles in flows is often modeled in a combined Eulerian–Lagrangian framework. The flow is evaluated on an Eulerian grid, while particles are modeled as Lagrangian points whose positions and velocities are evolved in time, resulting in particle trajectories embedded in the time-dependent flow field. The method essentially resolves the flow field in complex geometries in detail but uses a closure model for the particle dynamics aimed at including the essential particle–fluid interactions at the cost of detailed small-scale physics. Rarefaction effects are usually included through the phenomenological Cunningham correction on the drag force experienced by the particles. In this Lagrangian point-particle approach, any explicit reference to the finite size and the shape of the particles, and their local orientation in the flow field, is typically ignored. In this work we aim to address this gap by deriving, from fully-resolved Direct Simulation Monte Carlo (DSMC) studies, heuristic or closure models for the drag force acting on prolate and oblate spheroidal particles with different aspect ratios, and a fixed orientation, in uniform ambient rarefied flows covering the transition regime between the continuum and free-molecular limits. These closure models predict the drag in the transition regime for all considered parameter settings (validated with DSMC data). The continuum limit is enforced a priori and we retrieve the free-molecular limit with reasonable accuracy (based on comparisons with literature data). We also include in the models the capability to predict effects related to basic gas-surface interactions via the tangential momentum accommodation coefficient. We furthermore assess the validity of the proposed closure model for particle dynamics in proximity to solid walls. This investigation extends our previous work, which focused on small aspect ratio spheroids with exclusively diffusive gas-surface interactions [see Livi et al. (2022)]. The derived models are obtained for isothermal, subsonic flows relevant for particle contamination control in semiconductor manufacturing.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0997754624000967/pdfft?md5=02cbba47811cf272f4b0537ed161d5ab&pid=1-s2.0-S0997754624000967-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141693336","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A comprehensive comparison of passive flow controls on the wind turbine blade lift and drag performances: A CFD approach","authors":"","doi":"10.1016/j.euromechflu.2024.07.009","DOIUrl":"10.1016/j.euromechflu.2024.07.009","url":null,"abstract":"<div><p>Flow control techniques used on wind turbines have been shown to significantly increase energy generation when compared to traditional wind turbines. Although various flow control methods have been introduced in the last two decades, the comparison between these methods is still the least conducted by researchers. Therefore, the present study aims to evaluate the performance of an airfoil utilizing both single and dual passive flow control methods, such as droop, flap, microcylinder, slot, and spoiler with optimal parameters. In this study, a numerical model was developed and applied with the same boundary conditions as those in the experiment. The results of the developed numerical simulation were then validated with experimental and other numerical studies. Mosaic mesh was utilized and the results were compared with conventional mesh types. Even though the mosaic mesh requires a lower number of computational elements, it demonstrated higher computational accuracy when compared to hexcore, polyhedra and tetrahedral type meshes. After obtaining an accurate numerical model, parametric studies were then conducted. The findings mainly highlighted that the airfoil with a microcylinder consistently generated higher performance than droop, flap, spoiler, slot and conventional airfoil. The mean relative improvement was about 2.6%. In an extensive study, eight combinations of flow controls were proposed and evaluated. The highest performances were achieved with the combination of microcylinder and flap, up to 27.9% and the combination of microcylinder and slot, reaching up to 50.2%, for low and high AOAs, respectively.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141693985","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Dynamics and applications of finite-size fibre-like objects in turbulent flows","authors":"","doi":"10.1016/j.euromechflu.2024.07.007","DOIUrl":"10.1016/j.euromechflu.2024.07.007","url":null,"abstract":"<div><p>This review delves into the dynamics of fibre-laden turbulent flows, a field that has garnered substantial attention due to its relevance in both natural and engineering contexts. The focus here is mainly on finite-size fibres, those exceeding the Kolmogorov scale, diverging from the commonly studied smaller ones. The study synthesises current understanding of the behaviour and organisation of both rigid and flexible finite-size fibres within turbulent flows, underscoring the added complexity these anisotropic particles introduce compared to their spherical counterparts. The influence of the length, the curvature and the inertia on the dynamics of rigid and flexible fibres is addressed. Fibre-based novel experimental methods, such as Fibre Tracking Velocimetry, are highlighted. Ultimately, this paper seeks to provide a clearer picture of the intricate dynamics at play in fibre-laden turbulent flows and their practical implications in various fields.</p></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-07-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0997754624000955/pdfft?md5=aa7ebb504cf4ff04056cad591a8a3955&pid=1-s2.0-S0997754624000955-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141700730","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}