Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized Method of Cells”

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 325 records · Page 18

Direct measurement of storage and loss behavior in AFM force–distance experiments using the modified Fourier transformation

Force–distance curve experiments are commonly performed in atomic force microscopy (AFM) to obtain the viscoelastic characteristics of materials, such as the storage and loss moduli or compliances. The classic methods used to obtain these characteristics consist of fitting a viscoelastic material model to the experimentally obtained AFM data. Here, we demonstrate a new method that utilizes the modified discrete Fourier transform to approximate the storage and loss behavior of a material directly from the data, without the need for a fit. Additionally, one may still fit a model to the resulting storage and loss behavior if a parameterized description of the material is desired. In contrast to fitting the data to a model chosen a priori, departing from a model-free description of the material's frequency behavior guides the selection of the model, such that the user may choose the one that is most appropriate for the particular material under study. To this end, we also include modified Fourier domain descriptions of commonly used viscoelastic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Engineering of Lead Selenide Quantum Dot Based Devices and Core/Shell Heterostructures

Near infrared (NIR) emitting colloidal quantum dots (QDs) such as PbSe are interesting materials for implementation in various optoelectronic devices such as solar cells, photodetectors, and radiation detectors. Material properties like size tunable emission wavelengths, facile solution processability, the possibility of carrier multiplication and Auger assisted up-conversion, and high Z number provide PbSe QDs with unique advantages over many currently available commercial materials used in these devices. In this dissertation, after an introduction to QDs in general, we first focus on synthesis of PbSe-based heterostructured QDs for study via optical spectroscopy. PbSe/CdSe QDs synthesized via cation exchange reaction were used as seeds to study the proper conditions for further CdSe or CdS shell growth. Optical spectroscopy studies show that growth of thicker CdSe shells leads to tunable visible emission, while growth of thick CdS shells can greatly increase the photoluminescence quantum yield (PLQY) observed from the intermediate CdSe shell. NIR PL lifetimes can be greatly enhanced by the growth of CdS shells, due to strong delocalization of the hole wavefunction within the PbSe core, while the electron is delocalized throughout the QD including in the CdS shell, resulting in radiative lifetimes in the tens of microseconds. Two pulse delay measurements show that Auger assisted up conversion is also taking place in these engineered QD dot heterostructures. In the sections that follow, we look at QD devices, first with a focus on novel devices fabricated from PbSe QDs, using an amine synthesis method that allows for facile in-solution ligand exchange. QD film mobilities are measured for PbSe QDs capped with various short ionic ligands. Carrier densities in PbSe QD films are calculated from C-V measurements, showing for the first time that carrier densities can be modulated by simply changing the capping ligand. vii Radiation detectors utilizing PbSe QDs capped with KI and NH4I ligands show tangible response to incident alpha radiation. Current output response to alpha radiation is dependent on strength of radiation source and amount of bias voltage applied. In the final section a critical review of hybrid layered 2D-QD photodetectors is presented. In a typical example, PbSe and PbS QDs act as the absorbing layer in phototransistor devices which utilize layered 2D materials such as graphene and transition metal dichalcogenides as transport layers. Resulting devices display high gains and improved directivities compared to QD only and 2D only devices.

36 MATERIALS SCIENCE↗

Discovery of Energy Storage Molecular Materials Using Quantum Chemistry-Guided Multiobjective Bayesian Optimization

Redox flow batteries (RFBs) are a promising technology for stationary energy storage applications due to their flexible design, scalability, and low cost. In RFBs, energy is carried in flowable redoxactive materials (redoxmers) which are stored externally and pumped to the cell during operation. Further improvements in the energy density of RFBs necessitates redoxmer designs with wider redox potential windows and higher solubility. Additionally, designing redoxmers with a fluorescence-enabled self-reporting functionality allows monitoring of the state of health of RFBs. To accelerate the discovery of redoxmers with desired properties, state-of-the- art machine learning ( ML) methods, such as multiobjective Bayesian optimization (MBO), are useful. Here, we first employed density functional theory calculations to generate a database of reduction potentials, solvation free energies, and absorption wavelengths for 1400 redoxmer molecules based on a 2,1,3-benzothiadiazole (BzNSN) core structure. From the computed properties, we identified 22 Pareto-optimal molecules that represent best trade-off among all of the desired properties. We further utilized these data to develop and benchmark an MBO approach to identify candidates quickly and efficiently with multiple targeted properties. With MBO, optimal candidates from the 1400-molecule data set can be identified at least 15 times more efficiently compared to the brute force or random selection approach. Importantly, we utilized this approach for discovering promising redoxmers from an unseen database of 1 million BzNSN-based molecules, where we discovered 16 new Pareto-optimal molecules with significant improvements in properties over the initial 1400 molecules. Furthermore, we anticipate that this active learning technique is general and can be utilized for the discovery of any class of functional materials that satisfies multiple desired property criteria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhancing Post-Irradiation Examination (PIE) Efficiency with High-Throughput Experiments

Nuclear materials testing is essential for understanding changes in reactor environments, but traditional methods are costly and time-consuming. High-Throughput Experimentation (HTE) enables rapid screening and optimization of material properties by significantly increasing the number of specimens in the experiment. In this work, HTE methods were employed to irradiate novel materials at the Advanced Test Reactor (ATR) at Idaho National Laboratory. A total of 367 sub-sized tensile specimens were irradiated using the standard capsule design. To facilitate HTE, automated systems such as the Automated Mechanical Testing System were implemented for mechanical testing within hot cell environments, thereby reducing the need for manual interaction with hot cell manipulators and significantly lowering both the cost and duration of experiments. Digital image correlation (DIC) was used to measure displacement during mechanical testing, providing precise, non-contact measurement capabilities. The methodologies developed in this research will support the advancement of HTE techniques for future nuclear materials testing experiments at ATR. This approach allows for more efficient and cost-effective evaluation of new materials, ultimately leading to improved reactor performance, enhanced safety, and extended longevity within the nuclear industry.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Persistent memory as an effective alternative to random access memory in metagenome assembly

Abstract Background The assembly of metagenomes decomposes members of complex microbe communities and allows the characterization of these genomes without laborious cultivation or single-cell metagenomics. Metagenome assembly is a process that is memory intensive and time consuming. Multi-terabyte sequences can become too large to be assembled on a single computer node, and there is no reliable method to predict the memory requirement due to data-specific memory consumption pattern. Currently, out-of-memory (OOM) is one of the most prevalent factors that causes metagenome assembly failures. Results In this study, we explored the possibility of using Persistent Memory (PMem) as a less expensive substitute for dynamic random access memory (DRAM) to reduce OOM and increase the scalability of metagenome assemblers. We evaluated the execution time and memory usage of three popular metagenome assemblers (MetaSPAdes, MEGAHIT, and MetaHipMer2) in datasets up to one terabase. We found that PMem can enable metagenome assemblers on terabyte-sized datasets by partially or fully substituting DRAM. Depending on the configured DRAM/PMEM ratio, running metagenome assemblies with PMem can achieve a similar speed as DRAM, while in the worst case it showed a roughly two-fold slowdown. In addition, different assemblers displayed distinct memory/speed trade-offs in the same hardware/software environment. Conclusions We demonstrated that PMem is capable of expanding the capacity of DRAM to allow larger metagenome assembly with a potential tradeoff in speed. Because PMem can be used directly without any application-specific code modification, these findings are likely to be generalized to other memory-intensive bioinformatics applications.

59 BASIC BIOLOGICAL SCIENCES↗

Non-Solvent Induced Phase Separation Enables Designer Redox Flow Battery Electrodes

Porous carbonaceous electrodes are performance-defining components in redox flow batteries (RFBs), where their properties impact the efficiency, cost, and durability of the system. The overarching challenge is to simultaneously fulfill multiple seemingly contradictory requirements—i.e., high surface area, low pressure drop, and facile mass transport—without sacrificing scalability or manufacturability. Here, non-solvent induced phase separation (NIPS) is proposed as a versatile method to synthesize tunable porous structures suitable for use as RFB electrodes. The variation of the relative concentration of scaffold-forming polyacrylonitrile to pore-forming poly(vinylpyrrolidone) is demonstrated to result in electrodes with distinct microstructure and porosity. Tomographic microscopy, porosimetry, and spectroscopy are used to characterize the 3D structure and surface chemistry. Flow cell studies with two common redox species (i.e., all-vanadium and Fe2+/3+) reveal that the novel electrodes can outperform traditional carbon fiber electrodes. It is posited that the bimodal porous structure, with interconnected large (>50 µm) macrovoids in the through-plane direction and smaller (<5 µm) pores throughout, provides a favorable balance between offsetting traits. Although nascent, the NIPS synthesis approach has the potential to serve as a technology platform for the development of porous electrodes specifically designed to enable electrochemical flow technologies.

36 MATERIALS SCIENCE↗

The crystal structure of feitknechtite (β-MnOOH) and a new MnOOH polymorph

Studies suggest that feitknechtite (β-MnOOH) is a prevalent, and perhaps necessary, intermediate phase during the synthesis of birnessite-like phases, the abiotic oxidation of Mn2+, and the transformation of biogenic hexagonal phyllomanganates to more complex Mn oxides in laboratory and natural systems. Researchers have generally described feitknechtite as consisting of pyrochroite-like (or cadmium iodide-like) Mn-O octahedral layers, but a detailed crystal structure has not been reported. For this work, we used TEM/SAED and powder XRD and Rietveld refinements to derive the unit cell and, for the first time, report a complete structure description for feitknechtite (β-MnOOH). Rietveld refinements were also completed for three natural feitknechtite/hausmannite samples, and time-resolved synchrotron XRD experiments were used to follow the thermal transformation of feitknechtite to hausmannite. Additionally, we identified and report the structure for a second, and perhaps novel, MnOOH polymorph (proposed designation ε-MnOOH), mixed with the synthetic feitknechtite, that is similar to β-MnOOH but with a different layer stacking.

58 GEOSCIENCES↗

Hydrogen Education for a Decarbonized Global Economy (H 2 EDGE) (Final Technical Report)

The energy sector is undergoing rapid growth and transformation, creating urgent demand for a skilled workforce to support emerging technologies and resilient systems. Hydrogen is an emerging solution for “hard to abate” sectors, including heavy industry and transportation that also provides optionality for the electric system and enables long‑duration energy storage. As hydrogen scales from traditional industrial uses to new applications, success depends on knowledgeable engineers, technicians, planners, and operators safely deploying hydrogen technologies across the full value chain: production, delivery, storage, and end use. The H 2 EDGE initiative, funded by the U.S. Department of Energy’s Hydrogen and Fuel Cell Technologies Office from 2020 to 2025, was launched to prepare people for careers hydrogen and related technologies. Through its interactive approach, H 2 EDGE equipped stakeholders with knowledge and tools to build a reliable and sustainable energy future together. H 2 EDGE built a national network of academic and industry partners to deliver modular training, share open-access curricula, and engage diverse talent pipelines. The program leveraged EPRI’s extensive industry membership and expanded on the train-the-trainer model of the GridEd program to connect educators, employers, and community organizations across all regions of the U.S. New methods were introduced to more systematically map hydrogen competencies and assess curricula using traditional and AI-assisted approaches. Outcomes included delivery of professional short courses, developing university curriculum, sponsoring faculty and student projects, and engaging stakeholders through workshops, site tours, and webinars.

08 HYDROGEN↗

Direct evidence of the effect of a moderate external magnetic field on stimulated Raman scattering in the kinetic regime

We present results from an experiment carried out at the OMEGA-EP laser facility that investigated the effect of a perpendicular magnetic field on stimulated Raman scattering (SRS) and report the first direct measurement of magnetic mitigation of SRS. A 13-T magnetic field generated by pulsed-power coils was imposed on a gas jet plasma, and a novel three-picket interaction beam was used to explore SRS reflectivity for several plasma conditions within single shots. The time-resolved backscattered light shows that SRS was mitigated by the external 13-T magnetic field in the kinetic regime (kλ D ∼ 0.3, where k is the electron plasma wave's wavenumber and λ D is electron Debye length) and at a density of n e /n cr ∼ 0.10. On the other hand, we also measured an enhancement of SRS reflectivity at lower density (n e /n cr < 0.08). We discuss experimental results in the context of magnetohydrodynamic and particle-in-cell simulations that scan SRS dynamics for a variety of plasma conditions. While the experimental evidence of SRS mitigation validates prior work on the kinetic simulation of SRS in an external magnetic field, we also find that other mechanisms such as SRS rescatter can lead to an enhancement of measured SRS reflectivity in simulations of parameters relevant to our experiment.

Gas jet↗

Corrosion testing needs and considerations for additively manufactured materials in nuclear reactors

Metal additive manufacturing holds significant promise as an enabling technology for the 21st century nuclear energy industry. Metal additive manufacturing (MAM) can allow the fabrication of novel materials and innovative component designs that are not achievable through conventional manufacturing. Due to its very different fabrication methods, as-fabricated MAM components are characteristically different from conventionally manufactured components. MAM materials exhibit very different microstructures from conventional cast or wrought materials. For example, austenitic stainless steels fabricated by laser powder bed fusion exhibit columnar grain structures, dislocation cell structures, and melt pool fingerprints. In addition, MAM fabrication may result in defects such as porosity, incomplete processing of the feedstock (e.g., lack of fusion in melt-based methods), and oxide inclusions. Heat treatments may further evolve the microstructure, microsegregation, and stresses within the component. Furthermore, MAM components have a rough surface with feature sizes on the order of the feedstock material, as opposed to smooth surfaces resulting from conventional machining and forming operations. As a result, the corrosion behavior of MAM components will likely be significantly different from that of conventionally formed components. Corrosive environments for structural materials within advanced reactor environments include molten fluoride and chloride salts, liquid sodium and lead-bismuth, and high-temperature helium. The Advanced Materials and Manufacturing Technologies program within the Department of Nuclear Energy in the United States Department of Energy is assessing the unique concerns of MAM component corrosion and testing methodologies in advanced nuclear reactor environments. We discuss these concerns and testing strategies in this paper.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Capturing thin structures in VOF simulations with two-plane reconstruction

A novel interface reconstruction strategy for volume of fluid (VOF) methods is introduced that represents the liquid-gas interface as two planes that co-exist within a single computational cell. In comparison to the piecewise linear interface calculation (PLIC), this new algorithm greatly improves the accuracy of the reconstruction, in particular when dealing with thin structures such as films. The placement of the two planes requires the solution of a non-linear optimization problem in six dimensions, which has the potential to be overly expensive. Further, an efficient solution to this optimization problem is presented here that exploits two key ideas: an algorithm for extracting multiple plane orientations from transported surface data, and an efficient and mass-conserving distance-finding algorithm that accounts for two planes with arbitrary orientation. Additionally, a simple and robust strategy is presented to accurately represent the surface tension forces produced at the interface of subgrid-thickness films. The performance of this new VOF reconstruction is demonstrated on several test cases that illustrate the capability to handle arbitrarily thin films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗

Dirac Node Pinning from Dzyaloshinskii-Moriya Interactions in a Kagome Spin Liquid

Recent experiments on the Kagome spin liquid candidate material YCu 3 ⁢(OH) 6 ⁢Br 2 ⁢[Br 1−𝑦 ⁢(OH) 𝑦 ] suggest the presence of Dirac fermionic spinons near the magnetization plateau at 1/9. Theories suggest that the spinons are charge neutral spin-1/2 excitations, in a 2⁢𝜋/3 flux, which triples the unit cell. Generally a gap is expected, and there is no symmetry protection for the Dirac nodes in this system. The question arises as to what causes the nodes and stabilizes them. In this work, we propose a node-creation and node-pinning mechanism driven by Dzyaloshinskii-Moriya (DM) interactions. Employing Gutzwiller-projected variational Monte Carlo calculations, we demonstrate that DM interactions induce a band closing phase transition in the spinon spectrum. There is a change in the Chern number when the bands are inverted. Together with the DM-generated internal gauge flux, the coupling to the spinon orbital magnetization counteracts the band reopening. Furthermore, this interplay energetically pins the Dirac nodes over a range of parameters, resulting in a pinning mechanism distinct from the usual one from symmetry protection.

Kagome lattice↗

A Novel use of Direct Simulation Monte-Carlo to Model Dynamics of COVID-19 Pandemic Spread

In this report, we evaluate a novel method for modeling the spread of COVID-19 pandemic. In this new approach we leverage methods and algorithms developed for fully-kinetic plasma physics simulations using Particle-In-Cell (PIC) Direct Simulation Monte-Carlo (DSMC) models. This approach then leverages Sandia-unique simulation capabilities, and High-Performance Computer (HPC) resources and expertise in particle-particle interactions using stochastic processes. Our hypothesis is that this approach would provide a more efficient platform with assumptions based on physical data that would then enable the user to assess the impact of mitigation strategies and forecast different phases of infection. This work addresses key scientific questions related to the assumptions this new approach must make to model the interactions of people using algorithms typically used for modeling particle interactions in physics codes (kinetic plasma, gas dynamics). The model developed uses rational/physical inputs while also providing critical insight; the results could serve as inputs to, or alternatives for, existing models. The model work presented was developed over a four-week time frame, thus far showing promising results and many ways in which this model/approach could be improved. This work is aimed at providing a proof-of-concept for this new pandemic modeling approach, which could have an immediate impact on the COVID-19 pandemic modeling, while laying a basis to model future pandemic scenarios in a manner that is timely and efficient. Additionally, this new approach provides new visualization tools to help epidemiologists comprehend and articulate the spread of this and other pandemics as well as a more general tool to determine key parameters needed in order to better predict pandemic modeling in the future. In the report we describe our model for pandemic modeling, apply this model to COVID-19 data for New York City (NYC), assess model sensitivities to different inputs and parameters and , finally, propagate the model forward under different conditions to assess the effects of mitigation and associated timing. Finally, our approach will help understand the role of asymptomatic cases, and could be extended to elucidate the role of recovered individuals in the second round of the infection, which is currently being ignored.

59 BASIC BIOLOGICAL SCIENCES↗

Can Cooling and Heating Functions Be Modeled with Homogeneous Radiation Fields?

Abstract Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between 5 ≤ z ≤ 10 in simulated galaxies from the Cosmic Reionization On Computers project. We compare three functions: (1) the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature, (2) the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM), and (3) the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as CLOUDY.

79 ASTRONOMY AND ASTROPHYSICS↗

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

Here we present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature thereby integrating over the thermal seed photon flux from the disk surface. The method either assumes the ion and electron temperatures are equal (1T) or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. This new approach leads to radiative efficiency values sime65% above Novikov–Thorne, with a larger fraction of total cooling in the corona as compared to simulations performed using the original target-temperature cooling function. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach, which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

79 ASTRONOMY AND ASTROPHYSICS↗

Positioning Accuracy in a Concurrent Robot-CNC Hybrid Manufacturing System

Abstract Additive manufacturing (AM) has gained notoriety for offering advantages over traditional manufacturing methods, such as increased design complexity and flexibility. However, it has not found widespread use beyond rapid prototyping. One hindrance to the acceptance of AM processes in industry is the time and cost of fabrication per component. While metal AM by itself can be inexpensive, extra manufacturing steps in the form of subtractive manufacturing (SM) may need to be performed to reach final part tolerances, leading to hybrid additive-subtractive manufacturing (HASM) of a part, which increases time and cost. A potential area to reduce cost is through increasing the efficiency of the HASM process by conducting additive and subtractive manufacturing simultaneously. Usually, HASM is performed in a process where AM is completed in one machine or cell and transferred to another machine or cell for SM in a sequential assembly line process. This efficiency decreases part cost, but high aspect ratio parts or parts with internal geometry that require interleaved additive deposition and machining cannot be produced. One unexplored solution to simultaneous HASM that allows for interleaved operations is to operate the deposition head and machining spindle concurrently within the same machine envelope, known as concurrent HASM (CHASM). In this type of process, both AM and SM occur simultaneously on a batch of small parts or a single large part, maintaining a high efficiency without sacrificing the full range of complex geometries that AM allows for. A potential approach to the single-machine method could be to combine a robot and mill within the same envelope. A challenge to this approach, however, is control of both systems. Most machine controllers have limited external communication or, if a robot has been integrated, only offer movement of either the robot or mill at any given time. As a result, systems must pause either the AM or SM process to switch between them rather than working simultaneously. The present work investigates the positional accuracy of such a CHASM system comprised of a robotic arm and a 3-axis mill. Open-loop tests with limited communication between machines are performed on the system to verify positional error during concurrent robot-mill movements. Under certain conditions, it is demonstrated that position error can stay within 2 mm for the duration of a single layer; however, these tests show that, generally, the open-loop positioning performance of the system is inadequate for CHASM without part-specific hand-tuning of parameters. Based on these results, a set of requirements for successful robot-CNC CHASM is proposed for future integrations.

Goodwin, Jesse↗

Numerical study of magnetic island coalescence using magnetohydrodynamics with adaptively embedded particle-in-cell model

Collisionless magnetic reconnection typically requires kinetic treatment that is, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with an adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where the kinetic physics is prominent. We compare the simulation results among three cases: (1) MHD with adaptively embedded PIC regions, (2) MHD with statically (or fixed) embedded PIC regions, and (3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations as well as the ion pressure tensor elements and ion agyrotropy. In order to reach good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature, where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

79 ASTRONOMY AND ASTROPHYSICS↗