Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Laws”

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 19 records

Application of an Empirical Density Law via Python for Aqueous Plutonium Chloride Systems in MCNP6

Current aqueous plutonium processing models for criticality safety often contain significant bias due to material modeling assumptions. These solutions include plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution. Furthermore, there is no current predictive capability for modeling plutonium metal-water mixtures that is approved for use at Los Alamos National Laboratory (LANL). Recent density measurements for aqueous plutonium chloride systems (PuCl 3 -HCl-H 2 O) now allow for the development and application of a more realistic density law. This work develops a Python-based density law for this ternary solution using an empirical method. This code can be used in conjunction with an MCNP6 input to determine the density and composition of a solution based upon user inputs of plutonium concentration, hydrochloric acid concentration, and temperature. The tool allows users to input plutonium and acid content of a solution in terms of molality, molarity, or concentration, and predicts density within the current data range within 1.4% of experimental data. Current preliminary MCNP6 calculations utilizing this tool have demonstrated a minimum decrease in system reactivity of 5% in comparison to the current modeling conventions. Thus, this tool enables more accurate criticality safety operational limits by better crediting chorine content while still maintaining necessary conservatism.

Density Law↗

A Semi-Empirical Density Law for Ternary, Homogeneous PuCl 3 /HCl/H 2 O Solutions

Nuclear material operations pose unique hazards that are not encountered in other chemical, energy, or manufacturing industries. One of these hazards is the potential for a nuclear criticality accident when handling fissile isotopes such as 235 U and 239 Pu. These hazards are particularly high when fissile material is dissolved in solution as the neutron behaviors of the system can change rapidly with the physical and chemical changes accessible in solution. Current estimates of solution density used for criticality safety are outdated and hinder fissile material handling. Developing new estimates for these safety calculations requires experimental characterization and the derivation of empirical density models. We have derived a density law describing PuCl 3 /HCl/H 2 O solutions from experimental data characterizing solution density. Density data was treated using a Pitzer-derived eight-parameter equation, defining density as a function of analyte concentrations, temperature, and interactions between these variables. The model is predictive across the concentration and temperature ranges from which it was derived. The potential effects of varying oxidation states of plutonium, which are easily accessible in aqueous media, on the bulk solution density of the ternary system were also investigated. The resulting Pitzer-derived density law was applied to a nuclear criticality safety model, and the impact of the experimental characterization of solution density relative to previous estimates was demonstrated to be significant and suggest that the current approach to estimating density in nuclear criticality safety calculations may lead to overly conservative controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of an Empirical Density Law via Python for Aqueous Plutonium Chloride Systems for MCNP6

Criticality safety models for aqueous plutonium chloride systems often contain a significant bias due to assumptions in material compositions. These systems are currently modeled as a fictitious metal-water mixture because little is known about the true solution density. Furthermore, no predictive density tools or capabilities for modeling aqueous plutonium chloride systems are approved for use at Los Alamos National Laboratory. Recent density measurements of this ternary system (PuCl 3 -HCl-H 2 O) have allowed for the development of a more realistic density law, which is applied in this work via an empirical method based in Python. This tool, entitled PuCS (Plutonium Chloride Solution tool) may be used to determine solution density and composition based on the plutonium content, acid content, and temperature for MCNP6 inputs. PuCS has been found to predict density within 2% of experimental data. In conclusion, MCNP6 calculations have found that crediting minimal amounts of free acid (0.5 M) may correspond to a ~12% decrease in peak reactivity in comparison to current modeling methods.

42 ENGINEERING↗

Application of a Density Law via Python for Aqueous Plutonium Nitrate

A predictive density tool has been developed in Python to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer Method and an empirical method were implemented into the tool, allowing for plutonium nitrate density calculations. Additionally, the Python tool can generate atom densities for a MCNP6.2 material card using the density from the selected method and directly the densities into a prepared MCNP6 input text file. The material card and density are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, plutonium isotope weight percentages and impurity concentrations. The Python tool has been validated and verified against the International Handbook of Evaluated Criticality Safety Benchmark Experiments to predict densities within a root mean square error of 1.0% for the Pitzer method and 1.8% for the empirical method. These errors in density were shown to lead to a ±0.5% error in MCNP6.2 calculated k effective for the Pitzer method and a ±1.7% error for the Empirical method. Simultaneous work is also being done at the University of New Mexico and Los Alamos National Laboratory to create a similar tool for plutonium chloride solutions, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also lead the way towards a better understanding of solution systems and potential relaxation in the conservatism of the current aqueous plutonium processing limits.

97 MATHEMATICS AND COMPUTING↗

Experimental Steps toward a Density Law for Chlorine-Crediting Criticality Models of Aqueous Plutonium solutions

Criticality accident prevention is an essential safety consideration for all operations involving fissionable and fissile materials. Predictive criticality calculations using advanced neutron transport codes, such as Monte Carlo N-Particle Transport (MCNP), are invaluable tools for designing and implementing operational limits. While indispensable, these tools are limited by the quality and accuracy of the inputs that the user provides to define the modeled system. Parameters such as atomic composition, shape, and material density must be accurately defined to obtain a meaningful result. In the case of material property details for a plutonium model, accurate material characterization data is sometimes sparse. Thus limitations in calculation accuracy can be reduced by improving our knowledge of material properties in the targeted fissile systems. One of the biggest remaining challenges to accurately defining fissile systems is the description of aqueous fissile solutions. Even simple properties, such as density, are not well-known for fissile solutions relevant to nuclear energy and security. Described here are initial efforts undertaken to improve criticality calculation inputs for fissile plutonium chloride solutions in water. This effort is focused on experimentally determining accurate solution characteristics for ternary plutonium chloride/hydrochloric acid/water systems, by measurement of water activity and solution density. The effect of inputting experimental densities for these solutions into MCNP criticality calculations is compared to the traditional approach of modeling an idealized (and fictitious) plutonium-water mixture. Expansion of these efforts to a working density law for aqueous plutonium chloride solutions is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of Chlorine Capture and a Proposed Density Law on the Reactivity of Plutonium Solution Systems

During fissionable material processing, all normal and credible abnormal conditions must remain safely subcritical. Nuclear Criticality Safety (NCS) uses a number of methods to determine subcriticality, one of which is the use of neutron transport codes such as MCNP6. In order to create models for use with MCNP6, both the geometry and materials in fissionable material processes must be known, or assumptions must be made and quantified for the impact to bias. One of the systems with a significant amount of bias due to material modeling assumptions is in the area of aqueous plutonium processing. These solutions are typically plutonium nitrate solutions or plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution and there is no current predictive capability approved for use at Los Alamos National Laboratory (LANL) for modeling them. This research is currently underway to fill the gap and develop an algorithm for use with MCNP6 to model the density of plutonium chloride solutions. The method is to be validated with experimental data for density, and also validated with critical experiments using MCNP6. Note that the Chlorine Worth Study (CWS) was performed in December 2021 to help bridge the gap in chlorine data for critical experiments, and is currently awaiting International Criticality Safety Benchmark Evaluation Project (ICSBEP) review. This study was performed by LANL at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS). Additional information regarding this experiment may currently be found in LA-UR- 22-29180. Additionally, the Chemistry-Actinide Analytical Chemistry (C-AAC) at LANL has performed a number of solution density measurements for PuCl 3 -HC 1 -H 2 O, allowing for such data be used to create a semi-empirical density via the Pitzer method. The published dataset for the measurements is documented in LA-UR-22-25454. This method has already been tested successfully for aqueous plutonium nitrate solutions in SCALE. Current solution density measurements exist of plutonium concentrations of 0-~142g/L, all at 2M HC1, at temperatures 20-40°C. Additional data was taken for HC1-corrected density values, which essentially mimics the data for a pure PuCl x -water solution. The calculations in this report aim to support the current research by demonstrating the difference in system reactivity for the current modeling method when compared to the new proposed modeling with a density law implementation, which is being written as a Python tool to be used with MCNP6.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Universal scaling laws and density slopes for dark matter haloes

Smalls scale challenges suggest some missing pieces in our current understanding of dark matter. In this work, a cascade theory for dark matter is proposed to provide extra insights, similar to the cascade phenomenon in hydrodynamic turbulence. The kinetic energy is cascaded in dark matter from small to large scales involves a constant rate $ε_u$ (≈–4.6×10 –7 m 2 / s 3 ). Confirmed by N-body simulations, the energy cascade leads to a two-thirds law for kinetic energy $v^2_r$ on scale $\textit{r}$ such that $v^2_r∝(ε_ur)^{2/3}$. Equivalently, a four-thirds law can be established for mean halo density $ρ_s$ enclosed in the scale radius $r_s$ such that $ρ_s ∝ε^{2/3}_uG^{–1}r^{–4/3}_s$, which was confirmed by galaxy rotation curves. Critical properties of dark matter might be obtained by identifying key constants on relevant scales. First, the largest halo scale $r_l$ can be determined by $–u^3_0/ε_u$, where u0 is the velocity dispersion. Second, the smallest scale rη is dependent on the nature of dark matter. For collisionless dark matter, $r_η∝(–Gℏ/ε_u)^{1/3}$ ≈ 10 –13 m , where ℏ is the Planck constant. An uncertainty principle for momentum and acceleration fluctuations is also postulated. For self-interacting dark matter, $r_η∝ε^2_uG^{–3}(σ/m)^3$, where σ/m is the cross-section of interaction. On halo scale, the energy cascade leads to an asymptotic density slope γ = –4/3 for fully virialized haloes with a vanishing radial flow, which might explain the nearly universal halo density. Based on the continuity equation, halo density is analytically shown to be closely dependent on the radial flow and mass accretion, such that simulated haloes can have different limiting slopes. A modified Einasto density profile is proposed accordingly.

79 ASTRONOMY AND ASTROPHYSICS↗

Correcting for Selection Biases in the Determination of the Hubble Constant from Time-Delay Cosmography

The time delay between multiple images of strongly lensed quasars has been used to infer the Hubble constant. The primary systematic uncertainty for time-delay cosmography is the mass-sheet transform (MST), which preserves the lensing observables while altering the inferred ⁠H 0 . The TDCOSMO collaboration used velocity dispersion measurements of lensed quasars and lensed galaxies to infer that mass sheets are present, which decrease the inferred H 0 by 8 per cent. Here, we test the assumption that the density profiles of galaxy–galaxy and galaxy–quasar lenses are the same. We use a composite star-plus-dark-matter mass profile for the parent deflector population and model the selection function for galaxy–galaxy and galaxy–quasar lenses. We find that a power-law density profile with an MST is a good approximation to a two-component mass profile around the Einstein radius, but we find that galaxy–galaxy lenses have systematically higher mass-sheet components than galaxy–quasar lenses. For individual systems, λ int correlates with the ratio of the half-light radius and Einstein radius of the lens. By propagating these results through the TDCOSMO hierarchical inference code, we find that H 0 is lowered by a further 3 per cent. Using a more recent measurement of velocity dispersions and our fiducial model for selection biases, we infer H 0 = 66 ± 4 (stat) ± 1 (model sys) ± 2 (measurement sys) km s -1 Mpc -1 for the TDCOSMO plus SLACS data set. The first residual systematic error is due to plausible alternative choices in modelling the selection function, and the second is an estimate of the remaining systematic error in the measurement of velocity dispersions for SLACS lenses. Accurate time-delay cosmography requires precise velocity dispersion measurements and accurate calibration of selection biases.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling the Milky Way – I. Method and first results fitting the thick disc and halo with DES-Y3 data

ABSTRACT We present a technique to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) generated from trilegal models to real data. We apply this technique, which we call mwfitting, to photometric data from the first 3 yr of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud, and the Sagittarius Stream, our main sample spans a total area of ∼2300 deg2. We further explore a smaller subset (∼1300 deg2) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3 per cent. We fit the DES data with an exponential thick disc model and an oblate double power-law halo model. We find that the best-fitting thick disc model has radial and vertical scale heights of 2.67 ± 0.09 kpc and 925 ± 40 pc, respectively. The stellar halo is fit with a broken power-law density profile with an oblateness of 0.75 ± 0.01, an inner index of 1.82 ± 0.08, an outer index of 4.14 ± 0.05, and a break at 18.52 ± 0.27 kpc from the Galactic centre. Several previously discovered stellar overdensities are recovered in the residual stellar density map, showing the reliability of mwfitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict Milky Way star counts for surveys such as the LSST and Euclid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Disentanglement of density and rotation dependences of the field penetration threshold on the J-TEXT tokamak

The field penetration threshold of magnetic perturbations has been observed to vary non-monotonically with an increase of density in ohmic plasmas on the J-TEXT tokamak. This observation appears contradicting the previous empirical density scaling law. Disentanglement of plasma density and rotation dependences of the field penetration threshold has been carried out. It shows that the field penetration threshold depends only weakly on the density but linearly on the plasma rotation. This result is not only important for the prediction of error field tolerance in fusion devices, but also opens a question on the role of density in the forced magnetic reconnection process in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bottom-Up Meets Top-Down: The Crossroads of Multiscale Chromatin Modeling

Chromatin can be viewed as a hierarchically structured fiber that regulates gene expression. It consists of a complex network of DNA and proteins whose characteristic dynamical modes facilitate compaction and rearrangement in the cell nucleus. These modes stem from chromatin's fundamental unit, the nucleosome, and their effects are propagated across length scales. Understanding the effects of nucleosome dynamics on the chromatin fiber, primarily through post-translational modifications that occur on the histones, is of central importance to epigenetics. Within the last decade, imaging and chromosome conformation capture techniques have revealed a number of structural and statistical features of the packaged chromatin fiber at a hitherto unavailable level of resolution. Such experiments have led to increased efforts to develop polymer models that aim to reproduce, explain, and predict the contact probability scaling and density heterogeneity. At nanometer scales, available models have focused on the role of the nucleosome and epigenetic marks on local chromatin structure. At micrometer scales, existing models have sought to explain scaling laws and density heterogeneity. Less work, however, has been done to reconcile these two approaches: bottom-up and top-down models of chromatin. Here, we highlight the multiscale simulation models that are driving toward an understanding of chromatin structure and function, from the nanometer to the micron scale, and we highlight areas of opportunity and some of the prospects for new frameworks that bridge these two scales. Taken together, experimental and modeling advances over the last few years have established a robust platform for the study of chromatin fiber structure and dynamics, which will be of considerable use to the chromatin community in developing an understanding of the interplay between epigenomic regulation and molecular structure.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Measuring the thermal and ionization state of the low- z IGM using likelihood free inference

ABSTRACT We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho/ \bar{\rho })^{\gamma -1}$ and the UV background photoionization rate $\Gamma _{{{{\rm H\, {\small I}}}}{}}$ of the intergalactic medium (IGM) based on the Voigt profile decomposition of the Ly α forest into a set of discrete absorption lines with Doppler parameter b and the neutral hydrogen column density $N_{\rm H\, {\small I}}$. Previous work demonstrated that the shape of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution is sensitive to the IGM thermal parameters T0 and γ, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density dN/dz, and we demonstrate that precise constraints can also be obtained on $\Gamma _{{{{\rm H\, {\small I}}}}{}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b-N_{{{{\rm H\, {\small I}}}}{}}$ distribution on IGM parameters trained on an ensemble of 624 nyx hydrodynamical simulations at z = 0.1, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS data sets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly α forest spectra at z ≃ 0.1, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma _{\log T_0} \sim 0.08$ dex, σγ ∼ 0.06, and $\sigma _{\log \Gamma _{{{{\rm H\, {\small I}}}}{}}} \sim 0.07$ dex).

79 ASTRONOMY AND ASTROPHYSICS↗

3-3-1 Self interacting dark matter and the galaxy core-cusp problem

The core–cusp problem remains a challenging discrepancy between observations and simulations in the standard Λ CDM model for the formation of galaxies. The problem is that Λ CDM simulations predict a steep power-law mass density profile at the center of galactic dark matter halos. However, observations of dwarf galaxies in the Local Group reveal a density profile consistent with a nearly flat distribution of dark matter near the center. A number of solutions to this dilemma have been proposed. In this paper we summarize investigations into the possibility that the dark matter particles themselves self interact and scatter. Such self-interacting dark matter (SIDM) particles can smooth out the dark-matter profile in high-density regions. We also review the theoretical proposal that self-interacting dark matter may arise as an additional Higgs scalar in the 3–3–1 extension of the Standard Model (SM). We present new simulations of galaxy formation and evolution for this formulation of self-interacting dark matter. Current constraints on this self-interacting dark matter are then summarized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Compressive Response and Energy Absorption of Additively Manufactured Elastomers with Varied Simple Cubic Architectures

Additive manufacturing, and particularly the vat photopolymerization process, enables the fabrication of complex geometries at high resolution and small length scales, making it well-suited for fabricating cellular structures (e.g., foams and lattices). Among these, elastomeric cellular structures are of growing interest due to their tunable compliance and energy dissipation. However, comprehensive data on the compressive behavior of these structures remains limited, especially for investigating the structure-property effects from changing the density and distribution of material within the cellular structure. This study explores how the mechanical response of polyurethane-based simple cubic structures changes when varying volume fraction, unit cell length, and unit cell patterning, which have not been systematically investigated previously in additively manufactured elastomers. Increasing volume fraction from 10% to 50% yielded significant changes in compressive stress–strain performance (decreasing strain at 0.5 MPa by 41.6% and increasing energy absorption density by 3962.5%). Although changing the unit cell length between 2.5 and 7 mm in ~30 mm parts did not result in statistically different stress–strain responses, modifying the configuration of struts of different thicknesses across designs with 30% volume fraction altered the stress–strain behavior (differences of 12.5% in strain at 0.5 MPa and 109.4% for energy absorption density). Power law relationships were developed to understand the interactions between volume fraction, unit cell length, and elastic modulus, and experimental data showed strong fits (R 2 > 0.91). These findings enhance the understanding of how multiple structural design aspects influence the performance of elastomeric cellular materials, providing a foundation for informing strategic design of tailorable materials for diverse mechanical applications.

36 MATERIALS SCIENCE↗

Thermal conductivity of two-dimensional disordered fibrous materials defined by interfiber thermal contact conductance and intrinsic conductivity of fibers

A general framework for the theoretical analysis and numerical calculations of the effective thermal conductivity of two-dimensional homogeneous and isotropic disordered fibrous materials is developed in this work based on the model of random contacts between soft-core spherocylinders. The analysis accounts for the interfiber contact conductance and intrinsic conductivity of fibers and is performed in a wide space of governing parameters that includes the fiber aspect ratio, Biot number calculated for a single thermal contact between fibers, and material density ranging from values corresponding to the percolation threshold up to those characteristic of dense fiber networks. For dense networks, exact theoretical equations for the thermal conductivity of materials composed of spherocylinders with an arbitrary aspect ratio and Biot number are derived. The effect of the intrinsic conductivity of fibers on the thermal transport in fibrous materials is found to depend on the density of contacts and can be significant in sufficiently dense fiber networks even if the Biot number for a single thermal contact is small. Semiempirical corrections to the theoretical equations are derived for small and moderate fiber densities. The power law exponent describing the approximate dependence of the conductivity on the density parameter is found to vary from values close to 1 up to values exceeding 2 when evaluated within different finite ranges of the density parameter. This finding explains the variability of scaling laws for thermal conductivity of fibrous materials suggested in the literature based on numerical simulations performed in different regions of the space of material parameters.

36 MATERIALS SCIENCE↗

Stress due to electric charge density distribution in a dielectric slab

The spatial distribution of electric field due to an imposed electric charge density profile in an infinite slab of dielectric material is derived analytically by integrating Gauss’s law. Various charge density distributions are considered, including exponential and power-law forms. Here, the Maxwell stress tensor is used to compute a notional static stress in the material due to the charge density and its electric field. Characteristics of the electric field and stress distributions are computed for example cases in polyethylene, showing that field magnitudes exceeding the dielectric strength would be required in order to achieve a stress exceeding the ultimate tensile strength.

42 ENGINEERING↗