Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STATISTICAL MECHANICS”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Application of the target decomposition theorem to a polarimetric random media model

With advances in polarimetric radar measurements of land surfaces, the need for understanding the underlying scattering mechanisms and dominant target features has become the focus of many studies. In particular, the maximum use of the polarimetric information to identify and/or separate parameters related to the surface features such as vegetation thickness, structure, water content, and soil surface characteristics will enhance the possibility of using polarimetric radars for monitoring the earth's surface from space. In this paper, Cloude's decomposition theorem is applied to a polarimetric random media model to simulate the radar measurements of vegetated canopies. The vegetated canopies are modeled as a three layer discrete random medium with leaves and branches in the first layer, tree trunks in the second layer, and a half space of homogeneous ground with rough interface as the bottom layer. The distorted born approximation (DBA) has been used to compute full Mueller matrix of the canopy, using canonical dielectric objects such as thin discs and cylinders as leaves, branches, and trunks, respectively. The Mueller matrix and the derived covariance matrix contain information on the second order statistics of radar signals at various polarizations from the canopy. To decompose the covariance matrix to its constituent targets, the eigenvalues and eigenvectors of the covariance matrix are computed in terms of the physical parameters of the canopy. In addition, each eigenvector explicitly shows the scattering mechanisms such as odd and even reflections in the canopy. Cloude's decomposition theorem is applied using the Pauli spin matrices as a basis and an expression for the degree of disorder or the entropy for the vegetated surface is found. Then, the physical parameters estimated from in situ measurements are used in the random media to obtain realistic covariance matrices. As a result, the sensitivity of the eigenvalue spectrums and the coefficients resulting from the target decomposition theorem to the physical parameters of the canopy are examined and the possible use of Cloude's theorem to estimate vegetation and soil parameters is discussed.

Saatchi, Sasan S.↗
NASA NTRSDate not supplied

A Method for Validating Causal Diagrams of Human Health Risk in Space Flight

The complexity of cause-and-effect relationships between spaceflight hazards and resulting health conditions clouds understanding of the totality of human system risk in space. In response, NASA has introduced Directed Acyclic Graphs (causal diagrams) into the human systems risk management process. These diagrams allow for a common understanding of the mechanisms that lead from unique hazards of spaceflight to the health outcomes important to agencies and astronauts. However, the paucity of available biomedical data from spaceflight creates a need for methods of validating causal models that can accommodate data from spaceflight model analogs. Here we outline one approach utilizing open-access rodent bone datasets from the Ames Life Sciences Data Archive. The properties of directed acyclic graphs themselves can provide an epistemological and statistical framework for validation of a priori causal representations of human system risk in space flight. The assumed causal connections on the graph creates sets of logical implications: variables that – if the causal diagram is correct – should be correlated, as well as sets that should be conditionally independent. By testing these implied correlations and conditional independencies both statistically and heuristically, we can provide evidence for or against specific causal pathways on the causal diagram. In addition to validation of expert-generated causal diagrams, machine learning techniques can learn the most likely structure of a causal diagram from a given dataset. Comparison with and reconciliation between machine-learned causal diagrams and expert-generated diagrams is another technique for challenging assumptions and improving our understanding of causal mechanisms. Accurately representing complex causation is essential to systemic understanding of human health risks in space travel. Having a robust system of validating causal diagrams helps us arrive at more accurate representations of causal systems. This process will be integral to developing the countermeasures necessary for extended exploration of the moon and Mars.

Robert Reynolds↗
NASA NTRS2008

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗
DOE OSTI2022

Pauling’s rules for oxide-based minerals: A re-examination based on quantum mechanical constraints and modern applications of bond-valence theory to Earth materials

Since their introduction in 1929, Pauling’s five rules have been used by scientists from many disciplines to rationalize and predict stable arrangements of atoms and coordination polyhedra in crystalline solids; amorphous materials such as silicate glasses and melts; nanomaterials, poorly crystalline solids; aqueous cation and anion complexes; and sorption complexes at mineral-aqueous solution interfaces. The predictive power of these simple yet powerful rules was challenged recently by George et al. (2020), who performed a statistical analysis of the performance of Pauling’s five rules for about 5000 oxide crystal structures. They concluded that only 13% of the oxides satisfy the last four rules simultaneously and that the second rule has the most exceptions. They also found that Pauling’s first rule is satisfied for only 66% of the coordination environments tested and concluded that no simple rule linking ionic radius to coordination environment will be predictive due to the variable quality of univalent radii. We address these concerns and discuss quantum mechanical calculations that complement Pauling’s rules, particularly his first (radius sum and radius ratio rule) and second (electrostatic valence rule) rules. We also present a more realistic view of the bonded radii of atoms, derived by determining the local minimum in the electron density distribution measured along trajectories between bonded atoms known as bond paths, i.e., the bond critical point (rc). Electron density at the bond critical point is a quantum mechanical observable that correlates well with Pauling bond strength. Moreover, a metal atom in a polyhedron has as many bonded radii as it has bonded interactions, resulting in metal and O atoms that may not be spherical. O atoms, for example, are not spherical in many oxide-based crystal structures. Instead, the electron density of a bonded oxygen is often highly distorted or polarized, with its bonded radius decreasing systematically from ~1.38 Å when bonded to highly electropositive atoms like sodium to 0.64 Å when bonded to highly electronegative atoms like nitrogen. Bonded radii determined for metal atoms match the Shannon (1976) radii for more electropositive atoms, but the match decreases systematically as the electronegativities of the M atoms increase. As a result, significant departures from the radius ratio rule in the analysis by George et al. (2020) is not surprising. We offer a modified, more fundamental version of Pauling’s first rule and demonstrate that the second rule has a one-to-one connection between the electron density accumulated between the bonded atoms at the bond critical point and the Pauling bond strength of the bonded interaction. Pauling’s second rule implicitly assumes that bond strength is invariant with bond length for a given pair of bonded atoms. Many studies have since shown that this is not the case, and Brown and Shannon (1973) developed an equation and a set of parameters to describe the relation between bond length and bond strength, now redefined as bond valence to avoid confusion with Pauling bond-strength. Brown (1980) used the valence-sum rule, together with the path rule and the valence-matching principle, as the three axioms of bond-valence theory (BVT), a powerful method for understanding many otherwise elusive aspects of crystals and also their participation in dynamic processes. We show how a priori bond-valence calculations can predict unstrained bond-lengths and how bond-valence mapping can locate low-Z atoms in a crystal structure (e.g., Li) or examine possible diffusion pathways for atoms through crystal structures. In addition, we briefly discuss Pauling’s third, fourth, and fifth rules, the first two of which concern the sharing of polyhedron elements (edges and faces) and the common instability associated with structures in which a polyhedron shares an edge or face with another polyhedron and contains high-valence cations. The olivine [α-(MgxFe1–x)2SiO4] crystal structure is used to illustrate the distortions from hexagonal close-packing of O atoms caused by metal-metal repulsion across shared polyhedron edges. We conclude by discussing several applications of BVT to Earth materials, including the use of BVT to: (1) locate H+ ions in crystal structures, including the location of protons in the crystal structures of nominally anhydrous minerals in Earth’s mantle; (2) determine how strongly bonded (usually anionic) structural units interact with weakly bonded (usually cationic) interstitial complexes in complex uranyl-oxide and uranyl-oxysalt minerals using the valence-matching principle; (3) calculate Lewis acid strengths of cations and Lewis base strengths of anions; (4) determine how (H2O) groups can function as bond-valence transformers by dividing one bond into two bonds of half the bond valence; (5) help characterize products of sorption reactions of aqueous cations (e.g., Co2+ and Pb2+) and oxyanions [e.g., selenate (Se6+O4)2- and selenite (Se4+O3)2-] at mineral-aqueous solution interfaces and the important role of protons in these reactions; and (6) help characterize the local coordination environments of highly charged cations (e.g., Zr4+, Ti4+, U4+, U5+, and U6+) in silicate glasses and melts.

Geochemistry & Geophysics↗
DOE OSTI2021

Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study

Development of an effective antiviral drug for Coronavirus Disease 2019 (COVID-19) is a global health priority. Although several candidate drugs have been identified through in vitro and in vivo models, consistent and compelling evidence from clinical studies is limited. The lack of evidence from clinical trials may stem in part from the imperfect design of the trials. We investigated how clinical trials for antivirals need to be designed, especially focusing on the sample size in randomized controlled trials. A modeling study was conducted to help understand the reasons behind inconsistent clinical trial findings and to design better clinical trials. We first analyzed longitudinal viral load data for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) without antiviral treatment by use of a within-host virus dynamics model. The fitted viral load was categorized into 3 different groups by a clustering approach. Comparison of the estimated parameters showed that the 3 distinct groups were characterized by different virus decay rates (p-value < 0.001). The mean decay rates were 1.17 d -1 (95% CI: 1.06 to 1.27 d -1 ), 0.777 d -1 (0.716 to 0.838 d -1 ), and 0.450 d -1 (0.378 to 0.522 d -1 ) for the 3 groups, respectively. Such heterogeneity in virus dynamics could be a confounding variable if it is associated with treatment allocation in compassionate use programs (i.e., observational studies). Subsequently, we mimicked randomized controlled trials of antivirals by simulation. An antiviral effect causing a 95% to 99% reduction in viral replication was added to the model. To be realistic, we assumed that randomization and treatment are initiated with some time lag after symptom onset. Using the duration of virus shedding as an outcome, the sample size to detect a statistically significant mean difference between the treatment and placebo groups (1:1 allocation) was 13,603 and 11,670 (when the antiviral effect was 95% and 99%, respectively) per group if all patients are enrolled regardless of timing of randomization. The sample size was reduced to 584 and 458 (when the antiviral effect was 95% and 99%, respectively) if only patients who are treated within 1 day of symptom onset are enrolled. We confirmed the sample size was similarly reduced when using cumulative viral load in log scale as an outcome. We used a conventional virus dynamics model, which may not fully reflect the detailed mechanisms of viral dynamics of SARS-CoV-2. The model needs to be calibrated in terms of both parameter settings and model structure, which would yield more reliable sample size calculation. In this study, we found that estimated association in observational studies can be biased due to large heterogeneity in viral dynamics among infected individuals, and statistically significant effect in randomized controlled trials may be difficult to be detected due to small sample size. The sample size can be dramatically reduced by recruiting patients immediately after developing symptoms. We believe this is the first study investigated the study design of clinical trials for antiviral treatment using the viral dynamics model.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Reconstructing neutrinoless double beta decay event kinematics in a xenon gas detector with vertex tagging

If neutrinoless double beta decay is discovered, the next natural step would be understanding the lepton number violating physics responsible for it. Several alternatives exist beyond the exchange of light neutrinos. Some of these mechanisms can be distinguished by measuring phase-space observables, namely the opening angle cos θ among the two decay electrons, and the electron energy spectra, T 1 and T 2 . In this work, we study the statistical accuracy and precision in measuring these kinematic observables in a future xenon gas detector with the added capability to precisely locate the decay vertex. For realistic detector conditions (a gas pressure of 10 bar and spatial resolution of 4 mm), we find that the average $\overline{cos θ}$ and $\overline{T_1}$ values can be reconstructed with a precision of 0.19 and 110 keV, respectively, assuming that only 10 neutrinoless double beta decay events are detected.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Lithium and boron isotopic compositions of olivine in chondrules from carbonaceous and ordinary chondrite meteorites: Implications for the origin of solar 11 B/ 10 B ratio

We report the origin of solar 11 B/ 10 B ratio ~4 remains an open question. It has been thought that a significant portion of boron in the Solar System was derived from continuous spallation nucleosynthesis during interactions between Galactic Cosmic Rays and C—N—O nuclei in the interstellar medium. However, because GCR-produced boron is characterized by 11 B/ 10 B ~2.5, an endmember with 11 B/ 10 B > 4 is required to account for the solar 11 B/ 10 B ratio. Two leading hypotheses for the sources of 11 B-rich components include low energy spallation in the Sun’s parental molecular cloud and the neutrino process during the supernova explosions. In this study, lithium and boron elemental and isotopic compositions of seven porphyritic olivine chondrules and one isolated olivine crystal from four meteorites (Allende, Yamato 81020, Asuka 12236, and QUE 97008) were determined in-situ to help constrain which of the two nucleosynthetic mechanisms has most likely supplied the forming Solar System with extra 11 B. Apparent isotopic variations in Li/Si, B/Si, δ 7 Li and δ 11 B were found in chondrule olivine crystals, but only three chondrules exhibit statistically resolved δ 11 B heterogeneities. Using these three chondrules as a basis for discussion, we evaluated the processes that could potentially cause elemental and isotopic variations of Li and B. We argue that the data can be best understood in the context of condensation of lithium-boron-rich material onto chondrule precursors in (an) initially heterogeneous gaseous reservoir(s), which could be realized if the Solar System derived 11 B-rich components from a supernova or supernovae.

58 GEOSCIENCES↗
DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Performance of the spectropolarimeter for the Space Telescope faint object spectrograph

The design and preliminary test results for the spectropolarimeter for the Faint Object Spectrograph (FOS) for the Space Telescope are presented. The mechanical design and optical specifications of the spectropolarimeter are described noting that a Wollaston prism with an internal wedge angle of 20 deg is fixed behind each of two rotatable waveplate retarders of different retardations. Either waveplate/prism combination can be positioned at either of the two FOS entrance ports. Magnesium fluoride is chosen as the birefringent crystal for the polarizing elements to allow linear and circular polarization measurements down to Lyman-alpha at 1216 A. Mechanical stability and repeatability were determined by operational testing to give polarization-position angles of + or - 0.5 deg, corresponding to degree-of-polarization measurements of + or - 0.1 percent. Faint-object accuracy, dependent on photoelectron statistics and hence on observation time, is calculated to be one percent in each 100-A-wide spectral band for a 20-min observation of an AO star with V = 15th magnitude.

Allen, R. G.↗
NASA NTRS1992

Microstructural design of fiber composites

The optimum performance design of composite microstructures is discussed. The forces driving progress in fiber composites are examined, and recent developments in the mechanics of laminated composites are surveyed, emphasizing thick laminates, hygrothermal effects, and thermal transient effects. The strength of continuous-fiber composites is discussed, presenting analyses of local load redistribution due to fiber breakages and treatments of statistical tensile strength theories. Modes of failure of laminated composites are examined. Elastic, physical, and viscoelastic properties as well as the strength and fracture behavior of short-fiber composites are studied, and it is shown how the performance of composites can be controlled by selecting material systems and their geometric distributions. 2D textile structural composites based on woven, knitted, and braided preforms are considered, and techniques for analyzing and modeling the thermomechanical behavior of 2D textile composites are presented. Recent developments in the processing of 3D textile preforms are introduced and the processing-microstructure relationship is demonstrated. Finite elastic deformation of flexible composites is addressed.

Chou, Tsu-Wei↗
NASA NTRS1995

Heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(s)

The heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) (eq 1) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(S) (eq 2) were investigated over the temperature range 223-296 K in a flow-tube reactor coupled to a quadrupole mass spectrometer. Either a chemical ionization mass spectrometer (CIMS) or an electron-impact ionization mass spectrometer (EIMS) was used to provide suitable detection sensitivity and selectivity. In order to mimic atmospheric conditions, partial pressures of HNO3 and N2O5 in the range 6 x 10(exp -8) - 2 x 10(exp -6) Torr were used. Granule sizes and surface roughness of the solid NaCl substrates were determined by using a scanning electron microscope. For dry NaCl substrates, decay rates of HNO3 were used to obtain gamma(1) = 0.013 +/- 0.004 (1sigma) at 296 K and > 0.008 at 223 K, respectively. The error quoted is the statistical error. After all corrections were made, the overall error, including systematic error, was estimated to be about a factor of 2. HCl was found to be the sole gas-phase product of reaction 1. The mechanism changed from heterogeneous reaction to predominantly physical adsorption when the reactor was cooled from 296 to 223 K. For reaction 2 using dry salts, gamma(2) was found to be less than 1.0 x 10(exp -4) at both 223 and 296 K. The gas-phase reaction product was identified as ClNO2 in previous studies using an infrared spectrometer. An enhancement in reaction probability was observed if water was not completely removed from salt surfaces, probably due to the reaction of N2O5(g) + H2O(s) yields 2HNO3(g). Our results are compared with previous literature values obtained using different experimental techniques and conditions. The implications of the present results for the enhancement of the hydrogen chloride column density in the lower stratosphere after the El Chichon volcanic eruption and for the chemistry of HCl and HNO3 in the marine troposphere are discussed.

Leu, Ming-Taun↗
NASA NTRS2000

Statistical Study of Auroral Kilometric Radiation Fine Structure Striations Observed by Polar

We have conducted a statistical survey of a semirandom sample of the auroral kilometric radiation (AKR) data observed by the plasma wave instrument wideband receiver on board the Polar spacecraft. We have determined that AKR fine structure patterns with very narrowband, negative drifting striations occur in approximately 6% of the high-resolution wideband spectrograms when AKR is present. Positive sloping striations are also observed, but at a much lower rate. More than 8200 AKR stripes have been scaled. The stripes are predominantly found in the 40 to 215-kHz frequency range and have a frequency extent of about 4 kHz and a duration of usually less than 2 s. The majority of the stripes have drift rates between -8 and -2 kHz/s, with a peak in the distribution between -6 and -4 kHz/s. There is also a much smaller group of striations with positive drift rates of up to about 5 or 6 kHz/s. We have further investigated the change of drift rate with frequency. Almost all striations are observed in the lowest two frequency bands of the wideband receiver (f < 215 kHz). There is an increase in the statistical drift rate with increasing frequency. The statistical slope of the striations increases with frequency from about -4.4 kHz/s at 75 kHz to about -5.7 kHz/s at 170 kHz. This frequency dependence of the drift rate is consistent, under certain conditions, with a production mechanism stimulated by an upward propagating electromagnetic ion cyclotron wave, as had been suggested earlier. However, such a changing drift rate is also compatible with a stimulated source region that propagates upward along the magnetic field line at the velocity of an ion beam accelerated by a local, upward directed electric field, as is typically observed in the auroral region. An explanation for this association is not apparent at this time.

Menietti, J. D.↗
NASA NTRS2018

Exoplanet Biosignatures: Future Directions

We introduce a Bayesian method for guiding future directions for detection of life on exoplanets. We describe empirical and theoretical work necessary to place constraints on the relevant likelihoods, including those emerging from better understanding stellar environment, planetary climate and geophysics, geochemical cycling, the universalities of physics and chemistry, the contingencies of evolutionary history, the properties of life as an emergent complex system, and the mechanisms driving the emergence of life. We provide examples for how the Bayesian formalism could guide future search strategies, including determining observations to prioritize or deciding between targeted searches or larger lower resolution surveys to generate ensemble statistics and address how a Bayesian methodology could constrain the prior probability of life with or without a positive detection.

Bayesian analysis↗
DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

60 - APPLIED LIFE SCIENCES↗
NASA NTRS2018

Effect of Ex Vivo Ionizing Radiation on Static and Fatigue Properties of Mouse Vertebral Bodies

For a variety of medical and scientific reasons, human bones can be exposed to a wide range of ionizing radiation levels. In vivo radiation therapy (0.05 kGy) is used in cancer treatment, and ex vivo irradiation (25-35 kGy) is used to sterilize bone allografts. Ionizing radiation in these applications has been shown to increase risk of fracture, decrease bone quality and degrade collagen integrity. Past studies have investigated the deleterious effects of radiation on cortical or trabecular bone specimens individually, but to date no studies have examined whole bones containing both cortical and trabecular tissue. Furthermore, a clear relationship between the dose and the mechanical and biochemical response of bone's extracellular matrix has yet to be established for doses ranging from cancer therapy to allograft sterilization (0.05-35 kGy). To gain insight into these issues, we conducted an ex vivo radiation study to investigate non-cellular (i.e. matrix) effects of ionizing radiation dose on vertebral whole bone mechanical properties, over a range of radiation doses (0.05-35 kGy), with a focus on any radiation-induced changes in collagen. With underlying mechanisms of action in mind, we hypothesized that any induced reductions in mechanical properties would be associated with changes in collagen integrity. METHODS: 20-week old female mice were euthanized and the lumbar spine was dissected using IACUC approved protocols. The lumbar vertebrae (L1- S1) were extracted from the spine via cuts through adjacent intervertebral discs, and the endplates, posterior processes, surrounding musculature, and soft tissues were removed (approx. 1.5mm diameter, approx. 2mm height). Specimens were randomly assigned to one of five groups for ex vivo radiation exposure: x-ray irradiation at 0.05, 1, 17, or 35 kGy, or a 0 kGy control. Following irradiation, the vertebrae were imaged using microcomputed tomography (micro-CT) and then subjected to either monotonic compressive loading to failure or uniform cyclic compressive loading. During cyclic testing, samples were loaded in force control to a force level that corresponded to a strain of 0.46%, as determined in advance by a linearly elastic micro-CT-based finite element analysis for each specimen. Tests were stopped at imminent fracture, defined as a rapid increase in strain. The main outcome for the monotonic test was the strength (maximum force); for cyclic testing it was the fatigue life (log of the number of cycles of loading at imminent failure). A fluorometric assay was used on the S1 vertebrae to measure the number of non-enzymatic collagen crosslinks[4]. A one-way ANOVA was performed on mechanical properties and collagen crosslinks; means were compared with controls using Dunnett's method, with a Tukey-Kramer post-hoc analysis when significance was found (p < or = 0.05). RESULTS: Compared to the unirradiated control group, the concentration of non-enzymatic collagen crosslinks was significantly increased for all irradiated groups (p < 0.0001), and being higher by at least 50% (Figure 1a). By contrast, the radiation effects on the collagen were only evident at the higher doses. For irradiation exposures of 17 kGy or more, strength decreased substantially as the radiation level was increased, but no effect was evident below 17 kGy (Figure 1b). There was no significant change in the stiffness or maximum displacement for any radiation dose (p>0.05). The finite element analysis prescribed force level for cyclic loading exceeded the measured (monotonic) strength of the 17 and 35 kGy irradiated groups (mean +/- SD, 20.6 +/- 5.6 N; 13.2 +/- 3.7 N, respectively) and therefore these groups were eliminated from the fatigue study. The fatigue life for the 0.05 and 1 kGy groups were similar to each other and were not statistically significantly different from the control group (Figure 1c).

radiation↗
DOE OSTI2021

A statistical analysis of developing knock intensity in a mixture with temperature inhomogeneities

Knock formation and its intensity for a stoichiometric ethanol/air mixture under a representative end-gas auto-ignition condition in IC engines with temperature inhomogeneities are investigated using multi-dimensional direct numerical simulations (DNS) with a 40-species skeletal mechanism of ethanol. Two- and three-dimensional simulations are performed by systematically varying temperature fluctuations and its most energetic length scale, l T . The volumetric fraction of the mixture regions that have the propensity to detonation development, F D , is proposed as a metric to predict the amplitude of knock intensity. It is found that with increasing l T , F D shows a good agreement with the heat release fraction of the mixture regions with pressure greater than equilibrium pressure, F H . The detonation peninsula is well captured by F D and F H when plotting them as a function of the volume-averaged ξ, ξ, (ξ= a/S sp is the ratio of the acoustic speed, a to the ignition front speed, S sp ) Decreasing l T is found to significantly reduce the super-knock intensity. The results suggest that decreasing l T , as in engines with tumble designs resulting in a smaller turbulence scale, will be effective in mitigating the the super-knock development.

42 ENGINEERING↗
DOE OSTI2022

Calabi-Yau CFTs and random matrices

Using numerical methods for finding Ricci-flat metrics, we explore the spectrum of local operators in two-dimensional conformal field theories defined by sigma models on Calabi-Yau targets at large volume. Focusing on the examples of K3 and the quintic, we show that the spectrum, averaged over a region in complex structure moduli space, possesses the same statistical properties as the Gaussian orthogonal ensemble of random matrix theory.

superstring vacua↗