Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximate Newton”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

A Newton-Raphson Method Approach to Adjusting Multi-Source Solar Simulators

NASA Glenn Research Center has been using an in house designed X25 based multi-source solar simulator since 2003. The simulator is set up for triple junction solar cells prior to measurements b y adjusting the three sources to produce the correct short circuit current, lsc, in each of three AM0 calibrated sub-cells. The past practice has been to adjust one source on one sub-cell at a time, iterating until all the sub-cells have the calibrated Isc. The new approach is to create a matrix of measured lsc for small source changes on each sub-cell. A matrix, A, is produced. This is normalized to unit changes in the sources so that Ax(delta)s = (delta)isc. This matrix can now be inverted and used with the known Isc differences from the AM0 calibrated values to indicate changes in the source settings, (delta)s = A ·'x.(delta)isc This approach is still an iterative one, but all sources are changed during each iteration step. It typically takes four to six steps to converge on the calibrated lsc values. Even though the source lamps may degrade over time, the initial matrix evaluation i s not performed each time, since measurement matrix needs to be only approximate. Because an iterative approach is used the method will still continue to be valid. This method may become more important as state-of-the-art solar cell junction responses overlap the sources of the simulator. Also, as the number of cell junctions and sources increase, this method should remain applicable.

Snyder, David B.↗
NASA NTRS2016

First Search for an X-Ray-Optical Reverberation Signal in an Ultraluminous X-Ray Source

Using simultaneous optical (VLT/FORS2) and X-ray (XMM-Newton) data of NGC 5408, we present the first ever attempt to search for a reverberation signal in an ultraluminous X-ray source (NGC 5408 X-1). The idea is similar to active galactic nucleus broad line reverberation mapping where a lag measurement between the X-ray and the optical flux combined with a Keplerian velocity estimate should enable us to weigh the central compact object. We find that although NGC 5408 X-1ʼs X-rays are variable on a timescale of a few hundred seconds (rms of 9.0 +/- 0.5%), the optical emission does not show any statistically significant variations. We set a 3σ upper limit on the rms optical variability of 3.3%. The ratio of the X-ray to the optical variability is an indicator of X-ray reprocessing efficiency. In X-ray binaries, this ratio is roughly 5. Assuming a similar ratio for NGC 5408 X-1, the expected rms optical variability is approximately equal to 2%, which is still a factor of roughly two lower than what was possible with the VLT observations in this study. We find marginal evidence (3 sigma) for optical variability on an approximately 24 hr timescale. Our results demonstrate that such measurements can be made, but photometric conditions, low sky background levels, and longer simultaneous observations will be required to reach optical variability levels similar to those of X-ray binaries.

: accretion↗
NASA NTRS2017

Nano Icy Moons Propellant Harvester

As one of just a few bodies identified in the solar system with a liquid ocean, Europa has become a top priority in the search for life outside of Earth. However, cost estimates for exploring Europa have been prohibitively expensive, with estimates of a NASA Flagship class orbiter and lander approaching $5 billion. ExoTerra's NIMPH offers an affordable solution that can not only land, but return a sample from the surface to Earth. NIMPH combines solar electric propulsion (SEP) technologies being developed for the asteroid redirect mission and microsatellite electronics to reduce the cost of a full sample return mission below $500 million. A key to achieving this order-of-magnitude cost reduction is minimizing the initial mass of the system. The cost of any mission is directly proportional to its mass. By keeping the mission within the constraints of an Atlas V 551 launch vehicle versus an SLS, we can significantly reduce launch costs. To achieve this we reduce the landed mass of the sample return lander, which is the largest multiplier of mission mass, and shrink propellant mass through high efficiency SEP and gravity assists. The NIMPH projects first step in reducing landed mass focuses on development of a micro-In Situ Resource Utilization (micro-ISRU) system. ISRU allows us to minimize landed mass of a sample return mission by converting local ice into propellants. The project reduces the ISRU system to a CubeSat-scale package that weighs just 1.74 kg and consumes just 242 W of power. We estimate that use of this ISRU vs. an identical micro-lander without ISRU reduces fuel mass by 45 kg. As the dry mass of the lander grows for larger missions, these savings scale exponentially. Taking full advantage of the micro-ISRU system requires the development of a micro-liquid oxygen-liquid hydrogen engine. The micro-liquid oxygen-liquid hydrogen engine is tailored for the mission by scaling it to match the scale of the micro-lander and the low gravity of the target moon. We also tailor the engine for a near stoichiometric mixture ratio of 7.5. Most high-performance liquid oxygen-liquid hydrogen engines inject extra liquid hydrogen to lower the average molecular weight of the exhaust, which improves specific impulse. However, this extra liquid hydroden requires additional power and processing time on the surface for the ISRU to create. This increases mission cost, and on missions within high radiation environments such as Europa, increases radiation shielding mass. The resulting engine weighs just 1.36 kg and produces 71.5 newton of thrust at 364 s specific impulse. Finally, the mission reduces landed mass by taking advantage of the SEP modules solar power to beam energy to the surface using a collimated laser. This allows us to replace an 45 kg MMRTG with a 2.5 kg resonant array. By using the combination of ISRU, a liquid oxygen-liquid hydrogen engine, and beamed power, we reduce the initial mass of the lander to just 51.5 kg. When combined with an SEP module to ferry the lander to Europa the initial mission mass is just 6397 kg - low enough to be placed on an Earth escape trajectory using an Atlas V 551 launch vehicle. By comparison, we estimate a duplicate lander using an MMRTG and semi-storable propellants such as liquid oxygen-methane would result in an order of magnitude increase in initial lander mass to 445 kg. Attempting to perform the trajectory with a 450 s liquid oxygen-liquid hydrogen engine would increase initial mass to approximately 135,000 kg. Using an Atlas V 1 U.S. Dollar per kg rate to Earth escape value of $27.7k per kg, just the launch savings are over $3.5 billion.

liquid↗
NASA NTRS2016

The XXL Survey I. Scientific Motivations - Xmm-Newton Observing Plan - Follow-up Observations and Simulation Programme

The quest for the cosmological parameters that describe our universe continues to motivate the scientific community to undertake very large survey initiatives across the electromagnetic spectrum. Over the past two decades, the Chandra and XMM-Newton observatories have supported numerous studies of X-ray-selected clusters of galaxies, active galactic nuclei (AGNs), and the X-ray background. The present paper is the first in a series reporting results of the XXL-XMM survey; it comes at a time when the Planck mission results are being finalized. Aims. We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of approx. 5 x 10(exp 15) erg/s/sq cm in the [0.5-2] keV band (completeness limit). The surveys main goals are to provide constraints on the dark energy equation of state from the space-time-distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programs. Methods. We describe the 542 XMM observations along with the associated multi- and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. Results. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. Conclusions. The XXL multi- data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve asa calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps.

X-rays: general – large-scale structure of Unive↗
NASA NTRS2008

High Power Amplifier and Power Supply

A document discusses the creation of a high-voltage power supply (HVPS) that is able to contain voltages up to -20 kV, keep electrical field strengths to below 200 V/mil (approximately equal to 7.87 kV/mm), and can provide a 200-nanosecond rise/fall time focus modulator swinging between cathode potential of 16.3 kV and -19.3 kV. This HVPS can protect the 95-GHz, pulsed extended interaction klystron (EIK) from arcs/discharges from all sources, including those from within the EIK fs vacuum envelope. This innovation has a multi-winding pulse transformer design, which uses new winding techniques to provide the same delays and rise/fall times (less than 10 nanoseconds) at different potential levels ranging from -20 kV to -16 kV. Another feature involves a high-voltage printed-wiring board that was corona-free at -20 kV DC with a 3- kV AC swing. The corona-free multilayer high-voltage board is used to simulate fields of less than 200 V/mil (approximately equal to 7.87 kV/mm) at 20 kV DC. Drive techniques for the modulator FETs (field-effect transistors) (four to 10 in a series) were created to change states (3,000-V swing) without abrupt steps, while still maintaining required delays and transition times. The packing scheme includes a potting mold to house a ten-stage modulator in the space that, in the past, only housed a four-stage modulator. Problems keeping heat down were solved using aluminum oxide substrate in the high-voltage section to limit temperature rise to less than 10 while withstanding -20 kV DC voltage and remaining corona-free.

Duong, Johnny↗
NASA NTRS2006

The Evolution of the Cluster X-ray Scaling Relations in the Wide Angle ROSAT Pointed Survey Sample at 0.6 < z < 1.0

The X-ray properties of a sample of 11 high-redshift (0.6 < z < 1 .O) clusters observed with Chardm and/or XMM-Newton are used to investigate the evolution of the cluster scaling relations. The observed evolution in the normalization of the L-T, M-T, M(sub 2)-T and M-L relations is consistent with simple self-similar predictions, in which the properties of clusters reflect the properties of the Universe at their redshift of observation. Under the assumption that the model of self-similar evolution is correct and that the local systems formed via a single spherical collapse, the high-redshift L-T relation is consistent with the high-z clusters having virialized at a significantly higher redshift than the local systems. The data are also consistent with the more realistic scenario of clusters forming via the continuous accretion of material. The slope of the L-T relation at high redshift (B = 3.32 +/- 0.37) is consistent with the local relation, and significantly steeper than the self-similar prediction of B = 2. This suggests that the same non-gravitational processes are responsible for steepening the local and high-z relations, possibly occurring universally at z is approximately greater than 1 or in the early stages of the cluster formation, prior to their observation. The properties of the intracluster medium at high redshift are found to be similar to those in the local Universe. The mean surface-brightness profile slope for the sample is Beta = 0.66 +/- 0.05, the mean gas mass fractions within R(sub 2500(z)) and R(200(z)) are 0.069 +/- 0.012 and 0.11 +/- 0.02, respectively, and the mean metallicity of the sample is 0.28 +/- 0.11 Z(sub solar).

Maughan, B. J.↗
NASA NTRS2007

Evidence for a Broad Relativistic Iron Line from the Neutron Star LMXB Ser X-1

We report on an analysis of XMM-Newton data from the neutron star low mass X-ray binary (LMXB) Serpens X-1 (Ser X-1). Spectral analysis of EPIC PN data indicates that the previously known broad iron Ka emission line in this source has a significantly skewed structure with a moderately extended red wing. The asymmetric shape of the line is well described with the laor and diskline models in XSPEC, which strongly supports an inner accretion disk origin of the line. To our knowledge this is the first strong evidence for a relativistic line in a neutron star LMXB. This finding suggests that the broad lines seen in other neutron star LMXBs likely originate from the inner disk as well. Detailed study of such lines opens up a new way to probe neutron star parameters and their strong gravitational fields. The laor model describes the line from Ser X-1 somewhat better than diskline, and suggests that the inner accretion disk radius is less than 6GM/c(exp 2). This is consistent with the weak magnetic fields of LMXBs, and may point towards a high compactness and rapid spin of the neutron star. Finally, the inferred source inclination angle in the approximate range 50-60 deg is consistent with the lack of dipping from Ser X-1.

Bhattacharyya, Sudip↗
NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗
DOE OSTI2024

First characterization of the emission behavior of Mrk 421 from radio to very high-energy gamma rays with simultaneous X-ray polarization measurements

Aims.We have performed the first broadband study of Mrk 421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. Methods.The data were collected as part of an extensive multiwavelength campaign carried out between May and June 2022 using MAGIC,Fermi-LAT,NuSTAR,XMM-Newton,Swift, and several optical and radio telescopes to complement IXPE data. Results.During the IXPE exposures, the measured 0.2–1 TeV flux was close to the quiescent state and ranged from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the very high-energy (VHE) and X-ray emission are positively correlated at a 4σsignificance level. The IXPE measurements reveal an X-ray polarization degree that is a factor of 2–5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, theSwift-XRT monitoring reveals an X-ray flux increase with a clear spectral hardening. This suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation,NuSTARdata reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counterclockwise), implying important changes in the particle acceleration efficiency on approximately hour timescales.

Astronomy & Astrophysics↗
NASA NTRS2004

An Extreme, Blueshifted Iron Line in the Narrow Line Seyfert 1 PG 1402+261

We report on a short, XMM-Newton observation of the radio-quiet Narrow Line Seyfert 1 PG 1402+261. The EPIC X-ray spectrum of PG 1402+261 shows a strong excess of counts between 6 - 9 keV in the rest frame. This feature can be modeled by an unusually strong (equivalent width 2 keV) and very broad energy at 7.3 keV appears blue-shifted with respect to the iron Kalpha emission band between 6.4 - 6.97 keV, whilst the blue-wing of the line extends to 9 keV in the quasar rest frame. The line profile can be fitted by reflection from the inner accretion disk, but an inclination angle of greater than 60 degrees is required to model the extreme blue-wing of the line. Furthermore the extreme strength of the line requires a geometry whereby the hard X-ray emission from PG1402+261 above 2 keV is dominated by the pure-reflection component from the disk, whilst little or none of the direct hard power-law is observed. Alternatively the spectrum above 2 keV may instead be explained by an ionized absorber, if the column density is sufficiently high (NH greater than 3 x 10(exp 23) per square centimeter) and if the matter is ionized enough to produce a deep (tau approximately equal to 1) iron K-shell absorption edge at 9 keV. This absorber could originate in a large column density, high velocity outflow, perhaps similar to those which appear to be observed in several other high accretion rate AGN. Further observations, especially at higher spectral resolution, are required to distinguish between the accretion disk reflection or outflow scenarios.

Reeves, J. N.↗
DOE OSTI2022

Iterative subspace algorithms for finite-temperature solution of Dyson equation

One-particle Green’s functions obtained from the self-consistent solution of the Dyson equation can be employed in the evaluation of spectroscopic and thermodynamic properties for both molecules and solids. However, typical acceleration techniques used in the traditional quantum chemistry self-consistent algorithms cannot be easily deployed for the Green’s function methods because of a non-convex grand potential functional and a non-idempotent density matrix. Moreover, the optimization problem can become more challenging due to the inclusion of correlation effects, changing chemical potential, and fluctuations of the number of particles. In this paper, we study acceleration techniques to target the self-consistent solution of the Dyson equation directly. We use the direct inversion in the iterative subspace (DIIS), the least-squared commutator in the iterative subspace (LCIIS), and the Krylov space accelerated inexact Newton method (KAIN). We observe that the definition of the residual has a significant impact on the convergence of the iterative procedure. Based on the Dyson equation, we generalize the concept of the commutator residual used in DIIS and LCIIS and compare it with the difference residual used in DIIS and KAIN. The commutator residuals outperform the difference residuals for all considered molecular and solid systems within both GW and GF2. For a number of bond-breaking problems, we found that an easily obtained high-temperature solution with effectively suppressed correlations is a very effective starting point for reaching convergence of the problematic low-temperature solutions through a sequential reduction of temperature during calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

A Stellar-mass Black Hole in the Ultra-luminous X-ray Source M82 X-1

We have analyzed the archival XMM-Newton data of the archetypal Ultra-Luminous X-ray Source (ULX) M82 X-1 with an LO5 ksec exposure when the source was in the steady state. Thanks to the high photon statistics from the large effective area and long exposure, we were able to discriminate different X-ray continuum spectral models. Neither the standard accretion disk model (where the radial dependency of the disk effective temperature is T(r) proportional to r(sup -3/4)) nor a power-law model gives a satisfactory fit. In fact, observed curvature of the M82 X-1 spectrum was just between those of the two models. When the exponent of the radial dependence (p in T(r) proportional to r(sup -P)) of the disk temperature is allowed to be free, we obtained p = 0.61 (sup +0.03)(sub -0.02). Such a reduction of p from the standard value 3/4 under extremely high mass accretion rates is predicted from the accretion disk theory as a consequence of the radial energy advection. Thus, the accretion disk in M82 X-1 is considered to be in the Slim disk state, where an optically thick Advection Dominant Accretion Flow (ADAF) is taking place. We have applied a theoretical slim disk spectral model to M82 X-1, and estimated the black hole mass approximately equal to 19 - 32 solar mass. We conclude that M82 X-1 is a stellar black hole which has been produced through evolution of an extremely massive star, shining at a several times the super-Eddington luminosity.

Okajima, Takashi↗
NASA NTRS1982

Application of strong fluctuation random medium theory to scattering of electromagnetic waves from a half-space of dielectric mixture

The strong fluctuation random medium theory is applied to calculate scattering from a half-space of dielectric mixture. The first and second moments of the fields are calculated, respectively, by using the bilocal and the distorted Born approximations, and the low frequency limit is taken. The singularity of the dyadic Green's function is taken into account. Expressions for the effective permittivity for the full space case are derived. It is shown that the derived result of the effect permittivity is identical to that of the Polder and van Santern mixing formula. The correlation function of the random medium is obtained by using simple physical arguments and is expressed in terms of the fractional volumes and particle sizes of the constituents of the mixture. Backscattering coefficients of a half-space dielectric mixture are also calculated. Numerical results of the effective permittivity and backscattering coefficients are illustrated using typical parameters encountered in microwave remote sensing of dry and wet snow. It is also shown that experimental data can be matched with the theory by using physical parameters of the medium as obtained from ground truth measurements.

Tsang, L.↗
NASA NTRS2004

Over-Luminous Elliptical Galaxies

The first paper from our work has been completed and accepted for publication. Another paper presents a study of the ESO 30601 70 galaxy group, combining Chandra, XMM-Newton, and optical observations. We find that the system is a true fossil galaxy group - a group whose optical light is dominated by a single galaxy. The group X-ray emission is composed of a central, dense, cool core (10 kpc in radius) and an isothermal medium beyond the central 10 kpc. The region between 10 and 50 kpc (the cooling radius) has the same temperature as the gas from 50 to 400 kpc, although the gas cooling time between 10 and 50 kpc (2-6 Gyr) is shorter than the Hubble time. Thus, the ESO 3060170 group does not have a group-sized cooling core. We suggest that the group cooling core may have been heated by a central active galactic nucleus (AGN) outburst in the past and that the small, dense, cool core is the truncated relic of a previous cooling core. The Chandra observations also reveal a variety of X-ray features in the central region, including a finger, an edge-like feature, and a small tail, all aligned along a north-south axis, as are the galaxy light and group galaxy distribution. The proposed AGN outburst may cause gas to slosh around the center and produce these asymmetric features. The observed flat temperature profile to 1/3rvir is not consistent with the predicted temperature profile in recent numerical simulations. We compare the entropy profile of the ESO 3060170 group with those of three other groups and find a flatter relation than that predicted by simulations involving only shock heating, S approximately r approximately 0.85. This is direct evidence of the importance of non-gravitational processes in group centers. We derive the mass profiles within 1/3rvir and find that the ESO 3060170 group is the most massive fossil group known.

Forman, William↗
DOE OSTI2021

Contingency Analysis Based on Partitioned and Parallel Holomorphic Embedding

In the steady-state contingency analysis, the traditional Newton-Raphson method suffers from non-convergence issues when solving post-outage power flow problems, which hinders the integrity and accuracy of security assessment. In this paper, we propose a novel robust contingency analysis approach based on holomorphic embedding (HE). Here, the HE-based simulator provides theoretical convergence guarantee, which is desirable because it avoids the influence of numerical issues and provides a credible security assessment conclusion. In addition, based on the multi-area characteristics of real-world power systems, a partitioned HE (PHE) method is proposed with an interfacebased partitioning of HE formulation. The PHE method does not undermine the numerical robustness of HE and significantly reduces the computation burden in large-scale contingency analysis. The PHE method is further enhanced by parallel or distributed computation to become parallel PHE (P2HE). Tests on a 458-bus system, a synthetic 419-bus system and a large-scale 21447-bus system demonstrate the advantages of the proposed methods in robustness and efficiency.

42 ENGINEERING↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Is Ice-Rafted Sediment in a North Pole Marine Record Evidence for Perennial Sea-ice Cover?

Ice-rafted sediments of Eurasian and North American origin are found consistently in the upper part (13 Ma BP to present) of the Arctic Coring Expedition (ACEX) ocean core from the Lomonosov Ridge, near the North Pole (approximately 88 degrees N). Based on modern sea-ice drift trajectories and speeds, this has been taken as evidence of the presence of a perennial sea-ice cover in the Arctic Ocean from the middle Miocene onwards. However, other high latitude land and marine records indicate a long-term trend towards cooling broken by periods of extensive warming suggestive of a seasonally ice-free Arctic between the Miocene and the present. We use a coupled sea-ice slab-ocean model including sediment transport tracers to map the spatial distribution of ice-rafted deposits in the Arctic Ocean. We use 6 hourly wind forcing and surface heat fluxes for two different climates: one with a perennial sea-ice cover similar to that of the present day and one with seasonally ice-free conditions, similar to that simulated in future projections. Model results confirm that in the present-day climate, sea ice takes more than 1 year to transport sediment from all its peripheral seas to the North Pole. However, in a warmer climate, sea-ice speeds are significantly faster (for the same wind forcing) and can deposit sediments of Laptev, East Siberian and perhaps also Beaufort Sea origin at the North Pole. This is primarily because of the fact that sea-ice interactions are much weaker with a thinner ice cover and there is less resistance to drift. We conclude that the presence of ice-rafted sediment of Eurasian and North American origin at the North Pole does not imply a perennial sea-ice cover in the Arctic Ocean, reconciling the ACEX ocean core data with other land and marine records.

Palaeoclimate↗