Pub Date : 2022-11-10DOI: 10.1080/23324309.2022.2150645
A. Z. Bozkır, D. C. Sahni, R. G. Türeci
Abstract In this study, the albedo problem is investigated with three different methods, method, method, and SVD method. The first two methods are used for the comparison of the SVD method results. Therefore, the main aim of this study is to study of the SVD method for albedo problem. The recently improved method is based on usage of the transformation of the integral part to a sum in the one-speed, source free, and homogeneous medium neutron transport equation. This sum term is written with two different ways. One of them is to use midpoint approximation and the second is to use linear approximation. The numerical results of the SVD method are compared with both the method and the method results. Another different situation of this study is to use a series of the Legendre polynomials for the outgoing flux over the surface.
{"title":"Half-Space Albedo Problem with Pure-Triplet Scattering and Legendre Polynomial Outgoing Flux","authors":"A. Z. Bozkır, D. C. Sahni, R. G. Türeci","doi":"10.1080/23324309.2022.2150645","DOIUrl":"https://doi.org/10.1080/23324309.2022.2150645","url":null,"abstract":"Abstract In this study, the albedo problem is investigated with three different methods, method, method, and SVD method. The first two methods are used for the comparison of the SVD method results. Therefore, the main aim of this study is to study of the SVD method for albedo problem. The recently improved method is based on usage of the transformation of the integral part to a sum in the one-speed, source free, and homogeneous medium neutron transport equation. This sum term is written with two different ways. One of them is to use midpoint approximation and the second is to use linear approximation. The numerical results of the SVD method are compared with both the method and the method results. Another different situation of this study is to use a series of the Legendre polynomials for the outgoing flux over the surface.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"393 - 417"},"PeriodicalIF":0.7,"publicationDate":"2022-11-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44410331","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}
Pub Date : 2022-09-19DOI: 10.1080/23324309.2022.2141780
R. G. Türeci
Abstract Recently developed the Anlı-Güngör scattering function is applied to the Milne problem in this study. It is the function of Legendre polynomials with t parameter which can be called as scattering parameter. The extrapolation distance values are calculated for varying t parameters and varying secondary neutron numbers with the HN method. The numerical results are calculated by 40th approximation in order to get high precision results. The calculated extrapolation distance results are very convergence. Since there is no a similar study for the Milne problem with the Anlı-Güngör scattering in the literature, the interpolation technique is applied to the HN method results. Thus, the isotropic scattering results could be determined by t → 0 limit.
{"title":"The Milne Problem with the Anlı-Güngör Scattering","authors":"R. G. Türeci","doi":"10.1080/23324309.2022.2141780","DOIUrl":"https://doi.org/10.1080/23324309.2022.2141780","url":null,"abstract":"Abstract Recently developed the Anlı-Güngör scattering function is applied to the Milne problem in this study. It is the function of Legendre polynomials with t parameter which can be called as scattering parameter. The extrapolation distance values are calculated for varying t parameters and varying secondary neutron numbers with the HN method. The numerical results are calculated by 40th approximation in order to get high precision results. The calculated extrapolation distance results are very convergence. Since there is no a similar study for the Milne problem with the Anlı-Güngör scattering in the literature, the interpolation technique is applied to the HN method results. Thus, the isotropic scattering results could be determined by t → 0 limit.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"354 - 371"},"PeriodicalIF":0.7,"publicationDate":"2022-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46781000","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}
Pub Date : 2022-09-19DOI: 10.1080/23324309.2022.2141779
D. Gülderen, D. C. Sahni, R. G. Türeci, A. Aydιn
Abstract The Milne problem is investigated for the linear-triplet anisotropic scattering with HN method in this study. The scattering function is the linear combination of linear anisotropic scattering and triplet anisotropic scattering in Mika’s scattering. The positivity condition is needed to find physical results since the scattering function defines the scattering probabilities. It defines the relationship between the scattering coefficients. HN method is based on using the Case method. Therefore, properties of the Case method should be derived for the scattering function. The linear anisotropic scattering is the dominant scattering according to the calculated results.
{"title":"The Milne Problem for Linear-Triplet Anisotropic Scattering with HN Method","authors":"D. Gülderen, D. C. Sahni, R. G. Türeci, A. Aydιn","doi":"10.1080/23324309.2022.2141779","DOIUrl":"https://doi.org/10.1080/23324309.2022.2141779","url":null,"abstract":"Abstract The Milne problem is investigated for the linear-triplet anisotropic scattering with HN method in this study. The scattering function is the linear combination of linear anisotropic scattering and triplet anisotropic scattering in Mika’s scattering. The positivity condition is needed to find physical results since the scattering function defines the scattering probabilities. It defines the relationship between the scattering coefficients. HN method is based on using the Case method. Therefore, properties of the Case method should be derived for the scattering function. The linear anisotropic scattering is the dominant scattering according to the calculated results.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"329 - 353"},"PeriodicalIF":0.7,"publicationDate":"2022-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42137859","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}
Pub Date : 2022-09-19DOI: 10.1080/23324309.2022.2145311
A. Tiwari, Akshay K. Singh
Abstract In this article, we have presented a mathematical analysis to study the peristaltic flow through vertical annulus filled with a porous material bounded by two concentric uniform tubes. This analysis can serve as a model which may help in understanding the mechanism of physiological flows and heat transfer in a vertical annulus subject to differentially heating in the presence of a temperature dependent sink. The inner tube is uniform and rigid, while the outer tube has a sinusoidal wave traveling down its wall. It is analyzed in a wave frame of reference moving with velocity of the wave under a zero Reynolds number and long wavelength approximation in the presence of a temperature dependent sink. The analytical solutions are obtained for temperature, axial velocity, stream function and axial pressure-gradient. We also present numerical integration in order to analyze the pressure rise and frictional forces on the inner and outer tubes. In order to have an estimate of the quantitative effects of various emerging physical parameters on flow characteristics which are involved in the solutions of the considered analysis, we have used the MATLAB software for plotting the contour graphs and discussed in details. We have observed that the intensity of heat sink increases when amplitude ratio is increased. The trapped bolus appears when Darcy number is small but the important observation is that for the large values of Darcy number, the trapped bolus disappears and fluid moves like a block, which shows some sort of rigidity.
{"title":"Effect of Temperature Dependent Sink on Peristaltic Transport in a Differentially Heated Vertical Annulus Filled with a Porous Material","authors":"A. Tiwari, Akshay K. Singh","doi":"10.1080/23324309.2022.2145311","DOIUrl":"https://doi.org/10.1080/23324309.2022.2145311","url":null,"abstract":"Abstract In this article, we have presented a mathematical analysis to study the peristaltic flow through vertical annulus filled with a porous material bounded by two concentric uniform tubes. This analysis can serve as a model which may help in understanding the mechanism of physiological flows and heat transfer in a vertical annulus subject to differentially heating in the presence of a temperature dependent sink. The inner tube is uniform and rigid, while the outer tube has a sinusoidal wave traveling down its wall. It is analyzed in a wave frame of reference moving with velocity of the wave under a zero Reynolds number and long wavelength approximation in the presence of a temperature dependent sink. The analytical solutions are obtained for temperature, axial velocity, stream function and axial pressure-gradient. We also present numerical integration in order to analyze the pressure rise and frictional forces on the inner and outer tubes. In order to have an estimate of the quantitative effects of various emerging physical parameters on flow characteristics which are involved in the solutions of the considered analysis, we have used the MATLAB software for plotting the contour graphs and discussed in details. We have observed that the intensity of heat sink increases when amplitude ratio is increased. The trapped bolus appears when Darcy number is small but the important observation is that for the large values of Darcy number, the trapped bolus disappears and fluid moves like a block, which shows some sort of rigidity.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"372 - 392"},"PeriodicalIF":0.7,"publicationDate":"2022-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46099739","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}
Pub Date : 2022-08-25DOI: 10.1080/23324309.2023.2166077
T. Burlacu, D. Lathouwers, Z. Perkó
Abstract In this paper we propose a solution to the need for a fast particle transport algorithm in Online Adaptive Proton Therapy capable of cheaply, but accurately computing the changes in patient dose metrics as a result of changes in the system parameters. We obtain the proton phase-space density through the product of the numerical solution to the one-dimensional Fokker-Planck equation and the analytical solution to the Fermi-Eyges equation. Moreover, a corresponding adjoint system was derived and solved for the adjoint flux. The proton phase-space density together with the adjoint flux and the metric (chosen as the energy deposited by the beam in a variable region of interest) allowed assessing the accuracy of our algorithm to different perturbation ranges in the system parameters and regions of interest. The algorithm achieved negligible errors ( ) for small Hounsfield unit (HU) perturbation ranges (–40 HU to 40 HU) and small to moderate errors (3% to 17%) – in line with the well-known limitation of adjoint approaches – for large perturbation ranges (–400 HU to 400 HU) in the case of most clinical interest where the region of interest surrounds the Bragg peak. Given these results coupled with the capability of further improving the timing performance it can be concluded that our algorithm presents a viable solution for the specific purpose of Online Adaptive Proton Therapy.
{"title":"A Deterministic Adjoint-Based Semi-Analytical Algorithm for Fast Response Change Computations in Proton Therapy","authors":"T. Burlacu, D. Lathouwers, Z. Perkó","doi":"10.1080/23324309.2023.2166077","DOIUrl":"https://doi.org/10.1080/23324309.2023.2166077","url":null,"abstract":"Abstract In this paper we propose a solution to the need for a fast particle transport algorithm in Online Adaptive Proton Therapy capable of cheaply, but accurately computing the changes in patient dose metrics as a result of changes in the system parameters. We obtain the proton phase-space density through the product of the numerical solution to the one-dimensional Fokker-Planck equation and the analytical solution to the Fermi-Eyges equation. Moreover, a corresponding adjoint system was derived and solved for the adjoint flux. The proton phase-space density together with the adjoint flux and the metric (chosen as the energy deposited by the beam in a variable region of interest) allowed assessing the accuracy of our algorithm to different perturbation ranges in the system parameters and regions of interest. The algorithm achieved negligible errors ( ) for small Hounsfield unit (HU) perturbation ranges (–40 HU to 40 HU) and small to moderate errors (3% to 17%) – in line with the well-known limitation of adjoint approaches – for large perturbation ranges (–400 HU to 400 HU) in the case of most clinical interest where the region of interest surrounds the Bragg peak. Given these results coupled with the capability of further improving the timing performance it can be concluded that our algorithm presents a viable solution for the specific purpose of Online Adaptive Proton Therapy.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"52 1","pages":"1 - 41"},"PeriodicalIF":0.7,"publicationDate":"2022-08-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49620860","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}
Pub Date : 2022-07-29DOI: 10.1080/23324309.2022.2110898
R. D. Garcia
Abstract Analytical discrete-ordinates (ADO) solutions are developed for two improved one-dimensional (1D) models of particle transport in ducts that include wall migration. One of the studied models is based on an approximation of the transverse and azimuthal dependencies of the angular flux in terms of two basis functions, while the other uses three. Particle migration in the wall is modeled by an exponential displacement kernel. Numerical results for duct reflection and transmission probabilities are reported for various test cases, including some cases of thermal-neutron transport through iron, concrete, and graphite ducts that have been defined and studied by other authors. A comparison is performed with a set of published results comprising realistic results obtained with the Monte Carlo code MCNP and results from a numerical implementation of the discrete-ordinates method for the model based on three basis functions. To resolve some issues raised during the comparison process, a numerical discrete-ordinates solution of the problem has also been implemented in the course of this work.
{"title":"Analytical Discrete-Ordinates Solutions for Improved 1D Models of Particle Transport in Ducts with Wall Migration","authors":"R. D. Garcia","doi":"10.1080/23324309.2022.2110898","DOIUrl":"https://doi.org/10.1080/23324309.2022.2110898","url":null,"abstract":"Abstract Analytical discrete-ordinates (ADO) solutions are developed for two improved one-dimensional (1D) models of particle transport in ducts that include wall migration. One of the studied models is based on an approximation of the transverse and azimuthal dependencies of the angular flux in terms of two basis functions, while the other uses three. Particle migration in the wall is modeled by an exponential displacement kernel. Numerical results for duct reflection and transmission probabilities are reported for various test cases, including some cases of thermal-neutron transport through iron, concrete, and graphite ducts that have been defined and studied by other authors. A comparison is performed with a set of published results comprising realistic results obtained with the Monte Carlo code MCNP and results from a numerical implementation of the discrete-ordinates method for the model based on three basis functions. To resolve some issues raised during the comparison process, a numerical discrete-ordinates solution of the problem has also been implemented in the course of this work.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"265 - 304"},"PeriodicalIF":0.7,"publicationDate":"2022-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43243419","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}
Pub Date : 2022-07-29DOI: 10.1080/23324309.2022.2110897
P. Picca, R. Furfaro, S. Dulla, P. Ravetto
Abstract The article presents the foundation of a novel methodology developed for the solution of the neutron transport equation, named the transport driven-diffusion approach, which can be considered as an evolution of the classic multiple collision method. The idea behind this method is based on the expansion of the full solution in terms of the contributions of the particles emitted by successive collisions plus a residual term, accounting for particles which have undergone more than a predefined number of collisions. In order to determine the contribution at each collision order, a transport equation with a source term is solved, while the estimation of the residue is based on a diffusion theory model. The physical rationale for the choice of the diffusion model for the residue is discussed and justified, as physics suggests that the diffusion assumptions become more applicable for the description of the particles having suffered a certain number of collisions rather than to the original transport problem. Some results are presented for a set of steady-state and time-dependent test cases. Their analysis shows the remarkable advantage of the method proposed in terms of accuracy and computational time, when compared to standard diffusion and multiple collision at the same order.
{"title":"“On the Foundation of Transport-Driven Diffusion for Neutron Transport Problems”","authors":"P. Picca, R. Furfaro, S. Dulla, P. Ravetto","doi":"10.1080/23324309.2022.2110897","DOIUrl":"https://doi.org/10.1080/23324309.2022.2110897","url":null,"abstract":"Abstract The article presents the foundation of a novel methodology developed for the solution of the neutron transport equation, named the transport driven-diffusion approach, which can be considered as an evolution of the classic multiple collision method. The idea behind this method is based on the expansion of the full solution in terms of the contributions of the particles emitted by successive collisions plus a residual term, accounting for particles which have undergone more than a predefined number of collisions. In order to determine the contribution at each collision order, a transport equation with a source term is solved, while the estimation of the residue is based on a diffusion theory model. The physical rationale for the choice of the diffusion model for the residue is discussed and justified, as physics suggests that the diffusion assumptions become more applicable for the description of the particles having suffered a certain number of collisions rather than to the original transport problem. Some results are presented for a set of steady-state and time-dependent test cases. Their analysis shows the remarkable advantage of the method proposed in terms of accuracy and computational time, when compared to standard diffusion and multiple collision at the same order.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"305 - 328"},"PeriodicalIF":0.7,"publicationDate":"2022-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46840519","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}
Pub Date : 2022-07-29DOI: 10.1080/23324309.2022.2111443
B. Ganapol
Abstract A new, highly precise benchmark for the monoenergetic 1 D neutron transport equation with isotropic scattering, based on Case’s singular eigenfunctions, is presented. Because of the nature of singular distributions, Case’s analytical solution notoriously resists straightforward numerical computation. To overcome this difficulty, two complementary Lagrange interpolation schemes are constructed to achieve extreme precision (∼8 to 10 places) for number of secondaries in the range 0.001 < c < 0.99999 and slab thicknesses in the range 1 < Δ <100.
摘要基于Case奇异特征函数，提出了具有各向同性散射的单能一维中子输运方程的高精度基准。由于奇异分布的性质，凯斯的解析解拒绝直接的数值计算。为了克服这一困难，构建了两种互补的拉格朗日插值方案，以在0.001 < c < 0.99999范围内的次级数和1 < Δ <100范围内的板料厚度实现极高的精度(~ 8到10位)。
{"title":"The Ln/Ln0 Method for 1D Neutron Transport in a Slab Medium","authors":"B. Ganapol","doi":"10.1080/23324309.2022.2111443","DOIUrl":"https://doi.org/10.1080/23324309.2022.2111443","url":null,"abstract":"Abstract A new, highly precise benchmark for the monoenergetic 1 D neutron transport equation with isotropic scattering, based on Case’s singular eigenfunctions, is presented. Because of the nature of singular distributions, Case’s analytical solution notoriously resists straightforward numerical computation. To overcome this difficulty, two complementary Lagrange interpolation schemes are constructed to achieve extreme precision (∼8 to 10 places) for number of secondaries in the range 0.001 < c < 0.99999 and slab thicknesses in the range 1 < Δ <100.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"239 - 264"},"PeriodicalIF":0.7,"publicationDate":"2022-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43679712","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}
Abstract Understanding, predicting and controlling the motion of the solid particles in a confined cavity is significant. The motion of a neutrally buoyant circular particle in a lid-driven square cavity is studied with the lattice Boltzmann method, where the effects of the initial position, particle size and Reynolds number are investigated. The obvious characteristic of the motion of the circular particle is the existence of the limit cycle, which is the competitive result of the inertia, wall-repulsion force and vortex behavior. The limit cycle is insensitive to the initial position of the circular particle, namely, no matter where the circular particle is placed initially, the limit cycle is the same. With the increase of the particle size, the wall-repulsion force becomes stronger, which is dominant over the centrifugal force, and the limit cycle shrinks toward the center of the square cavity. With the increase of the Reynolds number, a new secondary vortex develops at the top left corner of the square cavity, and the limit cycle is pushed toward the bottom right corner of the square cavity.
{"title":"Motion of a Neutrally Buoyant Circular Particle in a Lid-Driven Square Cavity: A Numerical Study","authors":"Junjie Hu, Dongke Sun, Shaohua Mao, Hongmei Wu, Songyang Yu, Maosen Xu","doi":"10.1080/23324309.2022.2100424","DOIUrl":"https://doi.org/10.1080/23324309.2022.2100424","url":null,"abstract":"Abstract Understanding, predicting and controlling the motion of the solid particles in a confined cavity is significant. The motion of a neutrally buoyant circular particle in a lid-driven square cavity is studied with the lattice Boltzmann method, where the effects of the initial position, particle size and Reynolds number are investigated. The obvious characteristic of the motion of the circular particle is the existence of the limit cycle, which is the competitive result of the inertia, wall-repulsion force and vortex behavior. The limit cycle is insensitive to the initial position of the circular particle, namely, no matter where the circular particle is placed initially, the limit cycle is the same. With the increase of the particle size, the wall-repulsion force becomes stronger, which is dominant over the centrifugal force, and the limit cycle shrinks toward the center of the square cavity. With the increase of the Reynolds number, a new secondary vortex develops at the top left corner of the square cavity, and the limit cycle is pushed toward the bottom right corner of the square cavity.","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"51 1","pages":"222 - 237"},"PeriodicalIF":0.7,"publicationDate":"2022-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44219679","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}
Pub Date : 2022-06-20DOI: 10.1080/23324309.2022.2086263
V. D. Camiola, G. Nastasi
Abstract
A hydrodynamical model for charge transport in narrow strips of graphene is here presented. The model takes into account the interactions with the well-known lattice vibrations and with the edge of the strip. The remarkable result is the modulation of the charge current due to the swapping of charge carriers between the conduction and the valence bands, controlled by the Fermi energy variation and by the thickness of the ribbon. The numerical test shows a behavior comparable with that one obtained by solving directly the Boltzmann equation but with a considerable reduction of the computational time.
{"title":"Bipolar Hydrodynamical Model for Charge Transport in Graphene Nanoribbons","authors":"V. D. Camiola, G. Nastasi","doi":"10.1080/23324309.2022.2086263","DOIUrl":"https://doi.org/10.1080/23324309.2022.2086263","url":null,"abstract":"<p><b>Abstract</b></p><p>A hydrodynamical model for charge transport in narrow strips of graphene is here presented. The model takes into account the interactions with the well-known lattice vibrations and with the edge of the strip. The remarkable result is the modulation of the charge current due to the swapping of charge carriers between the conduction and the valence bands, controlled by the Fermi energy variation and by the thickness of the ribbon. The numerical test shows a behavior comparable with that one obtained by solving directly the Boltzmann equation but with a considerable reduction of the computational time.</p>","PeriodicalId":54305,"journal":{"name":"Journal of Computational and Theoretical Transport","volume":"19 9","pages":""},"PeriodicalIF":0.7,"publicationDate":"2022-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138526233","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}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：