SEARCH · Search NASA

Results for “High Energy Density Physics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Dark Matter in the Universe and in the Galaxy

During the past four years, Prof. Kamionkowski and collaborators have made progress in research on the nature and distribution of dark-matter in the Universe and in the Galaxy, and on related topics in astrophysics and cosmology. We have made progress on research on the cosmic microwave background, large-scale structure, issues related to particle dark matter, and the gamma-ray-burst enigma. A significant fraction of the research supported by this ATP has been on the cosmic microwave background (CMB). Prof. Kamionkowski and collaborators showed how the polarization of the CMB could be used to detect long-wavelength gravitational waves, such as those produced by inflation. With Kosowsky, Prof. Kamionkowski calculated the amplitude of a stochastic gravitational-wave background that could be detected for a satellite experiment of a given sensitivity and angular resolution. They showed that polarization should improve the sensitivity oa MAP to these gravity waves, and that the Planck Surveyor should do even better. Prof. Kamionkowski, Caldwell, and a student calculated and illustrated the CMB temperature/polarization pattern produced by a single plane-wave gravitational wave. They calculated the amplitude of such a wave that would be detectable with MAP and Planck, and compared that with the sensitivity of traditional gravitational-wave detectors like LIGO and LISA. With Lue and Wang, the PI showed how parity violation from new high-energy physics could conceivably give rise to an observable signature in the CMB polarization. With Loeb, Prof. Kamionkowski showed how measurement of the polarization of CMB photons scattered by hot gas in a cluster could be used to determine the quadrupole moment of the CMB incident on that cluster. Prof. Kamionkowski and Jaffe calculated the amplitude of secondary anisotropies produced by scattering of CMB photons from reionized regions. Research has also been carried out on probing the large-scale distribution of mass in the Universe today, and on structure-formation theories. They investigated the possibility of determining the large-scale distribution of mass in the Universe via measurement of ellipticity-ellipticity correlations in the FIRST radio survey induced by weak gravitational lensing due to mass inhomogeneities along the line of sight. Dr. Summers, Prof. Kamionkowski, and a student investigated the distribution of protogalactic masses and angular momenta in an effort to understand how the luminosity function and angular-momentum distribution of disk galaxies arises from an initial power spectrum of density perturbations. Several projects related to the distribution and possible detection of dark matter in our Galactic halo were studied. The PI was involved in several projects involving the calculation of nuclear-reaction rates needed for stellar evolution and for predictions of solar-neutrino fluxes.

Kamionkowski, Marc↗

NASA NTRS2013

High-Strength, Superelastic Compounds

In a previous disclosure, the use of 60- NiTiNOL, an ordered intermetallic compound composed of 60 weight percent nickel and 40 weight percent titanium, was investigated as a material for advanced aerospace bearings due to its unique combination of physical properties. Lessons learned during the development of applications for this material have led to the discovery that, with the addition of a ternary element, the resulting material can be thermally processed at a lower temperature to attain the same desirable hardness level as the original material. Processing at a lower temperature is beneficial, not only because it reduces processing costs from energy consumption, but because it also significantly reduces the possibility of quench cracking and thermal distortion, which have been problematic with the original material. A family of ternary substitutions has been identified, including Hf and Zr in various atomic percentages with varying concentrations of Ni and Ti. In the present innovation, a ternary intermetallic compound consisting of 57.6 weight percent Ni, 39.2 weight percent Ti, and 3.2 weight percent Hf (54Ni-45Ti-1Hf atomic percent) was prepared by casting. In this material, Hf substitutes for some of the Ti atoms in the material. In an alternate embodiment of the innovation, Zr, which is close in chemical behavior to Hf, is used as the substitutional element. With either substitution, the solvus temperature of the material is reduced, and lower temperatures can be used to obtain the necessary hardness values. The advantages of this innovation include the ability to solution-treat the material at a lower temperature and still achieve the required hardness for bearings (at least 50 Rockwell C) and superelastic behavior with recoverable strains greater than 2%. Most structural alloys will not return to their original shape after being deformed as little as 0.2% (a tenth of that possible with superelastic materials like 60 NiTiNOL). Because lower temperatures can be used in the heat treatment process, less energy will be consumed, and there will be less dimensional distortion and quench cracking. This results in fewer scrap parts, less material waste from large amounts of material removal, and fewer machining steps to rework parts that are out of specification. This material has a combination of properties that have been previously unobtainable. The material has a Young s modulus of approximately 95 GPa (about half that of conventional steels), moderate density (10 to 15% lower than conventional steels), excellent corrosion resistance, and high hardness (58 to 62 HRC). These properties make this material uniquely suited for advanced bearings.

Stanford, Malcolm↗

NASA NTRS1990

Direct Numerical Simulations of Turbulent Convection: Variable Gravity and Uniform Rotation - Part 1

Turbulent convection may have played a major role in determining the structure and evolution of the primordial solar nebula, but current, incomplete models of convection and turbulence give very different results and remain largely untested in the absence of detailed astronomical observations. Numerical simulations provide an 'experimental' database for comparison with these models, and, to this end, direct numerical simulations of turbulent convection were performed with modifications intended to mimic some of the unique physical features of thin accretion disks, such as the primordial solar nebula: internal heating, a gravitational acceleration that is linearly proportional to the distance from midplane of the nebula, and rapid rotation. Peclet numbers in the simulations are comparable to those in solar nebula models; Rossby numbers in the simulations are an order of magnitude larger than those in solar nebula models because of the unrealistically high Prandtl and low, Reynolds numbers required to resolve all scales of the convective flow. We find that, despite the loss of buoyancy at midplane, turbulent motions easily penetrate the midplane region with little loss of intensity, providing efficient transport of heat and turbulent kinetic energy throughout the interior. A simple mixing length model modified to include rotation is found to give convective heat fluxes for the interior flow in rough agreement with the numerical simulations. We discuss these preliminary results with regard to assumptions about heating distributions and convective heat fluxes made in standard solar nebula models. More definitive comparisons with solar nebula modelling will become possible when more realistic effects of shear, density, stratification, and compression are included.

Cabot, W.↗

NASA NTRS2019

Photoionization Emission Models for the Cyg X-3 X-Ray Spectrum

We present model fits to the X-ray line spectrum of the well-known high-mass X-ray binary Cyg X-3. The primary observational data set is a spectrum taken with the Chandra X-ray Observatory High Energy Transmission Grating in 2006, though we compare it to all the other observations of this source taken so far by this instrument. We show that the density must be 1012 cm−3 in the region responsible for most of the emission. We discuss the influence of the dust scattering halo on the broadband spectrum, and we argue that dust scattering and extinction is not the most likely origin for the narrow feature seen near the Si K edge. We identify the features of a wind in the profiles of the strong resonance lines and show that the wind is more apparent in the lines from the lighter elements. We argue that this wind is most likely associated with the companion star. We show that the intensities of most lines can be fitted, crudely, by a single-component photoionized model. However, the iron K lines do not fit with this model. We show that the iron K line variability as a function of orbital phase is different from the lower-energy lines, which indicates that the lines arise in physically distinct regions. We discuss the interpretation of these results in the context of what is known about the system and similar systems.

Kallman, T.↗

NASA NTRS2019

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗

NASA NTRS2016

A Model for Jet-Surface Interaction Noise Using Physically Realizable Upstream Turbulence Conditions

This paper is a continuation of previous work in which a generalized Rapid Distortion Theory (RDT) formulation was used to model low-frequency trailing-edge noise. The research was motivated by proposed next-generation aircraft configurations where the exhaust system is tightly integrated with the airframe. Data from recent experiments at NASA on the interaction between high-Reynolds-number subsonic jet flows and an external flat plate showed that the power spectral density (PSD) of the far-field pressure underwent considerable amplification at low frequencies. For example, at the 90deg observation angle, the low-frequency noise could be as much as 10 dB greater than the jet noise itself. In this paper, we present predictions of the noise generated by the interaction of a rectangular jet with the trailing edge of a semi-infinite flat plate. The calculations are based on a formula for the acoustic spectrum of this noise source derived from an exact formal solution of the linearized Euler equations involving (in this case) one arbitrary convected scalar quantity and a Rayleigh equation Green's function. A low-frequency asymptotic approximation for the Green's function based on a two-dimensional mean flow is used in the calculations along with a physically realizable upstream turbulence spectrum, which includes a finite decorrelation region. Numerical predictions of the sound field, based on three-dimensional RANS solutions to determine the mean flow, turbulent kinetic energy and turbulence length and time scales, for a range of subsonic acoustic Mach number jets and nozzle aspect ratios are compared with experimental data. Comparisons of the RANS results with flow data are also presented for selected cases. We find that a finite decorrelation region in the turbulence spectrum increases the low-frequency algebraic decay (the low frequency "roll-off") of the acoustic spectrum with angular frequency thereby producing much closer agreement with noise data for Strouhal numbers less than 0.1. Secondly, the large-aspect-ratio theory is able to predict the low-frequency amplification due to the jet-edge interaction reasonably well, even for moderate aspect ratio nozzles. We show also that the noise predictions for smaller aspect ratio jets can be fine-tuned using the appropriate RANS-based mean flow and turbulence properties.

Jet↗

NASA NTRS1995

GRO: Black hole models for gamma-ray bursts

The Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory (CGRO) has established that the distribution of gamma-ray bursts (GRB's) is isotropic but is bound radially. This finding suggests that the bursts are either cosmological or they originate from an extended Galactic halo. The implied luminosities and the observed variability of the GRB's on time scales as short as one millisecond suggest that they originate from compact objects. We are presently studying black hole models for GRB's. Any such model must produce a non-thermal photon spectrum to agree with the observed properties. For a wide range of burst parameters the assumed bursting source consists of a non-thermal electron-positron-photon plasma of very high density. It seems possible to produce such a plasma in accretion onto black holes. In our on-going work, we are developing the kinetic theory for a non-equilibrium pair plasma. The main new features of our work are as follows: (1) We do not assume the presence of a thermal electron bath. (2) Non-thermal, high-energy pairs are allowed to have an arbitrary concentration and energy distribution. (3) There is no soft photon source in our model; initially all the photons in the plasma are either energetic X-rays or gamma-rays. (4) The initial energy distribution of the pairs as well as photons is arbitrary. (5) We collect the analytical expressions for the kinetic kernels for all relevant processes. And (6) we present a different approach to finding the time-evolution of pair and photon spectra, which is a combination of the kinetic-theory and the non-linear Monte-Carlo schemes. We have developed many Monte-Carlo programs to model various process, to take into account the time evolution, and to incorporate various physical effects which are unique to non-thermal plasmas. The hydrodynamics of fireballs in GRB's was studied before. Applying results from kinetic theory will improve our understanding of these systems.

Ruderman, Malvin↗

NASA NTRS2000

Temperature Distribution Within a Defect-Free Silicon Carbide Diode Predicted by a Computational Model

Most solid-state electronic devices diodes, transistors, and integrated circuits are based on silicon. Although this material works well for many applications, its properties limit its ability to function under extreme high-temperature or high-power operating conditions. Silicon carbide (SiC), with its desirable physical properties, could someday replace silicon for these types of applications. A major roadblock to realizing this potential is the quality of SiC material that can currently be produced. Semiconductors require very uniform, high-quality material, and commercially available SiC tends to suffer from defects in the crystalline structure that have largely been eliminated in silicon. In some power circuits, these defects can focus energy into an extremely small area, leading to overheating that can damage the device. In an effort to better understand the way that these defects affect the electrical performance and reliability of an SiC device in a power circuit, the NASA Glenn Research Center at Lewis Field began an in-house three-dimensional computational modeling effort. The goal is to predict the temperature distributions within a SiC diode structure subjected to the various transient overvoltage breakdown stresses that occur in power management circuits. A commercial computational fluid dynamics computer program (FLUENT-Fluent, Inc., Lebanon, New Hampshire) was used to build a model of a defect-free SiC diode and generate a computational mesh. A typical breakdown power density was applied over 0.5 msec in a heated layer at the junction between the p-type SiC and n-type SiC, and the temperature distribution throughout the diode was then calculated. The peak temperature extracted from the computational model agreed well (within 6 percent) with previous first-order calculations of the maximum expected temperature at the end of the breakdown pulse. This level of agreement is excellent for a model of this type and indicates that three-dimensional computational modeling can provide useful predictions for this class of problem. The model is now being extended to include the effects of crystal defects. The model will provide unique insights into how high the temperature rises in the vicinity of the defects in a diode at various power densities and pulse durations. This information also will help researchers in understanding and designing SiC devices for safe and reliable operation in high-power circuits.

Kuczmarski, Maria A.↗

NASA NTRSDate not supplied

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗

NASA NTRS2016

Modeling the Ionosphere-Thermosphere Response to a Geomagnetic Storm Using Physics-based Magnetospheric Energy Input: OpenGGCM-CTIM Results

The magnetosphere is a major source of energy for the Earth's ionosphere and thermosphere (IT) system. Current IT models drive the upper atmosphere using empirically calculated magnetospheric energy input. Thus, they do not sufficiently capture the storm-time dynamics, particularly at high latitudes. To improve the prediction capability of IT models, a physics-based magnetospheric input is necessary. Here, we use the Open Global General Circulation Model (OpenGGCM) coupled with the Coupled Thermosphere Ionosphere Model (CTIM). OpenGGCM calculates a three-dimensional global magnetosphere and a two-dimensional high-latitude ionosphere by solving resistive magnetohydrodynamic (MHD) equations with solar wind input. CTIM calculates a global thermosphere and a high-latitude ionosphere in three dimensions using realistic magnetospheric inputs from the OpenGGCM. We investigate whether the coupled model improves the storm-time IT responses by simulating a geomagnetic storm that is preceded by a strong solar wind pressure front on August 24, 2005. We compare the OpenGGCM-CTIM results with low-earth-orbit satellite observations and with the model results of Coupled Thermosphere-Ionosphere-Plasmasphere electrodynamics (CTIPe). CTIPe is an up-to-date version of CTIM that incorporates more IT dynamics such as a low-latitude ionosphere and a plasmasphere, but uses empirical magnetospheric input. OpenGGCMCTIM reproduces localized neutral density peaks at approx. 400 km altitude in the high-latitude dayside regions in agreement with in situ observations during the pressure shock and the early phase of the storm. Although CTIPe is in some sense a much superior model than CTIM, it misses these localized enhancements. Unlike the CTIPe empirical input models, OpenGGCM-CTIM more faithfully produces localized increases of both auroral precipitation and ionospheric electric fields near the high-latitude dayside region after the pressure shock and after the storm onset, which in turn effectively heats the thermosphere and causes the neutral density increase at 400 km altitude.

Connor, Hyunju K.↗

NASA NTRS2018

Imaging Plasma Density Structures in the Soft X-Rays Generated by Solar Wind Charge Exchange with Neutrals

Both heliophysics and planetary physics seek to understand the complex nature of the solar wind's interaction with solar system obstacles like Earth's magnetosphere, the ionospheres of Venus and Mars, and comets. Studies with this objective are frequently conducted with the help of single or multipoint in situ electromagnetic field and particle observations, guided by the predictions of both local and global numerical simulations, and placed in context by observations from far and extreme ultraviolet (FUV, EUV), hard X-ray, and energetic neutral atom imagers (ENA). Each proposed interaction mechanism (e.g., steady or transient magnetic reconnection, local or global magnetic reconnection, ion pick-up, or the Kelvin- Helmholtz instability) generates diagnostic plasma density structures. The significance of each mechanism to the overall interaction (as measured in terms of atmospheric/ionospheric loss at comets, Venus, and Mars or global magnetospheric/ionospheric convection at Earth) remains to be determined but can be evaluated on the basis of how often the density signatures that it generates are observed as a function of solar wind conditions. This paper reviews efforts to image the diagnostic plasma density structures in the soft (low energy, 0.1-2.0 keV) X-rays produced when high charge state solar wind ions exchange electrons with the exospheric neutrals surrounding solar system obstacles.

Sibeck, David G.↗

NASA NTRS1989

Secondary electron emission characteristics of untreated and ion-textured titanium

Experimentally determined values of true secondary electron emission and relative values of reflected primary electron yield are presented for untreated (simply machined) and ion-textured, high-purity titanium over ranges of primary electron beam energies and beam impingement angles. The purpose of the investigation was to explore the feasibility of using titanium as electrode material in the multistage depressed collectors (MDC's) used in microwave amplifier traveling wave tubes (TWT's) for space communications and aircraft applications. Because of its relatively low density and thermal expansion characteristics and relatively high strength, thermal emissivity, and melting temperature, titanium presents itself as a possible candidate for the MDC electrode application. A detailed description of the method of ion texturing the titanium is included. Although the ion-treated surface considered in this study is not presented as being optimum from the standpoint of secondary electron emission suppression, it nevertheless serves to demonstrate that the surface can be modified by this procedure to significantly reduce these emission characteristics relative to those of the untreated surface. Further studies can reasonably be expected to produce surfaces with even lower secondary emission characteristics. The titanium surface were tested at primary electron beam energies of 200 to 2000 eV and at direct (0 deg) to near-grazing (85 deg) beam impingement angles. True secondary electron emission and relative reflected primary electron yield characteristics of the surfaces were compared with each other and with textured titanium surface exhibited secondary electron emission characteristics sharply lower than those exhibited by untreated titanium or copper. Clearly, then, in consideration of the secondary electron emission suppression of ion-textured titanium along with its other favorable physical properties, it must be included as a potential candidate for use as MDC electrode material in some applications.

Curren, Arthur N.↗

NASA NTRS1995

Extreme ultraviolet spectroscopy and photometry of VV Puppis during a high accretion state

We determine the physical properties of the accretion region of the AM Her-type binary VV Puppis using extreme ultraviolet (EUV) medium-resolution spectroscopy and photometry obtained with the Extreme Ultraviolet Explorer (EUVE) observatory. The EUV continuum from VV Pup was detected in the wavelength range from 75 to 135 A and was simultaneously recorded with the Deep Survey/Spectrometer (DS/S) imaging telescope, allowing for the extraction of an accurate light curve. VV Pup appeared to have entered a high-accretion state just prior to the pointed EUVE observations. We use the EUV light curve to infer the diameter of the accretion region (d = 220 km) assuming a hemispherical geometry and a radius of 9000 km for the white dwarf. We perform a model atmosphere analysis and, based on the light curve properties and assuming a distance of 145 pc, we derive an effective temperature of the accretion region in the range 270,000 is less than T(sub eff) is less than 360,000 K and a neutral hydrogen column density in the local interstellar medium of n(sub H) = 1.9 - 3.7 x 10(exp 19)/sq cm. The total EUV/soft X-ray energy radiated by the accretion region is approximately 3.5 x 10(exp 32) ergs/s. Our results provide a first verification of past suggestions that deep heating of the white dwarf surface produces the soft X-ray flux from the polars. We present a possible detection of O VI absortion features, and we suggest that extensive EUVE observations targeting high-accretion events may result in oxygen and heavier element abundance determination in the accretion region.

Vennes, Stephane↗

NASA NTRS2013

Galactic Cosmic Ray Event-Based Risk Model (GERM) Code

This software describes the transport and energy deposition of the passage of galactic cosmic rays in astronaut tissues during space travel, or heavy ion beams in patients in cancer therapy. Space radiation risk is a probability distribution, and time-dependent biological events must be accounted for physical description of space radiation transport in tissues and cells. A stochastic model can calculate the probability density directly without unverified assumptions about shape of probability density function. The prior art of transport codes calculates the average flux and dose of particles behind spacecraft and tissue shielding. Because of the signaling times for activation and relaxation in the cell and tissue, transport code must describe temporal and microspatial density of functions to correlate DNA and oxidative damage with non-targeted effects of signals, bystander, etc. These are absolutely ignored or impossible in the prior art. The GERM code provides scientists data interpretation of experiments; modeling of beam line, shielding of target samples, and sample holders; and estimation of basic physical and biological outputs of their experiments. For mono-energetic ion beams, basic physical and biological properties are calculated for a selected ion type, such as kinetic energy, mass, charge number, absorbed dose, or fluence. Evaluated quantities are linear energy transfer (LET), range (R), absorption and fragmentation cross-sections, and the probability of nuclear interactions after 1 or 5 cm of water equivalent material. In addition, a set of biophysical properties is evaluated, such as the Poisson distribution for a specified cellular area, cell survival curves, and DNA damage yields per cell. Also, the GERM code calculates the radiation transport of the beam line for either a fixed number of user-specified depths or at multiple positions along the Bragg curve of the particle in a selected material. The GERM code makes the numerical estimates of basic physical and biophysical quantities of high-energy protons and heavy ions that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the first option, properties of monoenergetic beams are treated. In the second option, the transport of beams in different materials is treated. Similar biophysical properties as in the first option are evaluated for the primary ion and its secondary particles. Additional properties related to the nuclear fragmentation of the beam are evaluated. The GERM code is a computationally efficient Monte-Carlo heavy-ion-beam model. It includes accurate models of LET, range, residual energy, and straggling, and the quantum multiple scattering fragmentation (QMSGRG) nuclear database.

Cucinotta, Francis A.↗

NASA NTRS1994

The 1992/93 eclipse of 31 Cygni

Extensive new ultraviolet and optical spectra of an atmospheric eclipse define the physical properties throughout the wind and chromosphere of 31 Cyg. These data require mass loss of approximately 3 x 10(exp -8) solar masses/yr in a wind that may merge smoothly into the chromosphere. Considerations of how energy is injected into the wind, however, suggest that the chromosphere and wind are separate structures. Most, if not all, of the velocity structure in metallic lines, which we have heretofore simulated with Doppler widths in the range 15-25 km/s, results from differential expansion of the atmosphere. Electron densities in the inner R(sub star) of the chromosphere are in the range 1.5 x 10(exp 9) - 2 x 10(exp 8)/cu cm, which implies clumping of the gas. The ionization of oxygen and nitrogen is consistent with clumping by factors of 3-30 in the outer chromosphere, roughly the amount required to give enough gas pressure to support the chromosphere. Chromospheric gas in 31 Cyg becomes hotter with increasing height, thus with decreasing optical depth, in a way that seems similar for all the zeta Aur binaries. Excitation temperature for Fe II in 1992/93 rises from about 5000 K at the deepest points sampled to about 12,500 K high in the wind. Strengths of violet Balmer lines give an excitation temperature for hydrogen of 6200-6500 K above a radial mass column density of approximately 0.01 g/sq cm. This amount of excitation implies that Lyman-alpha is thermalized beneath about 4 x 10(exp -3) g/sq cm. The outer atmosphere was symmetrical to within a factor of 2 in 1992/93, although it was clearly variable at this level, and it had similar mass column densities as in 1982. One manifestation of the variability was a flow toward the B star at phases 0.013-0.022 spanning velocities 45-100 km/s. Several lines of evidence point to a complicated and variable ionization in the wind: At large distance from the K star, measured mass column densities are less by up to a factor of 3 than required by a smoothly flowing wind. Also, much of the gas beyond r = 350 solar radii = 1.75 R(sub star) has a negative radial velocity. Radial velocities of the shell lines imply the outer atmosphere rotates in the direction of orbital motion, possibly through the deflection of wind flow lines in this direction.

Eaton, Joel A.↗

NASA NTRSDate not supplied

Initial analysis of “stone” size Ryugu samples: current status

As a part of the initial analysis of the Ryugu samples, we perform a variety of analyses of millimeter-sized "stones". Our goals are to elucidate the entire formation process of C-type asteroid Ryugu from the viewpoint of petrology and mineralogy ande obtain necessary information by sample analysis, and then simulate the formation of Ryugu based on the evidencef obtained rom sample analysis. Eighteen stones (8 from Room A and 10 from Room C) were received from the ISAS curation facility on June 1, 2021, and brought into a fully nitrogen-displaced glove box at Tohoku University. At the same time, we also received the powder samples from Room A and Room C. To completely block the atmosphere from leaking into the container, all samples were put in the sample transport containers prepared by ISAS, transferred from the main chamber to the glovebox at ISAS, and then all containers were completely sealed in plastic bags with moisture and oxygen absorbers. No moisture or oxygen was detected when the bags were opened in the glove box at Tohoku University, so it was confirmed that there was no exposure to the atmosphere during transport. To date, a number of analyses have been carried out successfully and almost on schedule. The analysis started with the measurement of reflectance spectra, which are sensitive to atmospheric oxidation, hydroxylation, and adsorbed water. The visible, near-infrared, and mid-infrared reflectance spectra were measured while the samples were kept airtight. The spectra of powder samples and stone samples (as aggregates and as single stone) were successfully obtained. A major feature of the stone team's analysis is the use of synchrotron radiation facilities around the world. Since this analysis is non-destructive, stone samples whose reflectance spectra were measured were sent to KEK, SPring-8, ESRF (France), SOLEIL (France), DESY (Germany), and APS (USA). Using these synchrotron radiation facilities, high spatial resolution and sensitivity XRD, STXM, XANES, CT [1], IR-CT, FT-IR [2, 3], XRF, and Mössbauer analyses were performed. Most of the analyses were carried out under air-tight conditions on the stone samples and the particulates separated from the stone samples. These analyses allowed us to determine the three-dimensional distribution of minerals and elements, redox state, density and porosity of the stone samples. Furthermore, as a characteristic analysis of the stone team, light elemental analysis using negative Muon was performed at the MLF facility of J-PARC with an exceptionally long allocation of machine time. This is the only non-destructive method to measure the concentration of light elements in the whole (not the surface) of stone samples. Because the characteristic X-rays produced by muon irradiation are much higher in energy than the fluorescent X-rays produced by X-ray irradiation, there is little effect of self-absorption by the sample, and therefore, the concentration of light elements such as carbon, oxygen, and Na in the entire "stone" sample can be determined. Some stone samples are currently being measured for heat and strength physical properties in order to understand the physical properties of asteroid Ryugu. The data obtained from these measurements are useful for interpreting the remote-sensing data data taken from the surface layer of the asteroid Ryugu [4-7]. It is also important for understanding the behavior of the Ryugu material during impact events. The surfaces of many stone samples were observed by electron microscopy and other techniques, especially on natural “flat” surfaces formed on 5 stomes. As a result, characteristric mineral aggregates formed by the reactions with water and characteiristic impact features were observed on some of the samples. Based on the observations of surfaces and the synchrotron measurements of the whole stones, important objects such as characteristic structures and specific crystal aggregates for understanding the formation history of the asteroid were identified, and these parts were separated from the stone samples using an Xe beam (pFIB) and analyzed by various methods including transmission electron microscopy and synchrotron radiation analysis. Many stone samples, from which important objects have been separated, are embedded in epoxy resin and cut to produce many polished sections. Electron microscopy and spectroscopic measurements of the polished surfaces are being carried out to reveal the detailed mineralogical properties and elemental distribution inside the stone samples. In the fall, machine time for synchrotron radiation will begin, and we plan to analyze single crystals and characteristic objects separated from the stone samples.

Tomoki Nakamura↗

NASA NTRS1999

Engineering of Novel Biocolloid Suspensions

Colloidal suspensions are materials with a variety of uses from cleaners and lubricants to food, cosmetics, and coatings. In addition, they can be used as a tool for testing the fundamental tenets of statistical physics. Colloidal suspensions can be synthesized from a wide variety of materials, and in the form of monodisperse particles, which can self-assemble into highly ordered colloidal crystal structures. As such they can also be used as templates for the construction of highly ordered materials. Materials design of colloids has, to date, relied on entropic self-assembly, where crystals form as result of lower free energy due to a transition to order. Here, our goal is to develop a completely new method for materials fabrication using colloidal precursors, in which the self-assembly of the ordered colloidal structures is driven by a highly controllable, attractive interaction. This will greatly increase the range of potential structures that can be fabricated with colloidal particles. In this work, we demonstrate that colloidal suspensions can be crosslinked through highly specific biological crosslinking reactions. In particular, the molecules we use are protein-carbohydrate interactions derived from the immune system. This different driving force for self-assembly will yield different and novel suspensions structures. Because the biological interactions are heterotypic (A binding to B), this chemical system can be used to make binary alloys in which the two colloid subpopulations vary in some property - size, density, volume fraction, magnetic susceptibility, etc. An additional feature of these molecules which is unique - even within the realm of biological recognition - is that the molecules bind reversibly on reasonable time-scales, which will enable the suspension to sample different configurations, and allow us to manipulate and measure the size of the suspension dynamically. Because of the wide variety of structures that can be made from these novel colloids, and because the suspension structure can be altered dynamically, we believe this biocolloid system will yield a novel set of materials with many technological applications, including sensors (both biological and non-biological), optical filters and separation media.

Hammer, D. A.↗

NASA NTRS2016

Data-Driven Radiative Hydrodynamic Modeling of the 2014 March 29 X1.0 Solar Flare

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in Mg II hk, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DSTIBIS) in H(alpha) 6563A and Ca II 8542A, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multithreaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated Mg II line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of approx. 800 km below the transition region can turn the simulated Mg II line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

Da Costa, Fatima Rubio↗