Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Breakdown of electroneutrality in nanopores

Ion transport in extremely narrow nanochannels has gained increasing interest in recent years due to unique physical properties at the nanoscale and the technological advances that allow us to study them. It is tempting to approach this confined regime with the theoretical tools and knowledge developed for membranes and microfluidic devices, and naively apply continuum models, such as the Poisson-Nernst-Planck and Navier-Stokes equations. However, it turns out that some of the most basic principles we take for granted in larger systems, such as the complete screening of surface charge by counter-ions, can break down under extreme confinement. Additionally, we show that in a truly one-dimensional system of ions interacting with three-dimensional electrostatic interactions, the screening length is exponentially large, and can easily exceed the macroscopic length of a nanotube. Without screening, electroneutrality breaks down within the nanotube, with fundamental consequences for ion transport and electrokinetic phenomena. In this work, we build a general theoretical framework for electroneutrality breakdown in nanopores, focusing on the most interesting case of a one-dimensional nanotube, and show how it provides an elegant interpretation for the peculiar scaling observed in experimental measurements of ionic conductance in carbon nanotubes.

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

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2024

DECOVALEX-2023: Task G, SAFENET Final Report

DECOVALEX Task G deals with fracture mechanics at several scales using a combined approach of experimental work, related model development, benchmarking and experimental analysis. Fig. 0.1 provides a graphical abstract for Task G. The experimental basis for the related mechanical (M), hydro-mechanical (HM), and thermo-mechanical (TM) processes comes from the rock mechanics laboratories of Universities of Freiberg (TUBAF) and Edinburgh, as well as the Korea Institute of Civil Engineering and Building Technology (KICT). The experimental work in the rock laboratories is closely linked to the underground research laboratories (URLs) Reiche Zeche (Germany), KURT (Korea), and Mont Terri (Switzerland). The scientific key questions are related to fracture permeability evolution under THM conditions, anisotropy effects on fracturing processes, and thermal fracture slip, which are being addressed in specific steps of the task. A large variety of numerical methods have been developed and applied for experimental analysis, ranging from continuum to discontinuum approaches. A detailed comparison of mechanical and hydro-mechanical processes is given in the section 3 via the benchmarking exercises. As a result of DECOVALEX-2023 Task G we further improved our understanding and predictability of fracturing processes under THM conditions. This was based on robust numerical simulation methods and in-depth experimental analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Ultraviolet Imaging Telescope ultraviolet images - Large-scale structure, H II regions, and extinction in M81

The study employs UV images of M81 obtained by the Ultraviolet Imaging Telescope (UIT) during the December 1990 Astro-1 spacelab mission to determine 2490- and 1520-A fluxes from 46 H II regions and global surface brightness profiles. Comparison photometry in the V band is obtained from a ground-based CCD image. UV radial profiles show bulge and exponential disk components, with a local decrease in disk surface brightness inside the inner Lindblad Resonance about 4 arcmin from the nucleus. The V profile shows typical bulge plus exponential disk structure, with no local maximum in the disk. There is little change of UV color across the disk, although there is a strong gradient in the bulge. Observed m152-V colors of the H II regions are consistent with model spectra for young clusters, after dereddening using Av determined from m249-V and the Galactic extinction curve. The value of Av, so determined, is 0.4 mag greater on the average than Av derived from radio continuum and H-alpha fluxes.

Hill, Jesse K.↗
DOE OSTI2020

Tensor network formulation of two-dimensional gravity

We show how to formulate a lattice gauge theory whose naive continuum limit corresponds to two-dimensional (Euclidean) quantum gravity including a positive cosmological constant. More precisely the resultant continuum theory corresponds to gravity in a first-order formalism in which the local frame and spin connection are treated as independent fields. Recasting this lattice theory as a tensor network allows us to study the theory at strong coupling without encountering a sign problem. In two dimensions this tensor network is exactly soluble and we show that the system has a series of critical points that occur for pure imaginary coupling and are associated with first order phase transitions. We then augment the action with a Yang-Mills term which allows us to control the lattice spacing and show how to apply the tensor renormalization group to compute the free energy and look for critical behavior. Finally we perform an analytic continuation in the gravity coupling in this extended model and show that its critical behavior in a certain scaling limit depends only on the topology of the underlying lattice. We also show how the lattice gauge theory can be naturally generalized to generate the Polyakov or Liouville action for two dimensional quantum gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2010

On the Nature of the Compact Object in SS 433. Observational Evidence of X-Ray Photon Index Saturation

We present an analysis of the X-ray spectral properties observed from black hole , candidate (BHC) binary SS 433. We have analyzed Rossi X-ray Time Explorer (RXTE) data from this source, coordinated with Green Bank Interferometer/RATAN-600. We show that SS 433 undergoes a X-ray spectral transition from the low hard state (LHS) to the intermediate state (IS). We show that the X-ray broad-band energy spectra during all spectral states are well fit by a sum of so called "Bulk Motion Comptonization (BMC) component" and by two (broad and narrow) Gaussians for the continuum and line emissions respectively. In addition to these spectral model components we also find a strong feature that we identify as a" blackbody-like (BB)" component which color temperature is in the range of 4-5 keV in 24 IS spectra during the radio outburst decay in SS 433. Our observational results on the "high temperature BB" bump leads us to suggest the presence of gravitationally redshifted annihilation line emission in this source. In fact this spectral feature has been recently reproduced in Monte Carlo simulations by Laurent and Titarchuk. We have also established the photon index saturation at about 2.3 in index vs mass accretion correlation. This index-mass accretion correlation allows us to evaluate the low limit of black hole (BH) mass of compact object in SS 433, M(sub bh) approximately > 2 solar masses, using the scaling method using BHC GX 339-4 as a reference source. Our estimate of the BH mass in SS 433 is consistent with recent BH mass measurement using the radial-velocity measurements of the binary system by Hillwig & Gies who find that M(sub x)( = (4.3 +/- 0.8) solar masses. This is the smallest BH mass found up to now among all BH sources. Moreover, the index saturation effect versus mass accretion rate revealed in SS 433, like in a number of other BH candidates, is the strong observational evidence for the presence of a BH in SS 433.

Seifina, Elena↗
NASA NTRS2018

Morphology and Kinematics of Filaments in the Serpens and Perseus Molecular Clouds

We present H(exp 13)CO(exp +) (J= 1-0) and HNC (J= 1-0) maps of regions in Serpens South, Serpens Main, and NGC 1333 containing filaments. We also observe the Serpens regions using H(exp 13)CN (J= 1-0). These dense gas tracer molecular line observations carried out with CARMA have an angular resolution of approximately 7", a spectral resolution of approximately 0.16 km s (exp -1), and a sensitivity of 50-100 mJy beam (exp -1). Although the large-scale structure compares well with the Herschel dust continuum maps, we resolve finer structure within the filaments identified by Herschel. The H(exp 13)CO(exp +) emission distribution agrees with the existing CARMA diazenylium(exp +) (J= 1-0) maps, so they trace the same morphology and kinematics of the filaments. The H(exp 13)CO(exp +) maps additionally reveal that many regions have multiple structures partially overlapping in the line of sight. In two regions, the velocity differences are as high as 1.4 km s(exp -1). We identify eight filamentary structures having typical widths of 0.03 - 0.08 pc in these tracers. At least 50% of the filamentary structures have distinct velocity gradients perpendicular to their major axis, with average values in the range of 4-10 km s(exp -1) pc (exp -1). These findings are in support of the theoretical models of filament formation by 2D inflow in the shock layer created by colliding turbulent cells. We also find evidence of velocity gradients along the length of two filamentary structures; the gradients suggest that these filaments are inflowing toward the cloud core.

Dhabal, Arnab↗
NASA NTRS2022

An ALMA Study of the Massive Molecular Clump N159W-North in the Large Magellanic Cloud: A Possible Gas Flow Penetrating One of the Most Massive Protocluster Systems in the Local Group

Massive dense clumps in the Large Magellanic Cloud can be an important laboratory to explore the formation ofpopulous clusters. We report multiscale ALMA observations of the N159W-North clump, which is the most CO-intense region in the galaxy. High-resolution CO isotope and 1.3 mm continuum observations with an angularresolution of∼025(∼0.07 pc)revealed more thanfive protostellar sources with CO outflows within the mainridge clump. One of the thermal continuum sources, MMS-2, shows an especially massive/dense nature whosetotal H2mass and peak column density are∼104Meand∼1024cm−2, respectively, and harbors massive(∼100Me)starless core candidates identified as its internal substructures. The main ridge containing this sourcecan be categorized as one of the most massive protocluster systems in the Local Group. The CO high-resolutionobservations found several distinctfilamentary clouds extending southward from the star-forming spots. The CO(1–0)data set with a largerfield of view reveals a conical,∼30 pc long complex extending toward the northerndirection. These features indicate that a large-scale gas compression event may have produced the massive star-forming complex. Based on the striking similarity between the N159W-North complex and the other twopreviously reported high-mass star-forming clouds in the nearby regions, we propose a“teardrops inflow model”that explains the synchronized, extreme star formation across>50 pc, including one of the most massiveprotocluster clumps in the Local Group.

Kazuki Tokuda↗
DOE OSTI2021

A Minimal Information Set To Enable Verifiable Theoretical Battery Research

Batteries are an enabling technology for addressing sustainability through the electrification of various forms of transportation (1) and grid storage. (2) Batteries are truly multi-scale, multi-physics devices, and accordingly various theoretical descriptions exist to understand their behavior (3-5) ranging from atomistic details to techno-economic trends. As we explore advanced battery chemistries (6,7) or previously inaccessible aspects of existing ones, (8-10) new theories are required to drive decisions. (11-13) The decisions are influenced by the limitations of the underlying theory. Advanced theories used to understand battery phenomena are complicated and require substantial effort to reproduce. However, such constraints should not limit the insights from these theories. We can strive to make the theoretical research verifiable such that any battery stakeholder can assess the veracity of new theories, sophisticated simulations or elaborate analyses. We distinguish verifiability, which amounts to “Can I trust the results, conclusions and insights and identify the context where they are relevant?”, from reproducibility, which ensures “Would I get the same results if I followed the same steps?” With this motivation, we propose a checklist to guide future reports of theoretical battery research in Table 1. We hereafter discuss our thoughts leading to this and how it helps to consistently document necessary details while allowing complete freedom for creativity of individual researchers. Given the differences between experimental and theoretical studies, the proposed checklist differs from its experimental counterparts. (14,15) This checklist covers all flavors of theoretical battery research, ranging from atomic/molecular calculations (16-19) to mesoscale (20,21) and continuum-scale interactions, (9,22) and techno-economic analysis. (23,24) Finally, as more and more experimental studies analyze raw data, (25) we feel this checklist would be broadly relevant.

25 ENERGY STORAGE↗
DOE OSTI2024

Co-design for Particle Applications at Exascale

Co-design across the Exascale Computing Project (ECP) has been critical for both enabling science applications and bringing disparate communities together. Developing and porting applications to the various high-performance computing (HPC) architectures on pre-exascale and exascale computers has been quite challenging due to the diversity of hardware features and software stacks. The Co-design Center for Particle Applications (CoPA) has developed and enhanced the Cabana and PROGRESS/BML libraries to facilitate the creation of new particle applications, make existing particle applications exascale capable, and allow teams to explore new capabilities. Particle methods from atomistic, mesoscale, continuum, through cosmological scales have been built with Cabana, along with new possibilities for application coupling. Similarly, the PROGRESS/BML library has enabled quantum particle applications with linear algebra solvers to use advanced hardware. Across these CoPA-developed libraries, the co-design abstraction layer combines performance portability with math library support to facilitate separation of concerns and directly support science runs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

ISO Spectroscopy of Proto-Planetary Nebulae

The goal of this program was to determine the chemical properties of the dust shells around protoplanetary nebulae (PPNs) through a study of their short-wavelength (6-45 micron) infrared spectra. PPNs are evolved stars in transition from the asymptotic giant branch to the planetary nebula stages. Spectral features in the 10 to 20 gm region indicate the chemical nature (oxygen- or carbon-rich), and the strengths of the features relate to the physical properties of the shells. A few bright carbon-rich PPNs have been observed to show PAH features and an unidentified 21 micron emission feature. We used the Infrared Space Observatory (ISO) to observe a sample of IRAS sources that have the expected properties of PPNs and for which we have accurate positions. Some of these have optical counterparts (proposal SWSPPN01) and some do not (SWSPPN02). We had previously observed these from the ground with near-infrared photometry and, for those with visible counterparts, visible photometry and spectroscopy, which we have combined with these new ISO data in the interpretation of the spectra. We have completed a study of the unidentified emission feature at 21 micron in eight sources. We find the shape of the feature to be the same in all of the sources, with no evidence of any substructure. The ratio of the emission peak to continuum ranges from 0.13 to 1.30. We have completed a study of seven PPNs and two other carbon-rich objects for which we had obtained ISO 2-45 micron observations. The unidentified emission features at 21 and 30 micron were detected in six sources, including four new detections of the 30 micron feature. This previously unresolved 30 micron feature was resolved and found to consist of a broad feature peaking at 27.2 micron (the "30 micron" feature) and a narrower feature peaking at 25.5 micron (the "26 micron" feature). This new 26 micron feature is detected in eight sources and is particularly strong in IRAS Z02229+6208 and 16594-4656. The unidentified features at 3.3, 6,2, 7.7, and 11.3 micron, which are commonly observed in planetary nebulae and HII regions, are also seen in these PPNs. However, their strengths relative to the continuum plateaus at 8 and 12 micron are weaker than in planetary nebulae. The 6.9 micron feature, seen almost exclusively in PPNs, is strong. The spectral energy distributions of these PPNs were fitted with a radiative-transfer model, taking into account the emission features at 21, 26, and 30 micron. A significant fraction of the total energy output is emitted in these features: as high as 20% in the 30 micron feature and 8% in the 21 micron feature. The fact that so much energy is carried in these features suggests that the material responsible for this feature must be made of abundant elements, and most likely involves carbon. The change in the in feature strengths from stronger aliphatic bonds in PPNs to stronger aromatic bonds in PNs suggests a chemical and physical evolution in the carbonaceous circumstellar dust during this transition time scale of a few thousand years.

Hrivnak, Bruce J.↗
NASA NTRS2016

The Nature of the Torus in the Heavily Obscured AGN Markarian 3: an X-Ray Study

In this paper, we report the results of an X-ray monitoring campaign on the heavily obscured Seyfert galaxy, Markarian 3, carried out between the fall of 2014 and the spring of 2015 with NuSTAR, Suzaku and XMMNewton. The hard X-ray spectrum of Markarian 3 is variable on all the time-scales probed by our campaign, down to a few days. The observed continuum variability is due to an intrinsically variable primary continuum seen in transmission through a large, but still Compton-thin column density (N(sub H) approx. 0.8-1.1 x 10(exp 24)/sq cm). If arranged in a spherical-toroidal geometry, the Compton scattering matter has an opening angle approx. 66deg, and is seen at a grazing angle through its upper rim (inclination angle approx. 70deg). We report a possible occultation event during the 2014 campaign. If the torus is constituted by a system of clouds sharing the same column density, this event allows us to constrain their number (17 +/- 5) and individual column density, [approx. (4.9 +/- 1.5) x 10(exp 22)/ sq cm]. The comparison of IR and X-ray spectroscopic results with state-of-the art torus models suggests that at least two-thirds of the X-ray obscuring gas volume might be located within the dust sublimation radius. We report also the discovery of an ionized absorber, characterized by variable resonant absorption lines due to He- and H-like iron. This discovery lends support to the idea that moderate column density absorbers could be due to clouds evaporated at the outer surface of the torus, possibly accelerated by the radiation pressure due to the central AGN emission leaking through the patchy absorber.

X-rays: individual: Markarian 3↗
NASA NTRS1986

The structure of high-temperature flare plasma. II Comparison with observations

In Emslie, (1985) scaling law's which relate peak loop temperature to injected nonthermal electron flux in electron-heated models of the flare corona are derived. These predicted relationships are compared with recent coordinated observations in hard X-rays and soft X-rays. Satisfactory agreement is found for large events, while for smaller events the plasma electron temperature determined by soft X-ray spectral fitting is too high. A possible resolution of this apparent discrepancy through careful examination of the method used to determine the temperature of the soft X-ray emitting plasma is discussed. It is concluded that temperatures determined by spectral fitting over a series of lines are not necessarily representative of the true temperature of the plasma, since the fitting technique is plagued by the same difficulties as a more straightforward technique involving a single spectral line or portion of continuum (Craig and Brown, 1976). The differences between actual and derived temperatures are sufficiently large to remove the above discrepancy between observations and modeling of electron-heated coronae.

Emslie, A. G.↗
DOE OSTI2022

Diagnosing the DARHT Electron Beam-Target Interaction and Hydrodynamic Expansion [Dissertation]

Quantitative electron beam-target interaction studies on the Dual-Axis Radiographic Hydrodynamic Test (DARHT) electron linear induction accelerators have only been performed recently. This includes characterization of the temperature, density, pressure, extent, and expansion velocity of the plasma plume. The results presented in this dissertation present a detailed and unified overview of the accelerator systems, target heating physics, beam transport, and diagnostic tools. Additional information includes calibration sources, radiation hydrodynamics, spectroscopic-quality radiation transport modeling, experimental measurements, and analyses of electron beam driven aluminum experiments. The first set of spatially and temporally resolved spectroscopic measurements of electron beam driven aluminum are presented. Contamination quantification analyses are used to understand the origin of the strong Na-I 3p-3s lines that are observed in absorption within the aluminum plasma continuum. These results inform the creation of the first spectroscopic-quality radiation transport model that links several atomic physics codes to interpret the conditions from which the Na-I lines originate. A good agreement is found between the surface analysis results and the model which confirms the concentration of the sodium present within the aluminum alloy foil material. It also demonstrates, for the first time on electron beam-driven target experiments, the ability to interpret plasma conditions from measured absorption lines. In a second experimental campaign which focuses on pure aluminum, the Al-I 3p-4s and 3p3d doublets are both measured in emission. A detailed analysis of the Al-I 3p-4s doublet reveals that the lines undergo moderate self-absorption. A simple model of the self-absorption effect is successfully used to match measured spectra at various temperature/density/plasma scale length combinations. These measurements led to the realization of the minimum density that can be resolved by the spectrometers for the Al-I 3p-4s lines due to the large slit width required to observe a signal on aluminum. These measurements also demonstrate the sensitivity of visible and long wave UV spectroscopy to minor changes in both temperature and density. The simple self-absorption model will be useful for analysis of other beam-target interaction experiments with spectra exhibiting either self-absorption or full self-reversal. Substantial headway has been made on the modeling front by linking together several codes needed to model both the energy deposition, hydrodynamic motion, and atomic kinetics to produce synthetic spectral calculations that are compared with experimental measurements. There exists ample space for improvement, especially with benchmarking the hydrodynamics codes and equation-of-state tables with experimental measurements. The X-ray diagnostics required to make these new measurements along with the simulation capabilities required to interpret the results are under development and will be the subject of future studies.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Variability of compact radio sources at a wavelength of 1 millimeter

Broad-band 1 mm continuum observations of 3C 273, 3C 279, BL Lac, 3C 84, OJ 287, and 3C 345 made between 1977 March and 1981 January are reported. With the possible exception of 3C 279, all of the objects exhibited evidence for variability in the 1 mm flux density on the time scale of a few months. Flux variations at wavelengths of 1 mm and 2 cm are well correlated for the 'blazars' (BL Lac objects and optically violent variable quasars) observed; emission outbursts occur simultaneously and have similar amplitude at the two wavelengths. In contrast, an outburst in the flux from the Seyfert galaxy 3C 84 was observed at a wavelength of 1 mm in early 1980 but has not appeared at 2 cm as of 1 yr later. The observed relationship between 'blazar' emission variability at 1 mm and 2 cm is inconsistent with the canonical expanding-source model of radio variability or any of its modifications. Instead, models in which the number of synchrotron-radiating electrons is time variable (due to particle injection or acceleration behind a relativistic shock) are required.

Ennis, D. J.↗
DOE OSTI2025

Toward scalable bound-to-resonance extrapolations for few- and many-body systems

In open quantum many-body systems, the theoretical description of resonant states of many particles strongly coupled to the continuum can be challenging. Such states are commonplace in, for example, exotic nuclei and hadrons, and can reveal important information about the underlying forces at play in these systems. In this work, we demonstrate that the complex-augmented eigenvector continuation (CA-EC) method, originally formulated for the two-body problem with uniform complex scaling, can reliably perform bound-to-resonance extrapolations for genuine three-body resonances having no bound subsystems. Here, we first establish that three-body bound-to-resonance extrapolations are possible by benchmarking different few-body approaches, and we provide arguments to explain how the extrapolation works in the many-body case. We furthermore pave the way towards scalable resonance extrapolations in many-body systems by showing that the CA-EC method also works in the Berggren basis, studying a realistic application using the Gamow shell model.

Ab initio calculations↗
NASA NTRSDate not supplied

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield that used NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performance for protecting probes during planetary entry. Data from the Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. A space-grade silicone-based coating commercially known as NuSil CV-1144-0 [3] was applied to the entire MSL heatshield, including the MEDLI plugs, to mitigate the spread of dust from PICA. Modeling the thermal response of PICA-NuSil (PICA-N) system is still an open challenge. Ground testing of PICA-N models exhibited surface temperature jumps of the order of 150 K due to oxide scale formation and sub-sequent NuSil burn-off. It is therefore critical to include a validated model for the material response of the coating in engineering codes. A test campaign has been conducted at the NASA’s Langley HyMETS [4] facility to screen the response of PICA-N and gather detailed data on its behavior [5]. A first model of PICA-N thermal response has been developed using the Hy-METS experiments [6]. The objective of this work is to analyze the material response of the latest PICA-N model compared to the engineering model used to simulate the entry of MSL. The environment and material response around the MSL aeroshell during Mars atmospheric entry is simulated using a collection of tools. The Direct Simulation Monte Carlo SPARTA code [7] is used in the rarefied regime, the Data Parallel Line Relaxation (DPLR) code [8] is used in the continuum regime and radiative heating conditions are provided by the Nonequilibrium air radiation (NEQAIR) code [9] to estimate the environmental conditions. The thermal response inside the material is computed using the Porous material Analysis Toolbox based on Open-FOAM (PATO) [10,11,12]. Thermodynamic and chemistry properties are estimated using the Mutation++ library [13]. The approach implemented in PATO as a first cut PICA-N thermal response model is outlined in Figure 1. While the recession is less than the coating thickness, the Surface mass and energy balance Boundary Condition (SBC) uses the NuSil B’ tables. Once the recession removes the coating, the usual PICA B’ tables are used for the SBC. The B’ tables are computed using an equilibrium solver implemented in Mutation++, given the temperature, pressure, blowing rate, composition of the pyrolysis and environment gases, and the condensed species at the surface. Preliminary results of the 3D material response of the MSL heat-shield at the peak heating (80 sec after Entry Interface) are shown in Figure 2. Current NASA’s mission to Mars, Mars 2020, used the spare heatshield of MSL for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the PATO high-fidelity material response capability was benchmarked against flight data from MEDLI. This effort represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes. This bench-marking is awaiting the final release of the MEDLI-2 data.

Aerospace↗