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 181 records · Page 10

Lessons from the IEC Durability of Adhesion Accelerated Test Sequence

The IEC 62788-1-1 and IEC 63209-2 standards use aging sequences for durability of adhesion in photovoltaic (PV) modules, which may be evaluated using the single cantilever beam (SCB) test. Because the encapsulant forms critical interfaces with the front glass and solar cells, degradation at those interfaces under ultraviolet (UV) exposure, elevated temperature, and humidity can lead to interfacial delamination - compromising the long-term reliability. In this work, adhesion durability of UV-transmitting poly(ethylene-co-vinyl acetate) (EVA) encapsulant to glass and to silicon solar cells is evaluated after sequenced UV and damp heat aging (85C/85%RH). Laminates were prepared using StarPhire solar front glass with thin glass or PERC cells, and two EVA formulations with different concentrations of siloxane coupling agent. Adhesion was quantified by measuring critical debond energy using the SCB method. Both formulations exhibit similar qualitative trends, while different adhesion is observed at the periphery despite the use of low-shrink manufacturing. The results show that while glass/EVA adhesion remains stable or increases after UV exposure and shows only moderate changes after damp heat, the EVA/cell interface exhibits an irreversible loss of adhesion following UV and then damp heat exposure. Although glass/EVA interfaces generally exhibit lower debond energies, the EVA/cell interface is significantly more vulnerable to UV-driven degradation, identifying it as the dominant reliability risk location through early- and intermediate-module life. These results demonstrate that accelerated aging sequences can expose large, interface-specific losses in adhesion durability and underscore the importance of interface engineering for long-term PV module reliability.

14 SOLAR ENERGY↗

A robust method of nuclei isolation for single-cell RNA sequencing of solid tissues from the plant genus Populus

Single-cell transcriptome analysis has been extensively applied in humans and animal models to uncover gene expression heterogeneity between the different cell types of a tissue or an organ. It demonstrated its capability to discover key regulatory elements that determine cell fate during developmental programs. Single-cell analysis requires the isolation and labeling of the messenger RNA (mRNA) derived from each cell. These challenges were primarily addressed in mammals by developing microfluidic-based approaches. For plant species whose cells contain cell walls, these approaches have generally required the generation of isolated protoplasts. Many plant tissues’ secondary cell wall hinders enzymatic digestion required for individual protoplast isolation, resulting in an unequal representation of cell types in a protoplast population. This limitation is especially critical for cell types located in the inner layers of a tissue or the inner tissues of an organ. Consequently, single-cell RNA sequencing (scRNA-seq) studies using microfluidic approaches in plants have mainly been restricted to Arabidopsis roots, for which well-established procedures of protoplast isolation are available. Here we present a simple alternative approach to generating high-quality protoplasts from plant tissue by characterizing the mRNA extracted from individual nuclei instead of whole cells. We developed the protocol using two different plant materials with varying cellular complexity levels and cell wall structure, Populus shoot apices, and more lignified stems. Using the 10× Genomics Chromium technology, we show that this procedure results in intact mRNA isolation and limited leakage, with a broad representation of individual cell transcriptomes.

54 ENVIRONMENTAL SCIENCES↗

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗

Syntheses and Crystal Structures of Rare-Earth Oxyapatites Ca 2 RE 8 (SiO 4 ) 6 O 2 (RE = Pr, Tb, Ho, Tm)

Four different rare-earth oxyapatites of Ca 2 RE 8 (SiO 4 ) 6 O 2 (RE = Pr, Tb, Ho, Tm) were synthesized using a solution-based method followed by drying, calcination, and high-temperature sintering in air. X-ray powder diffraction and Raman spectroscopy were performed on the synthesized oxyapatites. Here, the RE oxyapatites crystallize in the hexagonal space group P6 3 /m with similar unit cell parameters, increasing linearly with larger RE cations. The unit cell volumes increase linearly whereas the densities decrease nonlinearly with larger RE cations. Raman spectra showed intense bands of the symmetric bending and stretching modes of SiO 4 at ~ 400 and 860 cm -1 regions, respectively. The bands generally shifted to higher frequencies with smaller RE cations in the structures.

36 MATERIALS SCIENCE↗

Security Constrained Distributed Transaction Model for Multiple Prosumers

Massive access of renewable energy has prompted demand-side distributed resources to participate in regulation and improve flexibility of power systems. With large-scale access of massive, decentralized, and diverse distributed resources, demand-side market members have transformed from traditional “consumers” to “prosumers”. To explore the distributed transaction model of prosumers, in this paper, a multi-prosumer distributed transaction model is proposed, and the Conditional Value-at-Risk (CVaR) theory is applied to quantify potential risks caused by the stochastic characteristics inherited from renewable energy. First, a prosumer model under constraints of the distribution network including photovoltaic units, fuel cells, energy storage system, central air conditioning and flexible loads is established, and a multi-prosumer distributed transaction strategy is proposed to achieve power sharing among multiple prosumers. Second, a prosumer transaction model based on CVaR is constructed to measure risks inherited from the uncertainty of PV output within the prosumer and ensure safety of system operation in extreme PV output scenarios. Then, the alternating direction multiplier method (ADMM) is utilized to solve the constructed model efficiently. Finally, distributed transaction costs of prosumers are distributed fairly based on the generalized Nash equilibrium to maximize social benefits. Simulation results show the multi-prosumer distributed transaction mechanism established under the proposed generalized Nash equilibrium method can encourage power sharing among prosumers, increasing their own income and social benefits. Also, the CVaR can assist decision making of prosumers in weighting the risks and benefits, improving system resilience through energy management of prosumers.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Landau level of fragile topology

Here, we study the Hofstadter butterfly and Landau levels of the twisted bilayer graphene (TBG). We show that the nontrivial fragile topology of the lowest two bands near the charge neutral point makes their Hofstadter butterfly generically connected with higher bands, closing the gap between the first and second conduction (valence) bands at a certain magnetic flux per unit cell. We also develop a momentum space method for calculating the TBG Hofstadter butterfly, from which we identify three phases where the Hofstadter butterflies of the lowest two bands and the higher bands are connected in different ways. We show this leads to a crossing between the $\textit{ν}$ = 4 Landau fan from the charge neutral point and the zero-field band gap at one flux per moiré unit cell, which corresponds to a magnetic field 25$θ^2$ T (twist angle $\textit{θ}$ in degrees). This provides an experimentally testable feature of the fragile topology. In general, we expect it to be a generic feature that the Hofstadter butterfly of the topological bands are connected with the Hofstadter spectra of other bands. We further show the TBG band theory with Zeeman splitting being the most sizable splitting could result in Landau fans at the charge neutral point and half fillings near the magic angle, and we predict their variations under an in-plane magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Annealing of ion tracks in apatite under pressure characterized in situ by small angle x-ray scattering

Fission track thermochronology is routinely used to investigate the thermal history of sedimentary basins, as well as tectonic uplift and denudation rates. While the effect of temperature on fission track annealing has been studied extensively to calibrate the application of the technique, the effect of pressure during annealing is generally considered to be negligible. However, a previous study suggested elevated pressure results in a significantly different annealing behaviour that was previously unknown. Here, we present a method to study track annealing in situ under high pressure by using synchrotron- based small angle x-ray scattering (SAXS). To simulate fission tracks in a controlled environment, ion tracks were created in apatite from Durango, Mexico using 2 GeV Au or Bi ions provided by an ion accelerator facility. Samples were annealed at 250°C at approximately 1 GPa pressure using diamond anvil cells (DAcs) with heating capabilities. Additional in situ annealing experiments at ambient pressure and temperatures between 320 and 390°C were performed for comparison. At elevated pressure a significantly accelerated annealing rate of the tracks was observed compared with annealing at ambient pressure. However, when extrapolated to geologically relevant temperatures and pressures, the effects become very small. The measurement methodology presented provides a new avenue to study materials behaviour in extreme environments.

58 GEOSCIENCES↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Colossal grain growth in Cd(Se,Te) thin films and their subsequent use in CdTe epitaxy by close-spaced sublimation

Many technologies deposit thin films on inexpensive substrates, resulting in small grains due to classic nucleation and grain growth theory. For example, state-of-the-art solar cells are made by depositing CdSeTe and CdTe layers on inexpensive glass coated with nanocrystalline transparent conducting oxides (TCOs), like SnO 2 . Characteristically, the grain size of these films is on the order of the film thickness, i.e. a few microns. CdTe small-grain films have poor electro-optical properties and require CdCl 2 passivation which fails to fully passivate grain boundaries, causes carrier compensation, and prevents implementing other II–VI alloys and materials to improve performance. Here, we present a method to increase grain size to 1 mm in CdSe x Te 1-x thin films deposited on glass/TCO substrates without CdCl 2 treatment. The colossal grain growth is driven by mechanisms distinct from classic nucleation, grain growth, and Ostwald ripening and only occurs at low selenium content (x ~ 0.1). We also demonstrate how these films can serve as templates for subsequent large-grain epitaxy of other compositions like CdTe, again without exposure to CdCl 2 . The results open new paths for thin film solar cell technology, and thin film devices in general.

36 MATERIALS SCIENCE↗

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Interplay of coulomb and exciton–phonon coupling controls singlet fission dynamics in two pentacene polymorphs

Pentacene is an important model organic semiconductor in both the singlet exciton fission (SF) and organic electronics communities. We have investigated the effect of changing crystal structure on the SF process, generating multiple triplet excitons from an initial singlet exciton, and subsequent triplet recombination. Unlike for similar organic semiconductors that have strong SF sensitive to polymorphism, we find almost no quantitative difference between the kinetics of triplet pair (TT) formation in the two dominant polymorphs of pentacene. Both pairwise dimer coupling and momentum-space crystal models predict much faster TT formation from the bright singlet excited state of the Bulk vs ThinFilm polymorph, contrasting with the experiment. GW and Bethe–Salpeter equation calculations, including exciton–phonon coupling, reveal that ultrafast phonon-driven transitions in the ThinFilm polymorph compensate the intrinsically slower purely Coulomb-mediated TT formation channel, rationalizing the similarity in observed rates. Taking into account the influence of subtle structural distinctions on both the direct and phonon-mediated SF pathways reveals a predictive capability to these methods, expected to be applicable to a wide variety of molecular crystals.

14 SOLAR ENERGY↗

Salt Sampling FY21 Technical Report

The goal of the salt sampling program at Argonne is to develop and deploy automated molten salt sampling approaches to enable high-precision in-process salt sample analysis to improve the timeliness of sampling-based accountancy measurements. Tools currently under development in support of this goal include (1) a modular vacuum sampler with an accompanying sample handling method for coupling vacuum sampling with high-precision at-line sample analysis, (2) a pneumatic sample generator that enables high-throughput sample analysis to improve the precision of existing analytical techniques, and (3) a windowless flow cell to enable on-line optical analysis of molten salt in a sampling loop. Compared to point sampling approaches (i.e., dip probes), vacuum sampling systems and on-line sampling loops facilitate access to a larger cross-section of a process fluid. This is known to improve the characterization of the process fluid by producing more representative samples and by enabling the analysis of a larger cross section of the fluid. A vacuum sampling approach for molten salts eliminates the risk of dross contamination of samples and avoids the use of moving parts in the salt. In FY21, two methods for integrating a vacuum sampler with a pneumatic sample generator were tested. These included direct fluidic coupling and coupling using a solid salt transfer mechanism. Solid salt transfer was ultimately selected over fluidic coupling, primarily to enable the transport of samples over longer distances to support automated at-line integration with high-precision techniques (such as microcalorimetry) that cannot withstand the extreme conditions near an electrorefining process. To facilitate rapid solid salt coupling, new mechanisms were developed for rapidly charging and discharging salt sample tubes at the vacuum sampler and pneumatic sample generator, respectively. While the charging mechanism will be deployed in FY22, the tube transfer method and discharge mechanism were tested in FY21. These were deployed at one of Argonne’s engineering-scale electrorefiners to implement at-line high-throughput pneumatic micro-sample generation capabilities. The method was used to generate precise uranium- and lanthanide-bearing electrorefiner micro-samples with the specific dimensions requested by researchers at Los Alamos National Laboratory for use in testing their novel microcalorimeter x-ray techniques. The solid salt transfer mechanism proved not only to be an effective means of integrating the precision sample generator with vacuum sampling, but also improved the performance of the sampler generator. Because the modular sampling approach described here eliminates the need for new high-radiation sample handling capabilities, salt-wetted seals, salt-wetted moving parts, and heated transfer lines outside the electrorefiner, it will address most of the remaining technical challenges for the at-line deployment of high-precision analytical techniques. This will enable significant reductions in the time delay for sampling-based accountancy measurements by eliminating the need for manual off-line sample processing and analysis. On-line optical analysis of molten salt in a sampling loop would provide complementary information to at-line and in-situ techniques. In FY21, an open-aperture molten salt gravity flow cell with windowless optical access to flowing salt was successfully demonstrated. Future work should include the refinement and performance testing of the on-line and at-line sampling tools, integration of additional analysis techniques, stakeholder outreach and collaboration, evaluation of the integrated methods, and analyses to determine how the various tools might fit into an integrated safeguards monitoring system of unattended near real time monitoring tools.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparative mitogenomics of kingdom Fungi – evolutionary insights and metagenomic applications

Mitochondria are essential components of eukaryotic cells, responsible for ATP production through oxidative phosphorylation. Despite their biological importance, unique challenges have hindered the adoption of automated mitochondrial genome (mitogenome) annotation methods, obstructing mitochondrial comparative genomics in a broad evolutionary context. Using Fungi as a study system and a Joint Genome Institute (JGI) annotated high-quality reference set, we observed broad patterns of mitochondrial evolution across the kingdom. We found that the median fungal mitogenome size is 58 kb and identified exceptionally large examples over 1 Mb in Pezizomycetes. All 14 expected oxidative phosphorylation protein-coding genes, plus rps3, were generally conserved. We found evidence of major evolutionary transitions within the Ascomycota, including the transfer of mitochondrially encoded atp8 and atp9 to the nuclear genomes across the Pezizomycotina and shifts in mitogenome tRNA patterns across the kingdom. We found substantial concordance between mitochondrial and nuclear evolution, enabling us to document 3131 total fungal mitogenomes from JGI-derived metagenomic datasets. We also identified 6467 total undeclared mitogenomes embedded in Genbank fungal nuclear assemblies. We provide interactive tools for mitogenome analysis through the JGI MycoCosm platform. Collectively, this work generated nearly 10 000 new fungal mitogenome annotations, providing a foundation and resources for future exploration of comparative fungal mitogenomics.

Ahrendt, Steven R. [USDOE Joint Genome Institute (↗

Dynamic Species Reduction for Multi-Cycle CFD Simulations (Final Technical Report)

This project primarily sought to address some of the computational cost concerns of detailed simulations by developing improved methods of handling species transport, and chemical kinetics evaluation in a commercial 3D Computational Fluid Dynamics (CFD) environment. Secondary goals were to apply these techniques to fuels and conditions of interest to better understand the key species and reactions required to adequately capture cycle to cycle coupling. Improved modeling will lead to better understanding of these combustion modes and their dependence on fuel composition, which can then enable more clean and efficient engines, ultimately benefiting the consumer as well as the general public. Two approaches were used to address the computational cost. A “Dynamic Species Reduction” (DSR) procedure was developed to remove chemical species from the simulation during periods when chemical reactions were not expected to be important, particularly during gas exchange when temperatures are low and little fuel remains. This modelling procedure automatically detects relevant species to retain in the simulation domain based on their local concentration, removes species below a specified concentration threshold, and then adjusts the remaining species mass to conserve not only the number of H, C, and O atoms in each cell, but also the relative proportions between species and the heating value of the mixture in every cell. The second procedure was a “Product Directed Remapping” (PDR) updates the algorithm used to group individual computational cells for chemical kinetics evaluation to account for non-uniform temperature distributions and the CO to CO 2 ratio in the cells. The model techniques developed in this work successfully demonstrated computational performance improvements for a range of conditions relevant for Highly Dilute SI and HCCI engine operation. Runtime reductions of 10% were observed for small mechanisms, with further reductions of up to 36% observed for larger mechanisms. Improvements were primarily related to reducing the number of chemical species tracked during the gas exchange process using the Dynamic Species Reduction method. With smaller benefits observed from changes to the kinetics binning and evaluation strategy in the post combustion region using the PDR method. The results of this work show the potential for improved computational runtime for complicated simulations. They can and should be extended to additional conditions and new bio-derived and renewable fuels as they are developed and new kinetic mechanisms become available.

10 SYNTHETIC FUELS↗

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Non-Contact Eddy Current Method for Assessing Proton Conductivity in Nafion Membranes

Accurate measurement of proton conductivity in Proton Exchange Membranes (PEMs) is vital for fuel cell performance and durability. This study explores eddy current testing as a non-destructive, non-contact and high-rate quality control (QC) method for measuring Nafion membrane proton conductivity, with potential for high-volume manufacturing applications.

33 ADVANCED PROPULSION SYSTEMS↗

Symbol alphabets from plabic graphs III: n = 9

Symbol alphabets of n-particle amplitudes in N = 4 super-Yang-Mills theory are known to contain certain cluster variables of G(4, n) as well as certain algebraic functions of cluster variables. In this paper we solve the C Z = 0 matrix equations associated to several cells of the totally non-negative Grassmannian, combining methods of arXiv:2012.15812 for rational letters and arXiv:2007.00646 for algebraic letters. We identify sets of parameterizations of the top cell of G + (5, 9) for which the solutions produce all of (and only) the cluster variable letters of the 2-loop nine-particle NMHV amplitude, and identify plabic graphs from which all of its algebraic letters originate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Intersection-distribution-based remapping between arbitrary meshes for staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics

In this paper, we present a new intersection-distribution-based remapping method between arbitrary polygonal meshes for indirect staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics. All cell-centered material quantities are conservatively remapped using intersections between the Lagrangian (old, source) mesh and the rezoned (new, target) mesh. The new nodal masses are obtained by conservative distribution of all material masses in each new cell to the cell's corners and then collecting those corner masses at new nodes. This distribution is done using a local constrained optimization approach for each cell in the new mesh. In order to remap nodal momentum we first define cell-centered momentum for each cell in the old mesh, conservatively remap this to the new mesh and then conservatively distribute the new zonal momentum to each cell's bounding nodes, again using local constrained optimization. Our method also conserves total energy by applying a new nodal kinetic energy correction that relies on a process similar to that used for remapping nodal mass and momentum. Cell-centered kinetic energy is computed, conservatively remapped and then distributed to nodes. The discrepancy between this conservatively remapped and actual nodal kinetic energy is then conservatively distributed to the internal energies of the materials in the cells surrounding each node. Unlike conventional cell-based corrections of this type, this new nodal kinetic energy correction has not been observed to drive material internal energy negative in any of our testing. Unlike flux based remapping, our new intersection-distribution method can be applied to remapping between source and target meshes that are arbitrarily different, which provides superior flexibility in the rezoning strategy. Our method is accurate, essentially conservative and essentially bounds preserving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗