Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “small cell problem”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Synthetic Pathways to gamma-Graphyne and Related Allotropes of Carbon

Graphynes, two-dimensional carbon lattices combining sp 1 and sp 2 hybridized atoms, were predicted theoretically more than three decades ago, but few structures have been realized to date. These carbons are believed to possess remarkable mechanical and electronic properties, including high charge carrier mobilities comparable to those in graphene (10 4 to 10 5 cm 2 V -1 s -1 ). Unlike graphene, certain graphynes are predicted to be intrinsic semiconductors. Among these intriguing structures, γ-graphyne stands out as the structurally simplest and most symmetric sp 1 /sp 2 lattice. γ-Graphyne was first theorized in 1987. In contrast with graphene, γ-graphyne will be a semiconductor with a small band gap suitable for fabrication of electronic devices. This solves one of the fundamental problems of carbon-based electronics, the necessity for inducing a band gap in graphene. γ-Graphyne has the potential to form the basis for the next generation of carbon-based electronics operating at speeds unattainable by traditional silicon chips. Unlike silicon, γ-graphyne is a direct band gap semiconductor, and it will feature exceptional strength comparable to that of other 2D carbon allotropes. Such combination of properties may enable a new generation of highly efficient, ultra-light and flexible solar cells. Despite being a potentially “magical” material, γ-graphyne remained synthetically elusive for over three decades. The primary goals of this project were: (1) Synthesis of bulk γ-graphyne phases through solution-phase 2D polymerizations; (2) Experimental exploration of the physical and chemical properties of γ-graphyne; and (3) Mechanistic and theoretical studies of the novel chemical transformations developed in Goal 1. Common pyrolytic and vapor-deposition methodologies used for the synthesis of graphitic allotropes are unsuitable for graphyne and other sp 1 -contaning structures, as acetylenes readily convert to graphene and amorphous carbon at high temperatures. In contrast, this proposal is based on solution-based 2D polymerization. The major advantages of this approach over the traditional high temperature techniques are the potential to adjust the structure of the material with atomic precision, and the possibility of using structurally complex and relatively fragile repeat units. The outcomes of this research can revolutionize carbon nanotechnology, expanding the field’s structural toolbox beyond primarily graphitic and benzenoid structures. Understanding the chemistry of sp 1 carbon allotropes can lead to entirely new classes of structures with unique properties, including graphyne ribbons, nanotubes, quantum dots, and heterostructures with other 2D materials. Furthermore, the development of reliable and robust synthetic pathways towards periodic covalent molecular sheets with atomically precise structures shall have a profound impact on chemistry and materials science.

2D polymerization↗
DOE OSTI2021

Decomposing Loosely Coupled Mixed-Integer Programs for Optimal Microgrid Design

Microgrids are frequently employed in remote regions, in part because access to a larger electric grid is impossible, difficult, or compromises reliability and independence. Although small microgrids often employ spot generation, in which a diesel generator is attached directly to a load, microgrids that combine these individual loads and augment generators with photovoltaic cells and batteries as a distributed energy system are emerging as a safer, less costly alternative. In this work, we present a model that seeks the minimum-cost microgrid design and ideal dispatched power to support a small remote site for one year with hourly fidelity under a detailed battery model; this mixed-integer nonlinear program (MINLP) is intractable with commercial solvers but loosely coupled with respect to time. A mixed-integer linear program (MIP) approximates the model, and a partitioning scheme linearizes the bilinear terms. We introduce a novel policy for loosely coupled MIPs in which the system reverts to equivalent conditions at regular time intervals; this separates the problem into subproblems that we solve in parallel. We obtain solutions within 5% of optimality in at most six minutes across 14 MIP instances from the literature and solutions within 5% of optimality to the MINLP instances within 20 minutes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length

A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Neutrino fast flavor instability in three dimensions for a neutron star merger

The flavor evolution of neutrinos in core collapse supernovae and neutron star mergers is a critically important unsolved problem in astrophysics. Following the electron flavor evolution of the neutrino system is essential for calculating the thermodynamics of compact objects as well as the chemical elements they produce. Accurately accounting for flavor transformation in these environments is challenging for a number of reasons, including the large number of neutrinos involved, the small spatial scale of the oscillation, and the nonlinearity of the system. We take a step in addressing these issues by presenting a method which describes the neutrino fields in terms of angular moments. We apply our moment method to neutron star merger conditions and show it simulates fast flavor neutrino transformation in a region where this phenomenon is expected to occur. By comparing with particle-in-cell calculations we show that the moment method is able to capture the three phases of growth, saturation, and decoherence, and correctly predicts the lengthscale of the fastest growing fluctuations in the neutrino field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Second Target Station Project (CHESS Technical Report)

CHESS is a direct geometry neutron spectrometer designed to detect and analyze weak signals intrinsic to small cross-sections (e.g., small mass, small magnetic moments or neutron absorbing materials). This instrument is optimized to enable unprecedented characterization of quantum materials, spin liquids, thermoelectric and battery materials, liquids, and soft matter. The ability to simultaneously measure dynamic processes over a wide energy range for very small samples will make CHESS the spectrometer of choice for the initial exploration of new materials. The broad dynamic range will also be well matched to measurements of relaxation processes and excitations in soft and biological matter. The 15 Hz repetition rate of STS enables use of multiple incident energies within a single source pulse, greatly expanding the information gained in a single experiment. An essential feature of CHESS is the capability for polarization analysis to separate nuclear from magnetic scattering or coherent from incoherent scattering in hydrogenous materials, and better understanding spin-anisotropic correlations. This instrument will employ advanced sample environments such as high-pressure cells, dilution refrigerators, high field cryo-magnets and polarization devices, as well as combinations of these, to solve problems at the forefront of materials research. CHESS will be one of the flagship spectrometers of the Second Target Station (STS), providing world leading capabilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Small-scale signatures of primordial non-Gaussianity in k-nearest neighbour cumulative distribution functions

ABSTRACT Searches for primordial non-Gaussianity in cosmological perturbations are a key means of revealing novel primordial physics. However, robustly extracting signatures of primordial non-Gaussianity from non-linear scales of the late-time Universe is an open problem. In this paper, we apply k-Nearest Neighbour cumulative distribution functions, kNN-CDFs, to the quijote-png simulations to explore the sensitivity of kNN-CDFs to primordial non-Gaussianity. An interesting result is that for halo samples with $M_\mathrm{ h}\langle 10^{14}$ M$_\odot$ $h^{-1}$, the kNN-CDFs respond to equilateral PNG in a manner distinct from the other parameters. This persists in the galaxy catalogues in redshift space and can be differentiated from the impact of galaxy modelling, at least within the halo occupation distribution (HOD) framework considered here. kNN-CDFs are related to counts-in-cells and, through mapping a subset of the kNN-CDF measurements into the count-in-cells picture, we show that our results can be modelled analytically. A caveat of the analysis is that we only consider the HOD framework, including assembly bias. It will be interesting to validate these results with other techniques for modelling the galaxy–halo connection, e.g. (hybrid) effective field theory or semi-analytical methods.

Coulton, William R. (ORCID:0000000212973673)↗
DOE OSTI2024

MedNAS: Multiscale Training-Free Neural Architecture Search for Medical Image Analysis

Deep neural networks have demonstrated impressive results in medical image analysis, but designing suitable architectures for each specific task is expertise dependent and time consuming. Neural architecture search (NAS) offers an effective means of discovering architectures. It has been highly successful in numerous applications, particularly in natural image classification. Yet, medical images possess unique characteristics, such as small regions and a wide variety of lesion sizes, that differentiate them from natural images. Furthermore, most current NAS methods struggle with high computational costs, especially when dealing with high-resolution image datasets. In this article, we present a novel evolutionary NAS method called multiscale training-free neural architecture search (MSTF-NAS) to address these challenges. Specifically, to accommodate the broad range of lesion region sizes in disease diagnosis, we develop a new reduction cell search space that enables the search algorithm to explicitly identify the optimal scale combination for multiscale feature extraction. Further, to overcome the issue of high computational costs, we utilize training-free indicators as performance measures for candidate architectures, which allows us to search for the optimal architecture more efficiently. More specifically, by considering the capability and simplicity of various networks, we formulate a multiobjective optimization problem that involves two training-free indicators and model complexity for candidate architectures. Extensive experiments on a large medical image benchmark and a publicly available breast cancer detection dataset are conducted. The empirical results demonstrate that our MSTF-NAS outperforms both human-designed architectures and current state-of-the-art NAS algorithms on both datasets, indicating the effectiveness of our proposed method.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Predicting Wind Loading and Instability in Solar Tracking PV Arrays

Wind loading and the fluctuating pressure loads it creates on PV panel surfaces are associated with multiple degradation mechanisms and failures. Modest wind speeds create reversing loads that can initiate cell cracks and weather cracked cells. Stronger wind speeds and extreme weather events can lead to larger scale forces and the aerodynamic instability known as torsional galloping. All these effects are dependent on the complex coupling between wind speed, panel orientation, and a myriad of other hardware and site-specific factors. In this work, we present the latest developments from our work to build an open-source, high-performance computing (HPC) fluid dynamics solver to predict and mitigate these effects. This simulation package allows users to easily specify different array layouts, solar-tracking angles, panel geometries, and weather conditions before automatically generating a refined computational mesh and solving for the unsteady loading on each panel surface. Small domains (e.g., a single panel row in isolation) can be solved on a modern laptop, while larger domains or very high-fidelity studies can be solved on distributed or HPC resources with minimal modifications to the underlying problem specification. We present preliminary case studies obtained using this simulation package and highlight how increased wind speeds combined with sub-optimal tracking angles can exacerbate degradation drivers.

aerodynamics↗
DOE OSTI2021

Current status and future outlook for patient‐derived cancer models from a rare cancer research perspective

Abstract Rare cancers are a group of approximately 200 malignancies with extremely low incidences and with a wide variety of genotypes and phenotypes. Collectively, they are more common than any single malignancy. However, given the small numbers of individuals diagnosed with rare cancers, it is difficult to design clinical trials with sufficient patient numbers. Therefore, few effective anticancer drugs have been developed, and evidence‐based medicine is not always feasible for rare cancers. Consequently, their clinical outcomes are generally poorer. Cancer research requires adequate models that faithfully recapitulate molecular features and reproduce treatment responses of the original tumors. Such models allow us to focus on more efficacious drugs in the clinical studies. For rare cancers, patient‐derived cancer models are particularly important because the enrollment of sufficient patients is rarely attainable within a reasonable period of time. However, extremely few models are available for rare cancers. For example, cell lines and xenografts are available for only a limited number of histological subtypes of sarcomas; therefore, most sarcoma research is performed without such models, and a lack of adequate cancer models causes a lag in therapeutic development. The establishment of novel rare cancer models will dramatically facilitate rare cancer research and treatment development in the near future. This review focuses on the status of patient‐derived rare cancer models and discusses their pivotal problems and possibilities, using sarcomas as a representative rare cancer type. Multi‐institutional collaboration will help address the scarcity of patient‐derived rare cancer models.

Kondo, Tadashi↗
DOE OSTI2022

A simple technique to classify diffraction data from dynamic proteins according to individual polymorphs

One often observes small but measurable differences in the diffraction data measured from different crystals of a single protein. These differences might reflect structural differences in the protein and may reveal the natural dynamism of the molecule in solution. Partitioning these mixed-state data into single-state clusters is a critical step that could extract information about the dynamic behavior of proteins from hundreds or thousands of single-crystal data sets. Mixed-state data can be obtained deliberately (through intentional perturbation) or inadvertently (while attempting to measure highly redundant single-crystal data). To the extent that different states adopt different molecular structures, one expects to observe differences in the crystals; each of the polystates will create a polymorph of the crystals. After mixed-state diffraction data have been measured, deliberately or inadvertently, the challenge is to sort the data into clusters that may represent relevant biological polystates. Here, this problem is addressed using a simple multi-factor clustering approach that classifies each data set using independent observables, thereby assigning each data set to the correct location in conformational space. This procedure is illustrated using two independent observables, unit-cell parameters and intensities, to cluster mixed-state data from chymotrypsinogen (ChTg) crystals. It is observed that the data populate an arc of the reaction trajectory as ChTg is converted into chymotrypsin.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Statistical Learning for Nonlinear Model Reduction from Local Simulations of Stochastic and Particle- and Agent-Based Systems

Stochastic physical systems across the sciences that have very high-dimensional state spaces, with a large number of fast degrees of freedom that force direct simulators to proceed by integration steps that are orders of magnitude smaller than events of interests (e.g., particle collisions). Examples range from molecular motion to dynamics of large populations of cells. A grand challenge in the simulation and understanding of such systems is the systematic construction of accurate, interpretable, reduced models, enabling faster simulations, revealing fundamental properties of the dynamics, and predicting phenomena of interest that the original simulator could not reached with sufficient accuracy or within a given computational budget. In this projected we developed novel statistical estimation/machine learning techniques for analyzing and building empirical reduced models for important families of high-dimensional stochastic systems, in particular: - we developed techniques for estimating interaction kernels in interacting particle- and agent-based systems, which are ubiquitous in Physics, Biology and many other sciences, given observed trajectories of the system; - we developed techniques for nonlinear model reduction for high-dimensional stochastic systems that have a small number of unknown, nonlinear slow variables, and a large number of fast modes, that are possibly of large magnitude, given observed short trajectories of the system in the form of bursts of trajectories from different initial conditions; - we developed novel techniques for estimating linear dynamical systems on graphs when both the dynamics and the underlying graph are unknown, and we have a sparse set of space-time observations; - we considered the problem of estimating an unknown nonlinear observation function of a standard process (e.g. Brownian motion), so that we can recognized if an observed dynamics is "just" a nonlinear version of a known dynamics; we also developed benchmarks for learning algorithms aimed at learning and classifying diffusion processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Engineering osmolysis susceptibility in Cupriavidus necator and Escherichia coli for recovery of intracellular products

Intracellular biomacromolecules, such as industrial enzymes and biopolymers, represent an important class of bio-derived products obtained from bacterial hosts. A common key step in the downstream separation of these biomolecules is lysis of the bacterial cell wall to effect release of cytoplasmic contents. Cell lysis is typically achieved either through mechanical disruption or reagent-based methods, which introduce issues of energy demand, material needs, high costs, and scaling problems. Osmolysis, a cell lysis method that relies on hypoosmotic downshock upon resuspension of cells in distilled water, has been applied for bioseparation of intracellular products from extreme halophiles and mammalian cells. However, most industrial bacterial strains are non-halotolerant and relatively resistant to hypoosmotic cell lysis. To overcome this limitation, we developed two strategies to increase the susceptibility of non-halotolerant hosts to osmolysis using Cupriavidus necator, a strain often used in electromicrobial production, as a prototypical strain. In one strategy, C. necator was evolved to increase its halotolerance from 1.5% to 3.25% (w/v) NaCl through adaptive laboratory evolution, and genes potentially responsible for this phenotypic change were identified by whole genome sequencing. The evolved halotolerant strain experienced an osmolytic efficiency of 47% in distilled water following growth in 3% (w/v) NaCl. In a second strategy, the cells were made susceptible to osmolysis by knocking out the large-conductance mechanosensitive channel (mscL) gene in C. necator. When these strategies were combined by knocking out the mscL gene from the evolved halotolerant strain, greater than 90% osmolytic efficiency was observed upon osmotic downshock. A modified version of this strategy was applied to E. coli BL21 by deleting the mscL and mscS (small-conductance mechanosensitive channel) genes. When grown in medium with 4% NaCl and subsequently resuspended in distilled water, this engineered strain experienced 75% cell lysis, although decreases in cell growth rate due to higher salt concentrations were observed. Our strategy is shown to be a simple and effective way to lyse cells for the purification of intracellular biomacromolecules and may be applicable in many bacteria used for bioproduction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Evaluation of bio-inspired flow fields in a mediated Li-S flow battery for grid energy storage

Lithium-sulfur is a redox flow battery with high energy density for applications in safe, reliable, and lasting scaling of energy. However, lithium-based batteries often encounter platting as a problem thanks to poor Li-ions deposition after cycling. Aiming to reduce this impact, a uniform and continuous flow of ions is needed. On this work, novel bio-inspired flow fields in the electrochemical cell were tested to improve ions flowability and lithium platting control, ultimately enhancing battery performance and life. To secure Li-S efficient, low-cost, and secure energy storage capabilities, we chose a configuration with decamethylferrocene and cobaltocene acting as redox mediators, Li metal as anode and sulfur kept in a separate catholyte reservoir. Flow test and battery results insinuated a beneficial influence of bio-inspired designs in flowing electrolyte uniformly with less pressure and pump power in comparison to other conventional designs used in the industry, with an encouraging ability to approach a cheap, safe, and reliable Li-S grid energy storage.

25 ENERGY STORAGE↗
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↗