Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass density ratio”

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 325 records · Page 18

Optimal Design of Functionally Graded Metallic Foam Insulations

The focus of our work has been on developing an insight into the physics that govern the optimum design of thermal insulation for use in thermal protection systems of launch vehicle. Of particular interest was to obtain optimality criteria for designing foam insulations that have density (or porosity) distributions through the thickness for optimum thermal performance. We investigate the optimum design of functionally graded thermal insulation for steady state heat transfer through the foam. We showed that the heat transfer in the foam has competing modes, of radiation and conduction. The problem assumed a fixed inside temperature of 400 K and varied the aerodynamic surface heating on the outside surface from 0.2 to 1.0 MW/sq m. The thermal insulation develops a high temperature gradient through the thickness. Investigation of the model developed for heat conduction in foams showed that at high temperatures (as on outside wall) intracellular radiation dominates the heat transfer in the foam. Minimizing radiation requires reducing the pore size, which increases the density of the foam. At low temperatures (as on the inside wall), intracellular conduction (of the metal and air) dominates the heat transfer. Minimizing conduction requires increasing the pore size. This indicated that for every temperature there was an optimum value of density that minimized the heat transfer coefficient. Two optimization studies were performed. One was to minimize the heat transmitted though a fixed thickness insulation by varying density profiles. The second was to obtain the minimum mass insulation for specified thickness. Analytical optimality criteria were derived for the cases considered. The optimality condition for minimum heat transfer required that at each temperature we find the density that minimizes the heat transfer coefficient. Once a relationship between the optimum heat transfer coefficient and the temperature was found, the design problem reduced to the solution of a simple nonlinear differential equation. Preliminary results of this work were presented at the American Society of Composites meeting, and the final version was submitted for publication in the AIAA Journal. In addition to minimizing the transmitted heat, we investigated the optimum design for minimum weight given an acceptable level of heat transmission through the insulation. The optimality criterion developed was different from that obtained for minimizing beat transfer coefficient. For minimum mass design, we had to find for a given temperature the optimum density, which minimized the logarithmic derivative of the insulation thermal conductivity with respect to its density. The logarithmic derivative is defined as the ratio of relative change in the dependent response (thermal conductivity) to the relative change in the independent variable (density). The results have been documented as a conference paper that will be presented at the upcoming AIAA.

Haftka, Raphael T.↗

The Ca II emission lines of SN 1987A

We study the evolution of the forbidden Ca II 7300 A and Ca II 8600 A emission lines in the spectrum SN 1987A. We can explain the observations with a model in which (1) the lines come from about 1.7 x 10 exp -4 solar mass of calcium distributed with filling factor about 0.1 in a sphere expanding with velocity 2500 km/s; (2) the electron number density is given by 4 x 10 exp 9/cu cm/(t/100d)-cubed; (3) the gas temperature decreases from T roughly 6400 K at 200 days to T roughly 3900 K at 500 days; and (4) ultraviolet pumping in the Ca II 3934 A, 3969 A (K, H) lines is required to explain the line ratios for t greater than about 350 days. The mass of calcium in our model is consistent with the LMC abundance of calcium in about 5 solar mass of hydrogen in the inner envelope of SN 1987A. If more calcium is produced by supernova nucleosynthesis, it must be in dense clumps that are too cool to emit the forbidden Ca II 7300 A and Ca II 8600 A lines.

Li, Hongwei↗

Current constraints on hidden mass in the Coma Cluster

Various possible means of accounting for the dynamical mass of the Coma Cluster are reviewed. We intentionally choose to emphasize conventional forms of hidden mass, specifically gas and low-luminosity stars, while adopting conventional mass-to-luminosity ratios for the observed luminous material in galaxies. We have attempted to combine all previous work relating to assumed gas parameters in order to examine precisely what constraints current observations, together with theoretical and observational upper limits, place on any gas in the Coma Cluster. In particular, we make use of a density distribution function that is consistent with the dynamical data and express all constraints in a manner which explicitly displays their dependence on Hubble constant, current observational flux levels, distance from the dynamical center of the cluster, and the mass of gas postulated.

Tarter, J.↗

Heavy metal toxicity as a kill mechanism in impact caused mass extinctions

Heavy metals that are known to be toxic exist in carbonaceous chrondrites at abundances considerably in excess to that of the terrestrial crust. An impactor of relatively undifferentiated cosmic matter would inject into the terrestrial environment large quantities of toxic elements. The abundances of toxic metals found in the Allende CV carbonaceous chondrite and the ratio of meteoritic abundance to crustal abundance are: Cr, 3630 PPM, 30X; Co, 662 PPM, 23X; ni, 13300 PPm, 134X; se, 8.2 PPM, 164X; Os, 0.828 PPM, 166X. The resulting areal density for global dispersal of impactor derived heavy metals and their dilution with terrestrial ejecta are important factors in the determination of the significance of impactor heavy metal toxicity as a kill mechanism in impact caused mass extinctions. A 10 km-diameter asteroid having a density of 3 gram per cu cm would yield a global areal density of impact dispersed chondritic material of 3 kg per square meter. The present areal density of living matter on the terrestrial land surface is 1 kg per square meter. Dilution of impactor material with terrestrial ejecta is determined by energetics, with the mass of ejecta estimated to be in the range of 10 to 100 times that of the mass of the impactor. Because a pelagic impact would be the most likely case, the result would be a heavy metal rainout.

Wdowiak, T. J.↗

Injection and distribution of liquids in the flow fields of blunt shapes at hypersonic speeds

A possible method for reducing the attenuation of radio transmission from hypersonic vehicles is by injection of a liquid into the flow field. The hot plasma layer surrounding the vehicle could presumably be quenched if the latent heat of the liquid and the heat capacity of the vapor could be utilized efficiently. In the case of stagnation-point injection, the momentum of the liquid could be utilized to modify the bow shock and possibly reduce the initial ionization. To attain these objectives of efficient coolant utilization or modification of the bow shock a knowledge of the liquid injection, atomization, spray penetration, and evaporation processes is required. This report presents the results of a study of liquid injection into hypersonic flow fields. A simplified theory for evaporation distances is derived and is shown to be in qualitative agreement with experimental data. Correlation of experimental data for maximum cross-current penetration of the liquid spray for injection near the shoulder pf a hemisphere cylinder, a sphere 9° cone, and a Mercury model are given. These data were obtained in the Langley Mach 8 variable density tunnel with both liquid nitrogen and water injection to provide a large range of liquid parameters. Flow-field disturbances, as indicated by schlieren photographs and surface pressures, increased with the ratio of coolant mass flux to airstream mass flux. For stagnation-point injection, the maximum shock standoff distance was correlated by the ratio of the liquid jet momentum flux to the airstream momentum flux.

Ivan E. Beckwith↗

Gravity Field and Internal Structure of Mercury from MESSENGER

Radio tracking of the MESSENGER spacecraft has provided a model of Mercury's gravity field. In the northern hemisphere, several large gravity anomalies, including candidate mass concentrations (mascons), exceed 100 milli-Galileos (mgal). Mercury's northern hemisphere crust is thicker at low latitudes and thinner in the polar region and shows evidence for thinning beneath some impact basins. The low-degree gravity field, combined with planetary spin parameters, yields the moment of inertia C/M(R(exp 2) = 0.353 +/- 0.017, where M and R are Mercury's mass and radius, and a ratio of the moment of inertia of Mercury's solid outer shell to that of the planet of C(sub m)/C = 0.452 +/- 0.035. A model for Mercury s radial density distribution consistent with these results includes a solid silicate crust and mantle overlying a solid iron-sulfide layer and an iron-rich liquid outer core and perhaps a solid inner core.

Smith, David E.↗

Gravity Field and Internal Structure of Mercury From MESSENGER

Radio tracking of the MESSENGER spacecraft has provided a model of Mercury’s gravity field. In the northern hemisphere, several large gravity anomalies, including candidate mass concentrations (mascons), exceed 100 milli-Galileos (mgal). Mercury’s northern hemisphere crust is thicker at low latitudes and thinner in the polar region and shows evidence for thinning beneath some impact basins. The low-degree gravity field, combined with planetary spin parameters, yields the moment of inertia C/MR(exp 2) = 0.353 ± 0.017, where M and R are Mercury’s mass and radius, and a ratio of the moment of inertia of Mercury’s solid outer shell to that of the planet of C(sub m)/C = 0.452 ± 0.035. A model for Mercury’s radial density distribution consistent with these results includes a solid silicate crust and mantle overlying a solid iron-sulfide layer and an iron-rich liquid outer core and perhaps a solid inner core.

David E Smith↗

Planetary atmospheres

During the period from 1979 to 1982, Pioneer and Voyager missions to Venus, Jupiter, and Saturn combined with ongoing analysis and data collection from the Viking-Mars mission have contributed greatly to an understanding of planetary atmospheres. Research concerning the atmosphere of Venus is discussed, taking into account the atmospheric composition, the isotopic ratios of noble gases, data concerning temperature and pressure, information regarding various atmospheric layers, atmospheric density data, and models regarding the dynamical properties of Venus' atmosphere. Information concerning the atmosphere of Mars is also examined. Isotopic ratios of carbon and oxygen are considered along with altitude profiles, diurnal heating and cooling, and thermal tides and their relationship to dust storms. Attention is given to the temperature structure of the Jovian atmosphere, the morphology of Jovian cloud features, features in the Saturnian atmosphere, information with respect to the atmosphere of Titan, the rotation rates of Uranus, masses and densities of Uranus and Neptune, and the identification of Pluto as a small icy body.

Beebe, R.↗

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

Early Results from GLASS-JWST. XI. Stellar Masses and Mass-to-light Ratio of z > 7 Galaxies

We exploit James Webb Space Telescope (JWST) NIRCam observations from the GLASS-JWST-Early Release Science program to investigate galaxy stellar masses at z > 7. We first show that JWST observations reduce the uncertainties on the stellar mass by a factor of at least 5–10, when compared with the highest-quality data sets available to date. We then study the UV mass-to-light ratio, finding that galaxies exhibit a a two orders of magnitude range of M / L UV values for a given luminosity, indicative of a broad variety of physical conditions and star formation histories. As a consequence, previous estimates of the cosmic stellar-mass density—based on an average correlation between UV luminosity and stellar mass—can be biased by as much as a factor of ~6. Our first exploration demonstrates that JWST represents a new era in our understanding of stellar masses at z > 7 and, therefore, of the growth of galaxies prior to cosmic reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

Upcycling linear low-density polyethylene waste to turbostratic graphene for high mass loading supercapacitors

Linear low-density polyethylene (LLDPE) waste is difficult to upcycle into more valuable carbon materials because it tends to completely decompose into small molecules during thermal processing. In this work, LLDPE is upcycled into a high quality turbostratic graphene using a pre-treatment step to oxidatively crosslink the polymer with the assistance of solid additives (KCl and K 2 CO 3 ) that improve crosslinking by increasing the effective surface area of the polymer melt during processing. After this pretreatment step, the crosslinked polymer could then be carbonized and catalytically graphenized between 400–950°C without decomposition of the polymer feedstock. The LLDPE derived graphene (LLDPE-G) obtained from this process has a Brunauer–Emmett–Teller (BET) specific surface area, up to 1800m 2 g -1 and average Raman I D /I G and I 2D /I G ratios of 0.85 and 0.57, respectively, indicating high quality graphene. When used as an electrode material in symmetric supercapacitors, LLDPE-G possesses a specific capacitance up to 175Fg -1 at a mass loading of 20mgcm -2 , which is two times the commercial requirement, yielding an areal capacitance of 3.5 Fcm -2 . Moreover, LLDPE-G exhibits exceptional cycling stability with a capacitance retention of 95.8 % after 100,000 cycles at a current density of 4.0 Ag -1 . Additionally, the KCl and K 2 CO 3 solids are recycled and reused over 3 complete reaction cycles to make new LLDPE-G with the material quality and electrocapacitive performance retained and verified after each cycle. Our approach creates new opportunities for upcycling waste LLDPE and other varieties of polyethylene into a higher value graphene used for electrochemical energy storage applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-linear flow modes of identified particles in Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5.02 TeV

The pT-differential non-linear flow modes, v 4,22 , v 5,32 , v 6,33 and v 6,222 for π ± , K ± , K$^0_S$ , $p + \bar{p}, Λ + \bar{Λ}$ and ϕ-meson have been measured for the first time at $\sqrt{s_{NN}}$ = 5.02 TeV in Pb-Pb collisions with the ALICE detector at the Large Hadron Collider. The results were obtained with a multi-particle technique, correlating the identified hadrons with reference charged particles from a different pseudorapidity region. These non-linear observables probe the contribution from the second and third order initial spatial anisotropy coefficients to higher flow harmonics. All the characteristic features observed in previous p T -differential anisotropic flow measurements for various particle species are also present in the non-linear flow modes, i.e. increase of magnitude with increasing centrality percentile, mass ordering at low p T and particle type grouping in the intermediate p T range. Hydrodynamical calculations (iEBE-VISHNU) that use different initial conditions and values of shear and bulk viscosity to entropy density ratios are confronted with the data at low transverse momenta. These calculations exhibit a better agreement with the anisotropic flow coefficients than the non-linear flow modes. These observations indicate that non-linear flow modes can provide additional discriminatory power in the study of initial conditions as well as new stringent constraints to hydrodynamical calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sustainable recovery of critical metals from spent lithium-ion batteries through gluconic acid-based bioleaching: Techno-economic analysis, life cycle assessment and process optimization

Recycling spent lithium-ion batteries (LIB) could potentially bridge the ever increasing supply and demand gap for critical metals and simultaneously facilitate the management of hazardous battery waste. This study investigated the optimization of gluconic acid-based bioleaching technology through design of experiments (DOE), combined with techno-economic analysis (TEA), and life cycle assessment (LCA) with the aim of maximizing the net present value (NPV) and minimizing global warming impacts of the process. Biolixiviant containing predominantly gluconic acid produced by the genetically engineered (ΔpstS, P 112 :mgdh) Gluconobacter oxydans B58 through fermentation using non-recyclable paper as a growth substrate was used for the LIB leaching. At optimal bioleaching conditions of gluconic acid (160 mM), leaching time (2.5 h), reducing agent FeSO 4 to metal, i.e., cobalt (Co), nickel (Ni) and manganese (Mn), mole ratio (0.88), temperature (55 °C) and pulp density (2.5 %), the leaching efficiency was 87 % 72 %, 94 %, and 88 % for Co, Ni, Mn and lithium (Li), respectively. TEA analysis confirmed that bioleaching plant with an annual black mass processing capacity of 10,000 metric tons and plant life of 30 years would be economically viable with an NPV and profit margin of $136 million and 11 %, respectively. The predicted carbon footprint of gluconic acid-based bioleaching for recovering 1 kg of Co (13.2 kg of CO 2 eq.) is lower compared to that of most state-of-the-art leaching technologies. Moreover, gluconic acid-based bioleaching effectively recovered target metals when tested for different black mass chemistries.

Bioleaching↗

Observations of solar wind ion charge exchange in the comet Halley coma

Giotto Ion Mass Spectrometer/High Energy Range Spectrometer (IMS/HERS) observations of solar wind ions show charge exchange effects and solar wind compositional changes in the coma of comet Halley. As the comet was approached, the He(++) to proton density ratio increased until about 1 hour before closest approach after which time it decreased. Abrupt increases in this ratio were also observed in the beginning and near the end of the so-called Mystery Region (8.6 - 5.5(10)(exp 5) km from the comet along the spacecraft trajectory). These abrupt increases in the density ratio were well correlated with enhanced fluxes of keV electrons as measured by the Giotto plasma electron spectrometer. The general increase and then decrease of the He(++) to proton density ratio is quantitatively consistent with a combination of the addition of protons of cometary origin to the plasma and loss of plasma through charge exchange of protons and He(++). In general agreement with the solar wind proton and He(++) observations, solar wind oxygen and carbon ions were observed to charge exchange from higher to lower charge states with decreasing distance to the comet. The more abrupt increases in the He(++) to proton and the He(++) to O(6+) density ratios in the mystery region require a change in the solar wind ion composition in this region while the correlation with energetic electrons indicates processes associated with the comet.

Fuselier, S. A.↗

Coronagraphic imaging of pre-main-sequence stars: Remnant evvelopes of star formation seen in reflection

We have obtained R- and I-band coronagraphic images of the vicinities of 11 pre-main sequence (PMS) stars to search for faint, small-scale reflection nebulae. The inner radius of the search and the field of view are 1.9 arcsec and 1x1 arcmin, respectively. Reflection nebulae were imaged around RY Tau, T Tau,DG Tau, SU Aur, AB Aur, FU Ori, and Z CMa. No nebulae were detected around HBC 347, GG Tau, V773 Tau, and V830 Tau. Categorically speaking, most of the classical T Tauri program stars and all the FU Orionis-type program stars are associated with the reflection nebulae, while none of the weak-line T Tauri program stars are associated with nebulae. The detected nebulae range in size from 250 to 37 000 AU. From the brightness ratios of the stars and nebulae, we obtain a lower limit to the visual extinction of PMS star light through the nebulae of (A(sub V))(sub neb) = 0.1. The lower limits of masses and volume densities of the nebulae associated with the classical T Tauri stars are 10(exp-6) Solar mass and N(sub H) = 10(exp 5)/cu cm, respectively. Lower limits for the nebulae around FU Orionis stars are 10(exp -5) Solar mass and n(sub H) = 10 (exp 5)/cu cm, respectively. Some reflection nebulae may trace the illuminated surfaces of the optically thick dust nebulae, so these mass estimates are not stringent. All the PMS stars with associated nebulae are strong far-infrared emitters. Both the far-infrared emission and the reflection nebulae appear to originate from the remnant envelopes of star formation. The 100 micrometers emitting regions of SU Aur and FU Ori are likely to be cospatial with the reflection nebulae. A spatial discontinuity between FU Ori and its reflection nebula may explain the dip in the far-infrared spectral energy distribution at 60 micrometers. The warped, disk-like nebulae around T Tau and Z CMa are aligned with and embrace the inner star/circumstellar disk systems. The arc-shaped nebula around DG Tau may be in contact with the coaligned inner star/disk system. These three-reflection nebulae may trace the surfaces of pseudodisks from which matter accretes onto the stars or the inner circumstellar disks. 19 stellar objects brighter than I = 19 were detected around 9 program stars. Using a color-magnitude diagram, we have identified three new PMS candidates aroun Z CMa and one previously known PMS candidate, GG Tau/c.

Nakajima, Tadashi↗

Launch Vehicle Performance for Bipropellant Propulsion Using Atomic Propellants With Oxygen

Atomic propellants for bipropellant launch vehicles using atomic boron, carbon, and hydrogen were analyzed. The gross liftoff weights (GLOW) and dry masses of the vehicles were estimated, and the 'best' design points for atomic propellants were identified. Engine performance was estimated for a wide range of oxidizer to fuel (O/F) ratios, atom loadings in the solid hydrogen particles, and amounts of helium carrier fluid. Rocket vehicle GLOW was minimized by operating at an O/F ratio of 1.0 to 3.0 for the atomic boron and carbon cases. For the atomic hydrogen cases, a minimum GLOW occurred when using the fuel as a monopropellant (O/F = 0.0). The atomic vehicle dry masses are also presented, and these data exhibit minimum values at the same or similar O/F ratios as those for the vehicle GLOW. A technology assessment of atomic propellants has shown that atomic boron and carbon rocket analyses are considered to be much more near term options than the atomic hydrogen rockets. The technology for storing atomic boron and carbon has shown significant progress, while atomic hydrogen is not able to be stored at the high densities needed for effective propulsion. The GLOW and dry mass data can be used to estimate the cost of future vehicles and their atomic propellant production facilities. The lower the propellant's mass, the lower the overall investment for the specially manufactured atomic propellants.

Palaszewski, Bryan↗

Upcycling Low Linear Density Polyethylene Waste into Turbostratic Graphene for High Mass Loading Supercapacitors

In this work, LLDPE was upcycled into a high quality turbostratic graphene using a pre-treatment step to oxidatively crosslink the polymer with the assistance of solid additives (KCl and K2CO3) that improve crosslinking by increasing the effective surface area of the polymer melt during processing. After this pretreatment step, the crosslinked polymer could then be carbonized and catalytically graphenized between 400-950 °C without complete decomposition of the material. The LLDPE derived graphene (LLDPE-G) obtained from this process has a Brunauer–Emmett–Teller (BET) specific surface area, up to 1800 m2g-1 and average Raman ID/IG and I2D/IG ratios of 0.85 and 0.57, respectively, indicating high quality graphene. When used as an electrode material in symmetric supercapacitors, LLDPE-G possesses an outstanding specific capacitance up to 175 Fg-1 at a mass loading of 20 mgcm-2, which is two times the commercial requirement, yielding an excellent areal capacitance of 3.5 Fcm-2. Moreover, LLDPE-G exhibits exceptional cycling stability with a capacitance retention of 95.8% after 100,000 cycles at a current density of 4.0 Ag-1. Additionally, the KCl and K2CO3 were recycled and reused over 3 complete cycles to make new LLDPE-G with the material quality and electrocapacitive performance retained and verified after each cycle. Our approach creates new opportunities for upcycling not only waste LLDPE but also other varieties of PE to high value graphene materials.

Gao, Yuan [NETL Site Support Contractor, National ↗

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗