NUMERICAL INTEGRAL OF RESISTANCE COEFFICIENTS IN DIFFUSION
Abstract
The resistance coefficients in the screened Coulomb potential of stellar plasma are evaluated to high accuracy. I have analyzed the possible singularities in the integral of scattering angle. There are possible singularities in the case of an attractive potential. This may result in a problem for the numerical integral. In order to avoid the problem, I have used a proper scheme, e.g., splitting into many subintervals where the width of each subinterval is determined by the variation of the integrand, to calculate the scattering angle. The collision integrals are calculated by using Romberg’s method, therefore the accuracy is high (i.e., ∼10{sup −12}). The results of collision integrals and their derivatives for −7 ≤ ψ ≤ 5 are listed. By using Hermite polynomial interpolation from those data, the collision integrals can be obtained with an accuracy of 10{sup −10}. For very weakly coupled plasma ( ψ ≥ 4.5), analytical fittings for collision integrals are available with an accuracy of 10{sup −11}. I have compared the final results of resistance coefficients with other works and found that, for a repulsive potential, the results are basically the same as others’; for an attractive potential, the results in cases of intermediate and strong couplingmore »
 Authors:
 Yunnan Observatories, Chinese Academy of Sciences, 396 Yangfangwang, Guandu District, Kunming, 650216 (China)
 Publication Date:
 OSTI Identifier:
 22663976
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Astrophysical Journal; Journal Volume: 834; Journal Issue: 2; Other Information: Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ACCURACY; COLLISION INTEGRALS; COLLISIONS; COMPARATIVE EVALUATIONS; COULOMB FIELD; DIFFUSION; HERMITE POLYNOMIALS; INTERPOLATION; PLASMA; SCATTERING; STAR MODELS; STARS; STRONGCOUPLING MODEL
Citation Formats
Zhang, Q. S., Email: zqs@ynao.ac.cn. NUMERICAL INTEGRAL OF RESISTANCE COEFFICIENTS IN DIFFUSION. United States: N. p., 2017.
Web. doi:10.3847/15384357/834/2/132.
Zhang, Q. S., Email: zqs@ynao.ac.cn. NUMERICAL INTEGRAL OF RESISTANCE COEFFICIENTS IN DIFFUSION. United States. doi:10.3847/15384357/834/2/132.
Zhang, Q. S., Email: zqs@ynao.ac.cn. Tue .
"NUMERICAL INTEGRAL OF RESISTANCE COEFFICIENTS IN DIFFUSION". United States.
doi:10.3847/15384357/834/2/132.
@article{osti_22663976,
title = {NUMERICAL INTEGRAL OF RESISTANCE COEFFICIENTS IN DIFFUSION},
author = {Zhang, Q. S., Email: zqs@ynao.ac.cn},
abstractNote = {The resistance coefficients in the screened Coulomb potential of stellar plasma are evaluated to high accuracy. I have analyzed the possible singularities in the integral of scattering angle. There are possible singularities in the case of an attractive potential. This may result in a problem for the numerical integral. In order to avoid the problem, I have used a proper scheme, e.g., splitting into many subintervals where the width of each subinterval is determined by the variation of the integrand, to calculate the scattering angle. The collision integrals are calculated by using Romberg’s method, therefore the accuracy is high (i.e., ∼10{sup −12}). The results of collision integrals and their derivatives for −7 ≤ ψ ≤ 5 are listed. By using Hermite polynomial interpolation from those data, the collision integrals can be obtained with an accuracy of 10{sup −10}. For very weakly coupled plasma ( ψ ≥ 4.5), analytical fittings for collision integrals are available with an accuracy of 10{sup −11}. I have compared the final results of resistance coefficients with other works and found that, for a repulsive potential, the results are basically the same as others’; for an attractive potential, the results in cases of intermediate and strong coupling show significant differences. The resulting resistance coefficients are tested in the solar model. Comparing with the widely used models of Cox et al. and Thoul et al., the resistance coefficients in the screened Coulomb potential lead to a slightly weaker effect in the solar model, which is contrary to the expectation of attempts to solve the solar abundance problem.},
doi = {10.3847/15384357/834/2/132},
journal = {Astrophysical Journal},
number = 2,
volume = 834,
place = {United States},
year = {Tue Jan 10 00:00:00 EST 2017},
month = {Tue Jan 10 00:00:00 EST 2017}
}

A pathintegral solution to truly nonlinear FokkerPlanck equations is derived. Such equations exhibit in the drift and diffusion coefficients a functional dependence on the distribution function. This type of implicit time dependence is shown to introduce terms into the propagator function of the exact same order in the time step tau, as does an explicit time dependence if the functional dependence is sufficiently smooth. A standard discrete lattice formulation of the path integral is then used to reproduce the appropriate, truly nonlinear FokkerPlanck equation. This discrete formulation provides a basis for an efficient numerical algorithm and is applied with excellentmore »

Numerical solution of stochastic differential equations with constant diffusion coefficients
We present RungeKutta methods of high accuracy for stochastic differential equations with constant diffusion coefficients. We analyze L/sub 2/ convergence of these methods and present convergence proofs. For scalar equations a secondorder method is derived, and for systems a method of order oneandonehalf is derived. We further consider a variance reduction technique based on Hermite expansions for evaluating expectations of functions of sample solutions. Numerical examples in two dimensions are presented. 
Mutual diffusion coefficients, electrical conductances, osmotic coefficients, and ionic transport coefficients l/sub ij/ for aqueous CuSO/sub 4/ at 25/sup 0/C
Experimental thermodynamic and transport data are presented for aqueous CuSO/sub 4/ solutions at 25/sup 0/C from low concentrations to near saturation. Included are diffusion coefficients (Rayleigh interferometry), electrical conductances (DC method), and osmotic coefficients (isopiestic). These data and corresponding literature data are critically compared, as are literature density data. The resulting best values, together with previously compared transference number results, provide an extensive set of critically reviewed data for aqueous CuSO/sub 4/. Onsager coefficients l/sub ij/ have been calculated for this salt and are compared with data for other valence types.