Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “framing”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Scaffolding civic infrastructures: Examining the role of civic technoscience in public engagements with oil and gas pipelines

Oil and gas pipelines bring attention to the importance of public participation in the management of large-scale infrastructure projects. Participation transforms how communities understand benefits and risks, and can result in safer and more resilient projects. However, participation can be hindered by procedural injustices in planning processes and lack of industry transparency. Based on a survey of 103 civil society groups conducted in 2021, this paper investigates how groups mobilize technical projects and expert knowledge to overcome these. We frame these mobilizations as expressions of civic technoscience, and argue their pivotal role in building capacity for effective engagements with pipelines. To this extent, we furthermore suggest that civic technoscience is fundamental to building civic infrastructures for reforming participatory governance gaps. We conclude with recommendations for how civic technoscience can better support the intentions of civic infrastructuring by coordinating activities across coalitions, as well as by shifting the scales at which groups engage with the governance of pipelines.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Thermal and structural analysis of a new WEST reciprocating probe head with Titanium-Zirconium-Molybdenum armour material

The Tore Supra tokamak was transformed into an X-point divertor fusion device in the frame of the WEST (W-for tungsten-Environment in Steady-state Tokamak) project, launched in support of the ITER tungsten divertor strategy. WEST began operation in 2017. A new reciprocating probe head, equipped with a Titanium-Zirconium-Molybdenum (TZM) heat shield, was recently built and operated on WEST during the C9 Campaign in early 2024. TZM was chosen as armour material for its very good mechanical properties between 700 °C and 1400 °C. This material was already extensively used for probe heads and other plasma facing components on the Alcator C-Mod tokamak. Here, this paper presents the design of the new WEST probe head. It is equipped with tungsten filaments that can be heated to emit electrons and provide a direct measurement of the plasma potential in the scrape-off layer, and it is also equipped with a pair of Langmuir probes. The probe plunges vertically into the plasma, reaches its maximum dive depth, then returns to its protected resting position in about 0.25 s. Thermal loads are calculated using as input the thermal properties of TZM and radial profiles of the heat flux measured in the past by other reciprocating probes, and scaled upwards to the expected heat loads in the WEST scrape-off layer. An analytical estimation of Electromagnetic (EM) loads will also be presented, as well as a thermal and structural analysis of the probe for 15 cycles at maximum plunge depth. The probe is observed by an infrared camera. Comparisons between the heat flux deduced by the Langmuir probe and the infrared radiance of the probe housing during plasma exposure are finally presented.

Diagnostic↗
DOE OSTI2024

Size effect characteristics and influences on fatigue behavior of laser powder bed fusion of thin wall GRCop-42 copper alloy

Size effects, influencing a material's strength, elongation, fatigue limit, and longevity, depend on the operative and dominant deformation and failure mechanisms. This study explores the size effects in additive manufactured (AM) GRCop-42 (Cu-4at%Cr-2at%Nb) thin wall structures fabricated via laser-powder bed fusion (L-PBF) and their impact on fatigue life. The influence of internal defects and surface topography on the fatigue life of specimens in both as-built and hot isostatic pressed (HIP) conditions across different thicknesses is investigated. Where micro-computed tomography (μCT) was used to quantify the internal porosity of as-built, pristine HIP'd, and fatigued HIP'd specimens, and laser microscopy was employed to quantify the surface topography of specimens prior to fatigue. Additionally, quasi-static tests were used to establish baseline mechanical properties (i.e. yield strength (YS), ultimate tensile strength (UTS), and elongation) to frame fatigue testing conditions. Results indicate a significant enhancement in fatigue life for HIP'd specimens for both thicknesses, with internal defects depicting a greater impact than surface topography. Furthermore, fractographic analysis suggests that thicker specimens exhibit higher resistance to crack propagation during fatigue testing in the absence of substantial porosity. Thus, the size effects observed on the fatigue life of L-PBF GRCop-42 appears to be dominated by internal defects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the measurement of shape: With applications to lunar regolith

With the renewed commitment from NASA and other commercial entities for a presence on the Moon, the importance of understanding the characteristics of lunar regolith and how to utilize it have become the target of increasing scrutiny. Much of what is known about lunar regolith was collected during and immediately after the Apollo program, however, analytical techniques and instrumentation have advanced in leaps and bounds in the subsequent decades. Specifically, dynamic image analysis systems have advanced to the point that millions of particles can have morphological characteristics automatically determined in relatively short time frames. Particle morphological data was collected on several lunar samples and a pair of widely used regolith simulants to ascertain the accuracy of these simulants and to explore statistical analysis methods of these large datasets. It is found that these morphology datasets can vary widely depending on the particle size of the particles, and simple averaging of the data skews the results heavily towards the numerically abundant size ranges, the fines. Different reporting methods are suggested to ameliorate these problems. When applied to the lunar regolith, the particles are noted to be less morphologically complex than initially suspected. Compared to the lunar material, the simulants are found to contain some more morphological variability and angular grains. Such difference is likely due to the wildly different comminution processes that these different powder systems are subjected to.

2D shape↗
DOE OSTI2020

Single Stage PLL-less Decoupled Active and Reactive Power Control for Weak Grid Interactive Inverters

This paper presents a single stage phase locked loop-less (PLL-less) active and reactive power (PQ) control for single-phase weak grid interactive inverters. The absence of the PLL requirement in the proposed PQ control enhances the stability margin comparing to conventional current control approaches in both the stationary or the synchronous reference frame. Additionally, the proposed control scheme enables direct PQ control with a single loop control structure. This paper demonstrates that the proposed PLL-less PQ control scheme is asymptotically stable if the controller gains are positive. The necessary conditions for supplying the rated active power into the weak grid are derived. The provided analysis shows that certain amount of reactive power injection is necessary to avoid point of common coupling (PCC) voltage collapse when the rated active power is injected into the weak grid. Several case studies are provided to demonstrate the performance of the proposed PLL-less PQ control scheme under weak grid condition.

Khan, Ahmad↗
DOE OSTI2026

Optimal CO 2 storage management considering safety constraints in multi-stakeholder multi-site GCS projects: A Markov game perspective

Geological carbon storage (GCS) projects could involve a diverse array of stakeholders or players from public, private, and regulatory sectors, each with different objectives and responsibilities. Given the complexity, scale, and long-term nature of GCS operations, determining whether individual stakeholders can independently optimize their interests — or whether collaborative coalition agreements are needed — remains a central question for effective GCS project planning and management. To access large, high-quality storage resources, future GCS deployment may increasingly occur in geologically connected sites, where shared geological features such as pressure space and reservoir pore capacity can lead to competitive behavior among stakeholders. In this work, we propose a paradigm based on Markov games to quantitatively investigate how different coalition structures affect the goals of stakeholders. We frame this multi-stakeholder multi-site problem as a multi-agent reinforcement learning problem with safety constraints. Our approach enables agents to learn optimal strategies while complying with safety regulations. We present an example where multiple operators are injecting CO 2 into their respective project areas in a geologically connected basin. To address the high computational cost of repeated simulations of high fidelity models, a previously developed surrogate model based on the Embed-to-Control (E2C) framework is employed. Our results demonstrate the effectiveness of the proposed framework in addressing optimal management of CO 2 storage when multiple stakeholders with different objectives and goals are involved.

58 GEOSCIENCES↗
DOE OSTI2026

Evolution of storage monitoring – update in response to commercial and regulatory drivers

Carbon Capture and Storage (CCS) is in transition from first-of-a kind projects and research-orientated pilots to commercially-motivated applications. Monitoring results from many newly developed and planned large scale commercial projects are limited; however, it is worthwhile to assess their evolution and consider new strategies as part of an effort to assess and document best practices. Commercial monitoring is targeted to activities that comply with regulatory drivers and de-risk investments. Commercial monitoring also supports accounting that storage has occurred and is tied to project financing. It deals with long time frames and large volumes injected into multiple wells and multiple projects in favorable areas. We see developing trends toward reproducible workflows that systematically reduce risks and clarify expectations for oversight and long-term surveillance. Monitoring techniques showing increasing trends include injection zone pressure as a history-matching and compliance tool. To reduce cost and environmental impact of time-lapse seismic data collection, deploying new approaches and tools, such as use of fibre and installed sources are increasingly applied. Concern over the risk of induced seismicity by regulatory bodies and the general public has increased, which has also resulted in increased monitoring. Some techniques used in the early research phases have been sidelined or used only in restricted applications. For example, geochemical analyses in the injection zone as well as the environment are now being deployed less than it was in research-oriented programs, except in the US where it is required by the permitting process. Expectations of frequent area-wide near surface monitoring have also decreased.

25 ENERGY STORAGE↗
DOE OSTI2022

Flow structures of the cross-flow over a five-layer helically coiled steam generator geometry

The helically-coiled heat exchanger (HCSG) offers advantages over straight tubes such as compactness in geometry, increased heat transfer coefficients, and capability to absorb thermal expansion. In this study, the time-resolved velocity field data of shell-side crossflow over a five-layer helically coiled steam generator geometry was obtained using Particle Image Velocimetry (PIV) between adjacent rods in four regions with up to 5000 frames per second (5 kHz) at approximately Re u = 3600. The slant rod bundle created different flow patterns in different regions. Proper Orthogonal Decomposition (POD) revealed high-energy-mode flow structures. Due to the confined flow channel geometry, the large flow structures were not advected resulting in no strong POD mode paring. Continuous Wavelet Transform (CWT) visualized the characteristics of flow fluctuation in time and frequency domain simultaneously, which ranged in Strouhal number, S u , between 0.03 and 1.17. In the wake regions of the rods, flows changed in both pattern and magnitude over time. Once a pattern formed, it continued for relatively shorter than straight tube arrays, before changing the pattern again. The histograms of the POD time coefficients presented this multimodal flow characteristics. The current study focuses more on whole flow field analyses in the selected regions rather than local point-wise analyses. Finally, this study revealed that various Strouhal numbers can exist in the helically coiled heat exchanger geometry, which was not observed in straight tube arrays. In particular, Strouhal numbers smaller than 0.1 do not seem to be caused by vortex shedding but seem to be the result of the combination of complex geometry and possible multimodal trends. Based on the current observation, further works would involve intensive instantaneous point-wise analyses on the interactions between important locations such as separation points, stagnation points, shear layers, etc.

42 ENGINEERING↗
DOE OSTI2026

Local convection characteristics of inline arrangement of Kagome-shaped unit cells in a square duct

Transient liquid crystal thermography experiments have been conducted to determine detailed convective heat transfer coefficients at the endwalls of lattice-frame configurations based on Kagome-shaped unit cells. Kagome unit cells with porosity of 0.88 were arranged in an inline manner, where two such arrangements were studied. In the first arrangement, a total of ten unit cells were placed next to each other along the streamwise direction resulting in a continuous configuration. In the 2nd arrangement, alternate unit cells from the continuous configuration were dropped, resulting in a discrete configuration which featured a total of 5 Kagome unit cells. Due to the asymmetric nature of the strut connections within the unit cell, the convective heat transfer coefficient maps were determined for the two opposite walls where the struts meet the endwalls. Transient liquid crystal experiments were conducted for Reynolds number ranging between 10,000 and 30,000. Here, the study was focused on the developing nature of flow along the streamwise direction and the local convection characteristics for continuous and discrete arrangement of unit cells. For the continuous configuration, the Nusselt number ratios (N u /Nu 0 ) varied between 3.15–3.46 and 2.75–2.89 for the two walls A and B, respectively. For the discrete configuration, convective heat transfer coefficients varied between 2.56–2.71 and 2.09–2.40 for the two walls. Kagome unit cells have the potential to be fabricated through inexpensive manufacturing routes such as wire-woven method and these unit cells find their applications in the areas which require different heat transfer levels on opposite walls.

Kagome↗
DOE OSTI2022

A physics-informed multi-agents model to predict thermo-oxidative/hydrolytic aging of elastomers

This paper introduces a novel physics-informed multi-agents constitutive model to propose prediction in quasi-static constitutive behavior of cross-linked elastomer and the loss of mechanical performance during environmental aging. The presented model is used to simulate the effect of single-mechanism chemical aging (i.e. thermal-inducedor hydrolytic aging) on the behavior of the material in this hybrid framework. Those environmental single-mechanism damages change the polymer matrix over time due to massive chain scission, chain formations, and changing the arrangement of molecules in the polymer matrix. Here we propose a data-driven super-constrained machine-learned engine to represent damage in the polymer matrix and capture the changes in material behavior, including its inelastic features such as Mullins effect and permanent set in the course of aging. We have simplified the 3D stress–strain tensor mapping problem into a small number of super-constrained 1D mapping problems by means of a sequential order reduction. An assembly of multiple replicated conditional neural-network learning-agents (L-agents) is trained to systematically simplify the high-dimensional mapping problem into multiple 1D problems, each represented by a different type of agent. Our hybrid framework is designed to capture the effect of deformation history, aging time, and aging temperature. The model is validated with respect to a comprehensive set of experiments specifically designed to benchmark model capabilities and also against available data in the literature. Thermodynamic consistency and frame independency have been verified. Besides acceptable predictive abilities, a significant reduction of computational cost to predict behavior at multiple states of deformation is the most significant feature of this model.

42 ENGINEERING↗
DOE OSTI2025

Ultra-high temperature testing and performance of L-PBF C103

Additive Manufacturing (AM) of refractory alloys is gaining traction as a materials processing route for components subject to extreme temperature environments. Due to the low oxidation resistance of refractory alloys, novel methods for evaluating their elevated temperature performance must be developed. In this work, a Gleeble® 3800 thermomechanical load frame was modified to evaluate the mechanical properties of laser powder bed fusion (L-PBF) consolidated niobium alloy C103 ranging from room temperature (RT) to 1400 °C. The fixturing and sample geometry were designed to accommodate Joule heating and prevent damage to the test chamber. Oxidation of the samples was minimized via testing in vacuum level of 1E-5 Torr. Ultimate tensile strength (UTS), yield strength (YS), elongation, and strain-hardening behavior were determined as a function of temperature. L-PBF C103 presented an average UTS of ∼650 MPa and over 25 % elongation at RT. Above RT, the UTS and YS dropped then leveled off from 500 °C to 1000 °C with values ranging from ∼400 MPa to ∼460 MPa, which is consistent with dynamic strain aging observed in this class of alloys. The strength rapidly declined after 1200 °C to ∼150 MPa at 1400 °C. Fractography indicated ductile fractures for the C103 at all test temperatures, and Electron Backscatter Diffraction (EBSD) analysis revealed a textured microstructure and the presence of dynamic recrystallization within the necked region of the sample tested at 1400 °C.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Profiling B cell immunodominance after SARS-CoV-2infection reveals antibody evolution to non-neutralizing viral targets

Dissecting the evolution of memory B cells (MBCs) against SARS-CoV-2 is critical for understanding antibody recall upon secondary exposure. Here, we used single-cell sequencing to profile SARS-CoV-2-reactive B cells in 38 COVID-19 patients. Using oligo-tagged antigen baits, we isolated B cells specific to the SARS-CoV-2 spike, nucleoprotein (NP), open reading frame 8 (ORF8), and endemic human coronavirus (HCoV) spike proteins. SARS-CoV-2 spike-specific cells were enriched in the memory compartment of acutely infected and convalescent patients several months post symptom onset. With severe acute infection, substantial populations of endemic HCoV-reactive antibody-secreting cells were identified and possessed highly mutated variable genes, signifying preexisting immunity. Finally, MBCs exhibited pronounced maturation to NP and ORF8 over time, especially in older patients. Monoclonal antibodies against these targets were non-neutralizing and non-protective in vivo. These findings reveal antibody adaptation to non-neutralizing intracellular antigens during infection, emphasizing the importance of vaccination for inducing neutralizing spike-specific MBCs.

antibody↗
DOE OSTI2020

Critical evaluation and thermodynamic modeling of the Pd–Sn system

The single crystalline material PdSn 4 , a homologue of the Dirac nodal arc semimetal PtSn 4 , is a promising candidate to search for new topological states with fascinating quantum physical properties. In the present work, the thermochemistry and phase diagram information of the Pd–Sn system published in the literature was collected and critically reviewed, two thermodynamic reassessments of the Pd–Sn system were carried out in the frame of the CALPHAD approach to optimize the Gibbs free energy of each phase presented in the system. The liquid phase was described using respectively the Bragg-Williams random mixing model and the associate solution model. All the reliable thermodynamic property and phase equilibria data can be generally described using the presently obtained thermodynamic descriptions of the Pd–Sn system. Issues related to the further improvement of the thermodynamic descriptions of the Pd–Sn system are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cross-sectoral impact of emerging technologies on U.S. manufacturing resilience and competitiveness

Introducing new technologies in one energy-intensive industry can affect how other industries operate and stay resilient, yet these cross-sector interactions are often underappreciated in conventional technology roadmaps. In practice, industrial systems do not evolve in isolation. They are linked through shared upstream and downstream dependencies, such as electricity and fuel supply, critical materials, transportation networks, and enabling infrastructure. As a result, large-scale technology deployment in one sector can reshape resource availability, infrastructure demand, and operational risk in others. These interdependencies mean that technology deployment decisions in one sector can create unintended bottlenecks or cascading benefits in others. Here, this article argues that a cross-sector, system-of-systems perspective is essential for evaluating and scaling emerging technologies in energy-intensive industries. By framing industrial transformation as an interconnected systems challenge rather than a set of isolated sectoral decisions, the study highlights how interdependence shapes technology feasibility, adoption pathways, and resilience outcomes. The article illustrates how cross-sector linkages can amplify both risks and benefits, and it emphasizes the importance of integrated planning approaches that account for shared dependencies, cascading impacts, and co-optimization opportunities. Adopting this broader perspective can support more robust technology roadmaps, improve strategic coordination across industries, and strengthen the long-term resilience of the industrial sector as a whole.

Nain, Preeti [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2025

Establishing a silica gel zone in well annulus and evaluating its performance in blocking vertical water flow

Wells are often constructed for monitoring purposes with relatively long screen lengths (e.g., >10 m). Vertical water flows can occur within the artificial or natural filterpack annulus that surrounds the screened interval, bypassing packer assemblies installed inside the wellbore. Attempts to isolate discrete vertical zones during groundwater sampling are unsuccessful when annular flow occurs and lead to remedy decisions based on biased or incorrect interpretations. Blocking vertical annular water flow and contaminant transport will help obtain more accurate concentrations of contaminants from sampling in targeted depth intervals. The application of silica gels formed from the injected colloidal silica CS suspensions is a novel approach to minimize or prevent movement of vertical movement of groundwater in the surrounding filterpack annulus. In this work, we tested the feasibility of injecting CS suspensions to target locations and developed a modified CS formulation that is injectable and prevents gravity sinking. We studied the distribution and penetration of silica gel at laboratory scale in mock well annulus with surrounding formations. We evaluated the performance of the silica gel in blocking vertical water flow in the annulus and in minimizing chemical transport through the gel zone. CS suspension formulations have been defined that are ready for injection, stay in target locations, and form gel within desired time frames. Injection of CS suspensions achieved uniform distribution in a well annulus filter pack, fully occupied the annulus pore space, and penetrated the formation surrounding the filter packer with a sufficient distance to create a hydraulic annular seal when the injection was applied at a sufficient rate. The depth of penetration into the formation was dependent on the permeability contrast between the filter pack and the surrounding formation. Silica gel that formed in the annulus blocked vertical water flow and stopped the chemical transport through the gel zone. In conclusion, this research reveals that using CS suspension injection and sequential gelation (CS-GEL) is a promising technology for blocking vertical water flow and chemical transport through the filter pack in targeted zones within the annulus of long-screened well systems.

Colloidal silica suspension↗
DOE OSTI2020

Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form

We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Nonlinear sparse Bayesian learning for physics-based models

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗