Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial power processing”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Hybrid Automotive Engine Using Ethanol-Burning Miller Cycle

A proposed hybrid (internal-combustion/ electric) automotive engine system would include as its internal-combustion subsystem, a modified Miller-cycle engine with regenerative air preheating and with autoignition like that of a Diesel engine. The fuel would be ethanol and would be burned lean to ensure complete combustion. Although the proposed engine would have a relatively low power-to-weight ratio compared to most present engines, this would not be the problem encountered if this engine were used in a non-hybrid system since hybrid systems require significantly lower power and thus smaller engines than purely internal-combustion-engine-driven vehicles. The disadvantage would be offset by the advantages of high fuel efficiency, low emission of nitrogen oxides and particulate pollutants, and the fact that ethanol is a renewable fuel. The original Miller-cycle engine, named after its inventor, was patented in the 1940s and is the basis of engines used in some modern automobiles, but is not widely known. In somewhat oversimplified terms, the main difference between a Miller-cycle engine and a common (Otto-cycle) automobile engine is that the Miller-cycle engine has a longer expansion stroke while retaining the shorter compression stroke. This is accomplished by leaving the intake valve open for part of the compression stroke, whereas in the Otto cycle engine, the intake valve is kept closed during the entire compression stroke. This greater expansion ratio makes it possible to extract more energy from the combustion process without expending more energy for compression. The net result is greater efficiency. In the proposed engine, the regenerative preheating would be effected by running the intake air through a heat exchanger connected to the engine block. The regenerative preheating would offer two advantages: It would ensure reliable autoignition during operation at low ambient temperature and would help to cool the engine, thereby reducing the remainder of the power needed for cooling and thereby further contributing to efficiency. An electrical resistance air preheater might be needed to ensure autoignition at startup and during a short warmup period. Because of the autoignition, the engine could operate without either spark plugs or glow plugs. Ethanol burns relatively cleanly and has been used as a motor fuel since the invention of internal-combustion engines. However, the energy content of ethanol per unit weight of ethanol is less than that of Diesel fuel or gasoline, and ethanol has a higher heat of vaporization. Because the Miller cycle offers an efficiency close to that of the Diesel cycle, burning ethanol in a Miller-cycle engine gives about as much usable output energy per unit volume of fuel as does burning gasoline in a conventional gasoline automotive engine. Because of the combination of preheating, running lean, and the use of ethyl alcohol, the proposed engine would generate less power per unit volume than does a conventional automotive gasoline engine. Consequently, for a given power level, the main body of the proposed engine would be bulkier. However, because little or no exhaust cleanup would be needed, the increase in bulk of the engine could be partially offset by the decrease in bulk of the exhaust system. The regenerative preheating also greatly reduces the external engine cooling requirement, and would translate to reduced engine bulk. It may even be possible to accomplish the remaining cooling of the engine by use of air only, eliminating the bulk and power consumption of a water cooling system. The combination of a Miller-cycle engine with regenerative air preheating, ethyl alcohol fuel, and hybrid operation could result in an automotive engine system that satisfies the need for a low pollution, high efficiency, and simple engine with a totally renewable fuel.

Weinstein, Leonard↗
DOE OSTI2021

Tuning Surface Stoichiometry of SOFC Electrodes at the Molecular and Nano-scale for Enhanced Performance and Durability

This project achieved the following objectives. Different cation segregation behaviors of different common SOFC cathodes, including La 0.6 Sr 0.4 Co0.2Fe 0.8 O 3-δ (LSCF), Sr 0.5 Sm 0.5 O 3-δ (SSC), and PrBa 0.5 Sr 0.5 Co 1.5 Fe 0.5 O5 +δ (PBSCF), were determined. We showed how oxygen partial pressure and gas impurities impact the stability of SOFC cathodes. We developed atomic layer deposition (ALD) coating techniques for electrodes and showed that by introducing different elements and different ALD oxidizers, electrode surface chemistry can be altered, resulting in enhanced oxygen reduction kinetics. In addition, we developed solution infiltration technique to enhance performance and durability of electrodes. Different infiltrates as well as different thermal treatment processes were screened to identify the optimal surface modification process that yields both low impedance and high durability. The modified cathode shows excellent stability and has a low area specific resistance (ASR) of only 0.2 Ωcm 2 at 600 °C after over 2000 hours of operation. Further, we developed ceramic anodes, SrFe(Co,Mo)O 3 (SFCM) and SrFe(Ni,Mo)O 3 (SFNM), and enhanced anode oxidation kinetics by solution infiltration or in situ catalyst exsolution from the ceramic anode surface. The modified anodes showed improved performance in full SOFCs, and the optimized anode modification shows high stability for over 400 hours at 550 °C. Moreover, the modified anode also demonstrates high durability in methane. This work provides fundamental understanding of electrode surface chemistry and demonstrates a simple, facile, cost-effective approach to enhance catalytic activity and durability of SOFC electrodes.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS2012

Implementation of a Transition Model in a NASA Code and Validation Using Heat Transfer Data on a Turbine Blade

The purpose of this report is to summarize and document the work done to enable a NASA CFD code to model laminar-turbulent transition process on an isolated turbine blade. The ultimate purpose of the present work is to down-select a transition model that would allow the flow simulation of a variable speed power turbine to be accurately performed. The flow modeling in its final form will account for the blade row interactions and their effects on transition which would lead to accurate accounting for losses. The present work only concerns itself with steady flows of variable inlet turbulence. The low Reynolds number k- model of Wilcox and a modified version of the same model will be used for modeling of transition on experimentally measured blade pressure and heat transfer. It will be shown that the k- model and its modified variant fail to simulate the transition with any degree of accuracy. A case is thus made for the adoption of more accurate transition models. Three-equation models based on the work of Mayle on Laminar Kinetic Energy were explored. The three-equation model of Walters and Leylek was thought to be in a relatively mature state of development and was implemented in the Glenn-HT code. Two-dimensional heat transfer predictions of flat plate flow and two-dimensional and three-dimensional heat transfer predictions on a turbine blade were performed and reported herein. Surface heat transfer rate serves as sensitive indicator of transition. With the newly implemented model, it was shown that the simulation of transition process is much improved over the baseline k- model for the single Reynolds number and pressure ratio attempted; while agreement with heat transfer data became more satisfactory. Armed with the new transition model, total-pressure losses of computed three-dimensional flow of E3 tip section cascade were compared to the experimental data for a range of incidence angles. The results obtained, form a partial loss bucket for the chosen blade. In time the loss bucket will be populated with losses at additional incidences. Results obtained thus far will be discussed herein.

Ameri, Ali A.↗
NASA NTRS2017

Modeling for Battery Prognostics

For any battery-powered vehicles (be it unmanned aerial vehicles, small passenger aircraft, or assets in exoplanetary operations) to operate at maximum efficiency and reliability, it is critical to monitor battery health as well performance and to predict end of discharge (EOD) and end of useful life (EOL). To fulfil these needs, it is important to capture the battery's inherent characteristics as well as operational knowledge in the form of models that can be used by monitoring, diagnostic, and prognostic algorithms. Several battery modeling methodologies have been developed in last few years as the understanding of underlying electrochemical mechanics has been advancing. The models can generally be classified as empirical models, electrochemical engineering models, multi-physics models, and molecular/atomist. Empirical models are based on fitting certain functions to past experimental data, without making use of any physicochemical principles. Electrical circuit equivalent models are an example of such empirical models. Electrochemical engineering models are typically continuum models that include electrochemical kinetics and transport phenomena. Each model has its advantages and disadvantages. The former type of model has the advantage of being computationally efficient, but has limited accuracy and robustness, due to the approximations used in developed model, and as a result of such approximations, cannot represent aging well. The latter type of model has the advantage of being very accurate, but is often computationally inefficient, having to solve complex sets of partial differential equations, and thus not suited well for online prognostic applications. In addition both multi-physics and atomist models are computationally expensive hence are even less suited to online application An electrochemistry-based model of Li-ion batteries has been developed, that captures crucial electrochemical processes, captures effects of aging, is computationally efficient, and is of suitable accuracy for reliable EOD prediction in a variety of operational profiles. The model can be considered an electrochemical engineering model, but unlike most such models found in the literature, certain approximations are done that allow to retain computational efficiency for online implementation of the model. Although the focus here is on Li-ion batteries, the model is quite general and can be applied to different chemistries through a change of model parameter values. Progress on model development, providing model validation results and EOD prediction results is being presented.

Prognostics↗
DOE OSTI2020

Spatially-Resolved Characterization Techniques and Their Implications for Nuclear Debris Formation.

Debris from nuclear tests has a complex formation process and can inform multiple fields of study such as geology, atmospheric science, shock physics, and chemistry under extreme conditions. Macroscopic nuclear debris forms as the result of fireball interaction with surrounding materials (e.g., structural and environmental component), which rapidly undergo melting and vaporization followed by condensation, convective/diffusive mixing, and solidification over the course of seconds. In atmospheric events, some of this material may disperse over long distances, but much of the material is deposited close-in to ground zero, often in multicomponent, partially (or entirely) amorphous debris formations. During the U.S. nuclear testing program, this material was collected and analyzed, particularly for radionuclide composition. Fallout formation models were developed from nuclear test data based on such radionuclide compositional analyses. These models were not just used to understand fallout dispersion (i.e. the spread of radioactivity over geographical regions) but also guide radiochemical interpretations of historical nuclear tests.6 Developments in analytical techniques over the past several decades have made it possible to make new analyses on historical debris, some of which may be more than a half-century old. For example, inductively-coupled plasma mass spectrometry (ICP-MS) has been used to look at the trace elements in glassy fallout material from aboveground nuclear tests (including the Trinity test). Advanced analyses using X-ray absorption allowed for the measurement of oxidation state of fallout constituents, including actinides. The ability to measure trace actinides and their subsequent oxidation state is important to understanding how the surrounding environment may have influenced the resultant fallout composition and may have implications for why certain fractionation trends have been observed. Here we present data from historic nuclear test debris illustrating the power of spatially resolved methods to connect interaction of the near-field environment with the explosion and provide new insights into nuclear debris formation.

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

Achievable Performance and Effective Interrogator Design for SAW RFID Sensor Tags

For many NASA missions, remote sensing is a critical application that supports activities such as environmental monitoring, planetary science, structural shape and health monitoring, non-destructive evaluation, etc. The utility of the remote sensing devices themselves is greatly increased if they are passive that is, they do not require any on-board power supply such as batteries and if they can be identified uniquely during the sensor interrogation process. Additional passive sensor characteristics that enable greater utilization in space applications are small size and weight, long read ranges with low interrogator power, ruggedness, and operability in extreme environments (vacuum, extreme high/low temperature, high radiation, etc.) In this paper, we consider one very promising passive sensor technology, called surface acoustic wave (SAW) radio-frequency identification (RFID), that satisfies all of these criteria. Although SAW RFID tags have great potential for use in numerous space-based remote sensing applications, the limited collision resolution capability of current generation tags limits the performance in a cluttered sensing environment. That is, as more SAW-based sensors are added to the environment, numerous tag responses are superimposed at the receiver and decoding all or even a subset of the telemetry becomes increasingly difficult. Background clutter generated by reflectors other than the sensors themselves is also a problem, as is multipath interference and signal distortion, but the limiting factor in many remote sensing applications can be expected to be tag mutual interference. This problem may be greatly mitigated by proper design of the SAW tag waveform, but that remains an open research problem, and in the meantime, several other related questions remain to be answered including: What are the fundamental relationships between tag parameters such as bit-rate, time-bandwidth-product, SNR, and achievable collision resolution? What are the differences in optimal or near-optimal interrogator designs between noise-limited environments and interference-limited environments? What are the performance characteristics of different interrogator designs in term of parameters such as transmitter power level, range, and number of interfering tags? In this paper, we present the results of a research effort aimed at providing at least partial answers to all of these questions.

Barton, Richard J.↗
DOE OSTI2025

Employing Eye Trackers to Reduce Nuisance Alarms

When process operators anticipate an alarm prior to its annunciation, that alarm loses information value and becomes a nuisance. This study investigated using eye trackers to measure and adjust the salience of alarms with three methods of gaze-based acknowledgement (GBA) of alarms that estimate operator anticipation. When these methods detected possible alarm anticipation, the alarm’s audio and visual salience was reduced. A total of 24 engineering students (male = 14, female = 10) aged between 18 and 45 were recruited to predict alarms and control a process parameter in three scenario types (parameter near threshold, trending, or fluctuating). The study evaluated whether behaviors of the monitored parameter affected how frequently the three GBA methods were utilized and whether reducing alarm salience improved control task performance. The results did not show significant task improvement with any GBA methods (F(3,69) = 1.357, p = 0.263, partial η 2 = 0.056). However, the scenario type affected which GBA method was more utilized (X 2 (2, N = 432) = 30.147, p < 0.001). Alarm prediction hits with gaze-based acknowledgements coincided more frequently than alarm prediction hits without gaze-based acknowledgements (X 2 (1, N = 432) = 23.802, p < 0.001, OR = 3.877, 95% CI 2.25–6.68, p < 0.05). Participant ratings indicated an overall preference for the three GBA methods over a standard alarm design (F(3,63) = 3.745, p = 0.015, partial η 2 = 0.151). This study provides empirical evidence for the potential of eye tracking in alarm management but highlights the need for additional research to increase validity for inferring alarm anticipation.

99 - GENERAL AND MISCELLANEOUS↗
NASA NTRS2026

Investigation of Bipropellant Plume-Induced Contamination Effects on Coverglass Materials

Contamination and degradation of external spacecraft materials by unburned and partially combusted species from bipropellant thruster plumes has long been observed as a key component of the induced space environment. Space shuttle flight experiments and returned flight hardware from the International Space Station (ISS) have both experienced microscopic impact features induced by high-velocity thruster plume droplets. Analytical results have shown that droplet impingement angle relative to a receiving surface plays a key role in the surface damage. Although impacts with normal impingement angles contribute more severely to surface degradation than highly oblique angles, surface effects at higher impingement angles should not be dismissed. Thruster plume-induced materials degradation is a complex phenomenon that depends on a variety of parameters, including but not limited to material type, system temperature and pressure, plume composition, and thruster firing specifications such as number of pulses, pulse duration, and sample distance from the thruster. For space applications, attaining the vacuum pressure and temperature conditions necessary for flight-like plume expansion and exposure conditions is not a trivial task. The German Aerospace Center (Deutsches Zentrum für Luft- und Raumfahrt, DLR) is a facility uniquely capable of simulating such conditions. Test coupons were exposed to bipropellant thruster firings under high vacuum at the DLR facility. Percent area coverage (PAC) and droplet size distributions were evaluated for the uncoated and coated solar array coverglass materials over a range of impingement angles (0̊ to 75̊). A post-test imaging workflow was developed that aimed to quantify changes in sample surface morphology obtained from scanning electron microscopy (SEM) images using the Image Processing and Analysis in Java (ImageJ) tool; an opensource image processing software. The goal was to create a framework through which to evaluate the effect of bipropellant-induced PAC and droplet size distribution on solar array coverglass optical transmission losses. Understanding this relationship is important because optical transmission losses are known to lead to current reduction in solar power generation systems. In addition to the development of surface characterization workflows, valuable lessons learned as they pertain to future investigations and experiments will be discussed. The authors hope that sharing these lessons will facilitate more utilization of DLR’s unique capabilities as well as open the conversation for how best to address experimental characterization of flight-like plume expansion and its impacts on materials surface degradation effects.

Gateway↗
DOE OSTI2020

Using the Carbon Capture Simulation Initiative (CCSI) Tool to Design the Experiments in the Parametric Campaign of a Novel Compact Absorber for Carbon Capture

Gas absorption towers with structured packing and solvent have been used for Carbon Dioxide (CO 2 ) Capture for about many decades. To overcome process limitations and practical disadvantages for CO 2 capture from the stationary emitter (e.g. NG and coal power plant), many new designs have been proposed and explored at the various scales in the last decade with aim of either low energy penalty or low capital cost. To reduce the size of the absorption tower and hence the total cost of CO 2 capture, the University of Kentucky Center for Applied Energy Research Center (UK CAER) has designed and built a novel CO 2 capture absorption tower or Compact Absorber, integrated into an existing large-bench scale CO 2 capture unit. The Compact Absorber has three sections. The top of the column is a fogging section where the solvent is sprayed through a nozzle producing droplets flowing downward in a co-current fashion with the flue gas. The center of the column is a frothing section where the solvent and flue gas flow through regenerative frothing screens designed by Industrial Climate Solutions, Inc. The bottom of the column is a typical structured packing section were the flue gas and solvent flow in a counter-current fashion. The parametric campaign will be conducted in order to optimize the operating parameters for CO 2 capture including liquid/gas ratio, lean loading, and temperature, liquid residence time. A simulated flue gas with 14% CO 2 will be used along with a UK CAER developed proprietary solvent. The 100-hour parametric campaign is designed using a statistical approach of the Sequential Design of Experiments (sDOE). sDOE is one of the CCSI tools that provides an adaptive statistical approach for designing future experiments based on the results of previous experiments. Application of a typical DOE provides the user with the minimum number of experiments required to get the same data, but sDOE allows the user to make an informed choice of experiments based on the results of previous experiments. The complete absorption column has been constructed and has been partially commissioned. Initial data has been collected by operating using the fogging section and the frothing section. The fogging section produces solvent droplets of about 100 μm sauter mean diameter and as small as 25 μm using a hydraulic nozzle by BETE. The frothing section produces bubbles of about 5mm with high mixing of solvent promoting the higher mass transfer from gas to liquid. The absorber reaches the capture efficiency of about 50% with only two sections in operation. Based on the current results, it can be deduced that increasing the solvent feed temperature and including the packed section for absorption the capture efficiency will increase further. Initial data will be collected using all three sections of the absorber and will be used for sDOE. Non-Uniform Space Filling model of sDOE will be used to prioritize the input conditions resulting into maximum capture efficiency. sDOE is performed using the platform called Framework Optimization, Quantification of Uncertainty, and Surrogates (FOQUS). The method and results demonstrating the progress of the parametric campaign from the initial set of experiments to the final stage of obtaining optimized parameters using sDOE tool will be presented in detail.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRS2009

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗
NASA NTRS1990

Impact fracture experiments simulating interstellar grain-grain collisions

Oxide and silicate grains condensing during the early phases of the formation of the solar system or in the outflow of stars are exposed to high partial pressures of the low-z elements H, C, N and O and their simple gaseous compounds. Though refractory minerals are nominally anhydrous and non-carbonate, if they crystallize in the presence of H2O, N2 and CO or CO2 gases, they dissolve traces of the gaseous components. The question arises: How does the presence of dissolved gases or gas components manifest itself when grain-grain collisions occur. What are the gases emitted when grains are shattered during a collision event. Researchers report on fracture experiments in ultrahigh vacuum (UHV, approximately less than 10 to the -8th power mbar) designed to measure (by means of a quadrupole mass spectrometer, QMS, with microns to ms time resolution) the emission of gases and vapors during and after impact (up to 1.5 sec). Two terrestrial materials were chosen which represent structural and compositional extremes: olivine (San Carlos, AZ), a densely packed Mg-Fe(2+) silicate from the upper mantle, available as 6 to 12 mm single crystals, and obsidian (Oregon), a structurally open, alkaline-SiO2-rich volcanic glass. In the olivine crystals OH- groups have been identified spectroscopically, as well as H2 molecules. Obsidian is a water-rich glass containing OH- besides H2O molecules. Olivine from the mantle often contains CO2, either as CO2-rich fluid in fluid inclusions or structurally dissolved or both. By analogy to synthetic glasses CO2 in the obsidian may be present in form of CO2 molecules in voids of molecular dimensions, or as carbonate anions, CO3(2-). No organic molecules have been detected spectroscopically in either material. Results indicate that refractory oxide/silicates which contain dissolved traces of the H2O and CO/CO2 components but no spectroscopically detectable traces of organics may release complex H-C-O (possibly H-C-N-O) molecules upon fracture, plus metal vapor. This points: (1) at complex reaction mechanisms between dissolved H2O, CO/CO2 (and N2) components within the mineral structure or during fracture, and (2) at the possibility that similar emission processes occur following grain-grain collisions in interstellar dust clouds.

Freund, Friedemann↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2013

Domain-Specific Languages and Diagram Customization for a Concurrent Engineering Environment

A major open question for advocates of Model-Based Systems Engineering (MBSE) is the question of how system and subsystem engineers will work together. The Systems Modeling Language (SysML), like any language intended for a large audience, is in tension between the desires for simplicity and for expressiveness. In order to be more expressive, many specialized language elements may be introduced, which will unfortunately make a complete understanding of the language a more daunting task. While this may be acceptable for systems modelers, it will increase the challenge of including subsystem engineers in the modeling effort. One possible answer to this situation is the use of Domain-Specific Languages (DSL), which are fully supported by the Unified Modeling Language (UML). SysML is in fact a DSL for systems engineering. The expressive power of a DSL can be enhanced through the use of diagram customization. Various domains have already developed their own schematic vocabularies. Within the space engineering community, two excellent examples are the propulsion and telecommunication subsystems. A return to simple box-and-line diagrams (e.g., the SysML Internal Block Diagram) are in many ways a step backward. In order allow subsystem engineers to contribute directly to the model, it is necessary to make a system modeling tool at least approximate in accessibility to drawing tools like Microsoft PowerPoint and Visio. The challenge is made more extreme in a concurrent engineering environment, where designs must often be drafted in an hour or two. In the case of the Jet Propulsion Laboratory's Team X concurrent design team, a subsystem is specified using a combination of PowerPoint for drawing and Excel for calculation. A pilot has been undertaken in order to meld the drawing portion and the production of master equipment lists (MELs) via a SysML authoring tool, MagicDraw. Team X currently interacts with its customers in a process of sharing presentations. There are several inefficiencies that arise from this situation. The first is that a customer team must wait two weeks to a month (which is 2-4 times the duration of most Team X studies themselves) for a finalized, detailed design description. Another is that this information must be re-entered by hand into the set of engineering artifacts and design tools that the mission concept team uses after a study is complete. Further, there is no persistent connection to Team X or institutionally shared formulation design tools and data after a given study, again reducing the direct reuse of designs created in a Team X study. This paper presents the underpinnings of subsystem DSLs as they were developed for this pilot. This includes specialized semantics for different domains as well as the process by which major categories of objects were derived in support of defining the DSLs. The feedback given to us by the domain experts on usability, along with a pilot study with the partial inclusion of these tools is also discussed.

Cole, Bjorn↗
NASA NTRS2008

Improving the Optical Quality Factor of the WGM Resonator

Resonators usually are characterized with two partially dependent values: finesse (F) and quality factor (Q). The finesse of an empty Fabry-Perot (FP) resonator is defined solely by the quality of its mirrors and is calculated as F=piR(exp 1/2)/(1-R). The maximum up-to-date value of reflectivity R approximately equal to 1 - 1.6 x 10(exp -6) is achieved with dielectric mirrors. An FP resonator made with the mirrors has finesse F=1.9 x 10(exp 6). Further practical increase of the finesse of FP resonators is problematic because of the absorption and the scattering of light in the mirror material through fundamental limit on the reflection losses given by the internal material losses and by thermodynamic density fluctuations on the order of parts in 109. The quality factor of a resonator depends on both its finesse and its geometrical size. A one-dimensional FP resonator has Q=2 F L/lambda, where L is the distance between the mirrors and lambda is the wavelength. It is easy to see that the quality factor of the resonator is unlimited because L is unlimited. F and Q are equally important. In some cases, finesse is technically more valuable than the quality factor. For instance, buildup of the optical power inside the resonator, as well as the Purcell factor, is proportional to finesse. Sometimes, however, the quality factor is more valuable. For example, inverse threshold power of intracavity hyperparametric oscillation is proportional to Q(exp 2) and efficiency of parametric frequency mixing is proportional to Q(exp 3). Therefore, it is important to know both the maximally achievable finesse and quality factor values of a resonator. Whispering gallery mode (WGM) resonators are capable of achieving larger finesse compared to FP resonators. For instance, fused silica resonators with finesse 2.3 x 10(exp 6) and 2.8 x 10(exp 6) have been demonstrated. Crystalline WGM resonators reveal even larger finesse values, F=6.3 x 10(exp 6), because of low attenuation of light in the transparent optical crystals. The larger values of F and Q result in the enhancement of various nonlinear processes. Low-threshold Raman lasing, optomechanical oscillations, frequency doubling, and hyperparametric oscillations based on these resonators have been recently demonstrated. Theory predicts a possibility of nearly 10(exp 14) room-temperature optical Q-factors of optical crystalline WGM resonators, which correspond to finesse levels higher than 10(exp 9). Experiments have shown numbers a thousand times lower than that. The difference occurs due to media imperfections. To substantially reduce the optical losses caused by the imperfections, a specific, multi-step, asymptotic processing of the resonator is implemented. The technique has been initially developed to reduce microwave absorption in dielectric resonators. One step of the process consists of mechanical polishing performed after high temperature annealing. Several steps repeat one after another to lead to significant reduction in optical attenuation and, as a result, to the increase of Q-factor as well as finesse of the resonator which demonstrates a CaF2 WGM resonator with F greater than 10(exp 7) and Q greater than 10(exp 11).

Savchenkov, Anatoliy↗
DOE OSTI2022

Observation of breakdown wave mechanism in avalanche ionization produced atmospheric plasma generated by a picosecond CO 2 laser

Understanding the formation and long-timescale evolution of atmospheric plasmas produced by ultrashort, long-wavelength infrared (LWIR) pulses is an important but partially understood problem. Of particular interest are plasmas produced in air with a peak laser intensity ~10 12 W/cm 2 , the so-called clamping intensity observed in LWIR atmospheric guiding experiments where tunneling and multi-photon ionization operative at near-IR or shorter wavelengths are inoperative. We find that avalanche breakdown on the surface of aerosol (dust) particles can act to seed the breakdown of air observed above the 200 GW/cm 2 threshold when a train of 3 ps 10.6 μm laser pulses separated by 18 ps is used. The breakdown first appears at the best focus but propagates backward toward the focusing optic as the plasma density approaches critical density and makes forward propagation impossible. The velocity of the backward propagating breakdown can be as high as 10 9 cm/s, an order of magnitude greater than measured with ns pulse-produced breakdown, and can be explained rather well by the so-called breakdown wave mechanism. Transverse plasma expansion with a similar velocity is assisted by UV photoionization and is observed as a secondary longitudinal breakdown mechanism in roughly 10% of the shots. When a cm-size, TW power beam is propagated, interception of aerosol particles is guaranteed and several (40 cm -3 ) breakdown sites appear, each initially producing a near-critical density plasma. On a 10 ns–1 μs timescale, shockwaves from each site expand radially and coalesce to produce a large hot gas channel. The radial velocity of the expansion agrees well with the prediction of the blast wave theory developed for ultrafast atmospheric detonations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

The Nuclear Cryogenic Propulsion Stage

Nuclear Thermal Propulsion (NTP) development efforts in the United States have demonstrated the technical viability and performance potential of NTP systems. For example, Project Rover (1955 - 1973) completed 22 high power rocket reactor tests. Peak performances included operating at an average hydrogen exhaust temperature of 2550 K and a peak fuel power density of 5200 MW/m3 (Pewee test), operating at a thrust of 930 kN (Phoebus-2A test), and operating for 62.7 minutes in a single burn (NRX-A6 test). Results from Project Rover indicated that an NTP system with a high thrust-to-weight ratio and a specific impulse greater than 900 s would be feasible. Excellent results were also obtained by the former Soviet Union. Although historical programs had promising results, many factors would affect the development of a 21st century nuclear thermal rocket (NTR). Test facilities built in the US during Project Rover no longer exist. However, advances in analytical techniques, the ability to utilize or adapt existing facilities and infrastructure, and the ability to develop a limited number of new test facilities may enable affordable development, qualification, and utilization of a Nuclear Cryogenic Propulsion Stage (NCPS). Bead-loaded graphite fuel was utilized throughout the Rover/NERVA program, and coated graphite composite fuel (tested in the Nuclear Furnace) and cermet fuel both show potential for even higher performance than that demonstrated in the Rover/NERVA engine tests.. NASA's NCPS project was initiated in October, 2011, with the goal of assessing the affordability and viability of an NCPS. FY 2014 activities are focused on fabrication and test (non-nuclear) of both coated graphite composite fuel elements and cermet fuel elements. Additional activities include developing a pre-conceptual design of the NCPS stage and evaluating affordable strategies for NCPS development, qualification, and utilization. NCPS stage designs are focused on supporting human Mars missions. The NCPS is being designed to readily integrate with the Space Launch System (SLS). A wide range of strategies for enabling affordable NCPS development, qualification, and utilization should be considered. These include multiple test and demonstration strategies (both ground and in-space), multiple potential test sites, and multiple engine designs. Two potential NCPS fuels are currently under consideration - coated graphite composite fuel and tungsten cermet fuel. During 2014 a representative, partial length (approximately 16") coated graphite composite fuel element with prototypic depleted uranium loading is being fabricated at Oak Ridge National Laboratory (ORNL). In addition, a representative, partial length (approximately 16") cermet fuel element with prototypic depleted uranium loading is being fabricated at Marshall Space Flight Center (MSFC). During the development process small samples (approximately 3" length) will be tested in the Compact Fuel Element Environmental Tester (CFEET) at high temperature (approximately 2800 K) in a hydrogen environment to help ensure that basic fuel design and manufacturing process are adequate and have been performed correctly. Once designs and processes have been developed, longer fuel element segments will be fabricated and tested in the Nuclear Thermal Rocket Element Environmental Simulator (NTREE) at high temperature (approximately 2800 K) and in flowing hydrogen.

Houts, Michael G.↗
DOE OSTI2024

The Monolithic Heat Pipe Microreactor Reference Plant Model

This work introduces a reference plant model for a generic monolithic heat-pipe-cooled microreactor. The model will serve as a springboard to develop future evaluation models in the licensing process of similar microreactor designs at the U.S. Nuclear Regulatory Commission. This model has been developed with the Comprehensive Reactor Analysis Bundle and its specifications are based on open literature publications for the eVinci TM design. BlueCRAB is the U.S. Nu- clear Regulatory Commission non-light-water reactor analysis system based on the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications includes tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C. No moderator or burnable poison pins are used in the design. The fuel enrichment is reduced to control excess reactivity in the core. This core design is not optimized and only serves for testing purposes, since the primary objective of this work is to exercise the multiphysics coupling for this type of reactor system. A three dimensional (3D) core heterogeneous Griffin discrete ordinates (SN) transport model allows the precise calculation of the flux distribution and pin powers. Griffin transfers the power density distribution and obtains a temperature distribution to and from BISON. The BISON model com- putes the 3D core temperature distribution and is coupled to 876 Sockeye subapplications running a heat pipe model. This 3D conduction model is coupled to the various heat pipes via heat flux boundary conditions. The model includes a small gap between the heat pipe and the monolith. Convective heat transfer boundaries with either ambient temperature or condenser temperature as heat sinks are imposed at the model boundaries. The 2D Sockeye heat pipe model uses a vapor- only methodology, which provides the needed resolution for transient calculations and allows the determination of various heat pipe limits. This approach is superior to the superconductor model traditionally used in steady-state calculations. BlueCRAB computes steady-state power and temperature distributions that serve as the initial condition for a loss-of-heat-sink transient simulation. The steady-state results show significant peaking due to the position of the control drum, but this is a characteristic of the particular design used, which is not optimized at this stage. The transient results show the reactor power slowly stabilizing towards a 3% power level after the partial loss of secondary heat removal. Several recriticalities are observed due to cooling through the secondary system but the reactor is self-stabilizing and behaves as expected.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2012

Achievable Performance and Effective Interrogator Design for SAW RFID Sensor Tags

For many NASA missions, remote sensing is a critical application that supports activities such as environmental monitoring, planetary science, structural shape and health monitoring, non-destructive evaluation, etc. The utility of the remote sensing devices themselves is greatly increased if they are passive V that is, they do not require any on-board power supply such as batteries V and if they can be identified uniquely during the sensor interrogation process. Additional passive sensor characteristics that enable greater utilization in space applications are small size and weight, long read ranges with low interrogator power, ruggedness, and operability in extreme environments (vacuum, extreme high/low temperature, high radiation, etc.) In this paper, we consider one very promising passive sensor technology, called surface acoustic wave (SAW) radio-frequency identification (RFID), that satisfies all of these criteria. In general, RFID is a method of identifying items using radio waves to interrogate tags encoded with a unique identifier that are affixed to the items of interest. In the case of passive tags, only the interrogator, which transmits power to the tags in the form of radio-frequency electromagnetic radiation, requires access to a power supply. Passive RFID technologies are used today in many applications, including asset tracking and management, security and access control, and remote sensing. To date, most of the development and application in RFID technology has focused on either asset/inventory tracking and control or security and access control because these are the largest commercial application areas. Recently however, there has been growing interest in using passive RFID technology for remote sensing applications, and SAW devices are at the forefront of RFID sensing technology development. Although SAW RFID tags have great potential for use in numerous space-based remote sensing applications, the limited collision resolution capability of current generation tags limits the performance in a cluttered sensing environment. That is, as more SAW-based sensors are added to the environment, numerous tag responses are superimposed at the receiver and decoding all or even a subset of the telemetry becomes increasingly difficult. Background clutter generated by reflectors other than the sensors themselves is also a problem, as is multipath interference and signal distortion, but the limiting factor in many remote sensing applications can be expected to be tag mutual interference. This problem may be greatly mitigated by proper design of the SAW tag waveform, but that remains an open research problem, and in the meantime, several other related questions remain to be answered including: (1) What are the fundamental relationships between tag parameters such as bit-rate, time-bandwidth-product, SNR, and achievable collision resolution? (2) What are the differences in optimal or near-optimal interrogator designs between noise-limited environments and interference-limited environments? (3) What are the performance characteristics of different interrogator designs in term of parameters such as transmitter power level, range, and number of interfering tags? In this paper, we will present the results of a research effort aimed at providing at least partial answers to all of these questions.

Barton Richard J.↗