Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability generating function”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Genetic Algorithm for Optimization: Preprocessor and Algorithm

Genetic algorithm (GA) inspired by Darwin's theory of evolution and employed to solve optimization problems - unconstrained or constrained - uses an evolutionary process. A GA has several parameters such the population size, search space, crossover and mutation probabilities, and fitness criterion. These parameters are not universally known/determined a priori for all problems. Depending on the problem at hand, these parameters need to be decided such that the resulting GA performs the best. We present here a preprocessor that achieves just that, i.e., it determines, for a specified problem, the foregoing parameters so that the consequent GA is a best for the problem. We stress also the need for such a preprocessor both for quality (error) and for cost (complexity) to produce the solution. The preprocessor includes, as its first step, making use of all the information such as that of nature/character of the function/system, search space, physical/laboratory experimentation (if already done/available), and the physical environment. It also includes the information that can be generated through any means - deterministic/nondeterministic/graphics. Instead of attempting a solution of the problem straightway through a GA without having/using the information/knowledge of the character of the system, we would do consciously a much better job of producing a solution by using the information generated/created in the very first step of the preprocessor. We, therefore, unstintingly advocate the use of a preprocessor to solve a real-world optimization problem including NP-complete ones before using the statistically most appropriate GA. We also include such a GA for unconstrained function optimization problems.

Sen, S. K.↗
NASA NTRS1999

Statistical Perspectives on Stratospheric Transport

Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Sparling, L. C.↗
NASA NTRS1996

Shape deformation of the organ of Corti associated with length changes of outer hair cell

Cochlear outer hair cells (OHC) are commonly assumed to function as mechanical effectors as well as sensory receptors in the organ of Corti (OC) of the inner ear. OHC in vitro and in organ explants exhibit mechanical responses to electrical, chemical or mechanical stimulation which may represent an aspect of their effector process that is expected in vivo. A detailed description, however, of an OHC effector operation in situ is still missing. Specifically, little is known as to how OHC movements influence the geometry of the OC in situ. Previous work has demonstrated that the motility of isolated OHCs in response to electrical stimulation and to K(+)-gluconate is probably under voltage control and causes depolarisation (shortening) and hyperpolarization (elongation). This work was undertaken to investigate if the movements that were observed in isolated OHC, and which are induced by ionic stimulation, could change the geometry of the OC. A synchronized depolarization of OHC was induced in guinea pig cochleae by exposing the entire OC to artificial endolymph (K+). Subsequent morphometry of mid-modiolar sections from these cochleae revealed that the distance between the basilar membrane (BM) and the reticular lamina (RL) had decreased considerably. Furthermore, in the three upper turns OHC had significantly shortened in all rows. The results suggest that OHC can change their length in the organ of Corti (OC) thus deforming the geometry of the OC. The experiments reveal a tonic force generation within the OC that may change the position of RL and/or BM, contribute to damping, modulate the BM-RL-distance and control the operating points of RL and sensory hair bundles. Thus, the results suggest active self-adjustments of cochlear mechanics by slow OHC length changes. Such mechanical adjustments have recently been postulated to correspond to timing elements of animal communication, speech or music.

Non-NASA Center↗
NASA NTRS2005

Formation of Nanophase Iron in Lunar Soil Simulant for Use in ISRU Studies

For the prospective return of humans to the Moon and the extensive amount of premonitory studies necessary, large quantities of lunar soil simulants are required, for a myriad of purposes from construction/engineering purposes all the way to medical testing of its effects from ingestion by humans. And there is only a limited and precious quantity of lunar soil available on Earth (i.e., Apollo soils) - therefore, the immediate need for lunar soil simulants. Since the Apollo era, there have been several simulants; of these JSC-1 (Johnson Space Center) and MLS-1 (Minnesota Lunar Simulant) have been the most widely used. JSC-1 was produced from glassy volcanic tuff in order to approximate lunar soil geotechnical properties; whereas, MLS-1 approximates the chemistry of Apollo 11 high-Ti soil, 10084. Stocks of both simulants are depleted, but JSC-1 has recently gone back into production. The lunar soil simulant workshop, held at Marshall Space Flight Center in January 2005, identified the need to make new simulants for the special properties of lunar soil, such as nanophase iron (np-Fe(sup 0). Hill et al. (2005, this volume) showed the important role of microscale Fe(sup 0) in microwave processing of the lunar soil simulants JSC-1 and MLS-1. Lunar soil is formed by space weathering of lunar rocks (e.g., micrometeorite impact, cosmic particle bombardment). Glass generated during micrometeorite impact cements rock and mineral fragments together to form aggregates called agglutinates, and also produces vapor that is deposited and coats soil grains. Taylor et al. (2001) showed that the relative amount of impact glass in lunar soil increases with decreasing grain size and is the most abundant component in lunar dust (less than 20 micrometer fraction). Notably, the magnetic susceptibility of lunar soil also increases with the decreasing grain size, as a function of the amount of nanophase-sized Fe(sup 0) in impact-melt generated glass. Keller et al. (1997, 1999) also discovered the presence of abundant np-Fe(sup 0) particles in the glass patinas coating most soil particles. Therefore, the correlation of glass content and magnetic susceptibility can be explained by the presence of the np-Feo particles in glass: small particles contain relatively more np-Fe(sup 0) as glass coatings because the surface area versus mass ratio of the grain size is so increased. The magnetic properties of lunar soil are important in dust mitigation on the Moon (Taylor et al. 2005). Thus material simulating this property is important for testing mitigation methods using electromagnetic field. This np- Fe(sup 0) also produces a unique energy coupling to normal microwaves, such as present in kitchen microwave ovens. Effectively, a portion of lunar soil placed in a normal 2.45 GHz oven will melt at greater than 1200 C before your tea will boil at 100 C, a startling and new discovery reported by Taylor and Meek (2004, 2005). Several methods have been investigated in attempts to make nanophase-sized Feo dispersed within silicate glass; like in the lunar glass. We have been successful in synthesizing such a product and continue to improve on our recipe. We have performed extensive experimentation on this subject to date. Ultimately it will probably be necessary to add this np-Fe(sup 0) bearing silicate glass to lunar soil stimulant, like JSC-1, to actually produce the desired magnetic and microwave coupling properties for use in appropriate ISRU experimentation.

Liu, Yang↗
NASA NTRS1991

Self assembly properties of primitive organic compounds

A central event in the origin of life was the self-assembly of amphiphilic, lipid-like compounds into closed microenvironments. If a primitive macromolecular replicating system could be encapsulated within a vesicular membrane, the components of the system would share the same microenvironment, and the result would be a step toward true cellular function. The goal of our research has been to determine what amphiphilic molecules might plausibly have been available on the early Earth to participate in the formation of such boundary structures. To this end, we have investigated primitive organic mixtures present in carbonaceous meteorites such as the Murchison meteorite, which contains 1-2 percent of its mass in the form of organic carbon compounds. It is likely that such compounds contributed to the inventory of organic carbon on the prebiotic earth, and were available to participate in chemical evolution leading to the emergence of the first cellular life forms. We found that Murchison components extracted into non-polar solvent systems are surface active, a clear indication of amphiphilic character. One acidic fraction self-assembles into vesicular membranes that provide permeability barriers to polar solutes. Other evidence indicates that the membranes are bimolecular layers similar to those formed by contemporary membrane lipids. We conclude that bilayer membrane formation by primitive amphiphiles on the early Earth is feasible. However, only a minor fraction of acidic amphiphiles assembles into bilayers, and the resulting membranes require narrowly defined conditions of pH and ionic composition to be stable. It seems unlikely, therefore, that meteoritic infall was a direct source of membrane amphiphiles. Instead, the hydrocarbon components and their derivatives more probably would provide an organic stock available for chemical evolution. Our current research is directed at possible reactions which would generate substantial quantities of membranogenic amphiphiles. One possibility is photochemical oxidation of hydrocarbons.

Deamer, D. W.↗
DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Utilization of Machine Learning Techniques for Managing the Tracking and Data Relay Satellite Constellation

National Aeronautics and Space Administration’s (NASA) Goddard Space Flight Center (GSFC) operates a constellation of ten geosynchronous Tracking and Data Relay Satellites (TDRS). The TDRS constellation consists of multiple geosynchronous communication relay satellites located around the equator so they can provide continual coverage of any mission in low earth orbit. The TDRS are located primarily in three oceanic regions around the earth. NASA’s White Sands Complex provides the ground communication support for TDRS located over the Atlantic and Pacific Oceans. Another TDRS ground station in Guam supports the TDRS over the Indian Ocean. With these satellites the TDRS network can provide continuous coverage of satellites in low-earth orbit. The NASA Space Network (SN) project office at GSFC manages the constellation of spacecraft. Major customers of the TDRS constellation include, but are not limited to, the International Space Station and the Hubble Space Telescope. The TDRS constellation has three generations of satellites and has been active for over 30 years providing reliable communication links between customer satellites and corresponding ground stations. However, one of the major concerns for TDRS, and in any space mission, is to ensure the health and safety of the spacecraft. Generally, engineers use telemetry data to monitor and analyze the performance and state of health of the spacecraft. Telemetry data contains hundreds of parameters that monitor each important component in the spacecraft, which can be utilized to recognize and characterize the behavior of the spacecraft. Each parameter contains considerable information to represent time-dependent properties of each spacecraft subsystem and component. During the entire life of a TDRS spacecraft, thousands of gigabytes of telemetry data are transmitted in real-time from the spacecraft to the ground station at the White Sands Complex in Las Cruces, New Mexico, and recorded as historical data sets for engineers to process and analyze the events that occurred on-orbit. These parameters contain the function of multiple spacecraft subsystems, such as the attitude control system (ACS), Thermal, Electrical Power Subsystem (EPS), etc. . The first and second generations have exceeded their required lifetime and NASA is keen to manage these spacecrafts carefully in order to maximize the remaining life using the spacecraft telemetry. The challenge is to know when the risk of losing a spacecraft in geosynchronous orbit exceeds the benefit of continued operations for customer support. In the TDRS fleet, the EPS is the most critical subsystem related to spacecraft operations. Failure of the EPS would strand a spacecraft in geosynchronous orbit. Since EPS provides power to the spacecraft, component failures ultimately lead to the inability to support the spacecraft loads and the communications payload. For instance, TDRS-8 has several anomalies in EPS including the Bus Voltage Limiter (BVL) shunt current, solar array loss of circuits, and failed battery cells. Any of these anomalies can cause critical issues to the spacecraft. Therefore, developing a system to analyze and perform early detection of a potential anomaly is an important issue in telemetry data analysis. In recent years, Telemetry Mining (TM) has been proposed to process telemetry data by using Data Mining (DM) techniques such as classification, clustering, regression and anomaly detection. Anomaly detection, also known as outlier detection, has been widely used in many data mining areas such as remote sensing, medical data processing and digital image processing. The goal of anomaly detection is to detect abnormal data, which contains a relatively low probability of occurrence among the entire data set. Early detection of anomalies is one of the most significant issues in managing the spacecraft configuration. If anomalies can be detected early enough, then the redundant resources can be used to extend the life of the operational spacecraft. We present an unsupervised anomaly detection method to process the EPS data extracted from TDRS-8. This is different from traditional analytical methods, which use telemetry data to illustrate behavior and physical meaning of each spacecraft component. TM connects multiple parameters as a vector and then conducts data analysis on this high dimension telemetry vector. This method is looking at the properties of a high dimensional vector that is able to consider the relationship between different parameters in the anomaly detection problem. This kind of method performs much better than the traditional limit checking method. In addition, we propose a new approach of real-time anomaly detection to process telemetry data in real-time, which can then be applied to spacecraft monitoring with high reliability, low cost and high accuracy.

Machine Learning (ML)↗
DOE OSTI2022

Hot and dense matter equation of state probability distributions for astrophysical simulations

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. Furthermore, we compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1999

Drop Breakup in Fixed Bed Flows as Model Stochastic Flow Fields

We examine drop breakup in a class of stochastic flow fields as a model for the flow through fixed fiber beds and to elucidate the general mechanisms whereby drops breakup in disordered, Lagrangian unsteady flows. Our study consists of two parallel streams of investigation. First, large scale numerical simulations of drop breakup in a class of anisotropic Gaussian fields will be presented. These fields are generated spectrally and have been shown in a previous publication to be exact representations of the flow in a dilute disordered bed of fibers if close interactions between the fibers and the drops are dynamically unimportant. In these simulations the drop shape is represented by second and third order small deformation theories which have been shown to be excellent for the prediction of drop breakup in steady strong flows. We show via these simulations that the mechanisms of drop breakup in these flows are quite different than in steady flows. The predominant mechanism of breakup appears to be very short lived twist breakups. Moreover, the occurrence of breakup events is poorly predicted by either the strength of the local flow in which the drop finds itself at breakup, or the degree of deformation that the drop achieves prior to breakup. It is suggested that a correlation function of both is necessary to be predictive of breakup events. In the second part of our research experiments are presented where the drop deformation and breakup in PDMS/polyisobutylene emulsions is considered. We consider very dilute emulsions such that coalescence is unimportant. The flows considered are simple shear and the flow through fixed fiber beds. Turbidity, small angle light scattering, dichroism and microscopy are used to interrogate the drop deformation process in both flows. It is demonstrated that breakup at very low capillary numbers occurs in both flows but larger drop deformation occurs in the fixed bed flow. Moreover, it is witnessed that breakup in the bed occurs continuously during flow and apparently with uniform probability through the bed length. The drop deformations witnessed in our experiments are larger than those predicted by the numerical simulations, and future plans to investigate these differences are discussed.

Shaqfeh, Eric S. G.↗
NASA NTRS2001

Extracellular calcium sensing and extracellular calcium signaling

The cloning of a G protein-coupled extracellular Ca(2+) (Ca(o)(2+))-sensing receptor (CaR) has elucidated the molecular basis for many of the previously recognized effects of Ca(o)(2+) on tissues that maintain systemic Ca(o)(2+) homeostasis, especially parathyroid chief cells and several cells in the kidney. The availability of the cloned CaR enabled the development of DNA and antibody probes for identifying the CaR's mRNA and protein, respectively, within these and other tissues. It also permitted the identification of human diseases resulting from inactivating or activating mutations of the CaR gene and the subsequent generation of mice with targeted disruption of the CaR gene. The characteristic alterations in parathyroid and renal function in these patients and in the mice with "knockout" of the CaR gene have provided valuable information on the CaR's physiological roles in these tissues participating in mineral ion homeostasis. Nevertheless, relatively little is known about how the CaR regulates other tissues involved in systemic Ca(o)(2+) homeostasis, particularly bone and intestine. Moreover, there is evidence that additional Ca(o)(2+) sensors may exist in bone cells that mediate some or even all of the known effects of Ca(o)(2+) on these cells. Even more remains to be learned about the CaR's function in the rapidly growing list of cells that express it but are uninvolved in systemic Ca(o)(2+) metabolism. Available data suggest that the receptor serves numerous roles outside of systemic mineral ion homeostasis, ranging from the regulation of hormonal secretion and the activities of various ion channels to the longer term control of gene expression, programmed cell death (apoptosis), and cellular proliferation. In some cases, the CaR on these "nonhomeostatic" cells responds to local changes in Ca(o)(2+) taking place within compartments of the extracellular fluid (ECF) that communicate with the outside environment (e.g., the gastrointestinal tract). In others, localized changes in Ca(o)(2+) within the ECF can originate from several mechanisms, including fluxes of calcium ions into or out of cellular or extracellular stores or across epithelium that absorb or secrete Ca(2+). In any event, the CaR and other receptors/sensors for Ca(o)(2+) and probably for other extracellular ions represent versatile regulators of numerous cellular functions and may serve as important therapeutic targets.

Review↗
NASA NTRS2014

Smooth Pursuit of Flicker-Defined Motion

We examined the pursuit response to stimuli defined by space-variant flicker of a dense random dot carrier pattern. On each frame, every element of the pattern could change polarity, with a probability given by a two-dimensional Gaussian distribution. A normal distribution produces a circular region of twinkle, while inverting the distribution results in a spot of static texture in a twinkling surround. In this latter case, the carrier texture could be stationary, or could move with the twinkle modulator, thereby producing first-order motion in the region of the spot. While the twinkle-defined spot produces a strong sensation of motion, the complementary stimulus defined by the absence of twinkle does not, when viewed peripherally, it appears to move in steps even when the generating distribution moves smoothly. We examined pursuit responses to these stimuli using two techniques: 1) the eye movement correlogram, obtained by cross-correlating eye velocity with the velocity of a randomly-moving stimulus; and 2) delayed visual feedback, where transient stabilization of a target can produce spontaneous oscillations of the eye, with a period empirically observed to vary linearly with the applied delay. Both techniques provide an estimate of the internal processing time, which can be as short as 100 milliseconds for a first-order target. Assessed by the correlogram method, the response to flicker-defined motion is delayed by more than 100 milliseconds, and significantly weaker (especially in the vertical dimension). When initially presented in the delayed feedback condition, purely saccadic oscillation is observed. One subject eventually developed smooth oscillations (albeit with significant saccadic intrusions), showing a period-versus-delay slope similar to that observed for first-order targets. This result is somewhat surprising, given that we interpret the slope of the period-versus-delay-function as reflecting the balance between position- and velocity-sensitive inputs to pursuit.

eye movements↗
DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling

In the Monte Carlo many-body perturbation (MC-MP) method, the conventional correlation-correction formula, which is a long sum of products of low-dimensional integrals, is first recast into a short sum of high-dimensional integrals over electron-pair and imaginary-time coordinates. These high-dimensional integrals are then evaluated by the Monte Carlo method with random coordinates generated by the Metropolis–Hasting algorithm according to a suitable distribution. The latter algorithm, while advantageous in its ability to sample nearly any distribution, introduces autocorrelation in sampled coordinates, which in turn increases the statistical uncertainty of the integrals and thus the computational cost. It also involves wasteful rejected moves and an initial “burn-in” step as well as displays hysteresis. Here, an algorithm is proposed that directly produces a random sequence of electron-pair coordinates for the same distribution used in the MC-MP method, which is free from autocorrelation, rejected moves, a burn-in step, or hysteresis. Furthermore, this direct-sampling algorithm is shown to accelerate second- (MC-MP2) and third-order Monte Carlo many-body perturbation (MC-MP3) calculations by up to 222% and 38%, respectively.

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

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Enigma Version 12

Enigma Version 12 software combines model building, animation, and engineering visualization into one concise software package. Enigma employs a versatile user interface to allow average users access to even the most complex pieces of the application. Using Enigma eliminates the need to buy and learn several software packages to create an engineering visualization. Models can be created and/or modified within Enigma down to the polygon level. Textures and materials can be applied for additional realism. Within Enigma, these models can be combined to create systems of models that have a hierarchical relationship to one another, such as a robotic arm. Then these systems can be animated within the program or controlled by an external application programming interface (API). In addition, Enigma provides the ability to use plug-ins. Plugins allow the user to create custom code for a specific application and access the Enigma model and system data, but still use the Enigma drawing functionality. CAD files can be imported into Enigma and combined to create systems of computer graphics models that can be manipulated with constraints. An API is available so that an engineer can write a simulation and drive the computer graphics models with no knowledge of computer graphics. An animation editor allows an engineer to set up sequences of animations generated by simulations or by conceptual trajectories in order to record these to highquality media for presentation. Enigma Version 12 Lyndon B. Johnson Space Center, Houston, Texas 28 NASA Tech Briefs, September 2013 Planetary Protection Bioburden Analysis Program NASA's Jet Propulsion Laboratory, Pasadena, California This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware.

Shores, David↗
DOE OSTI2021

Overview of the fusion nuclear science facility, a credible break-in step on the path to fusion energy

The Fusion Nuclear Science Facility (FNSF) is examined here as part of a two step program from ITER to commercial power plants. This first step is considered mandatory to establish the materials and component database in the real fusion in-service environment before proceeding to larger electricity producing facilities. The FNSF can be shown to make tremendous advances beyond ITER, toward a power plant, particularly in plasma duration and fusion nuclear environment. A moderate FNSF is studied in detail, which does not generate net electricity, but does reach the power plant blanket operating temperatures. The full poloidal Dual Coolant Lead Lithium (DCLL) blanket is chosen, with alternates being the Helium Cooled Lead Lithium (HCLL) and Helium Cooled Ceramic Breeder/Pebble Bed (HCCB/PB). Several power plant relevant choices are made in order to follow the philosophy of targeted technologies. Any fusion core component must be qualified by fusion relevant neutron testing and highly integrated non-nuclear testing before it can be installed on the FNSF in order to avoid the high probability of constant failures in a plasma-vacuum system. A range of missions for the FNSF, or any fusion nuclear facility on the path toward fusion power plants, are established and characterized by several metrics. A conservative physics strategy is pursued to accommodate the transition to ultra-long plasma pulses, and parameters are chosen to represent the power plant regime to the extent possible. An operating space is identified, and from this, one point is chosen for further detailed analysis, with R = 4.8 m, a = 1.2 m, IP = 7.9 MA, BT = 7.5 T, βN Gr = 0.9, fBS = 0.52, q95 = 6.0, H98 ∼1.0, and Q = 4.0. The operating space is shown to be robust to parameter variations. A program is established for the FNSF to show how the missions for the facility are met, with a He/H, a DD and 5 DT phases. The facility requires ∼25 years to complete its DT operation, including 7.8 years of neutron production, and the remaining spent on inspections and maintenance. The DD phase is critical to establish the ultra-long plasma pulse lengths. The blanket testing strategy is examined, and shows that many sectors have penetrations for heating and current drive (H/CD), diagnostics, or Test Blanket Modules (TBMs). The hot cell is a critical facility element in order for the FNSF to perform its function of developing the in-service material and component database. The pre-FNSF R&D is laid out in terms of priority topics, with the FNSF phases driving the time-lines for R&D completion. A series of detailed technical assessments of the FNSF operating point are reported in this issue, showing the credibility of such a step, and more detailed emphasis on R&D items to pursue. These include nuclear analysis, thermo-mechanics and thermal-hydraulics, liquid metal thermal hydraulics, transient thermo-mechanics, tritium analysis, maintenance assessment, magnet specification and analysis, materials assessments, core and scrape-off layer (SOL)/divertor plasma examinations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Diverse Chemo-Dynamical Properties of Nitrogen-Rich Stars Identified from Low-Resolution Spectra

The second generation of stars in the globular clusters (GCs) of the Milky Way (MW) exhibit unusually high N, Na, or Al, compared to typical Galactic halo stars at similar metallicities. The halo field stars enhanced with such elements are believed to have originated in disrupted GCs or escaped from existing GCs. We identify such stars in the metallicity range –3.0 < [Fe/H] < 0.0 from a sample of ~ 36,800 giant stars observed in the Sloan Digital Sky Survey and Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, and present their dynamical properties. The N-rich population (NRP) and N-normal population (NNP) among our giant sample do not exhibit similarities in either in their metallicity distribution function (MDF) or dynamical properties. We find that, even though the MDF of the NRP looks similar to that of the MW’s GCs in the range of [Fe/H] < –1.0, our analysis of the dynamical properties does not indicate similarities between them in the same metallicity range, implying that the escaped members from existing GCs may account for a small fraction of our N-rich stars, or the orbits of the present GCs have been altered by the dynamical friction of the MW. We also find a significant increase in the fraction of N-rich stars in the halo field in the very metal-poor (VMP; [Fe/H] < –2.0) regime, comprising up to ~ 20% of the fraction of the N-rich stars below [Fe/H] = –2.5, hinting that partially or fully destroyed VMP GCs may have in some degree contributed to the Galactic halo. A more detailed dynamical analysis of the NRP reveals that our sample of N-rich stars do not share a single common origin. Although a substantial fraction of the N-rich stars seem to originate from the GCs formed in situ, more than 60% of them are not associated with those of typical Galactic populations, but probably have extragalactic origins associated with Gaia Sausage/Enceladus, Sequoia, and Sagittarius dwarf galaxies, as well as with presently unrecognized progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Mars Surface Ionizing Radiation Environment: Need for Validation

Protection against the hazards from exposure to ionizing radiation remains an unresolved issue in the Human Exploration and Development of Space (HEDS) enterprise [1]. The major uncertainty is the lack of data on biological response to galactic cosmic ray (GCR) exposures but even a full understanding of the physical interaction of GCR with shielding and body tissues is not yet available and has a potentially large impact on mission costs. "The general opinion is that the initial flights should be short-stay missions performed as fast as possible (so-called 'Sprint' missions) to minimize crew exposure to the zero-g and space radiation environment, to ease requirements on system reliability, and to enhance the probability of mission success." The short-stay missions tend to have long transit times and may not be the best option due to the relatively long exposure to zero-g and ionizing radiation. On the other hand the short-transit missions tend to have long stays on the surface requiring an adequate knowledge of the surface radiation environment to estimate risks and to design shield configurations. Our knowledge of the surface environment is theoretically based and suffers from an incomplete understanding of the physical interactions of GCR with the Martian atmosphere, Martian surface, and intervening shield materials. An important component of Mars surface robotic exploration is the opportunity to test our understanding of the Mars surface environment. The Mars surface environment is generated by the interaction of Galactic Cosmic Rays (GCR) and Solar Particle Events (SPEs) with the Mars atmosphere and Mars surface materials. In these interactions, multiple charged ions are reduced in size and secondary particles are generated, including neutrons. Upon impact with the Martian surface, the character of the interactions changes as a result of the differing nuclear constituents of the surface materials. Among the surface environment are many neutrons diffusing from the Martian surface and especially prominent are energetic neutrons with energies up to a few hundred MeV. Testing of these computational results is first supported by ongoing experiments at the Brookhaven National Laboratory but equally important is the validation to the extent possible by measurements on the Martian surface. Such measurements are limited by power and weight requirements of the specific mission and simplified instrumentation by necessity lacks the full discernment of particle type and spectra as is possible with laboratory experimental equipment. Yet, the surface measurements are precise and a necessary requisite to validate our understanding of the surface environment. At the very minimum the surface measurements need to provide some spectral information on the neutron environment. Of absolute necessity is the precise knowledge of the detector response functions for absolute comparisons between the computational model of the surface environment and the detector measurements on the surface.

Wilson, J. W.↗