Daoru Frank Han

A 3D immersed finite element method with non-homogeneous interface flux jump for applications in particle-in-cell simulations of plasma-lunar surface interactions

Motivated by the need to handle complex boundary conditions efficiently and accurately in particle-in-cell (PIC) simulations, this paper presents a three-dimensional (3D) linear immersed finite element (IFE) method with non-homogeneous flux jump conditions for solving electrostatic field involving complex boundary conditions using structured meshes independent of the interface. This method treats an object boundary as part of the simulation domain and solves the electric field at the boundary as an interface problem. In order to resolve charging on a dielectric surface, a new 3D linear IFE basis function is designed for each interface element to capture the electric field jump on the interface. Numerical experiments are provided to demonstrate the optimal convergence rates in L 2 and H 1 norms of the IFE solution. This new IFE method is integrated into a PIC method for simulations involving charging of a complex dielectric surface in a plasma. A numerical study of plasma-surface interactions at the lunar terminator is presented to demonstrate the applicability of the new method.

A Nonhomogeneous Immersed-Finite-Element Particle-in-Cell Method for Modeling Dielectric Surface Charging in Plasmas

We present a particle-in-cell (PIC) method using a nonhomogeneous immersed-finite-element (IFE) field solver for modeling dielectric surface charging of complex-shaped objects in plasmas. The IFE solver allows PIC codes using a Cartesian mesh applied to simulations involving arbitrarily shaped objects with a similar accuracy as that using a body-fitting mesh. The object surface is treated as an interface. Surface charging is calculated directly from charge deposition at the interface, and the electrostatic fields on both sides of the interface are resolved self-consistently. The capability of the nonhomogeneous IFE-PIC method is demonstrated by a simulation study of the charging of an irregular-shaped asteroid in the solar wind.

Kinetic Simulations of Plasma Plume Potential in a Vacuum Chamber

Direct simulation Monte Carlo and particle-in-cell simulations are carried out to study the potential of a mesothermal plasma plume in a vacuum chamber. The results show that the beam potential with respect to the ambient in a vacuum chamber is different from that in space because the facility plasma can prematurely terminate the plume expansion process. As a result, the plume potential measured in a vacuum chamber may be significantly lower than that under the in-space condition. This can lead to underestimation of the backflow of charge-exchange ions and ionized contaminants in plasma thruster plume modeling.

Simulations of Ion Thruster Plume Contamination with a Whole Grid Sputtered Mo Source Model

A particle simulation based source model is developed to calculate the density distribution of the sputtered Mo atoms for a whole ion optics grid. The source model is used in PIC simulation of ion thruster plume contamination for 3-grid and 2-grid ion thrusters. The results show that the commonly used point-source approximation for sputtered Mo atoms is oversimplified and would lead to over-prediction of contamination deposition.

Uncertainty Quantification Integrated to the CFD Modeling of Synthetic Jet Actuators

The Point Collocation Non-Intrusive Polynomial Chaos (NIPC) method has been applied to two stochastic synthetic jet actuator problems used as test cases in the CFDVAL2004 workshop to demonstrate the integration of computationally efficient uncertainty quantification to the high-fidelity CFD modeling of synthetic jet actuators. In Case1 where the synthetic jet is issued into quiescent air, the NIPC method is used to quantify the uncertainty in the long-time averaged u and v-velocities at several locations in the flow field, due to the uniformly distributed uncertainty introduced in the amplitude and frequency of the oscillation of the piezo-electric membrane. Fifth order NIPC expansions were used to obtain the uncertainty information, which showed that the variation in the v-velocity is high in the region directly above the jet slot and the variation in the u-velocity is maximum in the region immediately adjacent to the slot. Even with a ten percent variation in the amplitude and frequency, the long-time averaged u and v velocity profiles could not match the experimental measurements at y=0.1mm above the slot indicating that the discrepancy may be due to other uncertainty sources in CFD or measurement errors. In Case 2 which includes a cross flow, the free stream velocity is treated as an uncertain input variable. Fifth degree NIPC expansions were employed to quantify the uncertainty in phase averaged velocity profiles as well as long-time averaged wall pressure and skin friction coefficient distributions. The results of Case 2 show that the uncertainty in phase averaged velocity profiles gets larger when approaching the main stream. The size of a separation bubble observed in this case remains relatively insensitive to the uncertain free stream velocity within the tolerance range considered.