SEARCH · Engineering Papers
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.
PyDislocDyn [1] is a suite of python programs designed to perform various calculations for basic research in dislocation dynamics in metals with various crystal symmetries in the continuum limit.
Over the past decade, there has been enormous progress in the description of nuclear structure from first principles, using interactions from Chiral Effective Field Theory that are rooted in Quantum Chromodynamics, the fundamental theory of strong interactions, and many-body methods that solve the Schrödinger equation with systematically improvable approximations. Amongst those are the family of In-Medium Similarity Renormalization Group (IMSRG) framework developed by the PI and his co-workers. While these methods scale polynomially in the size N of the single-particle basis, computational efforts still grows dramatically as we increase the degrees of freedom for the nucleons by relaxing symmetries or introducing continuum couplings (see below), or as we push to improved truncations to provide precise inputs for experimental efforts, in particular in fundamental symmetry searches. In order to address this growing computational cost, one focus area of this award was the exploration of compression and factorization methods. The key to success or failure is the presence of low-rank structures within the matrix elements of NN and 3N interactions or the IMSRG evolution operator, and a means to reformulate the method that will let us exploit them.
Context. Galactic cosmic rays (CRs) play a crucial role in ionisation, dissociation, and excitation processes within dense cloud regions where UV radiation is absorbed by dust grains and gas species. CRs regulate the abundance of ions and radicals, leading to the formation of more and more complex molecular species, and determine the charge distribution on dust grains. A quantitative analysis of these effects is essential for understanding the dynamical and chemical evolution of star-forming regions. Aims. The CR-induced photon flux has a significant impact on the evolution of the dense molecular medium in its gas and dust components. This study evaluates the flux of UV photons generated by CRs to calculate the photon-induced dissociation and ionisation rates of a vast number of atomic and molecular species, as well as the integrated UV photon flux. Methods. To achieve these goals, we took advantage of recent developments in the determination of the spectra of secondary electrons, in the calculation of state-resolved excitation cross sections of H 2 by electron impact, and of photodissociation and photoionisation cross sections. Results. We calculated the H 2 level population of each rovibrational level of the X, B, C, B', D, B", D', and a states. We then computed the UV photon spectrum of H 2 in its line and continuum components between 72 and 700 nm, with unprecedented accuracy, as a function of the CR spectrum incident on a molecular cloud, the H 2 column density, the isomeric H 2 composition, and the dust properties. The resulting photodissociation and photoionisation rates are, on average, lower than previous determinations by a factor of about 2, with deviations of up to a factor of 5 for the photodissociation of species such as AlH, C 2 H 2 , C 2 H 3 , C 3 H 3 , LiH, N 2 , NaCl, NaH, O$^+_2$, S 2 , SiH, l-C 4 , and l-C 5 H. A special focus is given to the photoionisation rates of H 2 , HF, and N 2 , as well as to the photodissociation of H 2 , which we find to be orders of magnitude higher than previous estimates. We give parameterisations for both the photorates and the integrated UV photon flux as a function of the CR ionisation rate, which implicitly depends on the H 2 column density, as well as the dust properties.
Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.
Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.
The static aerodynamics for the Mars Exploration Rover (MER) aeroshell are presented. This aerodynamic database was an integral part of the end-to-end simulation used in pre- entry analysis for determining the MER entry design requirements for development of the MER entry system, as well as targeting the MER landing sites. The database was constructed using the same approach used for Mars Pathfinder (MPF). However, the MER aerodynamic database is of much higher fidelity and tailored to the MER entry trajectories. This set of data includes direct simulation Monte Carlo calculations covering the transitional regime of the entry trajectory and computational fluid dynamics calculations describing the aerodynamics in the hypersonic and supersonic continuum regimes. An overview of the methodology used to generate the data is given along with comparisons to important features in the MPF aerodynamics and related heritage data. The MER and MPF comparison indicates that trajectory specific data is required to properly model the flight characteristics of a.blunt entry capsule at Mars.
The height, size, directivity, and lifetime of decametric storm sources are investigated. The height of the continuum sources is determined from the rotation rate, and the gradient of electron density for the regions of storm origin is computed by assuming that the radiation originates at the plasma frequency. The mean angular size of the sources is found to be quite large (increasing with decreasing frequency), the storm continuum is strongly directive toward the disk center, and east-west asymmetry is observed at decameter wavelengths. Two distinct classes of type III bursts are observed during storms: 'off-fringe' (displaced in position from the continuum source) and 'on-fringe' (coinciding in position with the continuum). A model of the storm region is proposed in which the 'on-fringe' bursts originate in regions of diverging field lines above closed magnetic loops and the 'off-fringe' bursts are excited by energetic electron streams having access to open field lines at the base of the loops.
Modeling of microwave heating of multiphase systems requires resolving the electromagnetic field from the single particle to the entire cavity. We introduce a multiscale methodology for computationally affordable simulations of microwave heating of multiphase systems consisting of one phase dispersed in a continuum phase. The methodology homogenizes the original multiphase system, by taking advantage of the large separation of length scales, to calculate the effective permittivity, effective thermal conductivity, and volumetric power absorbed in each phase. The methodology is rigorously assessed against particle-resolved detailed numerical simulations of model systems. We demonstrate a significant reduction in computational cost while retaining the accuracy of detailed simulations. Lastly, the present methodology enables high throughput exploration of novel multiphase system designs utilizing microwave heating.
Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.
The project involved the study of ultraviolet (UV) and visible emission mechanisms in astrophysical and atmospheric environments. In many situations, the emission is a direct consequence of a charge transferring collision of an ion with a neutral with capture of an electron to an excited state of the product ion. The process is also important in establishing the ionization and thermal balance of an astrophysical plasma. As little of the necessary collision data are available, the main thrust of the project was the calculation of total and state-selective charge transfer cross sections and rate coefficients for a very large number of collision systems. The data was computed using modern explicit techniques including the molecular-orbital close-coupling (MOCC), classical trajectory Monte Carlo (CTMC), and continuum distorted wave (CDW) methods. Estimates were also made in some instances using the multichannel Landau-Zener (MCLZ) and classical over-the-barrier (COB) models. Much of the data which has been computed has been formatted for inclusion in a charge transfer database on the World Wide Web (cfadc.phy.ornl.gov/astro/ps/data/). A considerable amount of data has been generated during the lifetime of the grant. Some of it has not been analyzed, but it will be as soon as possible, the data placed on our website, and papers ultimately written.
Numerical models used to study transport phenomena in laser-powder bed fusion processes often rely on assumptions and simplifications to reduce their computational expense. One common simplification is in the description of latent heat evolution during the solid-liquid phase change (i.e., the solidification pathway), justified by the fact that the mushy zone thickness is similar to the numerical grid spacing used for continuum transport models. The lack of resolution of transport phenomena in the mushy zone motivates the use of computationally convenient solidification paths such as linear or sigmoidal relationships over pathways derived from fundamental solidification theory such as equilibrium or Scheil models. In the present work, an uncertainty quantification (UQ) framework is used to analyze the influence of solidification pathway selection on the solidification dynamics and melt pool geometries in laser based additive manufacturing (AM) of IN625. Results show the solidification pathway has a quantifiable influence on the cooling rate at the liquidus isotherm, mushy zone thickness, and solidification time. Due to similarities in the latent heat evolution at the beginning of solidification, the equilibrium and Scheil models predict similar cooling rates near the liquidus isotherm, however the wider freezing range of Scheil leads to a wider mushy zone compared to equilibrium. The non-physical latent heat release profiles of sigmoidal and linear paths lead to significant overpredictions of cooling rates at the liquidus isotherm compared to equilibrium and Scheil. Finally, these results indicate that careful consideration should be given to the choice of solidification pathway to ensure reliable model predictions.
This report documents a validation study on the standard Open Hole Compression (OHC) laminate test specimen conducted as part of the NASA Advanced Composites Project (ACP). Tests were conducted on OHC specimens with hard, quasi-isotropic, and soft layups using digital image correlation and X-Ray computed tomography to capture the structural response and damage evolution. Progressive damage models were constructed for use with the CompDam continuum damage mechanics code following the best practices established during the ACP. Detailed interrogation of the analysis results and comparison with experimental measurements provide a basis for assessing the capability of the modeling approach for OHC. The structural response is found to be captured well, with strengths predicted within 3\% of the experimental values for hard and quasi-isotropic laminates. In the soft laminate, the model predicts failure to be more brittle than the nonlinear, ductile response that was measured. Damage states extracted from the models at the same load level as test measurements are overlaid to show directly the similarities and differences between test and analysis results. Studying the damage evolution predicted by the analysis reveals that the failure process is a competition between fiber damage and delamination/sub-laminate buckling, with fiber damage dominating the collapse in the hard laminate and sub-laminate buckling governing in the soft laminate. Finally, a series of parametric studies varying numerical solution parameters (mesh size, mass scaling) and physical properties (fiber direction compressive strength and toughness) reveal sensitivities and deficiencies of the model. To the authors' knowledge, this study is the first for OHC specimens to include detailed evaluation of damage mode interactions, direct overlay of predicted and measured damage states, and sensitivity of the predicted results to difficult-to-measure fiber direction material properties.
A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.
Due to the failure of thermodynamics for low temperature near-extremal black holes, it has long been conjectured that a 'thermodynamic mass gap' exists between an extremal black hole and the lightest near-extremal state. For non-supersymmetric near-extremal black holes in Einstein gravity with an AdS2 throat, no such gap was found. Rather, at that energy scale, the spectrum exhibits a continuum of states, up to non-perturbative corrections. In this paper, we compute the partition function of near-BPS black holes in supergravity where the emergent, broken, symmetry is PSU(1, 1|2). To reliably compute this partition function, we show that the gravitational path integral can be reduced to that of a $\mathcal{N}=4$ supersymmetric extension of the Schwarzian theory, which we define and exactly quantize. In contrast to the non-supersymmetric case, we find that black holes in supergravity have a mass gap and a large extremal black hole degeneracy consistent with the Bekenstein–Hawking area. Furthermore, our results verify a plethora of string theory conjectures, concerning the scale of the mass gap and the counting of extremal micro-states.
We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.
The development of a global, 3-D magnetohydrodynamic computational model is described. This model is used to quantitatively describe the detailed continuum field and plasma interaction process of the solar wind with cometary and planetary magneto/ionopause shapes. The present solar wind/terrestrial planet interaction (Level 1) model (which is based on an axisymmetric gas dynamic plus frozen field approximation to the full MHD equations) is extended to a full 3-D gas dynamic (Level 2) approximation. A mass loading capability in the Level 1 interaction model is implemented.
The objective of this research project was to develop robust and efficient computational tools for deformation and life analysis techniques. This included: development of viscoplastic continuum models, associated numerical integration techniques, suitable micromechanics models/methods, and coupled deformation and damage algorithms. The work was involved with models/methods applicable to composite structures on both microscale and macroscale levels. The microscale analysis utilized Aboudi's generalized method of cells (GMS) micromechanical model as well as finite element models of similar composite architectures, i.e., square and hexagonal fiber pack models.
With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies