Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium reconstruction”

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 127 records · Page 7

Assessment of Pioneer Venus Entry Heating with Coupled Radiation and Ablation

The objective of this work is to reconstruct thermocouple temperature measurements taken at the stagnation point of the Pioneer Venus large probe using uncoupled and coupled ablation calculations. The time history of the large probe forebody entry heating was obtained from a trajectory-based analysis that employed a material response solver to model equilibrium and nonequilibrium ablation chemistry at the surface. In the uncoupled approach, equilibrium ablation chemistry was assumed in the material response calculations, but the material response outputs were not returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer due to ablation (e.g., convective blockage and radiative energy absorption) did not influence the flowfield solutions. In the coupled approach, nonequilibrium ablation chemistry was performed and the material response outputs were returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer were accounted for in the flowfield solutions. Both uncoupled and coupled approaches produced reconstructed temperatures that were within 10% of the flight measurements through the entry heat pulse, with the coupled results comparing more favorably with the measurements. Overall, the reconstructed thermocouple temperatures in this work significantly improved on the results of previous reconstructions for the large probe.

Jannuel V V Cabrera↗

Super Resolving Unrolled Neural Networks for Remote Sensing

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

47 OTHER INSTRUMENTATION↗

Non-linear enhancement of ultrafast X-ray diffraction through transient resonances

Diffraction-before-destruction imaging with ultrashort X-ray pulses can visualize non-equilibrium processes, such as chemical reactions, with sub-femtosecond precision in the native environment. Here, a nanospecimen diffracts a single X-ray flash before it disintegrates. The sample structure can be reconstructed from the coherent diffraction image (CDI). State-of-the-art X-ray snapshots lack high spatial resolution because of weak diffraction signal. Bleaching effects from photo-ionization significantly restrain image brightness scaling. We find that non-linear transient ion resonances can overcome this barrier if X-ray laser pulses are shorter than in most experiments. We compared snapshots from individual ≈ 100 nm Xe nanoparticles as a function of pulse duration and incoming X-ray fluence. Our experimental results and Monte Carlo simulations suggest that transient resonances can increase ionic scattering cross sections significantly beyond literature values. This provides a novel avenue towards substantial improvement of the spatial resolution in CDI in combination with sub-femtosecond temporal precision at the nanoscale.

47 OTHER INSTRUMENTATION↗

The Paradox of a Wet (High H2O) and Dry (Low H2O/Ce) Mantle: High Water Concentrations in Mantle Garnet Pyroxenites from Hawaii

Water dissolved as trace amounts in anhydrous minerals has a large influence on the melting behavior and physical properties of the mantle. The water concentration of the oceanic mantle is inferred from the analyses of Mid-Ocean Ridge Basalt (MORB) and Oceanic Island Basalt (OIB). but there is little data from actual mantle samples. Moreover, enriched mineralogies (pyroxenites, eclogites) are thought as important sources of heterogeneity in the mantle, but their water concentrations and their effect on the water budget and cycling in the mantle are virtually unknown. Here, we analyzed by FTIR water in garnet clinopyroxenite xenoliths from Salt Lake Crater, Oahu, Hawaii. These pyroxenites are high-pressure (>20kb) crystal fractionates from alkalic melts. The clinopyroxenes (cpx) have 260 to 576 ppm wt H2O, with the least differentiated samples (Mg#>0.8) in the 400-500 ppm range. Orthopyroxene (opx) contain 117-265 ppm H2O, about half of that of cpx, consistent with other natural sample studies, but lower than cpx/opx equilibrium from experimental data. The pyroxenite cpx and opx H2O concentrations are at the high-end of on-and off-craton peridotite xenolith concentrations and those of Hawaiian spinel peridotites. In contrast, garnet has extremely low water contents (<5ppm H2O). There is no correlation between H2O in cpx and lithophile element concentrations. Phlogopite is present in some samples, and its modal abundance shows a positive correlation in Mg# with cpx, implying equilibrium. However, there is no correlation between H2O concentrations and or the presence of phlogopite. These data imply that cpx and opx may be at water saturation, far lower than experimental data suggest. Reconstructed bulk rock pyroxenite H2O ranges from 200-460 ppm (average 331 +/- 75 ppm), 2 to 8 times higher than H2O estimates for the MORB source (50-200 ppm), but in the range of E-MORB, OIB and the source of rejuvenated Hawaiian magmas. The average bulk rock pyroxenite H2O/Ce is 69 +/-35, lower than estimates of the MORB source (approx 150) or FOZO, C (200-250) mantle component, but consistent with "dry" EM sources (<100). These data suggest that a metasomatized, refertilized oceanic lithosphere that contains pyroxenitic veins (e.g. the lower part of an oceanic plate, where ascending melts can become trapped and crystallize), will have both higher water concentrations and low H2O/Ce, and may contribute to EM-type OIB sources, like that of Samoa basalts. Therefore, a low H2O/Ce mantle source may not necessarily be "dry".

Peslier, Anne H.↗

Evidence for Dry Carbonatite Metasomatism in the Oceanic Lithosphere from Peridotite Xenoliths of Samoa and Lanzarote

Water in Earth’s mantle affects processes like magmatism and plate tectonics. Experiments show that CO2-rich fluids lower the water solubility in olivine, implying that CO2-rich melts/fluids may dehydrate the lithosphere during metasomatism. To test this hypothesis, we report water concentrations (by polarized FTIR) of olivines, orthopyroxenes (OPX) and clinopyroxenes (CPX) from Savai’i (Samoa) and Lanzarote (Canary Islands) peridotite xenoliths with evidence of carbonatite metasomatism. Savai’i peridotites are highly depleted harzburgites and dunites with spinel Cr# (Cr/(Cr+Al)) ranging from 0.4 to 0.76 (estimated degree of melting: 191.5%). Strong Light Rare Earth Element (LREE) enrichments with Ti and Zr depletions in OPX and CO2-rich fluid inclusions (via Raman spectroscopy) are consistent with carbonatite metasomatism. Olivine, OPX and reconstructed bulk rock water concentrations (0.67-3.8, 17-89 and 4-26 ppm H2O, respectively) are low and show no apparent relationship with extent of carbonatite metasomatism. Calculated water concentrations of melts in equilibrium with Savai’i OPX (OPX/melt partitioning of water 0.0063 to 0.011) are, on average (0.540.32 wt% H2O), lower than host Samoan lavas (0.63 to 1.5 wt% H2O), despite the LREE enrichments in OPX. Lanzarote peridotites are also highly depleted (degree of melting from spinel Cr#: 171.8%).Water concentrations are low in olivines (1.7-5.3 ppm H2O) and variable in pyroxenes (OPX: 42-103 ppm H2O; CPX: 105-301 ppm H2O), and show no apparent correlation with indicators of carbonatite metasomatism. Both Savai’i and Lanzarote peridotites show negative correlations between water and degree of melting (i.e. Mg/(Mg+Fe), Cr#), suggesting melt depletion rather than metasomatism may have influenced their water concentrations. Calculated water concentrations of melts in equilibrium with Lanzarote CPX (average 1.90.75 wt% H2O; CPX/melt partitioning of water 0.011 to 0.012) are similar to those for Western Canaries lavas (average 1.80.31 wt%; CPX/melt partitioning of water 0.016 to 0.021) inferred from their CPX phenocrysts. However, calculated Ce concentrations in such melts (352 to 378 ppm; CPX/melt partitioning of Ce 0.07) are an order of magnitude greater than the lavas, and similar to carbonatites. This leads to H2O/Ce to be an order of magnitude lower in the inferred melts (26 to 57) than estimates for Western Canary lavas (280150). These low H2O/Ce ratios may suggest H2O loss from CPX during ascent, but the lack of strong water diffusion gradients in Lanzarote minerals does not support this. Instead we hypothesize that carbonatite metasomatism resulted in greater enrichment of Ce over H2O. Assuming carbonatite magmas are water rich, this implies a lower partitioning of water between minerals and melts during metasomatism, as suggested by experiments. Our data suggests carbonatite metasomatism does not result in significant re-hydration of the lithosphere, in contrast to silicate metasomatism as previously observed in Hawaiian peridotites.

Ashley, Aaron↗

Multi-probe analysis of the galaxy cluster CL J1226.9+3332: Hydrostatic mass and hydrostatic-to-lensing bias

The precise estimation of the mass of galaxy clusters is a major issue for cosmology. Large galaxy cluster surveys rely on scaling laws that relate cluster observables to their masses. From the high-resolution observations of ~45 galaxy clusters with the NIKA2 and XMM-Newton instruments, the NIKA2 Sunyaev-Zel’dovich Large Program should provide an accurate scaling relation between the thermal Sunyaev-Zel’dovich effect and the hydrostatic mass. In this paper we present an exhaustive analysis of the hydrostatic mass of the well-known galaxy cluster CL J1226.9+3332, the highest-redshift cluster in the NIKA2 Sunyaev-Zel’dovich Large Program at z = 0.89. We combined the NIKA2 observations with thermal Sunyaev-Zel’dovich data from the NIKA, Bolocam, and MUSTANG instruments and XMM-Newton X-ray observations, and tested the impact of the systematic effects on the mass reconstruction. We conclude that slight differences in the shape of the mass profile can be crucial when defining the integrated mass at R 500 , which demonstrates the importance of the modelling in the mass determination. We prove the robustness of our hydrostatic mass estimates by showing the agreement with all the results found in the literature. Another key factor for cosmology is the bias of the masses estimated assuming the hydrostatic equilibrium hypothesis. Based on the lensing convergence maps from the Cluster Lensing And Supernova survey with Hubble (CLASH) data, we obtain the lensing mass estimate for CL J1226.9+3332. From this we are able to measure the hydrostatic-to-lensing mass bias for this cluster, which spans from 1–b HSE/lens ~0.7 to 1, presenting the impact of data sets and mass reconstruction models on the bias.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessing the Impact of Laurentide Ice-sheet Topography on Glacial Climate

Simulations of past climates require altered boundary conditions to account for known shifts in the Earth system. For the Last Glacial Maximum (LGM) and subsequent deglaciation, the existence of large Northern Hemisphere ice sheets caused profound changes in surface topography and albedo. While ice-sheet extent is fairly well known, numerous conflicting reconstructions of ice-sheet topography suggest that precision in this boundary condition is lacking. Here we use a high-resolution and oxygen-isotopeenabled fully coupled global circulation model (GCM) (GISS ModelE2-R), along with two different reconstructions of the Laurentide Ice Sheet (LIS) that provide maximum and minimum estimates of LIS elevation, to assess the range of climate variability in response to uncertainty in this boundary condition.We present this comparison at two equilibrium time slices: the LGM, when differences in ice-sheet topography are maximized, and 14 ka, when differences in maximum ice-sheet height are smaller but still exist. Overall, we find significant differences in the climate response to LIS topography, with the larger LIS resulting in enhanced Atlantic Meridional Overturning Circulation and warmer surface air temperatures, particularly over northeastern Asia and the North Pacific. These up- and downstream effects are associated with differences in the development of planetary waves in the upper atmosphere, with the larger LIS resulting in a weaker trough over northeastern Asia that leads to the warmer temperatures and decreased albedo from snow and sea-ice cover. Differences between the 14 ka simulations are similar in spatial extent but smaller in magnitude, suggesting that climate is responding primarily to the larger difference in maximum LIS elevation in the LGM simulations. These results suggest that such uncertainty in ice-sheet boundary conditions alone may significantly impact the results of paleoclimate simulations and their ability to successfully simulate past climates, with implications for estimating climate sensitivity to greenhouse gas forcing utilizing past climate states.

glaciology↗

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

Inference of response functions with the help of machine-learning algorithms

Response functions are a key quantity to describe the near-equilibrium dynamics of strongly interacting many-body systems. Recent techniques that attempt to overcome the challenges of calculating these ab initio have employed expansions in terms of orthogonal polynomials. We employ a neural network prediction algorithm to reconstruct a response function 𝑆⁡(𝜔) defined over a range in frequencies 𝜔. Here, we represent the calculated response function as a truncated Chebyshev series whose coefficients can be optimized to reduce the representation error. We compare the quality of response functions obtained using coefficients calculated using a neural network (NN) algorithm with those computed using the Gaussian integral transform (GIT) method. In the regime where only a small number of terms in the Chebyshev series are retained, we find that the NN scheme outperforms the GIT method.

Kurkcuoglu, Doga Murat [Fermi National Accelerator↗

Operando characterization and regulation of metal dissolution and redeposition dynamics near battery electrode surface

Mn dissolution has been a longstanding, ubiquitous issue that negatively impacts the performance of Mn-based battery materials. Mn dissolution involves complex chemical and structural transformations at the electrode-electrolyte interface. The continuously evolving electrodeelectrolyte interface has posed great challenges for characterizing the dynamic interfacial process and quantitatively establishing the correlation with battery performance. In this study, we visualize and quantify the temporally and spatially resolved Mn dissolution/redeposition (D/R) dynamics of electrochemically operating Mn-containing cathodes. The particle-level and electrode-level analyses reveal that the D/R dynamics is associated with distinct interfacial degradation mechanisms at different states of charge. Our results statistically differentiate the contributions of surface reconstruction and Jahn-Teller distortion to the Mn dissolution at different operating voltages. Introducing sulfonated polymers (Nafion) into composite electrodes can modulate the D/R dynamics through trapping the dissolved Mn species and rapidly establishing the local Mn D/R equilibrium. This work represents an inaugural effort to pinpoint the chemical and structural transformations responsible for Mn dissolution via an operando synchrotron study and develops an effective method to regulate Mn interfacial dynamics for improving battery performance.

25 ENERGY STORAGE↗

Measuring Galactic dark matter through unsupervised machine learning

ABSTRACT Measuring the density profile of dark matter in the Solar neighbourhood has important implications for both dark matter theory and experiment. In this work, we apply autoregressive flows to stars from a realistic simulation of a Milky Way-type galaxy to learn – in an unsupervised way – the stellar phase space density and its derivatives. With these as inputs, and under the assumption of dynamic equilibrium, the gravitational acceleration field and mass density can be calculated directly from the Boltzmann equation without the need to assume either cylindrical symmetry or specific functional forms for the galaxy’s mass density. We demonstrate our approach can accurately reconstruct the mass density and acceleration profiles of the simulated galaxy, even in the presence of Gaia-like errors in the kinematic measurements.

Astronomy & Astrophysics↗

Light element ($$Z=1,2$$) production from spontaneous ternary fission of $$^{252}$$Cf

Abstract The yields of light elements ( $$Z=1,2$$ Z = 1 , 2 ) obtained from spontaneous ternary fission of $$^{252}$$ 252 Cf are treated within a nonequilibrium approach, and the contribution of unstable nuclei and excited bound states is taken into account. These light cluster yields may be used to probe dense matter, and to infer in-medium corrections such as Pauli blocking which is determined by the nucleon density. Continuum correlations are calculated from scattering phase shifts using the Beth-Uhlenbeck formula, and the effect of medium modification is estimated. The relevant distribution is reconstructed from the measured yields of isotopes. This describes the state of the nucleon system at scission and cluster formation, using only three Lagrange parameters which are the nonequilibrium counterparts of the temperature and chemical potentials, as defined in thermodynamic equilibrium. We concluded that a simple nuclear statistical equilibrium model neglecting continuum correlations and medium effects is not able to describe the measured distribution of H and He isotopes. Moreover, the freeze-out concept may serve as an important ingredient to the nonequilibrium approach using the relevant statistical operator concept.

Röpke, G.↗

Does sedimentary organic delta 13C record variations in quaternary ocean [CO2(aq)]?

Ocean surface water [CO2(aq)] variations based on glacial/interglacial changes in sediment delta 13Corg are shown to compare favorably with reconstructions based on ice core [CO2]. In particular, an approximate 80 microatmospheres increase in atmospheric pCO2 during the last glacial-interglacial transition is calculated to correspond to a 3-4 micromolar increase in ocean surface water [CO2(aq)] at atmospheric equilibrium. A widespread marine delta 13Corg decrease of 1-2% accompanied this event and was not preceded by an equivalent isotopic change in surface water total dissolved inorganic carbon. These observations support the hypothesis that [CO2(aq)] influences photosynthetic isotope fractionation between marine inorganic and organic carbon pools, and therefore that plankton/sediment delta 13Corg may serve as a proxy for surface water [CO2(aq)].

NASA Program Exobiology↗

InSight Aerothermal Environment Assessment

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

Beck, Robin↗

InSight's Reconstructed Aerothermal Environments

The InSight spacecraft was proposed to be a build-to-print copy of the Phoenix vehicle due to the knowledge that the lander payload would be similar and the trajectory would be similar. However, the InSight aerothermal analysts, based on tests performed in CO2 during the Mars Science Laboratory mission (MSL) and completion of Russian databases, considered radiative heat flux to the aftbody from the wake for the first time for a US Mars mission. The combined convective and radiative heat flux was used to determine if the as-flown Phoenix thermal protection system (TPS) design would be sufficient for InSight. All analyses showed that the design would be adequate. Once the InSight lander was successfully delivered to Mars on November 26, 2018, work began to reconstruct the atmosphere and trajectory in order to evaluate the aerothermal environments that were actually encountered by the spacecraft and to compare them to the design environments.The best estimated trajectory (BET) reconstructed for the InSight atmospheric entry fell between the two trajectories considered for the design, when looking at the velocity versus altitude values. The maximum heat rate design trajectory (MHR) flew at a higher velocity and the maximum heat load design trajectory (MHL) flew at a lower velocity than the BET. For TPS sizing, the MHL trajectory drove the design. Reconstruction has shown that the BET flew for a shorter time than either of the design environments, hence total heat load on the vehicle should have been less than used in design. Utilizing the BET, both DPLR and LAURA were first run to analyze the convective heating on the vehicle with no angle of attack. Both codes were run with axisymmetric, laminar flow in radiative equilibrium and vibrational non-equilibrium with a surface emissivity of 0.8. Eight species Mitcheltree chemistry was assumed with CO2, CO, N2, O2, NO, C, N, and O. Both codes agreed within 1% on the forebody and had the expected differences on the aftbody. The NEQAIR and HARA codes were used to analyze the radiative heating on the vehicle using full spherical ray-tracing. The codes agreed within 5% on most aftbody points of interest.The LAURA code was then used to evaluate the conditions at angle of attack at the peak heating and peak pressure times. Boundary layer properties were investigated to confirm that the flow over the forebody was laminar for the flight.Comparisons of the aerothermal heating determined for the reconstructed trajectory to the design trajectories showed that the as-flown conditions were less severe than design

radiation↗

Elevated Electron Temperature Coincident with Observed Fusion Reactions in a Sheared-Flow-Stabilized Z Pinch

The sheared-flow-stabilized Z pinch concept has been studied extensively and is able to produce fusion-relevant plasma parameters along with neutron production over several microseconds. In this work, we present elevated electron temperature results spatially and temporally coincident with the plasma neutron source. An optical Thomson scattering apparatus designed for the FuZE device measures temperatures in the range of 1–3 keV on the axis of the device, 20 cm downstream of the nose cone. The 17-fiber system measures the radial profiles of the electron temperature. Scanning the laser time with respect to the neutron pulse time over a series of discharges allows the reconstruction of the T e temporal response, confirming that the electron temperature peaks simultaneously with the neutron output, as well as the pinch current and inductive voltage generated within the plasma. Comparison to spectroscopic ion temperature measurements suggests a plasma in thermal equilibrium. The elevated T e confirms the presence of a plasma assembled on axis, and indicates limited radiative losses, demonstrating a basis for scaling this device toward net gain fusion conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Myths of German Graphite in World War II, with Original Translations

We re-examine the common narrative that a 1941 experimental error by physicists Walther Bothe and Peter Jensen led Germany to abandon graphite as a reactor moderator during World War II. We first detail the history of both German and American graphite experiments, noting that the Americans faced similar setbacks but succeeded only at the end of a costly 18-month graphite purification program. We then use Monte Carlo N-Particle simulations to reconstruct Bothe's 1941 experiment. We find the thermal absorption cross section of Bothe's Siemens electrographite to be 12.2 mb, in contrast to his reported 7.9 mb. This discrepancy arises because the neutrons in Bothe's experiment did not reach thermal equilibrium, leading to an underestimation of neutron absorption. Additionally, despite misconceptions that the Germans were unaware of boron impurities, we share evidence that Wilhelm Hanle accurately measured boron and cadmium impurities in the electrographite. To support our findings, we provide 9 excerpted and 3 complete English translations of classified wartime reports by Heisenberg, Joos, Bothe, Jensen, Höcker, Hanle, and Kremer. Our work aims to illuminate the rational constraints behind Germany's decision to forgo graphite moderation.

Park, Patrick J.↗