Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial accelerator”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Improving the electrostatic design of the Jefferson Lab 300 kV DC photogun

The 300 kV DC high voltage photogun at Jefferson Lab was redesigned to deliver electron beams with a much higher bunch charge and improved beam properties. The original design provided only a modest longitudinal electric field ( E z ) at the photocathode, which limited the achievable extracted bunch charge. To reach the bunch charge goal of approximately few nC with 75 ps full-width at half-maximum Gaussian laser pulse width, the existing DC high voltage photogun electrodes and anode–cathode gap were modified to increase E z at the photocathode. In addition, the anode aperture was spatially shifted with respect to the beamline longitudinal axis to minimize the beam deflection introduced by the non-symmetric nature of the inverted insulator photogun design. We present the electrostatic design of the original photogun and the modified photogun and beam dynamics simulations that predict vastly improved performance. We also quantify the impact of the photocathode recess on beam quality, where recess describes the actual location of the photocathode inside the photogun cathode electrode relative to the intended location. A photocathode unintentionally recessed/misplaced by sub-millimeter distance can significantly impact the downstream beam size.

47 OTHER INSTRUMENTATION↗
NASA NTRS2014

A High-Resolution X-Ray and Optical Study of SN1006: Asymmetric Expansion and Small-Scale Structure in a Type Ia Supernova Remnant

We introduce a deep (670 ks) X-ray survey of the entire SN 1006 remnant from the Chandra X-Ray Observatory, together with a deep Hα image of SN 1006 from the 4 m Blanco telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of the X-ray proper motions around the entire periphery, carried out over a 9 yr baseline. We find that the expansion velocity varies significantly with azimuth. The highest velocity of approx.7400 km/s (almost 2.5 times that in the northwest (NW)) is found along the southeast (SE) periphery, where both the kinematics and the spectra indicate that most of the X-ray emission stems from ejecta that have been decelerated little, if at all. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Neon emission is strongest in a sharp filament just behind the primary shock along the NWrim, where the pre-shock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in Hα, features that we interpret as ejecta knots that have reached relatively dense regions of the surrounding interstellar medium, but that appear in the interior in projection. Recent three-dimensional hydrodynamic models for Type Ia supernovae display small-scale features that strongly resemble the ones seen in X-rays in SN 1006; an origin in the explosion itself or from subsequent hydrodynamic instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front, as expected from diffusive shock acceleration theory, to numerous regions along both the northeast and southwest rims of the shell. Our data require that a precursor be thinner than about 3, and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of seven or more in a narrow precursor region, promoting diffusive particle acceleration.

galaxies: active↗
DOE OSTI2023

Data-driven predictions of the time remaining until critical global warming thresholds are reached

Leveraging artificial neural networks (ANNs) trained on climate model output, we use the spatial pattern of historical temperature observations to predict the time until critical global warming thresholds are reached. Although no observations are used during the training, validation, or testing, the ANNs accurately predict the timing of historical global warming from maps of historical annual temperature. The central estimate for the 1.5 °C global warming threshold is between 2033 and 2035, including a ±1σ range of 2028 to 2039 in the Intermediate (SSP2-4.5) climate forcing scenario, consistent with previous assessments. However, our data-driven approach also suggests a substantial probability of exceeding the 2 °C threshold even in the Low (SSP1-2.6) climate forcing scenario. While there are limitations to our approach, our results suggest a higher likelihood of reaching 2 °C in the Low scenario than indicated in some previous assessments—though the possibility that 2 °C could be avoided is not ruled out. Explainable AI methods reveal that the ANNs focus on particular geographic regions to predict the time until the global threshold is reached. Our framework provides a unique, data-driven approach for quantifying the signal of climate change in historical observations and for constraining the uncertainty in climate model projections. Given the substantial existing evidence of accelerating risks to natural and human systems at 1.5 °C and 2 °C, our results provide further evidence for high-impact climate change over the next three decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Single-shot measurement of the nonlinear refractive index of air at 92 µm with a picosecond terawatt CO 2 laser

We developed a simple, accurate single-shot method to determine the nonlinear refractive index of air by measuring the evolution of the spatial shape of a laser beam propagating through the atmosphere. A distinctive feature of this new method, which relies on a modified Fresnel propagation model for data analysis, is the use of a hard aperture for producing a welldefined, high-quality beam from a comparatively nonuniform quasi-flat-top beam which is typical for highpeak-power lasers. The nonlinear refractive index of air for a very short (2 ps) long-wave infrared (LWIR) laser pulse was measured for the first time, yielding n 2 = 3.0 × 10 -23 m 2 /W at 9.2 µm. This result is 40% lower than a corresponding measurement with longer (200 ps) LWIR pulses at a similar wavelength.

43 PARTICLE ACCELERATORS↗
NASA NTRS1995

The spatial variation of the infrared-to-radio ratio in spiral galaxies

We have produced two-dimensional maps of the intensity ratio, Q(sub 60), of 60 micron infrared to 20 cm radio continuum emission, for a set of 25 nearby galaxies, mostly spirals. The ratio maps were obtained from infrared images made using IRAS data with the maximum correlation method, and radio images made using VLA data. Before taking the ratio, the radio images were processed so as to have the same resolution properties as the infrared images; the final spatial resolution in all cases is approximately 1 min, corresponding to 1 - 2 kpc for most galaxies. This resolution represents a significant improvement over previous studies. Our new high-resolution maps confirm the slow decrease of Q(sub 60) with increasing radial distance from the nucleus, but show additional structure which is probably associated with separate sites of active star formation in the spiral arms. The maps show Q(sub 60) to be more closely related to infrared surface brightness than to the radial distance r in the galaxy disk. We note also that the Q(sub 60) gradients are absent (or at least reduced) for the edge-on galaxies, a property which can be attributed to the dilution of contrast due to the averaging of the additional structure along the line of sight. The results are all in qualitative agreement with the suggestion that the radio image represents a smeared version of the infrared image, as would be expected on the basis of current models in which the infrared-radio correlation is driven by the formation of massive stars, and the intensity distribution of radio emission is smeared as a result of the propagation of energetic electrons accelerated during the supernova phase.

Marsh, K. A.↗
DOE OSTI2025

Introducing new resonant soft x-ray scattering capability in SSRL

Resonant soft x-ray scattering (RSXS) is a powerful technique for probing both spatial and electronic structures within solid-state systems. Here, we present a newly developed RSXS capability at beamline 13-3 of the Stanford Synchrotron Radiation Lightsource, designed to enhance materials science research. This advanced setup achieves a base sample temperature as low as 9.8 K combined with extensive angular motions (azimuthal ϕ and flipping χ), enabling comprehensive exploration of reciprocal space. Two types of detectors—an Au/GaAsP Schottky photodiode and a charge-coupled device detector with over 95% quantum efficiency—are integrated to effectively capture scattered photons. Extensive testing has confirmed the enhanced functionality of this RSXS setup, including its temperature and angular performance. The versatility and effectiveness of the system have been demonstrated through studies of various materials, including superlattice heterostructures and high-temperature superconductors.

Kuo, Cheng-Tai [SLAC National Accelerator Laborato↗
DOE OSTI2023

Coupled Target-Beam-Moderator Optimization for the Second Target Station

This report describes the results for a coupled target-beam-moderator optimization analysis for the Second Target Station (STS) at ORNL's Spallation Neutron Source. This study is a continuation of the optimization analysis for the moderators in the preliminary design of STS performed in 2022. In the 2022 analysis the dimensions of the moderators are parameterized, while the target and the proton beam profile are kept constant. In this analysis the target height and the proton beam profile are added as parameters. This allows to study the coupled effects of changing target, moderator and beam dimensions. Similar to the 2022 analysis, this work is performed with an automated optimization workflow that uses the optimization toolbox DAKOTA, parameterized geometries in CREO and SpaceClaim, the unstructured mesh generation in Attila4MC, and the particle transport code MCNP6.2©. This workflow enables an efficient optimization using high-fidelity geometries. The main conclusions of this analysis are the following: • Coupled beam-target-moderator optimization provides a few additional percent performance gain over stand-alone moderator optimization. • The moderator performance is not very sensitive to the target height (between ≈60 and ≈80 mm) as long as the beam profile is chosen adequately. • The moderator performance is sensitive to the choice of beam spatial standard deviations, even when the footprint is kept constant. • The optimal moderator radius is the same for a beam footprint of 30 cm 2 , 62.5 cm 2 , and 90 cm 2 . Also the slope of the super-gaussian beam profile does not significantly impact the optimal moderator radius. • The optimal parameters and sensitivities are very similar to the 2022 optimization analysis. These results only indicate a a difference in the optimal radius of the cylindrical moderator, however, this has been corrected in the final design moderator optimization. The main purpose of this report is to document the simulations, results and lessons learned. The most impactful results are summarized in. We also note that the target geometry used in this work is not the final design.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Experimental evaluation of localized creep deformation in grade 91 steel weldments

Spatially resolved measurement of localized creep deformation in heterogeneous creep resistant steel weldments is crucial but challenging for lifetime assessments of critical steam components in power plants. Here, experimental approaches were established to quantitatively evaluate commonly observed localized creep deformation in multi-pass Grade 91 steel weldments. An in-situ digital image correlation (DIC) system was utilized with a creep testing frame to monitor and measure both full-field strain and localized strain accumulation across the weldments during long-term creep testing at elevated temperatures. The in situ DIC method measured not only the creep deformation behavior of the weld metal, heat affected zone (HAZ), and base metal, but also creep strain evolution for each sub-region within the HAZ itself, including coarse-grained HAZ, fine-grained HAZ (FGHAZ), and intercritical HAZ (ICHAZ). The DIC results revealed that local creep strain in the ICHAZ reached up to 90% strain before final rupture, whereas nominal creep strain measured by the standard extensometer of the tested cross-weld specimen was below 10%, indicative of Type IV cracking of the Grade 91 weld. Microstructural analyses revealed that the faster creep degradation/deformation in the HAZ was caused mainly by accelerated matrix grain recrystallization/growth and a reduced pinning effect from the segregated and coarsened precipitates in the FGHAZ and ICHAZ. The ultimate creep rupture occurred in the ICHAZ owing to its lowest creep resistance induced by the largest recrystallized grain size, the lowest fraction of coincidence site lattice, and the lowest local strain energies/dislocation densities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tagging efficiency study of incoherent diffractive vector meson production at the second interaction region at the Electron-Ion Collider

The Electron-Ion Collider (EIC) is an upcoming accelerator facility aimed at exploring the properties of quarks and gluons in nucleons and nuclei, shedding light on their structure and dynamics. The inaugural experimental apparatus, ePIC (electron-Proton and Ion Collider), is designed as a general-purpose detector to address the National Academy of Sciences and the Nuclear Science Advisory Committee physics program at the EIC. The wider EIC community is strongly supporting a second interaction region and an associated second detector to enhance the full science program. In this study, we evaluate how this second interaction region and detector can be complementary to ePIC. The layout of an interaction region for the second detector offers a secondary focus that provides better forward detector acceptance at scattering angles near θ ~ 0 mrad, which can specifically enhance the exclusive, tagging, and diffractive physics program. Here, this article presents an analysis of a tagging program using the second interaction region layout with incoherent diffractive vector meson production. The current design of the second EIC interaction region is evaluated for its vetoing capabilities of incoherent events required for the study of coherent diffractive measurements. We find an increased vetoing performance compared to the ePIC interaction region, thus improving measurements which are important for the spatial imaging of nucleons and nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Cluster in the Auroral Acceleration Region

Due to a fortuitous evolution of the Cluster orbit, the Cluster spacecraft penetrated for the first time in its mission the heart of Earth's auroral acceleration region (AAR) in December 2009 and January 2010. During this time a special AAR campaign was carried out by the various Cluster instrument teams with special support from ESA and NASA facilities. We present some of the first multi-spacecraft observations of the waves, particles and fields made during that campaign. The Cluster spacecraft configuration during these AAR passages was such that it allowed us to explore the differences in the signatures of waves, particles, and fields on the various spacecraft in ways not possible with single spacecraft. For example, one spacecraft was more poleward than the other three (C2), one was at higher altitude (C1), and one of them (0) followed another (C4) through the AAR on approximately the same track but delayed by three minutes. Their separations were generally on the order of a few thousand km or less and occasionally two of them were lying along the same magnetic field line. We will show some of the first analyses of the data obtained during the AAR campaign, where upward and downward current regions, and the waves specifically associated with those regions, as well as the auroral cavities, were observed similarly and differently on the various spacecraft, helping us to explore the spatial, as well as the temporal, aspects of processes occurring in the AAR.

Pickett, Jolene S.↗
DOE OSTI2026

Composition and persistence of soil organic matter along eroding and depositional transects in buried vs. modern soil layers: A case of the Brady paleosol at Wauneta, Nebraska

Paleosols form when soils are buried through deposition by aeolian, colluvial, alluvial or other processes. Burial of former topsoil isolates soil organic matter (SOM) from surface conditions, allowing carbon to accumulate and potentially remain stable for millennia. In this study, SOM composition, distribution, and persistence were analyzed in the Brady Soil of Nebraska, USA to compare SOM spatial variability in modern and buried soils, as well as the impact of erosional exposure on SOM stability. The Brady Soil, formed as a surface soil during the Pleistocene-Holocene transition and now a paleosol buried up to 6 m deep (or more) by loess deposition during the Holocene, was sampled along burial (up to 5.8 m depth) and erosional (up to 1.8 m depth) transects to compare SOM dynamics in different geomorphic settings. Fourier Transform Infrared Spectroscopy (FTIR) and Fourier Transform ion cyclotron resonance mass spectrometry (FTICR-MS) were used to analyze SOM composition, while δ 13 C isotope analyses identified SOM sources and radiocarbon values were used to estimate turnover rates. Results confirmed a vegetation shift from C3 to C4 plants after Brady Soil formation, reflecting warming climatic conditions. Increasing SOM age and decreasing δ 13 C and δ 15 N values with depth indicated slowing of decomposition rate in buried soils. Higher pH in the Brady Soil suggested greater base cation content, supporting SOM stabilization through organo-mineral associations and aggregate formation. However, exposure of the Brady Soil due to surface erosion caused faster SOM turnover. This result suggested susceptibility of buried SOM to losses via decomposition upon erosional exposure, possibly accelerated by priming in response to modern SOM inputs. These findings highlight the potential loss of carbon stocks in buried soils under future climate change, as shifts in soil physicochemical properties may destabilize long-preserved SOM.

Environmental sciences↗
NASA NTRSDate not supplied

Sensorimotor Countermeasures

Exploration class missions will include multiple transitions between gravitational environments, sometimes after long periods in microgravity, which will impact the neurovestibular system and sensorimotor capabilities. New countermeasures and assessment tools are needed to enhance preflight disorientation training, maintain inflight physical performance, and accelerate recovery following transitions between gravity environments so that crew are enabled to perform critical exploration mission tasks. The Sensorimotor Countermeasures capability area encompasses three projects. 1) The Upright Proprioception Retention via Inflight Training and Evaluation (UPRITE) system is being developed as an in-flight countermeasure that will mitigate the degradation of balance control mechanisms, which will improve post-flight postural stability. The development of such a countermeasure is not simply a matter of adapting a ground-based training program for spaceflight. The way that we control balance in a gravitational environment cannot be duplicated in 0g. Therefore, the UPRITE system challenges proprioception and tactile function using a 0g configuration. The steps to develop the countermeasure have been broken into seven phases: training type, controllable board factors, protocol development, 0g hardware test, training efficacy, ISS hardware demonstration, and countermeasure verification. 2) Unobtrusive Monitoring Tools and Operational Assessments (aka Head/Body Assessments) are being developed to quantify crew member sensorimotor adaptation following gravitational transitions. Three areas are needed in this development 1) explore, develop, and validate unobtrusive monitoring tools, 2) define sensorimotor performance metrics, and 3) explore capabilities to enhance sensorimotor assessment tasks for operational use and validate assessment tasks via high-fidelity analogs. 3) The Spatial Disorientation Trainer is being developed as a portable training capability to simulate the performance of landing and recovery type tasks while experiencing vertigo due to post-flight vestibular alterations. These simulations can be used 1) across research studies investigating the impacts of vestibular disruption on operational performance/assessments, 2) to train astronauts for upcoming spaceflight missions, and 3) to educate space medicine and operational personnel prior to providing post-flight support.

Scott Jonathan Wood↗
NASA NTRS2019

Analysis Ready Data in Analytics Optimized Data Stores for Analysis of Big Earth Data in the Cloud

Cloud computing offers the possibility of making the analysis of Big Data approachable for a wider community due to affordable access to computing power, an ecosystem of usable tools for parallel processing, and migration of many large datasets to archives in the cloud, allowing data-proximal computing. Generally, data analysis acceleration in the cloud comes from running multiple nodes in a split-combine-apply strategy. Data systems such as the Earth Observing System Data and Information System are in a position to "pre-split" the data by storing them in a data store that is optimized for data parallel computing, i.e., an Analytics-Optimized Data Store (AODS). A variety of approaches to AODS are possible, from highly scalable databases to scalable filesystems to data formats optimized for cloud access (e.g., zarr and cloud-optimized datasets), with the optimal choice dependent on both the types of analysis and the geospatial structure of the data. A key question is how much preprocessing of the data to do, both before splitting and as the first part of the apply step. Again, the geospatial structure of the data and the analysis type influence the decision, with the added complexity of the user type. Trans-disciplinary users who are not well-versed in the nuances of quality-filtering and georeferencing of remote sensing orbit/swath/scene data tend to ask for more highly processed data, relying on the data provider to make sensible decisions on preprocessing parameters. (This accounts for the popularity of "Level 3" gridded data, despite the lower spatial resolution it provides.) In this case, data can be preprocessed before the split, resulting in higher performance in the rest of the "apply" step, which can be transformative for use cases such as interactive data exploration at scale. Discipline researchers who are experienced with remote sensing data often prefer more flexibility in customizing the preprocessing data into Analysis Ready Data, resulting in more need for on-the-fly preprocessing.

Lynnes, Christopher↗
NASA NTRS2011

The Changing California Coast: The Effect of a Variable Water Budget on Coastal Vegetation Succession

The land-ocean interface along the central coast of California is one of the most diverse biogeographic regions of the state. This area is composed of a species-rich mosaic of coastal grassland, shrubland, and forest vegetation types. An acceleration of conifer encroachment into shrublands and shrub encroachment into grasslands along the coast has been recently documented. These vegetation changes are believed to be driven primarily by fire suppression and changing grazing patterns. Climatic variables such as precipitation, fog, cloud cover, temperature, slope, and elevation also play an important role in vegetation succession. Our study area is located along the central California coast, which is characterized by a precipitation gradient from the relatively wetter and cooler north to the drier and warmer south. Some studies indicate changing fog patterns along this coast, which may greatly impact vegetation. A decrease in water availability could slow succession processes. The primary objective of this project is to determine if vegetation succession rates are changing for the study area and to identify climate and ecosystem variables which contribute to succession, specifically the transition among grassland, shrubland, and forest. To identify vegetation types and rates of succession, we classified two Landsat TM 5 scenes from 1985 to 2010 with a resulting overall accuracy of 82.4%. Vegetation succession was correlated to changes in maximum and minimum temperatures, precipitation, and elevation for each sub-region of the study area. Fog frequency was then compared between the northern and southern regions of the study area for determining the spatial relation between fog frequency and the percent of vegetation change.

Hsu, Wei-Chen↗
DOE OSTI2023

Patterns and rates of soil movement and shallow failures across several small watersheds on the Seward Peninsula, Alaska

Abstract. Thawing permafrost can alter topography, ecosystems, and sediment and carbon fluxes, but predicting landscape evolution of permafrost-influenced watersheds in response to warming and/or hydrological changes remains an unsolved challenge. Sediment flux and slope instability in sloping saturated soils have been commonly predicted from topographic metrics (e.g., slope, drainage area). In addition to topographic factors, cohesion imparted by soil and vegetation and melting ground ice may also control spatial trends in slope stability, but the distribution of ground ice is poorly constrained and hard to predict. To address whether slope stability and surface displacements follow topographic-based predictions, we document recent drivers of permafrost sediment flux present on a landscape in western Alaska that include creep, solifluction, gullying, and catastrophic hillslope failures ranging in size from a few meters to tens of meters, and we find evidence of rapid and substantial landscape change on an annual timescale. We quantify the timing and rate of surface movements using a multi-pronged, multi-scalar dataset including aerial surveys, interannual GPS surveys, synthetic aperture radar interferometry (InSAR), and climate data. Despite clear visual evidence of downslope soil transport of solifluction lobes, we find that the interannual downslope surface displacement of these features does not outpace downslope displacement of soil in locations where lobes are absent (downslope movement means: 7 cm yr−1 for lobes over 2 years vs. 10 cm yr−1 in landscape positions without lobes over 1 year). Annual displacements do not appear related to slope, drainage area, or modeled total solar radiation but are likely related to soil thickness, and volumetric sediment fluxes are high compared to temperate landscapes of comparable bedrock lithology. Time series of InSAR displacements show accelerated movement in late summer, associated with intense rainfall and/or deep thaw. While mapped slope failures do cluster at slope–area thresholds, a simple slope stability model driven with hydraulic conductivities representative of throughflow in mineral and organic soil drastically overpredicts the occurrence of slope failures. This mismatch implies permafrost hillslopes have unaccounted-for cohesion and/or throughflow pathways, perhaps modulated by vegetation, which stabilize slopes against high rainfall. Our results highlight the breadth and complexity of soil transport processes in Arctic landscapes and demonstrate the utility of using a range of synergistic data collection methods to observe multiple scales of landscape change, which can aid in predicting periglacial landscape evolution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Kinematics of the ion tail of comet P/Swift-Tuttle

We have obtained long-slit high resolution spectra of the H2O(+) 6199 A complex in the near tail of comet P/Swift-Tuttle. The observations were made using the Hamilton echelle spectrometer fed by the Lick Observatory 0.6 m coude auxiliary telescope. For most of our observations, the spectral slit was aligned along the Sun-tail axis and the cometary nucleus was placed at one end of the slit, giving us spectra having the spatial and spectral resolution needed to measure the radial velocity and velocity dispersion continuously down the cometary tail out to a distance of 4X10(exp 5) km. The radial velocities confirm the earlier more restricted observations by Rauer & Jockers (1993) and by Wyckoff & Lindholm (1994) showing that the tail motions are indeed bulk flows in the antisolar direction. Out to 3X10(exp 5) km in the tail typical bulk flows are at a speed of approximately 30 km/s. The velocity dispersion, (sigma(sub r)), of the H2O(+) lines follows a pattern that is quite systematic; sigma(sub r) is smallest near the cometary nucleus, and steadily increases down the tail. The highest velocity dispersions are found ahead of the nucleus and off the tail axis. These velocity dispersions are equivalent to ion temperatures ranging from 10(exp 5) to 10(exp 6) K. We note a clear anticorrelation between the H2O(+) line intensities (related to the ion density) and the bulk flow and dispersion velocities; direct mass loading of the solar wind by the observed water ions may be responsible. We discuss several approximate equipartition methods used to infer local magnetic fields induced by the interaction of the cometary ions with the solar wind particle/field stream. Typical fields derived are near 50 nT. The measured tailward accelerations are consistent with this order of magnitude B field.

Spinrad, Hyron↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗