Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical sensitivity analysis”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2005

Inference of Ice Cloud Properties from High-spectral Resolution Infrared Observations

The theoretical basis is explored for inferring the microphysical properties of ice crystal from high-spectral resolution infrared observations. A radiative transfer model is employed to simulate spectral radiances to address relevant issues. The extinction and absorption efficiencies of individual ice crystals, assumed as hexagonal columns for large particles and droxtals for small particles, are computed from a combination of the finite- difference time-domain (FDTD) technique and a composite method. The corresponding phase functions are computed from a combination of FDTD and an improved geometric optics method (IGOM). Bulk scattering properties are derived by averaging the single- scattering properties of individual particles for 30 particle size distributions developed from in situ measurements and for additional four analytical Gamma size distributions for small particles. The non-sphericity of ice crystals is shown to have a significant impact on the radiative signatures in the infrared (IR) spectrum; the spherical particle approximation for inferring ice cloud properties may result in an overest&ation of the optical thickness and an inaccurate retrieval of effective particle size. Furthermore, we show that the error associated with the use of the Henyey-Greenstein phase function can be as larger as 1 K in terms of brightness temperature for larger particle effective size at some strong scattering wavenumbers. For small particles, the difference between the two phase functions is much less, with brightness temperatures generally differing by less than 0.4 K. The simulations undertaken in this study show that the slope of the IR brightness temperature spectrum between 790-960/cm is sensitive to the effective particle size. Furthermore, a strong sensitivity of IR brightness temperature to cloud optical thickness is noted within the l050-1250/cm region. Based on this spectral feature, a technique is presented for the simultaneous retrieval of the visible optical thickness and effective particle size from high spectral resolution infrared data under ice cloudy con&tion. The error analysis shows that the uncertainty of the retrieved optical thickness and effective particle size has a small range of variation. The error for retrieving particle size in conjunction with an uncertainty of 5 K in cloud'temperature, or a surface temperature uncertainty of 2.5 K, is less than 15%. The corresponding e m r in the uncertainty of optical thickness is within 5-2096, depending on the value of cloud optical thickness. The applicability of the technique is demonstrated using the aircraft-based High- resolution Interferometer Sounder (HIS) data from the Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) in 1996 and the First ISCCP Regional Experiment - Arctic Clouds Experiment (FIRE-ACE) in 1998.

Huang, Hung-Lung↗
DOE OSTI2022

Prospects and Limitations of Predicting Fuel Ignition Properties from Low-Temperature Speciation Data

Using chemical kinetic modeling and statistical analysis, we investigate the possibility of correlating key chemical “markers”–typically small molecules–formed during very lean (φ ~ 0.001) oxidation experiments with near-stoichiometric (φ ~ 1) fuel ignition properties. One goal of this work is to evaluate the feasibility of designing a fuel-screening platform, based on small laboratory reactors that operate at low temperatures and use minimal fuel volume. Buras et al. [Combust. Flame2020,216, 472–484] have shown that convolutional neural net (CNN) fitting can be used to correlate first-stage ignition delay times (IDTs) with OH/HO2 measurements during very lean oxidation in low-T flow reactors with better than factor-of-2 accuracy. In this work, we test the limits of applying this correlation-based approach to predict the low-temperature heat release (LTHR) and total IDT, including the sensitivity of total IDT to the equivalence ratio, φ. We demonstrate that first-stage IDT can be reliably correlated with very lean oxidation measurements using compressed sensing (CS), which is simpler to implement than CNN fitting. LTHR can also be predicted via CS analysis, although the correlation quality is somewhat lower than for first-stage IDT. In contrast, the accuracy of total IDT prediction at φ = 1 is significantly lower (within a factor of 4 or worse). Furthermore, these results can be rationalized by the fact that the first-stage IDT and LTHR are primarily determined by low-temperature chemistry, whereas total IDT depends on low-, intermediate-, and high-temperature chemistry. Oxidation reactions are most important at low temperatures, and therefore, measurements of universal molecular markers of oxidation do not capture the full chemical complexity required to accurately predict the total IDT even at a single equivalence ratio. As a result, we find that φ-sensitivity of ignition delay cannot be predicted at all using solely correlation with lean low-T chemical speciation measurements.

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

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Shock Waves↗
DOE OSTI2023

Chemical bonding in americium oxides probed by X-ray spectroscopy

The electronic structure and the chemical state in Am binary oxides and Am-doped UO 2 were studied by means of X-ray absorption spectroscopy at shallow Am core (4d and 5d) edges. In particular, the Am 5f states were probed and the nature of their bonding to the oxygen states was analyzed. The interpretation of the experimental data was supported by the Anderson impurity model (AIM) calculations which took into account the full multiplet structure due to the interaction between 5f electrons as well as the interaction with the core hole. The sensitivity of the branching ratio of the Am 4d 3/2 and 4d 5/2 X-ray absorption lines to the chemical state of Am was shown using Am binary oxides as reference systems. The observed ratio for Am-doped UO 2 suggests that at least at low Am concentrations, americium is in the Am(III) state in the UO 2 lattice. To confirm the validity of the applied AIM approach, the analysis of the Am 4f X-ray photoelectron spectra of AmO 2 and Am 2 O 3 was also performed which revealed a good agreement between experiment and calculations. As a whole, AmO 2 can be classified as the charge-transfer compound with the 5f occupancy (n ƒ ) equal to 5.73 electrons, while Am 2 O 3 is rather a Mott–Hubbard system with n ƒ = 6.05.

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

Colloidal CdSe nanocrystals are inherently defective

Abstract Colloidal CdSe nanocrystals (NCs) have shown promise in applications ranging from LED displays to medical imaging. Their unique photophysics depend sensitively on the presence or absence of surface defects. Using simulations, we show that CdSe NCs are inherently defective; even for stoichiometric NCs with perfect ligand passivation and no vacancies or defects, we still observe that the low energy spectrum is dominated by dark, surface-associated excitations, which are more numerous in larger NCs. Surface structure analysis shows that the majority of these states involve holes that are localized on two-coordinate Se atoms. As chalcogenide atoms are not passivated by any Lewis base ligand, varying the ligand should not dramatically change the number of dark states, which we confirm by simulating three passivation schemes. Our results have significant implications for understanding CdSe NC photophysics, and suggest that photochemistry and short-range photoinduced charge transfer should be much more facile than previously anticipated.

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

Lessons from Astrobiological Planetary Analogue Exploration in Iceland: Biomarker Assay Performance and Downselection

Understanding the sensitivity of biomarker assays to the local physicochemical environment, and the underlying spatial distribution of the target biomarkers in 'homogeneous' environments, can increase mission science return. We have conducted four expeditions to Icelandic Mars analogue sites in which an increasingly refined battery of physicochemical measurements and biomarker assays were performed, staggered with scouting of further sites. Completed expeditions took place in 2012 (location scouting and field assay use testing), 2013 (sampling of two major sites with three assays and observational physicochemical measurements), 2015 (repeat sampling of prior sites and one new site, scouting of new sites, three assays and three instruments), and 2016 (preliminary sampling of new sites with analysis of returned samples). Target sites were geologically recent basaltic lava flows, and sample loci were arranged in hierarchically nested grids at 10 cm, 1 m, 10 m, 100 m, and >1 km order scales, subject to field constraints. Assays were intended to represent a diversity of potential biomarker types (cell counting via nucleic acid staining and fluorescence microscopy, ATP quantification via luciferase luminescence, and relative DNA quantification with simple domain-level primers) rather a specific mission science target, and were selected to reduce laboratory overhead, require limited consumables, and allow rapid turnaround. All analytical work was performed in situ or in a field laboratory within a day's travel of the field sites unless otherwise noted. We have demonstrated the feasibility of performing ATP quantification and qPCR analysis in a field-based laboratory with single-day turnaround. The ATP assay was generally robust and reliable and required minimal field equipment and training to produce a large amount of useful data. DNA was successfully extracted from all samples, but the serial-batch nature of qPCR significantly limited the number of primers (hence classifications) and replicates that could be run in a single day. Fluorescence microscopy did not prove feasible under the same constraints, primarily due to the large number of person-hours required to view, analyze, and record results from the images; however, this could be mitigated with higher-quality imaging instruments and appropriate image analysis software.

Planetary↗
DOE OSTI2022

Evolution of Bonding and Magnetism via Changes in Valence Electron Count in CuFe 2– x Co x Ge 2

A series of solid solutions, CuFe 2–x Co x Ge 2 (x = 0, 0.2, 0.4, 0.8, and 1.0), have been synthesized by arc-melting and characterized by powder X-ray and neutron diffraction, magnetic measurements, Mössbauer spectroscopy, and electronic band structure calculations. All compounds crystallize in the CuFe 2 Ge 2 structure type, which can be considered as a three-dimensional framework built of fused MGe6 octahedra and MGe 5 trigonal bipyramids (M = Fe and Co), with channels filled by rows of Cu atoms. As the Co content (x) increases, the unit cell volume decreases in an anisotropic fashion: the b and c lattice parameters decrease while the a parameter increases. The changes in all the parameters are nearly linear, thus following Vegard’s law. CuFe 2 Ge 2 exhibits two successive antiferromagnetic (AFM) orderings, corresponding to the formation of a commensurate AFM structure, followed by an incommensurate AFM structure observed at lower temperatures. Additionally, as the Co content increases, the AFM ordering temperature (TN) gradually decreases, and only one AFM transition is observed for x ≥ 0.2. The magnetic behavior of unsubstituted CuFe 2 Ge 2 was found to be sensitive to the preparation method. The temperature-dependent zero-field 57 Fe Mössbauer spectra reveal two hyperfine split components that evolve in agreement with the two consecutive AFM orderings observed in magnetic measurements. In contrast, the field-dependent spectra obtained for fields ≥ 2 T reveal a parallel arrangement of the moments associated with the two crystallographically unique metal sites. Electronic band structure calculations and chemical bonding analysis reveal a mix of strong M–M antibonding and non-bonding states at the Fermi level, in support of the overall AFM ordering observed in zero field. The substitution of Co for Fe reduces the population of the M–M antibonding states and the overall density of states at the Fermi level, thus suppressing the T N value.

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

The Role of Spin–Orbit Coupling in the Linear Absorption Spectrum and Intersystem Crossing Rate Coefficients of Ruthenium Polypyridyl Dyes

The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g., triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using ab initio electronic structure methods are invaluable to identify structure-property relations for dye sensitizers. The primary challenge for simulations to elucidate the electronic and nuclear origins of these properties is a spin-orbit interaction which drives transitions between electronic states. Here, in this work, we present a computational analysis of the spin-orbit corrected linear absorption cross sections and intersystem crossing rate coefficients for a derivative set of phosphonated tris(2,2'-bipyridine)ruthenium(2+) dye molecules. After sampling the ground state vibrational distributions, the predicted linear absorption cross sections indicate that the mixture between singlet and triplet states plays a crucial role in defining the line shape of the metal-to-ligand charge transfer bands in these derivatives. Additionally, an analysis of the intersystem crossing rate coefficients suggests that transitions from the singlet into the triplet manifolds are ultrafast with rate coefficients on the order of 10 13 s -1 for each dye molecule.

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

An experimental and modeling study on auto-ignition kinetics of ammonia/methanol mixtures at intermediate temperature and high pressure

A rapid compression machine (RCM) has been applied to measure the ignition delay times of NH 3 /CH 3 OH mixtures covering pressures of 20 and 40 bar, equivalence ratios of 0.5, 1.0 and 2.0, and temperatures between 845 and 1100 K. Here the measurements show that the NH 3 /CH 3 OH mixtures become more reactive with increasing methanol addition. Addition of merely 1% (molar basis) of CH 3 OH to NH 3 lowers the ignition temperature around 100 K at 40 bar in comparison to pure NH 3 . The ignition delay is a complex function of fuel mixture and stoichiometry. For the 1% CH 3 OH mixture, the leaner mixtures are more reactive, while the reverse trend is found for mixtures with 5%, 20% and pure CH 3 OH. Analysis of the pressure profiles shows three distinct ignition modes for NH 3/ CH 3 OH mixtures, facilitated by the pre-ignition heat release from NH 3 consumption. A detailed mechanism for ignition of NH 3 /CH 3 OH fuel blends has been developed, capable of reproducing the ignition behavior of mixtures with reasonable accuracy. A subset for amine / methanol interactions was established, with rate constants for the key reaction between NH 2 and CH 3 OH calculated from ab initio theory. A sensitivity analysis indicates that the critical reactions during the auto-ignition process vary with the CH 3 OH mole fraction in the fuel. The ammonia chemistry, namely NH 2 + NO, NH2 + NO 2 and NH 2 + HO 2 , is dominant for the mixture with 1% CH 3 OH, while the reactions related to CH 3 OH and H 2 O 2 are more important for the 20% CH 3 OH mixture. The interaction between ammonia and methanol shows a more prominent effect on auto-ignition for mixtures with 5% CH 3 OH in fuel as compared to those with 1% and 20% CH 3 OH. According to the modeling results, methanol addition is found to enrich the O/H radical pool, consuming ammonia and promoting auto-ignition through different reaction pathways.

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

One Electron Multiple Proton Transfer in Model Organic Donor–Acceptor Systems: Implications for High-Frequency EPR

EPR spectroscopy is an important spectroscopic method for identification and characterization of radical species involved in many biological reactions. The tyrosyl radical is one of the most studied amino acid radical intermediates in biology. Often in conjunction with histidine residues, it is involved in many fundamental biological electron and proton transfer processes, such as in the water oxidation in photosystem II. As biological processes are typically extremely complicated and hard to control, molecular bio-mimetic model complexes are often used to clarify the mechanisms of the biological reactions. In this work, we present theoretical calculations to investigate the sensitivity of magnetic resonance parameters to proton-coupled electron transfer events, as well as conformational substates of the molecular constructs which mimic the tyrosine-histidine (Tyr-His) pairs found in a large variety of proteins. Upon oxidation of the phenol, the Tyr analog, these complexes can perform not only one-electron one-proton transfer (EPT), but also one-electron two-proton transfers (E2PT). It is shown that in aprotic environment the g X -components of the electronicg-tensor are extremely sensitive to the first proton transfer from the phenoxyl oxygen to the imidazole nitrogen (EPT product), leading to a significant increase of the g X -value of up to 0.003, but are not sensitive to the second proton transfer (E2PT). In the latter case, the change of the g X -value is much smaller (ca. 0.0001), which is too small to be distinguished even by high-frequency EPR. The 14 N hyperfine values are also too similar to allow differentiation between the different protonation states in EPT and E2PT. The magnetic resonance parameters were also calculated as a function of the rotation angles around single bonds. It was demonstrated that rotation of the phenoxyl group results in large positive changes (> 0.001) in the g X -values. Analysis of the data reveals that the main source of these changes is related to the strength of the H-bond between phenoxyl oxygen and the proton(s) on N 1 and N 2 positions of the imidazole.

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

Preferential Direction of Electron Transfers at a Dye–Metal Oxide Interface with an Insulating Fluorinated Self-Assembled Monolayer and MgO

Slowing non-desirable electron transfer reactions at metal oxide-dye interfaces is important for many technologies. Specifically, after an interfacial photoinduced charge separation event at a metal oxide-dye interface, it is critically important to limit the rate of electron transfer reactions back to sensitizers and to limit electron transfer reactions between the electrolyte and the metal oxide. Ruthenium-based dyes at metal oxide interfaces are widely used in many fields; however, these dyes often have poor surface insulation resulting in fast recombination kinetics with transition metal-based redox shuttles in an electrolyte. This work explores two semiconductor surface modification strategies designed to minimize recombination events of electrons in TiO 2 with oxidized redox shuttles using a fluorinated siloxane insulator (PFTS) and a metal oxide insulator (MgO) with a well-known Ru dye, B11. Additionally, the influence of these treatments on the rate and duration of photoinduced interfacial charge separation at the TiO 2 -dye interface was examined. The TiO 2 -dye-redox shuttle systems were studied via dye-sensitized solar cell current-voltage curve, incident photon-to-current conversion efficiency, small-modulated photovoltage transient, time-correlated single photon counting, and transient absorption spectroscopy measurements. MgO was found to decrease the rate of the electron transfer reaction from the metal oxide to the electrolyte, decrease the rate of the electron transfer reaction to the oxidized dye from TiO 2 , increase the electron transfer reaction rate from a reduced redox shuttle to an oxidized dye, and decrease the electron injection rate from the photoexcited dye to TiO2. Interestingly, PFTS was found to desirably improve on these rates relative to MgO or untreated TiO 2 . Here, a model based on electrostatic interactions is presented to explain the exceptional behavior of PFTS with DFT computational analysis of PFTS supporting this model.

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

Application of Nuclear Technology to Art Identification Problems: First Annual Report

Throughout the centuries, as man has become increasingly affluent, his interest in antiquities and objects of art, and the price he has been willing to pay for them, has increased. Inevitably, one result has been to attract the unscrupulous to the production and sale of counterfeit paintings, sculptures, and other works of art. The skill of forgers varies but often is highly sophisticated. However, dealers, museum curators, and other experts have also become more skillful in the detection of forgeries. Scientific tools of increasing sensitivity and sophistication have gradually been applied to the examination of the materials of art and archaeology. Such tools frequently support and render more reliable the judgment of the experts. The quality of forgeries may improve further as their makers become acquainted with the new methods of detection and in turn learn to circumvent or confound the methods. Therefore, the development of still more advanced methods of examination has great utility in the art world. This is particularly true if the ultimate effect is to make circumvention of the methods of examination so costly that the economic incentives for producing forgeries may be substantially reduced, if not eliminated.

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

Nanoscale solid-state nuclear quadrupole resonance spectroscopy using depth-optimized nitrogen-vacancy ensembles in diamond

Nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) spectroscopy of bulk quantum materials have provided insight into phenomena, such as quantum phase criticality, magnetism, and superconductivity. With the emergence of nanoscale 2D materials with magnetic phenomena, inductively detected NMR and NQR spectroscopy are not sensitive enough to detect the smaller number of spins in nanomaterials. The nitrogen-vacancy (NV) center in diamond has shown promise in bringing the analytic power of NMR and NQR spectroscopy to the nanoscale. However, due to depth-dependent formation efficiency of the defect centers, noise from surface spins, band bending effects, and the depth dependence of the nuclear magnetic field, there is ambiguity regarding the ideal NV depth for surface NMR of statistically polarized spins. In this work, we prepared a range of shallow NV ensemble layer depths and determined the ideal NV depth by performing NMR spectroscopy on statistically polarized 19 F in Fomblin oil on the diamond surface. We found that the measurement time needed to achieve a signal-to-noise ratio of 3 using XY8-N noise spectroscopy has a minimum at an NV ensemble depth of 5.5 ± 1.5 nm for ensembles activated from 100 ppm nitrogen concentration. To demonstrate the sensing capabilities of NV ensembles, we perform NQR spectroscopy on the 11 B of hexagonal boron nitride flakes. We compare our best diamond to previous work with a single NV and find that this ensemble provides a shorter measurement time with excitation diameters as small as 4 μm. Furthermore, this analysis provides ideal conditions for further experiments involving NMR/NQR spectroscopy of 2D materials with magnetic properties.

2D materials↗
NASA NTRS2019

Experimental Demonstration and System Analysis for Plasmonic Force Propulsion

This project developed a higher-fidelity model of a recently envisioned small spacecraft propulsion system for precision pointing and proximity control. Plasmonic force propulsion harnesses solar light focused onto plasmon reactive subwavelength nanostructures to accelerate and expel nanoparticle propellant via strong optical forces. The goal of the project was to show that plasmonic space propulsion can provide the level of proximity and attitude control envisioned for future NASA nano/picosatellite missions, a level that is better than state-of-the-art approaches. We achieved this goal by showing that plasmonic force thrusters are feasible for a range of advanced mission concepts requiring swarm formations in a deep space environment. We performed three case studies that evaluated the performance of the plasmonic force propulsion thruster in a deep space, microsatellite swarm formation. These case studies assumed the propulsion system could generate thrust at the level predicted from our Phase 1 study (1.6 μN). Through these cases we were able to analyze the concept within a mission specific context through detailed orbital dynamics calculations. Results indicate that, with the Phase 1 estimated thrust level, the approach is promising for providing attitude control to swarm formation spacecraft. Further, we achieved goals related to technology development. Specifically, we experimentally demonstrated nanoparticle acceleration due to plasmonic forces with asymmetric nanostructures excited by focused laser light. Additionally, we investigated the thrust sensitivity and nanoparticle propellant injection dependencies upon thermal effects. As a result of our study, plasmonic force propulsion is at an early TRL 3. Active research and design has been conducted analytically and in the laboratory. Furthermore, practical applications such as the three case studies have been identified for the scientific basic principles that were observed. Future efforts related to fundamental understanding of these techniques should focus on 1) developing a standalone array of asymmetric nanostructures that can effectively interact with a stream or reservoir of particles or 2) experimentally evaluate a dielectrophoretic injector for nanoparticle propellant. The main limitation discovered about plasmonic propulsion regards performance estimates significantly below the Phase 1 estimations. Specifically, original assumptions in the Phase 1 project (notably, a linear array of asymmetric nanostructures) is not a viable approach to achieving significant acceleration, high exhaust velocity, of nanoparticles. More specifically, we assumed in Phase 1 that nanoparticles would be accelerated in series by a long linear array of asymmetric nanostructures. That is, the acceleration of the nanoparticle would build and increase with the kick received by each subsequent nanostructure. This is fundamentally flawed. The potential profile of a single nanostructure is such that it prohibits this phenomenon. The potential energy associated with the plasmon-generated dielectrophoretic force is a potential well, which is good for trapping nanoparticles, but cannot provide significant acceleration of particles to expel them out and away from the nanostructure. Further, a nanoparticle expelled from the first nanostructure would need to overcome the potential barrier for entry into the next nanostructure accelerating stage. Fundamentally, this effect means that a linear array of nanostructures is not a viable accelerating structure. Correspondingly then, acceleration can, or should, only be provided by one nanostructure, and the net acceleration and thrust force of a single nanostructure is small (~cm/s exhaust velocities, sub-nN level thrust vs. the 100’s m/s, μN originally envisioned). While our experiments demonstrated acceleration and manipulation of a nanoparticle using laser light in aqueous environment, the achievable energy and momentum addition to the nanoparticle from a single nanostructure stage is too low for useful propulsion. In terms of thrust prediction, the estimated thrust of 1.6 μN in Phase 1 is reduced to a few nN of thrust with this new insight and understanding of the concept. This thrust level is too small to achieve attitude control of swarms as originally envisioned.

Plasmonic Force↗
DOE OSTI2022

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

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

Simulation and analysis of the 3D double shock-induced recompression of spall using a hybrid PAGOSA with FLIP+MPM

In this work, an innovative hybrid algorithm PAGOSA with FLIP + MPM is first presented by coupling PAGOSA with the particle FLIP + MPM, and applied to systematically exploring the 3D recompression of spall induced by double-shock waves in ductile materials as the few existing studies are currently limited to the pure longitudinal recompression of spall, and ignore the 3D effects. This algorithm solves the difficulties the grid-based methods encounter when applied to capturing the fracture in material. We first validate the capabilities of PAGOSA with FLIP + MPM to predict different fracture/fragmentation cases by solving two benchmark problems and comparing with the analytical solutions or the experiment results. The convergences and the infinity-norm errors are also investigated. Subsequently, the 3D effects are demonstrated by simulating the spallation in material using PAGOSA with FLIP + MPM and comparing with those in the longitudinal cases. Additionally, some factors that affect the spallation with 3D effects are also discussed. Finally, the recompressions of spall in material with/without 3D effects are systematically simulated and analyzed. The recompression of spall driven by a high explosive is also simulated to show the ability of PAGOSA with FLIP + MPM to handle this kind of problem. The conditions that can lead to a recompression of spall are explored. Also, numerical results show that PAGOSA with FLIP + MPM can accurately predict different fracture cases, that the 3D effects play an important role in fracture/fragmentation, and the recompression of spall is very sensitive to the occurrence conditions, shows great differences when considering 3D effects in real applications. Moreover, the present hybrid method can be easily extended to other grid-based techniques employed for fracture in materials.

3D effects↗
NASA NTRS2012

International Space Station 2A Array Modal Analysis

On December 9th 2009, the International Space Station (ISS) 2A solar array mast experienced prolonged longeron shadowing during a Soyuz undocking. Analytical reconstruction of induced thermal and dynamic structural loads showed an exceedance of the mast buckling limit. Possible structural damage to the solar array mast could have occurred during this event. A Low fidelity video survey of the 2A mast showed no obvious damage of the mast longerons or battens. The decision was made to conduct an on-orbit dynamic test of the 2A array on December 18th, 2009. The test included thruster pluming on the array while photogrammetry data was recorded. The test was similar to other Dedicated Thruster Firings (DTFs) that were performed to measure structural frequency and damping of a solar array. Results of the DTF indicated lower frequency mast modes than model predictions, thus leading to speculation of mast damage. A detailed nonlinear analysis was performed on the 2A array model to assess possible solutions to modal differences. The setup of the parametric nonlinear trade study included the use of a detailed array model and the reduced mass and stiffness matrices of the entire ISS being applied to the array interface. The study revealed that the array attachment structure is nonlinear and thus was the source of error in the model prediction of mast modes. In addition, a detailed study was performed to determine mast mode sensitivity to mast longeron damage. This sensitivity study was performed to assess if the ISS program has sufficient instrumentation for mast damage detection.

Laible, Michael↗