Anomalous transport in magnetically confined plasmas is investigated using temporal fractional transport equations.The use of temporal fractional transport equations means that the order of the partial derivative with...Anomalous transport in magnetically confined plasmas is investigated using temporal fractional transport equations.The use of temporal fractional transport equations means that the order of the partial derivative with respect to time is a fraction. In this case, the Caputo fractional derivative relative to time is utilized, because it preserves the form of the initial conditions. A numerical calculation reveals that the fractional order of the temporal derivative α(α ∈(0, 1), sub-diffusive regime) controls the diffusion rate. The temporal fractional derivative is related to the fact that the evolution of a physical quantity is affected by its past history, depending on what are termed memory effects. The magnitude of α is a measure of such memory effects. When α decreases, so does the rate of particle diffusion due to memory effects. As a result,if a system initially has a density profile without a source, then the smaller the α is, the more slowly the density profile approaches zero. When a source is added, due to the balance of the diffusion and fueling processes, the system reaches a steady state and the density profile does not evolve. As α decreases, the time required for the system to reach a steady state increases. In magnetically confined plasmas, the temporal fractional transport model can be applied to off-axis heating processes. Moreover, it is found that the memory effects reduce the rate of energy conduction and hollow temperature profiles can be sustained for a longer time in sub-diffusion processes than in ordinary diffusion processes.展开更多
The existence of spiral solution for the two-dimensional transport equations is considered in the present paper. Based on the notion of generalized solutions in the sense of Lebesgue-stieltjes integral, the global wea...The existence of spiral solution for the two-dimensional transport equations is considered in the present paper. Based on the notion of generalized solutions in the sense of Lebesgue-stieltjes integral, the global weak solution of transport equations which includes δ-shocks and vacuum is constructed for some special initial data.展开更多
This paper deal with a nonlinear transport equation with delayed neutron andgeneral boundary conditions. We establish, via the nonlinear semigroups approach, the exis-tence and uniqueness of the mild solution, weak so...This paper deal with a nonlinear transport equation with delayed neutron andgeneral boundary conditions. We establish, via the nonlinear semigroups approach, the exis-tence and uniqueness of the mild solution, weak solution, strong solution and local solutionon LP-spaces (1 ≤ p 〈 +∞). Local and non local evolution problems are discussed.展开更多
In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the g...In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.展开更多
We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wav...We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wavenumber and is the mean free path of the electron). The Boltzmann transport equation and other relevant equations were solved analytically to obtain an expression for the AME current density, consisting of longitudinal and Hall components. Our numerical results indicate that both components of the AME current densities display oscillatory behaviour. Furthermore, geometric resonances and Weiss oscillations were each defined using the relationship between the current density and Surface Acoustic Wave (SAW) frequency and the inverse of the applied magnetic field, respectively. Our results show that the AME current density of bandgap graphene, which can be controlled to suit a particular electronic device application, is smaller than that of (gapless) graphene and is therefore, more suited for nanophotonic device applications.展开更多
This paper concerns the weak solutions of some Monge-Amp^re type equa- tions in the optimal transportation theory. The relationship between the Aleksandrov solutions and the viscosity solutions of the Monge-Ampere typ...This paper concerns the weak solutions of some Monge-Amp^re type equa- tions in the optimal transportation theory. The relationship between the Aleksandrov solutions and the viscosity solutions of the Monge-Ampere type equations is discussed. A uniform estimate for solution of the Dirichlet problem with homogeneous boundary value is obtained.展开更多
This paper presents two uniformly convergent numerical schemes for the two dimensional steady state discrete ordinates transport equation in the diffusive regime,which is valid up to the boundary and interface layers....This paper presents two uniformly convergent numerical schemes for the two dimensional steady state discrete ordinates transport equation in the diffusive regime,which is valid up to the boundary and interface layers.A five-point nodecentered and a four-point cell-centered tailored finite point schemes(TFPS)are introduced.The schemes first approximate the scattering coefficients and sources by piecewise constant functions and then use special solutions to the constant coefficient equation as local basis functions to formulate a discrete linear system.Numerically,both methods can not only capture the diffusion limit,but also exhibit uniform convergence in the diffusive regime,even with boundary layers.Numerical results show that the five-point scheme has first-order accuracy and the four-point scheme has second-order accuracy,uniformly with respect to the mean free path.Therefore a relatively coarse grid can be used to capture the two dimensional boundary and interface layers.展开更多
A B-spline Interpolation Transport Solver(BITS) based on a collocation method is developed. It solves transport equations as a generalized interpolation problem, taking the first-order accuracy in time and the second-...A B-spline Interpolation Transport Solver(BITS) based on a collocation method is developed. It solves transport equations as a generalized interpolation problem, taking the first-order accuracy in time and the second-order accuracy in space along with a predictor–corrector or under-relaxation iteration method. Numerical tests show that BITS can solve one-dimensional transport equations for tokamak plasma more accurately without additional computation cost, compared to the finite difference method transport solver which is widely used in existing tokamak transport codes.展开更多
We present an algorithm for numerical solution of transport equation in diffusive regimes, in which the transport equation is nearly singular and its solution becomes a solution of a diffusion equation. This algorithm...We present an algorithm for numerical solution of transport equation in diffusive regimes, in which the transport equation is nearly singular and its solution becomes a solution of a diffusion equation. This algorithm, which is based on the Least-squares FEM in combination with a scaling transformation, presents a good approximation of a diffusion operator in diffusive regimes and guarantees an accurate discrete solution. The numerical experiments in 2D and 3D case are given, and the numerical results show that this algorithm is correct and efficient.展开更多
In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. Th...In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. The convergence and error estimations of this scheme are proved. Its applications in the field of neutron logging show its effectiveness.展开更多
This paper analyzes a class of first order partial differential equations with delay (a model for the blood production system). The asymptotic behavior of solutions are studied.
The van der Waals(vdW)heterostructures of bilayer transition metal dichalcogenide obtained by vertically stacking have drawn increasing attention for their enormous potential applications in semiconductors and insulat...The van der Waals(vdW)heterostructures of bilayer transition metal dichalcogenide obtained by vertically stacking have drawn increasing attention for their enormous potential applications in semiconductors and insulators.Here,by using the first-principles calculations and the phonon Boltzmann transport equation(BTE),we studied the phonon transport properties of WS2/WSe2 bilayer heterostructures(WS2/WSe2-BHs).The lattice thermal conductivity of the ideal WS2/WSe2-BHs crystals at room temperature(RT)was 62.98 W/mK,which was clearly lower than the average lattice thermal conductivity of WS2 and WSe2 single layers.Another interesting finding is that the optical branches below 4.73 THz and acoustic branches have powerful coupling,mainly dominating the lattice thermal conductivity.Further,we also noticed that the phonon mean free path(MFP)of the WS2/WSe2-BHs(233 nm)was remarkably attenuated by the free-standing monolayer WS2(526 nm)and WSe2(1720 nm),leading to a small significant size effect of the WS2/WSe2-BHs.Our results systematically demonstrate the low optical and acoustic phonon modes-dominated phonon thermal transport in heterostructures and give a few important guidelines for the synthesis of van der Waals heterostructures with excellent phonon transport properties.展开更多
The well-posedness for the time-dependent neutron transport equationwith integral boundary conditions is established in L^1 space. Some spectral propertiesof the transport operator are discussed, the dominant eigenval...The well-posedness for the time-dependent neutron transport equationwith integral boundary conditions is established in L^1 space. Some spectral propertiesof the transport operator are discussed, the dominant eigenvalue is proved existing,and furthermore, the conservative law of migrating particle numbers is established.展开更多
It is well known that the boundary element method (BEM) is capable of converting a boundary- value equation into its discrete analog by a judicious application of the Green’s identity and complementary equation. Howe...It is well known that the boundary element method (BEM) is capable of converting a boundary- value equation into its discrete analog by a judicious application of the Green’s identity and complementary equation. However, for many challenging problems, the fundamental solution is either not available in a cheaply computable form or does not exist at all. Even when the fundamental solution does exist, it appears in a form that is highly non-local which inadvertently leads to a sys-tem of equations with a fully populated matrix. In this paper, fundamental solution of an auxiliary form of a governing partial differential equation coupled with the Green identity is used to discretize and localize an integro-partial differential transport equation by conversion into a boundary-domain form amenable to a hybrid boundary integral numerical formulation. It is observed that the numerical technique applied herein is able to accurately represent numerical and closed form solutions available in literature.展开更多
We show that by integrating out the electric field and incorporating proper boundary conditions,a Boltzmann equation can describe electron transport properties,continuously from the diffusive to ballistic regimes.Gene...We show that by integrating out the electric field and incorporating proper boundary conditions,a Boltzmann equation can describe electron transport properties,continuously from the diffusive to ballistic regimes.General analytical formulas of the conductance in D = 1,2,3 dimensions are obtained,which recover the Boltzmann–Drude formula and Landauer–B ¨uttiker formula in the diffusive and ballistic limits,respectively.This intuitive and efficient approach can be applied to investigate the interplay of system size and impurity scattering in various charge and spin transport phenomena,when the quantum interference effect is not important.展开更多
This paper deals with the solution of a parametric equation with generalized boundary condiiton in transport theory. It gives the distribution of parameter (so called delta-eigenvalue [1]) with which the equation has ...This paper deals with the solution of a parametric equation with generalized boundary condiiton in transport theory. It gives the distribution of parameter (so called delta-eigenvalue [1]) with which the equation has non-zero solution. A necessary and sufficient condition for the existence of; he control critical eigenvalue delta0 is established.展开更多
To evaluate the pollutant dispersion in background turbulent flows, most researches focus on statistical variation of concentrations or its fluctuations. However, those time-averaged quantities may be insufficient for...To evaluate the pollutant dispersion in background turbulent flows, most researches focus on statistical variation of concentrations or its fluctuations. However, those time-averaged quantities may be insufficient for risk assessment, because there emerge many high-intensity pollutant areas in the instantaneous concentration field. In this study, we tried to estimate the frequency of appearance of the high concentration areas in a turbulent flow based on the Probability Density Function (PDF) of concentration. The high concentration area was recognized by two conditions based on the concentration and the concentration gradient values. We considered that the estimation equation for the frequency of appearance of the recognized areas consisted of two terms based on each condition. In order to represent the two terms with physical quantities of velocity and concentration fields, simultaneous PIV (Particle Image Velocimetry) and PLIF (Planar Laser-Induced Fluorescence) measurement and PLIF time-serial measurement were performed in a quasi-homogeneous turbulent flow. According to the experimental results, one of the terms, related to the condition of the concentration, was found to be represented by the concentration PDF, while the other term, by the streamwise mean velocity and the integral length scale of the turbulent flow. Based on the results, we developed an estimation equation including the concentration PDF and the flow features of mean velocity and integral scale of turbulence. In the area where the concentration PDF was a Gaussian one, the difference between the frequencies of appearance estimated by the equation and calculated from the experimental data was within 25%, which showed good accuracy of our proposed estimation equation. Therefore, our proposed equation is feasible for estimating the frequency of appearance of high concentration areas in a limited area in turbulent mass diffusion.展开更多
We propose an optimized scheme to determine the smearing parameter in the Gaussian function that is used to replace the Dirac δ function in the first Brillouin zone sampling. The broadening width is derived by analyz...We propose an optimized scheme to determine the smearing parameter in the Gaussian function that is used to replace the Dirac δ function in the first Brillouin zone sampling. The broadening width is derived by analyzing the difference of the results from the phase-space method and Gaussian broadening method. As a demonstration, using the present approach,we investigate the phonon transport in a typical layered material, graphite. Our scheme is benchmarked by comparing with other zone sampling methods. Both the three-phonon phonon scattering rates and thermal conductivity are consistent with the prediction from the widely used tetrahedron method and adaptive broadening method. The computational efficiency of our scheme is more than one order of magnitude higher than the two other methods. Furthermore, the effect of fourphonon scattering in phonon transport in graphite is also investigated. It is found that four-phonon scattering reduces the through-plane thermal conductivity by 10%. Our methods could be a reference for the prediction of thermal conductivity of anisotropic material in the future.展开更多
In this work, we simulate numerically the thermal effects in nitrogen at atmospheric pressure caused by a negative corona DC discharge of low current. A mathematical function that simulates the injection of the therma...In this work, we simulate numerically the thermal effects in nitrogen at atmospheric pressure caused by a negative corona DC discharge of low current. A mathematical function that simulates the injection of the thermal transfer in gas is proposed. The simulated discharge is of a negative point to plane mass type, with an interelectrode distance of 12 mm and a symmetry about the axis of discharge. The spatial and temporal evolution of neutrals is analyzed based upon the equations of continuity, momentum and energy in a cylindrical geometry. For the geometry of the system, the FCT (flux corrected transport) procedure was adopted.展开更多
An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme (TSPAS), was generalized and implemented on a spherical icosahedral hexagonal grid (also referred to as a geodesic grid...An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme (TSPAS), was generalized and implemented on a spherical icosahedral hexagonal grid (also referred to as a geodesic grid) to solve the transport equation. The C grid discretization was used for the spatial discretization. To implement TSPAS on an unstructured grid, the original finite-difference scheme was further generalized. The two-step integration utilizes a combination of two separate schemes (a low-order monotone scheme and a high-order scheme that typically cannot ensure monotonicity) to calculate the fluxes at the cell walls (one scheme corresponds to one cell wall). The choice between these two schemes for each edge depends on a pre-updated scalar value using slightly increased fluxes. After the determination of an appropriate scheme, the final integration at a target cell is achieved by summing the fluxes that are computed by the different schemes. The conservative and shape-preserving properties of the generalized scheme are demonstrated. Numerical experiments are conducted at several horizontal resolutions. TSPAS is compared with the Flux Corrected Transport (FCT) approach to demonstrate the differences between the two methods, and several transport tests are performed to examine the accuracy, efficiency and robustness of the two schemes.展开更多
基金supported by the National Key R&D Program of China (Grant No. 2022YFE03090000)the National Natural Science Foundation of China (Grant No. 11925501)the Fundamental Research Fund for the Central Universities (Grant No. DUT22ZD215)。
文摘Anomalous transport in magnetically confined plasmas is investigated using temporal fractional transport equations.The use of temporal fractional transport equations means that the order of the partial derivative with respect to time is a fraction. In this case, the Caputo fractional derivative relative to time is utilized, because it preserves the form of the initial conditions. A numerical calculation reveals that the fractional order of the temporal derivative α(α ∈(0, 1), sub-diffusive regime) controls the diffusion rate. The temporal fractional derivative is related to the fact that the evolution of a physical quantity is affected by its past history, depending on what are termed memory effects. The magnitude of α is a measure of such memory effects. When α decreases, so does the rate of particle diffusion due to memory effects. As a result,if a system initially has a density profile without a source, then the smaller the α is, the more slowly the density profile approaches zero. When a source is added, due to the balance of the diffusion and fueling processes, the system reaches a steady state and the density profile does not evolve. As α decreases, the time required for the system to reach a steady state increases. In magnetically confined plasmas, the temporal fractional transport model can be applied to off-axis heating processes. Moreover, it is found that the memory effects reduce the rate of energy conduction and hollow temperature profiles can be sustained for a longer time in sub-diffusion processes than in ordinary diffusion processes.
基金supported by National Natural Science Foundation of China (10871199)one hundred talent project from the Chinese Academy of Sciences
文摘The existence of spiral solution for the two-dimensional transport equations is considered in the present paper. Based on the notion of generalized solutions in the sense of Lebesgue-stieltjes integral, the global weak solution of transport equations which includes δ-shocks and vacuum is constructed for some special initial data.
文摘This paper deal with a nonlinear transport equation with delayed neutron andgeneral boundary conditions. We establish, via the nonlinear semigroups approach, the exis-tence and uniqueness of the mild solution, weak solution, strong solution and local solutionon LP-spaces (1 ≤ p 〈 +∞). Local and non local evolution problems are discussed.
基金supported by the Simons Foundation:Collaboration Grantssupported by the AFOSR grant FA9550-18-1-0383.
文摘In this paper,we consider the high order method for solving the linear transport equations under diffusive scaling and with random inputs.To tackle the randomness in the problem,the stochastic Galerkin method of the generalized polynomial chaos approach has been employed.Besides,the high order implicit-explicit scheme under the micro-macro decomposition framework and the discontinuous Galerkin method have been employed.We provide several numerical experiments to validate the accuracy and the stochastic asymptotic-preserving property.
文摘We study the acoustomagnetoelectric (AME) effect in two-dimensional graphene with an energy bandgap using the semiclassical Boltzmann transport equation within the hypersound regime, (where represents the acoustic wavenumber and is the mean free path of the electron). The Boltzmann transport equation and other relevant equations were solved analytically to obtain an expression for the AME current density, consisting of longitudinal and Hall components. Our numerical results indicate that both components of the AME current densities display oscillatory behaviour. Furthermore, geometric resonances and Weiss oscillations were each defined using the relationship between the current density and Surface Acoustic Wave (SAW) frequency and the inverse of the applied magnetic field, respectively. Our results show that the AME current density of bandgap graphene, which can be controlled to suit a particular electronic device application, is smaller than that of (gapless) graphene and is therefore, more suited for nanophotonic device applications.
基金supported by National Natural Science Foundation of China(11071119)
文摘This paper concerns the weak solutions of some Monge-Amp^re type equa- tions in the optimal transportation theory. The relationship between the Aleksandrov solutions and the viscosity solutions of the Monge-Ampere type equations is discussed. A uniform estimate for solution of the Dirichlet problem with homogeneous boundary value is obtained.
基金supported by the NSFC Project No.10971116.M.Tang is supported by Natural Science Foundation of Shanghai under Grant No.12ZR1445400Shanghai Pujiang Program 13PJ1404700+1 种基金supported in part by the National Natural Science Foundation of China under Grant DMS-11101278the Young Thousand Talents Program of China.
文摘This paper presents two uniformly convergent numerical schemes for the two dimensional steady state discrete ordinates transport equation in the diffusive regime,which is valid up to the boundary and interface layers.A five-point nodecentered and a four-point cell-centered tailored finite point schemes(TFPS)are introduced.The schemes first approximate the scattering coefficients and sources by piecewise constant functions and then use special solutions to the constant coefficient equation as local basis functions to formulate a discrete linear system.Numerically,both methods can not only capture the diffusion limit,but also exhibit uniform convergence in the diffusive regime,even with boundary layers.Numerical results show that the five-point scheme has first-order accuracy and the four-point scheme has second-order accuracy,uniformly with respect to the mean free path.Therefore a relatively coarse grid can be used to capture the two dimensional boundary and interface layers.
基金the National MCF Energy R&D Program of China(No.2019YFE03040004)the Comprehensive Research Facility for Fusion Technology Program of China(No.2018-000052-73-01-001228)the National MCF Energy R&D Program of China(No.2019YFE03060000)。
文摘A B-spline Interpolation Transport Solver(BITS) based on a collocation method is developed. It solves transport equations as a generalized interpolation problem, taking the first-order accuracy in time and the second-order accuracy in space along with a predictor–corrector or under-relaxation iteration method. Numerical tests show that BITS can solve one-dimensional transport equations for tokamak plasma more accurately without additional computation cost, compared to the finite difference method transport solver which is widely used in existing tokamak transport codes.
基金This work was supported by National Natural Science Foundation of China(No.10371096)
文摘We present an algorithm for numerical solution of transport equation in diffusive regimes, in which the transport equation is nearly singular and its solution becomes a solution of a diffusion equation. This algorithm, which is based on the Least-squares FEM in combination with a scaling transformation, presents a good approximation of a diffusion operator in diffusive regimes and guarantees an accurate discrete solution. The numerical experiments in 2D and 3D case are given, and the numerical results show that this algorithm is correct and efficient.
文摘In this paper, a new numerical method, the coupling method of spherical harmonic function spectral and streamline diffusion finite element for unsteady Boltzmann equation in the neutron logging field, is discussed. The convergence and error estimations of this scheme are proved. Its applications in the field of neutron logging show its effectiveness.
文摘This paper analyzes a class of first order partial differential equations with delay (a model for the blood production system). The asymptotic behavior of solutions are studied.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51720105007,51806031,11602149,and GZ1257)the Fundamental Research Funds for the Central Universities,China(Grant Nos.DUT16RC(3)116 and DUT19RC(3)006)The computing resources from Supercomputer Center of Dalian University of Technology and ScGrid are greatly acknowledged。
文摘The van der Waals(vdW)heterostructures of bilayer transition metal dichalcogenide obtained by vertically stacking have drawn increasing attention for their enormous potential applications in semiconductors and insulators.Here,by using the first-principles calculations and the phonon Boltzmann transport equation(BTE),we studied the phonon transport properties of WS2/WSe2 bilayer heterostructures(WS2/WSe2-BHs).The lattice thermal conductivity of the ideal WS2/WSe2-BHs crystals at room temperature(RT)was 62.98 W/mK,which was clearly lower than the average lattice thermal conductivity of WS2 and WSe2 single layers.Another interesting finding is that the optical branches below 4.73 THz and acoustic branches have powerful coupling,mainly dominating the lattice thermal conductivity.Further,we also noticed that the phonon mean free path(MFP)of the WS2/WSe2-BHs(233 nm)was remarkably attenuated by the free-standing monolayer WS2(526 nm)and WSe2(1720 nm),leading to a small significant size effect of the WS2/WSe2-BHs.Our results systematically demonstrate the low optical and acoustic phonon modes-dominated phonon thermal transport in heterostructures and give a few important guidelines for the synthesis of van der Waals heterostructures with excellent phonon transport properties.
文摘The well-posedness for the time-dependent neutron transport equationwith integral boundary conditions is established in L^1 space. Some spectral propertiesof the transport operator are discussed, the dominant eigenvalue is proved existing,and furthermore, the conservative law of migrating particle numbers is established.
文摘It is well known that the boundary element method (BEM) is capable of converting a boundary- value equation into its discrete analog by a judicious application of the Green’s identity and complementary equation. However, for many challenging problems, the fundamental solution is either not available in a cheaply computable form or does not exist at all. Even when the fundamental solution does exist, it appears in a form that is highly non-local which inadvertently leads to a sys-tem of equations with a fully populated matrix. In this paper, fundamental solution of an auxiliary form of a governing partial differential equation coupled with the Green identity is used to discretize and localize an integro-partial differential transport equation by conversion into a boundary-domain form amenable to a hybrid boundary integral numerical formulation. It is observed that the numerical technique applied herein is able to accurately represent numerical and closed form solutions available in literature.
基金Project supported by the National Basic Research Program of China(Grant Nos.2015CB921202 and 2014CB921103)the National Natural Science Foundation of China(Grant No.11225420)
文摘We show that by integrating out the electric field and incorporating proper boundary conditions,a Boltzmann equation can describe electron transport properties,continuously from the diffusive to ballistic regimes.General analytical formulas of the conductance in D = 1,2,3 dimensions are obtained,which recover the Boltzmann–Drude formula and Landauer–B ¨uttiker formula in the diffusive and ballistic limits,respectively.This intuitive and efficient approach can be applied to investigate the interplay of system size and impurity scattering in various charge and spin transport phenomena,when the quantum interference effect is not important.
基金Project supported by the National Natural Science Foundation of China.
文摘This paper deals with the solution of a parametric equation with generalized boundary condiiton in transport theory. It gives the distribution of parameter (so called delta-eigenvalue [1]) with which the equation has non-zero solution. A necessary and sufficient condition for the existence of; he control critical eigenvalue delta0 is established.
文摘To evaluate the pollutant dispersion in background turbulent flows, most researches focus on statistical variation of concentrations or its fluctuations. However, those time-averaged quantities may be insufficient for risk assessment, because there emerge many high-intensity pollutant areas in the instantaneous concentration field. In this study, we tried to estimate the frequency of appearance of the high concentration areas in a turbulent flow based on the Probability Density Function (PDF) of concentration. The high concentration area was recognized by two conditions based on the concentration and the concentration gradient values. We considered that the estimation equation for the frequency of appearance of the recognized areas consisted of two terms based on each condition. In order to represent the two terms with physical quantities of velocity and concentration fields, simultaneous PIV (Particle Image Velocimetry) and PLIF (Planar Laser-Induced Fluorescence) measurement and PLIF time-serial measurement were performed in a quasi-homogeneous turbulent flow. According to the experimental results, one of the terms, related to the condition of the concentration, was found to be represented by the concentration PDF, while the other term, by the streamwise mean velocity and the integral length scale of the turbulent flow. Based on the results, we developed an estimation equation including the concentration PDF and the flow features of mean velocity and integral scale of turbulence. In the area where the concentration PDF was a Gaussian one, the difference between the frequencies of appearance estimated by the equation and calculated from the experimental data was within 25%, which showed good accuracy of our proposed estimation equation. Therefore, our proposed equation is feasible for estimating the frequency of appearance of high concentration areas in a limited area in turbulent mass diffusion.
基金support from the National Natural Science Foundation of China (Grant No. 51706134)。
文摘We propose an optimized scheme to determine the smearing parameter in the Gaussian function that is used to replace the Dirac δ function in the first Brillouin zone sampling. The broadening width is derived by analyzing the difference of the results from the phase-space method and Gaussian broadening method. As a demonstration, using the present approach,we investigate the phonon transport in a typical layered material, graphite. Our scheme is benchmarked by comparing with other zone sampling methods. Both the three-phonon phonon scattering rates and thermal conductivity are consistent with the prediction from the widely used tetrahedron method and adaptive broadening method. The computational efficiency of our scheme is more than one order of magnitude higher than the two other methods. Furthermore, the effect of fourphonon scattering in phonon transport in graphite is also investigated. It is found that four-phonon scattering reduces the through-plane thermal conductivity by 10%. Our methods could be a reference for the prediction of thermal conductivity of anisotropic material in the future.
文摘In this work, we simulate numerically the thermal effects in nitrogen at atmospheric pressure caused by a negative corona DC discharge of low current. A mathematical function that simulates the injection of the thermal transfer in gas is proposed. The simulated discharge is of a negative point to plane mass type, with an interelectrode distance of 12 mm and a symmetry about the axis of discharge. The spatial and temporal evolution of neutrals is analyzed based upon the equations of continuity, momentum and energy in a cylindrical geometry. For the geometry of the system, the FCT (flux corrected transport) procedure was adopted.
基金supported by the National Natural Science Foundation of China(Grant No.41505066)the Basic Scientific Research and Operation Foundation of Chinese Academy Meteorological Sciences(Grant Nos.2015Z002,2015Y005)the National Research and Development Key Program:Global Change and Mitigation Strategies(No.2016YFA0602101)
文摘An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme (TSPAS), was generalized and implemented on a spherical icosahedral hexagonal grid (also referred to as a geodesic grid) to solve the transport equation. The C grid discretization was used for the spatial discretization. To implement TSPAS on an unstructured grid, the original finite-difference scheme was further generalized. The two-step integration utilizes a combination of two separate schemes (a low-order monotone scheme and a high-order scheme that typically cannot ensure monotonicity) to calculate the fluxes at the cell walls (one scheme corresponds to one cell wall). The choice between these two schemes for each edge depends on a pre-updated scalar value using slightly increased fluxes. After the determination of an appropriate scheme, the final integration at a target cell is achieved by summing the fluxes that are computed by the different schemes. The conservative and shape-preserving properties of the generalized scheme are demonstrated. Numerical experiments are conducted at several horizontal resolutions. TSPAS is compared with the Flux Corrected Transport (FCT) approach to demonstrate the differences between the two methods, and several transport tests are performed to examine the accuracy, efficiency and robustness of the two schemes.