Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Expression systems”

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 OSTI2024

MultiSector Dynamics: 2023 Inaugural Workshop Report

Preface The MultiSector Dynamics (MSD) Community of Practice (CoP) hosted an inaugural workshop on October 3-5, 2023 at the University of California, Davis, to bring together members of the MSD community of practice to advance understanding of the co-evolution of human and natural systems, and to build the next generation of tools that bridge sectors, scales, and systems to realize a more resilient and equitable future. The theme of the workshop was "Advancing Complex Adaptive Human-Earth Systems Science in a World of Interconnected Risks". This document outlines the motivation for the workshop, its goals and objectives, the application process, the agenda, overviews of the training sessions offered to the workshop participants and a summary of each breakout session. The MSD workshop report further discusses the feedback from workshop participants and presents some reflections and next steps. The MSD Workshop organizers thank the DOE Office of Science, Earth and Environmental System Modeling, MultiSector Dynamics program area for financial support of its activities through the Integrated Multisector Multiscale Modeling (IM3) project. For more information related to the broader DOE MultiSector Dynamics Program please see https://climatemodeling.science.energy.gov/program-area/multisector-dynamics. D.L.M. and C.M.B. acknowledge support from the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory (ORNL), managed by UT-Battelle, LLC, for the US Department of Energy (DOE). Disclaimer This report was prepared as an account of work sponsored by an agency of the United States Government. Neither theUnited States Government nor any agency thereof, nor Battelle Memorial Institute, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, complete- ness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial products, process, or service by trade name,trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Pacific Northwest National Laboratory operated by Battelle for the United States Department of Energy Available from:Office of Scientific and Technical Information http://www.OSTI.gov multisectordynamics.org This work is made available under the terms of the Creative Commons Attribution- NonCommercial 4.0 International (CC BY-NC 4.0) https://creativecommons.org/licenses/by-nc/4 Suggested citation: Monier, E., Reed, P.M., Vernon, C.R., Hadjimichael, A., Brelsford, C.M., Burleyson, C.B., Dyreson, A.R., Fletcher, S.M., Giang, A., Gupta, R.S., Jackson, N.D., Jones, A.D., Lamontagne, J.R., McCollum, D.L., Morris, J.F., Moss, R.H., Peng, W., Saari, R.K., Srikrishnan, V., Szinai, J.K., Yoon, J. (2024) MultiSector Dynamics: 2023 Inaugural Workshop Report. MSD-LIVE Data Repository. doi:10.57931/2371710.

Monier, Erwan↗
DOE OSTI2021

Gas transport characteristics of supramolecular networks of metal-coordinated highly branched Poly(ethylene oxide)

Model systems are developed and investigated to better understand the effect of polyether-metal ion interactions on gas separation characteristics. These systems help answer current questions raised by the substantial body of research on metal-organic frameworks (MOFs) dispersed in polyethers to improve gas separation performance, where favorable interactions between the metal centers and polyethers are preferred to improve interfacial compatibility. Specifically, we investigate CO 2 /gas transport properties of supramolecular networks comprising cross-linked poly(ethylene oxide) (XLPEO) and dissociable salts, including LiClO 4 , Ni(BF 4 ) 2 , and Cu(BF 4 ) 2 . Increasing the salt content increases the glass transition temperature (T g ) and generally decreases gas diffusivity and permeability, which can be successfully described using a T g -integrated free volume model with an expression similar to the Vogel-Tammann-Fulcher (VTF) equation. Surprisingly, low loadings of LiClO 4 and Cu(BF 4 ) 2 (2 mass% or less) can increase gas permeability by 30%–70% without affecting the CO 2 /gas selectivity. This increase correlates with polyether-metal ion dynamics as measured by dielectric spectroscopy. Understanding how interaction-mediated dynamics affect gas transport will be instrumental to designing MOF-based mixed matrix materials for gas separations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Frequency Domain Multi-user Processing for Single Carrier Modulation with Cyclic Prefix

Frequency domain (FD) detection and precoding have been shown to be effective means of approaching the theoretical per-user capacity for single carrier modulation SCM) schemes. When a cyclic prefix is added to the SCM waveform, the circulant structures of the resulting convolutional channel matrix allows for relatively simple expressions for the FD processing, and achieves performance near the theoretical capacity even for high signal to noise ratio values. This paper details the signal processing for uplink and downlink multi-user operation with a multi-antenna base station.

5G and Beyond Communication↗
DOE OSTI2021

Nrf2 activation is involved in osteogenic differentiation of periodontal ligament stem cells under cyclic mechanical stretch

Highlights: • Cyclic mechanical stretch increases the nuclear accumulation of Nrf2 in PDLSCs. • T-BHQ promotes osteogenic differentiation under cyclic mechanical stretch. • Nrf2 activation may improve alveolar bone remodeling during orthodontic treatment. During orthodontic treatment, mechanical stretch serves a crucial function in osteogenic differentiation of periodontal ligament stem cells (PDLSCs). Up-regulated reactive oxygen species (ROS) level is a result of cyclic mechanical stretch in many cell types. Nuclear factor erythroid-2-related factor-2 (Nrf2) is a master regulator in various antioxidants expression. However, it is not known whether cyclic mechanical stretch could induce the ROS generation in PDLSCs and whether Nrf2 participated in this process. The present study was aimed to investigate the role of Nrf2 in PDLSCs under cyclic mechanical stretch. Our results showed that cyclic mechanical stretch increased ROS level and the nuclear accumulation of Nrf2 during osteoblast differentiation. Knocking down Nrf2 by siRNA transfection increased ROS formation and suppressed osteogenic differentiation in PDLSCs. T-BHQ, a Nrf2 activator, promoted the osteogenic differentiation in PDLSCs under cyclic mechanical stretch, and improved the microstructure of alveolar bone during orthodontic tooth movement in rats by employing micro-CT system. Taken together, Nrf2 activation was involved in osteogenic differentiation under cyclic mechanical stretch in PDLSCs. T-BHQ could promote the osteogenic differentiation in vitro and in vivo, suggesting a promising option for the remodeling of the alveolar bone during orthodontic tooth movement.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Assembly of ordered DNA-curli fibril complexes during Salmonella biofilm formation correlates with strengths of the type I interferon and autoimmune responses

Deposition of human amyloids is associated with complex human diseases such as Alzheimer’s and Parkinson’s. Amyloid proteins are also produced by bacteria. The bacterial amyloid curli, found in the extracellular matrix of both commensal and pathogenic enteric bacterial biofilms, forms complexes with extracellular DNA, and recognition of these complexes by the host immune system may initiate an autoimmune response. Here, we isolated early intermediate, intermediate, and mature curli fibrils that form throughout the biofilm development and investigated the structural and pathogenic properties of each. Early intermediate aggregates were smaller than intermediate and mature curli fibrils, and circular dichroism, tryptophan, and thioflavin T analyses confirmed the establishment of a beta-sheet secondary structure as the curli conformations matured. Intermediate and mature curli fibrils were more immune stimulatory than early intermediate fibrils in vitro . The intermediate curli was cytotoxic to macrophages independent of Toll-like receptor 2. Mature curli fibrils had the highest DNA content and induced the highest levels of Isg15 expression and TNFα production in macrophages. In mice, mature curli fibrils induced the highest levels of anti-double-stranded DNA autoantibodies. The levels of autoantibodies were higher in autoimmune-prone NZBWxF/1 mice than wild-type C57BL/6 mice. Chronic exposure to all curli forms led to significant histopathological changes and synovial proliferation in the joints of autoimmune-prone mice; mature curli was the most detrimental. In conclusion, curli fibrils, generated during biofilm formation, cause pathogenic autoimmune responses that are stronger when curli complexes contain higher levels of DNA and in mice predisposed to autoimmunity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Qubit Lattice Algorithms based on the Schrodinger-Dirac representation of Maxwell Equations and their Extensions

It is well known that Maxwell equations can be expressed in a unitary Schrodinger-Dirac representation for homogeneous media. However, difficulties arise when considering inhomoge- neous media. A Dyson map points to a unitary field qubit basis, but the standard qubit lattice algorithm of interleaved unitary collision-stream operators must be augmented by some sparse non-unitary potential operators that recover the derivatives on the refractive indices. The effect of the steepness of these derivatives on two dimensional scattering is examined with simulations showing quite complex wavefronts emitted due to transmissions/reflections within the dielectric objects. Maxwell equations are extended to handle dissipation using Kraus operators. Then, our theoretical algorithms are extended to these open quantum systems. A quantum circuit diagram is presented as well as estimates on the required number of quantum gates for implementation on a quantum computer.

Vahala, George↗
DOE OSTI2024

On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics

During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exact hopfion vortices in a 3D Heisenberg ferromagnet

Here, we find exact static soliton solutions for the unit spin vector field of an inhomogeneous, anisotropic three-dimensional Heisenberg ferromagnet. Each soliton is labeled by two integers n and m. It is a (modified) skyrmion in the z = 0 plane with winding number n, which twists out of the plane m times in the z-direction to become a 3D soliton. Here m arises due to the periodic boundary condition at the z-boundaries. We use Whitehead’s integral expression to find that the Hopf invariant of the soliton is an integer H = nm. It represents a hopfion vortex. Plots of the preimages of this topological soliton show that they are either unknots or nontrivial knots, depending on n and m. Any pair of preimage curves links H times, corroborating the interpretation of H as a linking number. We also calculate the exact energy of the hopfion vortex, and show that its topological lower bound has a sublinear dependence on H. Using Derrick’s scaling analysis, we demonstrate that the presence of a spatial inhomogeneity in the anisotropic interaction, which in turn introduces a characteristic length scale in the system, leads to the stability of the hopfion vortex.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Qubit Lattice Algorithms Based on the Schrödinger-Dirac Representation of Maxwell Equations and Their Extensions

It is well known that Maxwell equations can be expressed in a unitary Schrodinger-Dirac representation for homogeneous media. However, difficulties arise when considering inhomogeneous media. A Dyson map points to a unitary field qubit basis, but the standard qubit lattice algorithm of interleaved unitary collision-stream operators must be augmented by some sparse non-unitary potential operators that recover the derivatives on the refractive indices. Here, the effect of the steepness of these derivatives on two-dimensional scattering is examined with simulations showing quite complex wavefronts emitted due to transmissions/reflections within the dielectric objects. Maxwell equations are extended to handle dissipation using Kraus operators. Then, our theoretical algorithms are extended to these open quantum systems. A quantum circuit diagram is presented as well as estimates on the required number of quantum gates for implementation on a quantum computer.

2D electromagnetic scattering↗
DOE OSTI2025

Phonon-Spin Scattering from Unpaired f-electrons in U atoms and It’s Influence on Thermal Transport in Uranium-doped Thorium Dioxide Single Crystals

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO2) is investigated. ThO2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO2 system.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Thermal conductivity suppression in uranium-doped thorium dioxide due to phonon-spin interactions

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO 2 ) is investigated. ThO 2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO 2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO 2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO 2 system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

MODE: A Web Application for Interactive Visualization and Exploration of Omics Data

Studies generating transcriptomics, proteomics, lipidomics, and metabolomics (colloquially referred to as “omics”) data allow researchers to find biomarkers or molecular targets, or understand complex biological structures and functions by identifying changes in biomolecule abundance and expression between experimental conditions. Omics data is multi-dimensional and oftentimes summarization techniques such as principal component analysis (PCA) are used to identify high-level patterns in data. Though useful, these summaries don’t allow exploration of detailed patterns in omics data that may have biological relevance. The use of interactive HTML displays with plots allows researchers to interact with omics data at a detailed level, but building these displays requires significant coding expertise. To overcome this barrier, the software MODE was built to empower users to build their own interactive HTML displays to support scientific discovery. These displays are easily shareable, do not depend on a specific operating system, and allow users to effortlessly sort and filter plots by categorical or numerical variables. MODE allows users to build and share these displays with several options for plot design and meta selection. In conclusion, the MODE web application and its capabilities are presented and then demonstrated on lipidomics data from a leaf wounding study.

lipidomics↗
DOE OSTI2021

Diel Transcriptional Oscillations of a Plastid Antiporter Reflect Increased Resilience of Thalassiosira pseudonana in Elevated CO 2

Acidification of the ocean due to high atmospheric CO 2 levels may increase the resilience of diatoms causing dramatic shifts in abiotic and biotic cycles with lasting implications on marine ecosystems. Here, we report a potential bioindicator of a shift in the resilience of a coastal and centric model diatom Thalassiosira pseudonana under elevated CO 2 . Specifically, we have discovered, through EGFP-tagging, a plastid membrane localized putative Na + (K + )/H + antiporter that is significantly upregulated at >800 ppm CO 2 , with a potentially important role in maintaining pH homeostasis. Notably, transcript abundance of this antiporter gene was relatively low and constant over the diel cycle under contemporary CO 2 conditions. In future acidified oceanic conditions, dramatic oscillation with >10-fold change between nighttime (high) and daytime (low) transcript abundances of the antiporter was associated with increased resilience of T. pseudonana. By analyzing metatranscriptomic data from the Tara Oceans project, we demonstrate that phylogenetically diverse diatoms express homologs of this antiporter across the globe. We propose that the differential between night- and daytime transcript levels of the antiporter could serve as a bioindicator of a shift in the resilience of diatoms in response to high CO 2 conditions in marine environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Creating high-precision reference gas standards of 85Kr for groundwater age-dating

Absolute gas counting (AGC) was applied to two gas blends of 85Kr in argon-methane (P10) counting gas to establish a high-precision specific activity (Bq/cm3) reference value for characterizing 85Kr detection efficiency for groundwater age dating measurements. The AGC or length-compensated technique has been utilized by the metrology community for decades and is an accepted method for developing radioactive gas standards. The AGC capability at Pacific Northwest National Laboratory (PNNL) uses a set of nine unequal-length proportional counters with precisely-measured internal volumes, and a gas loading system with high-precision pressure and temperature sensors. A series of AGC measurements were collected at multiple pressures to determine the inverse pressure relationship (1/P) for 85Kr and define a wall-effect correction that accounts for events decaying into the detector wall and not depositing sufficient energy in the gas to be detected. In addition to the wall-effect, two additional corrections were evaluated and are discussed in detail. Specifically, the threshold effect which accounts for events deposited below the analysis threshold and a detection efficiency as a function of detector volume effect that was observed during analysis. A robust uncertainty model was developed using the Guide to the expression of Uncertainty in Measurements (GUM) approach. The combination of carefully scrutinized correction factors, precise measurements of pressure, temperature and detector volume, and robust counting statistics resulted in the determination of high-precision specific activity values with 0.50% or less total combined uncertainty for two Kr-in-P10 reference gas standards (KP10) that will enable new groundwater age-dating measurements at PNNL.

Kr-85↗
DOE OSTI2020

PcsR2 Is a LuxR-Type Regulator That Is Upregulated on Wheat Roots and Is Unique to Pseudomonas chlororaphis

LuxR solos are common in plant-associated bacteria and increasingly recognized for playing important roles in plant-microbe interkingdom signaling. Unlike the LuxR-type transcriptional regulators of prototype LuxR/LuxI quorum sensing systems, luxR solos do not have a LuxI-type autoinducer synthase gene associated with them. LuxR solos in plant-pathogenic bacteria are important for virulence and in plant endosymbionts contribute to symbiosis. In the present study, we characterized an atypical LuxR solo, PcsR2, in the biological control species Pseudomonas chlororaphis 30–84 that is highly conserved among sequenced P. chlororaphis strains. Unlike most LuxR solos in the plant-associated bacteria characterized to date, pcsR2 is not associated with a proline iminopeptidase gene and the protein has an atypical N-terminal binding domain. We created a pcsR2 deletion mutant and used quantitative RT-PCR to show that the expression of pcsR2 and genes in the operon immediately downstream was upregulated ~10-fold when the wild type strain was grown on wheat roots compared to planktonic culture. PcsR2 was involved in upregulation. Using a GFP transcriptional reporter, we found that expression of pcsR2 responded specifically to root-derived substrates as compared to leaf-derived substrates but not to endogenous AHLs. Compared to the wild type, the mutant was impaired in the ability to utilize root carbon and nitrogen sources in wheat root macerate and to colonize wheat roots. Phenazine production and most biofilm traits previously shown to be correlated with phenazine production also were diminished in the mutant. Gene expression of several of the proteins in the phenazine regulatory network including PhzR, Pip (phenazine inducing protein) and RpeA/RpeB were reduced in the mutant, and overexpression of these genes in trans restored phenazine production in the mutant to wild-type levels, indicating PcsR2 affects the activity of the these regulatory genes upstream of RpeA/RpeB via an undetermined mechanism. Our results indicate PcsR2 upregulates the expression of the adjacent operon in response to unknown wheat root-derived signals and belongs to a novel subfamily of LuxR-type transcriptional regulators found in sequenced P. chlororaphis strains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗