Search NASA⌕ Search

SEARCH · Search NASA

Results for “Debye Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simple phenomenological modeling of transition-region capacitance of forward-biased p-n junction diodes and transistor diodes

The derivation of a simple expression for the capacitance C(V) associated with the transition region of a p-n junction under a forward bias is derived by phenomenological reasoning. The treatment of C(V) is based on the conventional Shockley equations, and simpler expressions for C(V) result that are in general accord with the previous analytical and numerical results. C(V) consists of two components resulting from changes in majority carrier concentration and from free hole and electron accumulation in the space-charge region. The space-charge region is conceived as the intrinsic region of an n-i-p structure for a space-charge region markedly wider than the extrinsic Debye lengths at its edges. This region is excited in the sense that the forward bias creates hole and electron densities orders of magnitude larger than those in equilibrium. The recent Shirts-Gordon (1979) modeling of the space-charge region using a dielectric response function is contrasted with the more conventional Schottky-Shockley modeling.

Lindholm, F. A.↗

Microstructural analysis of undoped and moderately Sc-doped TiO2 anatase nanoparticles using Scherrer equation and Debye function analysis

Highlights: • UV light sensitive undoped and moderate Sc-doped anatase nanocrystallites were prepared without post-synthesis annealing. • Elongated nanocrystallites were found by Scherrer equation and Debye Function Analysis applied to XRD data and also HRTEM. • 4 at.% Sc doping reduced the amount of Ti vacancies by half. • Doping showed ~0.3% growth of lattice parameters. - Abstract: UV light sensitive undoped and moderate (4 at.%) Sc-doped anatase nanocrystallites were prepared by an efficient and environmentally benign method based on homogeneous hydrolysis of TiOSO{sub 4} and Sc{sub 2}(C{sub 2}O{sub 4}){sub 3.5}H{sub 2}O aqueous solutions using urea as a precipitation agent, without post-synthesis calcination or annealing. X-ray diffraction study of the obtained powders revealed single phase anatase in both samples. Two different methods were applied to estimate the anatase crystallite size and shape from the X-ray diffraction data. The first one is based on the more conventional Scherrer equation (as implanted in GSAS-II program) while the second one is based on Debye function analysis (as implanted in DUBUSSY suite version 2.2), which is considered to be more adequate for nanocrystallites. Although each program has its own shape models, both methods indicated similar results of elongated nanocrystallite geometry for Sc-doped TiO{sub 2}: (i) GSAS-II confirmed ellipsoid shape with equatorial size of 5.1(0) nm and axial size of 6.8(0), (ii) DEBUSSY suite 2.2 showed cylinder shape with diameter of equivalent circle in the ab-plane, D{sub ab} = 4.9(1.0) nm and crystallite length along the c-axis, L{sub c} = 5.4(3.2) nm. For the undoped TiO{sub 2} sample, equatorial size of 5.5(0) nm and axial size of 7.3(0) nm for the ellipsoid shape was determined according to GSAS-II, while D{sub ab} = 5.5(1.3) nm and L{sub c} = 6.0(2.2) nm for the cylinder shape was determined by DEBUSSY suite 2.2. These outcomes were found to corroborate with HRTEM analysis. Doping with 4 at.% Sc caused ~0.3% growth of lattice parameters from a = 3.7979 [4] Å, c = 9.4995 [9] Å for the undoped sample to a = 3.8110 [7] Å, c = 9.5274 [16] Å for the moderately Sc doped sample, while reduction in the Ti vacancies by half was confirmed by both methods. The introducing of Sc as electronically active secondary species into the crystal lattice of TiO{sub 2} can greatly alter its optical absorption. The incorporation of Sc into the crystal lattice of TiO{sub 2} favors the substitution of Ti, as both have approximately the same ionic size, and generates surface oxygen vacancies which can act as trapping center for the photogenerated electrons and reduce the recombination of electron–hole pairs. We can expect Sc doped TiO{sub 2} nanomaterials to have applications as effective photocatalysts for selective photocatalytic oxidation of various environmental organic pollutants.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Effects of laser radiation field on energies of hydrogen atom in plasmas

In this study, for the first time, the Schrödinger equation with more general exponential cosine screened Coulomb (MGECSC) potential is solved numerically in the presence of laser radiation field within the Ehlotzky approximation using the asymptotic iteration method. The MGECSC potential includes four different potential forms in consideration of different sets of the parameters in the potential. By applying laser field, the total interaction potential of hydrogen atom embedded in plasmas converts to double well-type potential. The plasma screening effects under the influence of laser field as well as confinement effects of laser field on hydrogen atom in Debye and quantum plasmas are investigated by solving the Schrödinger equation with the laser-dressed MGECSC potential. It is resulted that since applying a monochromatic laser field on hydrogen atom embedded in a Debye and quantum plasma causes to shift in the profile of the total interaction potential, the confinement effects of laser field on hydrogen atom in plasmas modeled by the MGECSC potential change localizations of energy states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First-principles study on structural mechanical and thermodynamic properties of HfMoTaTiZr

This study calculates the phase diagram of HfMoTaTiZr using the TCHEA1 database in Thermo-Calc software. Furthermore, it also investigates the structural and mechanical properties of HfMoTaTiZr by applying first-principles calculations. The elastic and other mechanical properties such as Young’s modulus, pugh ratio, hardness, bulk modulus was derived and discussed. These parameters were found to be in good harmony with previously available experimental reports. Vickers hardness of HfMoTaTiZr predicted from the Tians model agrees with the available experimental result. This study also investigates the thermal coefficient of linear expansion and vibrational heat capacity of the HfMoTaTiZr using the Debye–Grüneisen theory.

36 MATERIALS SCIENCE↗

Role of electrostatic fluctuations in doped semiconductors upon the transition from band to hopping conduction (by the example of p-Ge:Ga)

The electrostatic model of ionization equilibrium between hydrogen-like acceptors and v-band holes in crystalline covalent p-type semiconductors is developed. The range of applicability of the model is the entire insulator side of the insulator–metal (Mott) phase transition. The density of the spatial distribution of acceptor- and donor-impurity atoms and holes over a crystal was assumed to be Poissonian and the fluctuations of their electrostatic potential energy, to be Gaussian. The model takes into account the effect of a decrease in the energy of affinity of an ionized acceptor to a v-band hole due to Debye–Hückel ion screening by both free v-band holes and localized holes hopping over charge states (0) and (–1) of acceptors in the acceptor band. All donors are in charge state (+1) and are not directly involved in the screening, but ensure the total electroneutrality of a sample. In the quasiclassical approximation, analytical expressions for the root-mean-square fluctuation of the v-band hole energy W{sub p} and effective acceptor bandwidth W{sub a} are obtained. In calculating W{sub a}, only fluctuations caused by the Coulomb interaction between two nearest point charges (impurity ions and holes) are taken into account. It is shown that W{sub p} is lower than W{sub a}, since electrostatic fluctuations do not manifest themselves on scales smaller than the average de Broglie wavelength of a free hole. The delocalization threshold for v-band holes is determined as the sum of the diffusive-percolation threshold and exchange energy of holes. The concentration of free v-band holes is calculated at the temperature T{sub j} of the transition from dc band conductivity to conductivity implemented via hopping over acceptor states, which is determined from the virial theorem. The dependence of the differential energy of the thermal ionization of acceptors at the temperature 3T{sub j}/2 on their concentration N and degree of compensation K (the ratio between the donor and acceptor concentrations) is determined. Good quantitative agreement between the results of the calculation and data on the series of neutron transmutation doped p-Ge samples is obtained up to the Mott transition without using any fitting parameters.

36 MATERIALS SCIENCE↗

X-ray diffraction strains in laser-ablated aluminum, nickel, sodium, and invar: Pressures to 475 GPa

Dynamically compressed materials in longitudinal waves are described by two physical models: hydrostatic pressure, with equal, normal, principal stresses or material uniaxially strained in the wave propagation direction. These models are disparate, so experimental comparisons and evaluations are important. Polycrystalline material in a state of hydrostatic pressure will have no eccentricity of x-ray-diffracted Debye–Scherrer rings. A general three-dimensional solution of Bragg diffracted x rays based on principal crystallographic strains in the compression wave was found. The distortion of x-ray diffraction beams has been used for strain measurements; the analysis developed incorporates a strained reciprocal lattice and the incident x-ray beam. Strain-distorted Polanyi surfaces form an annulus of compression with an ellipsoid of revolution in reciprocal space, which is intersected by Ewald's sphere for Bragg diffraction. The in situ measurements for strain describe nanosecond diffraction evaluated using two planes, (hkl) and (h′k′l′), both in the same crystallographic phase. Diffraction from Al, Ni, Na, and Invar quantifies the compression axial strains in these materials: The compression axial ratios are 0.65, 1.05, 0.88, and 1.58 at pressures of 291, 402, 409, and 367 GPa for the respective materials. Crystal structure transformations with homogeneous pressurized stresses, mandating equal normal strains, should not be anticipated to agree with heterogeneous, uniaxially strained, and sheared crystalline phases. Measurements support in-plane strains increasing with pressure, p , in fcc and hcp aluminum as ε 22 = −0.21 p /TPa.

Alloys↗

Thermodynamics of Thomas-Fermi screened Coulomb systems

We obtain in closed analytic form, estimates for the thermodynamic properties of classical fluids with pair potentials of Yukawa type, with special reference to dense fully ionized plasmas with Thomas-Fermi or Debye-Hueckel screening. We further generalize the hard-sphere perturbative approach used for similarly screened two-component mixtures, and demonstrate phase separation in this simple model of a liquid mixture of metallic helium and hydrogen.

Firey, B.↗

Magnetic sheath adjacent to absorbing surfaces.

Major attention is given to the existence and structure of a 'magnetic sheath' adjacent to perfectly absorbing solid surfaces. The basic equations for a two-field continuum model are presented, and the quasi-neutral approximation and Bohm's criterion are briefly reviewed. The infinite conductivity MHD equations are recovered by invoking the small magnetic Debye length approximation. This approximation is relaxed, and the appropriate 'boundary layer' equation for the magnetic sheath is presented. The first magnetic sheath criterion is derived, which states that the conditions at the edge of the sheath, as characterized by a point in the Friedrichs diagram, must be located in certain restricted regions in order for a magnetic sheath to exist. The second magnetic sheath criterion is discussed.

Lam, S. H.↗

Giant Dipole Moments: Remarkable Effects Mono‐, Di‐, and Tri‐ Hydrated 5,6‐Diaminobenzene‐1,2,3,4‐Tetracarbonnitrile

The molecule 5,6-diaminobenzene-1,2,3,4-tetracarbonnitrile (MOI) was first synthesized by Müllen and coworkers in 2016 and boasts an ultrastrong dipole moment of $14.1\pm 0.7$ Debye in THF. Gas phase DFT computations do not fully reflect this ultrastrong dipole moment, demonstrating the role of solvent in increasing this dipole moment. Here, we investigate the effect of solvent molecule position on the dipole moment of this species, computationally examining systems with giant dipole moments. These systems are optimized in the gas phase with the B3LYP functional, employing the aug-cc-pVTZ and def2-TZVP basis sets, as well as the B3LYP-D3BJ/aug-cc-pVTZ functional in Orca. Single point DLPNO-CCSD/aug-cc-pVDZ results were obtained from Orca and Psi4, as well as DLPNO-CCSD(T)/CBS information from Psi4. Additionally, these are compared to the dipole moments of di- and tri-hydrated systems, and the SMD models for THF and water at the B3LYP/aug-cc-pVTZ level of theory. The dissociation energies, HOMO-LUMO energy gaps, and dipole moments are presented. These metrics show the nh1nh1′ THF system boasts the largest dissociation energy and dipole moment of the singly solvated systems, due to its strong hydrogen bonding. The importance of solvent placement is highlighted and may guide the synthesis of macromolecules or organic frameworks incorporating the MOI or MOI-like subunits. Remarkably, a single solvent molecule provides a good model for the difference between the gas phase and solvated species. The predicted gas phase dipole moments computed with B3LYP/aug-cc-pVTZ for the MOI, its monohydrated complex, dihydrated complex, and its trihydrated complex are 9.6, 14.2, 16.0, and 16.8 Debye, respectively.

dipole↗

A self-consistent two-fluid model of a magnetized plasma-wall transition

A self-consistent one-dimensional two-fluid model of the magnetized plasma-wall transition is presented. The model includes magnetic field, elastic collisions between ions and electrons, and creation/annihilation of charged particles. Two systems of differential equations are derived. The first system describes the whole magnetized plasma-wall transition region, which consists of the pre-sheath, the magnetized pre-sheath (Chodura layer), and the sheath, which is not neutral, but contains a positive space charge. The second system of equations describes only the neutral part of the plasma-wall transition region—this means only the pre-sheath and the Chodura layer, but not also the sheath. Both systems are solved numerically. The first system of equations has two singularities. The first occurs when ion velocity in the direction perpendicularly to the wall drops below the ion thermal velocity. The second occurs when the electron velocity in the direction perpendicularly to the wall exceeds the electron thermal velocity. The second system of differential equations only has one singularity, which has also been derived analytically. For finite electron to ion mass ratio, the integration of the second system always breaks down before the Bohm criterion is fulfilled. Some properties of the first system of equations are examined. It is shown that the increased collision frequency demagnetizes the plasma. On the other hand, if the magnetic field is so strong that the ion Larmor radius and the Debye length are comparable, the electron velocity in the direction perpendicularly to the wall reaches the electron thermal velocity before the ion velocity in the direction perpendicularly to the wall reaches the ion sound velocity. In this case, the integration of the model equations breaks down before the Bohm criterion is fulfilled and the sheath is formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamical Theory of X-Ray Diffraction for Restricted Beams: II. Diffuse Scattering by a Porous Crystal

A theory of X-ray diffuse scattering of restricted beams by a porous crystal is developed. Models of cylindrical and ellipsoidal pores, as well as the pores in the form of rectangular and triangular prisms, are considered. For these pores, the expressions for correlation eigenfunctions, static Debye—Weller factors, and correlation volumes are derived. The effect of the inclination and spatial correlation of pores on the equal-intensity diffuse-scattering contours is analyzed. The influence of the pore size fluctuation on the scattered intensity distribution is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure and disorder resulting in ultralow thermal conductivity in the defect chalcopyrite AgInSnSe 4

Quaternary chalcogenides are attractive for a variety of technological fields of interest due to the diverse physical properties these materials possess. Herein, the structural and thermal properties of the disordered chalcopyrite AgInSnSe 4 are reported revealing a distinct relation between its structure and physical properties. Modeling of the experimental temperature-dependent thermal properties revealed an ultralow thermal conductivity (0.47 Wm -1 K -1 at room temperature) due to lattice anharmonicity, a low speed of sound and a low Debye temperature. Moreover, first principles electronic structure calculations revealed that the Ag-Se tetrahedra within the crystal structure possess relatively weak bonds due to occupied antibonding states from p-d orbital hybridization, which suppress the thermal conductivity. Here, in order to further quantify our findings, we extended our analyses to include comparisons with other ternary and quaternary adamantine materials, all of which possess a fourfold tetrahedral coordination of atoms, in revealing the origin of the thermal properties in AgInSnSe 4 .

36 MATERIALS SCIENCE↗

Simulations of Solar Wind Plasma Flow Around a Simple Solar Sail

In recent years, a number of solar sail missions of various designs and sizes have been proposed (e.g., Geostorm). Of importance to these missions is the interaction between the ambient solar wind plasma environment and the sail. Assuming a typical 1 AU solar wind environment of 400 km/s velocity, 3.5 cu cm density, ion temperature of approx.10 eV, electron temperature of 40 eV, and an ambient magnetic field strength of 10(exp -4) G, a first order estimate of the plasma interaction with square solar sails on the order of the sizes being considered for a Geostorm mission (50 m x 50 m and 75 m x 75 m corresponding to approx.2 and approx.3 times the Debye length in the plasma) is carried out. First, a crude current balance for the sail surface immersed in the plasma environment and in sunlight was used to estimate the surface potential of the model sails. This gave surface potentials of approx.10 V positive relative to the solar wind plasma. A 3-D, Electrostatic Particle-in-Cell (PIC) code was then used to simulate the solar wind flowing around the solar sail. It is assumed in the code that the solar wind protons can be treated as particles while the electrons follow a Boltzmann distribution. Next, the electric field and particle trajectories are solved self-consistently to give the proton flow field, the electrostatic field around the sail, and the plasma density in 3-D. The model sail was found to be surrounded by a plasma sheath within which the potential is positive compared to the ambient plasma and followed by a separate plasma wake which is negative relative to the plasma. This structure departs dramatically from a negatively charged plate such as might be found in the Earth s ionosphere on the night side where both the plate and its negative wake are contiguous. The implications of these findings are discussed as they apply to the proposed Geostorm solar sail mission.

Garrett, Henry B.↗