SEARCH · Engineering Papers
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.
RENEW computer program simulates maintenance events and estimates data pertinent to maintenance demands. Developed in support of Space Station Freedom Program (SSFP) Work Package 4. Uses data on reliability and maintainability (R&M) as well as logistical data to estimate both average and time-dependent maintenance demands. Estimates failure and repair times by use of Monte Carlo simulations. Written in BASIC and Assembly language.
Empirical correlations for evaporation rates from single fuel droplets have existed since the 1930s. These correlations, which will be referred to in this article as Froessling/Ranz-Marshall types of correlations, are appropriate to the special cases of steady-state evaporation in the absence of chemical reaction. In a previous article by the authors, the quasi-steady evaporation and burning processes associated with a fuel drop in a convective environment are examined through a droplet model based on the boundary layer approach. For droplet Reynolds numbers of practical interest, this model produces very reasonable steady state as well as quasi-time-dependent droplet simulations, requiring relatively short computational times and yielding good agreement with the above-mentioned empirical correlations. The steady-state case, however, is usually relevant to practical combustor situations only when the drop has reached a nearly uniform temperature since the heating process of the drop cannot be considered to be quasi-steady. In the present study, the transient heating process of the droplet interior during evaporation and/or burning is taken into account, and thus calculations pertaining to the entire life-time of the droplet are carried out. It is of particular interest here to obtain simplified correlations to describe the transient behavior of evaporating and burning droplets; these may be incorporated with greater ease into spray calculations. Accordingly, we have chosen to use stagnation conditions in the present model in a modification of the Froessling/Ranz-Marshall correlations. These modified correlations, incorporating an effective transfer number, produce a fairly accurate representation of droplet evaporation and burning, while requiring only one tenth the computational effort used in a full boundary layer solution.
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.
The Two-dimensional Implicit Thermal Response and Ablation program, TITAN, was developed and integrated with a Navier-Stokes solver, GIANTS, for multidimensional ablation and shape change simulation of thermal protection systems in hypersonic flow environments. The governing equations in both codes are demoralized using the same finite-volume approximation with a general body-fitted coordinate system. Time-dependent solutions are achieved by an implicit time marching technique using Gauess-Siedel line relaxation with alternating sweeps. As the first part of a code validation study, this paper compares TITAN-GIANTS predictions with thermal response and recession data obtained from arc-jet tests recently conducted in the Interaction Heating Facility (IHF) at NASA Ames Research Center. The test models are graphite sphere-cones. Graphite was selected as a test material to minimize the uncertainties from material properties. Recession and thermal response data were obtained from two separate arc-jet test series. The first series was at a heat flux where graphite ablation is mainly due to sublimation, and the second series was at a relatively low heat flux where recession is the result of diffusion-controlled oxidation. Ablation and thermal response solutions for both sets of conditions, as calculated by TITAN-GIANTS, are presented and discussed in detail. Predicted shape change and temperature histories generally agree well with the data obtained from the arc-jet tests.
The United States Department of Energy is planning to develop the Advanced Stirling Radioisotope Generator (ASRG) for the National Aeronautics and Space Administration (NASA) for potential use on future space missions. The ASRG provides substantial efficiency and specific power improvements over radioisotope power systems of heritage designs. The ASRG would use General Purpose Heat Source modules as energy sources and the free-piston Advanced Stirling Convertor (ASC) to convert heat into electrical energy. Lockheed Martin Corporation of Valley Forge, Pennsylvania, is integrating the ASRG systems, and Sunpower, Inc., of Athens, Ohio, is designing and building the ASC. NASA Glenn Research Center of Cleveland, Ohio, manages the Sunpower contract and provides technology development in several areas for the ASC. One area is reliability assessment for the ASC heater head, a critical pressure vessel within which heat is converted into mechanical oscillation of a displacer piston. For high system efficiency, the ASC heater head operates at very high temperature (850 C) and therefore is fabricated from an advanced heat-resistant nickel-based superalloy Microcast MarM-247. Since use of MarM-247 in a thin-walled pressure vessel is atypical, much effort is required to assure that the system will operate reliably for its design life of 17 years. One life-limiting structural response for this application is creep; creep deformation is the accumulation of time-dependent inelastic strain under sustained loading over time. If allowed to progress, the deformation eventually results in creep rupture. Since creep material properties are not available in the open literature, a detailed creep life assessment of the ASC heater head effort is underway. This paper presents an overview of that creep life assessment approach, including the reliability-based creep criteria developed from coupon testing, and the associated heater head deterministic and probabilistic analyses. The approach also includes direct benchmark experimental creep assessment. This element provides high-fidelity creep testing of prototypical heater head test articles to investigate the relevant material issues and multiaxial stress state. Benchmark testing provides required data to evaluate the complex life assessment methodology and to validate that analysis. Results from current benchmark heater head tests and newly developed experimental methods are presented. In the concluding remarks, the test results are shown to compare favorably with the creep strain predictions and are the first experimental evidence for a robust ASC heater head creep life.
A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.
We present a description of the theoretical framework and best practice for using the paleo-climate model component of the Coupled Model Intercomparison Project (Phase 5) (CMIP5) to constrain future projections of climate using the same models. The constraints arise from measures of skill in hindcasting paleo-climate changes from the present over 3 periods: the Last Glacial Maximum (LGM) (21 thousand years before present, ka), the mid-Holocene (MH) (6 ka) and the Last Millennium (LM) (8501850 CE). The skill measures may be used to validate robust patterns of climate change across scenarios or to distinguish between models that have differing outcomes in future scenarios. We find that the multi-model ensemble of paleo-simulations is adequate for addressing at least some of these issues. For example, selected benchmarks for the LGM and MH are correlated to the rank of future projections of precipitationtemperature or sea ice extent to indicate that models that produce the best agreement with paleoclimate information give demonstrably different future results than the rest of the models. We also find that some comparisons, for instance associated with model variability, are strongly dependent on uncertain forcing timeseries, or show time dependent behaviour, making direct inferences for the future problematic. Overall, we demonstrate that there is a strong potential for the paleo-climate simulations to help inform the future projections and urge all the modeling groups to complete this subset of the CMIP5 runs.
Neutron noise techniques constitute several analysis methods applicable to non-destructive assay. One technique is the Rossi-α method to calculate the prompt neutron decay constant (α) or its inverse, the prompt neutron period (1/α), for assemblies of fissionable material. This work evaluates a high-data-throughput measurement at the National Criticality Experiments Research Center (NCERC) with organic scintillators measuring a subcritical assembly of highly enriched uranium (HEU) metal (93% 235 U). The assembly comprises hemi shells stacked together to form fully closed shells and is driven by a Cf source at the center. The assembly is a total of 59.85 kg HEU and a k eff of 0.99, according to MCNP6.2 KCODE simulations. Measurements were acquired with a three-by-four array of 5.08-cm-diameter by 5.08-cm-length trans-stilbene crystals 166 cm from the assembly center. Two types of coincident binning methods are used to build the Rossi-α distribution of coincident neutron detections: A) type 1 binning, also known as any-and-all forward time differences, and B) type 1 binning with a digital-delay technique that is analogous to Orndoff’s use of delay cabling in early Rossi-alpha measurements. Method B) disregards same-detector coincidences and once all nearest time coincidences are collected between all detectors after a trigger, a time delay of 0.75 µs is implemented. The measured prompt neutron decay constants for both techniques are calculated from single exponential fits and the two methods are compared. Method A) shows an apparent timing discontinuity near 500 ns time differences. This timing discontinuity interferes with the fitting method used to calculate the prompt neutron decay constant. Method B) mitigates the artificial timing discontinuity and removes the disagreement of the fit. Future work will model the time-dependent detector response to discern why this timing discontinuity occurs, discern how Method B) reduces the amplitude of this timing discontinuity, and apply this method to critical assembly measurements to ultimately confirm the recommendation to use Method B) for high-data-throughput measurements.
Understanding the emergent electronic structure in twisted atomically thin layers has led to the exciting field of twistronics. However, practical applications of such systems are challenging since the specific angular correlations between the layers must be precisely controlled and the layers have to be single crystalline with uniform atomic ordering. Here, in this work, an alternative, simple, and scalable approach is suggested, where nanocrystalline two-dimensional (2D) film on 3D substrates yields twisted-interface-dependent properties. Ultrawide-bandgap hexagonal boron nitride (h-BN) thin films are directly grown on high in-plane lattice mismatched wide-bandgap silicon carbide (4H-SiC) substrates to explore the twist-dependent structure-property correlations. Concurrently, nanocrystalline h-BN thin film shows strong non-linear second-harmonic generation and ultra-low cross-plane thermal conductivity at room temperature, which are attributed to the twisted domain edges between van der Waals stacked nanocrystals with random in-plane orientations. First-principles calculations based on time-dependent density functional theory manifest strong even-order optical nonlinearity in twisted h-BN layers. This work unveils that directly deposited 2D nanocrystalline thin film on 3D substrates could provide easily accessible twist-interfaces, therefore enabling a simple and scalable approach to utilize the 2D-twistronics integrated in 3D material devices for next-generation nanotechnology.
We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.
Here, we introduce high-order harmonic sideband spectroscopy and show that it can be a robust probe of charge migration (CM) in a halogenated carbon-chain molecule. We simulate both the CM and high-order harmonic generation (HHG) dynamics using ab initio time-dependent density-functional theory. We find that CM dynamics initiated along the molecular backbone induces sidebands in the HHG spectrum driven by a delayed laser pulse that is polarized perpendicular to the molecular axis. Monitoring the spectrum as either the HHG laser frequency or the relative delay is scanned allows for the extraction of detailed information about the time-domain characteristics of the CM process.
Abstract We present a new implementation of real‐time time‐dependent density functional theory (RT‐TDDFT) for calculating excited‐state dynamics of periodic systems in the open‐source Python‐based PySCF software package. Our implementation uses Gaussian basis functions in a velocity gauge formalism and can be applied to periodic surfaces, condensed‐phase, and molecular systems. As representative benchmark applications, we present optical absorption calculations of various molecular and bulk systems and a real‐time simulation of field‐induced dynamics of a (ZnO) 4 molecular cluster on a periodic graphene sheet. We present representative calculations on optical response of solids to infinitesimal external fields as well as real‐time charge‐transfer dynamics induced by strong pulsed laser fields. Due to the widespread use of the Python language, our RT‐TDDFT implementation can be easily modified and provides a new capability in the PySCF code for real‐time excited‐state calculations of chemical and material systems.
The decay of the pump-beam spatial structure under conditions of stimulated temperature scattering of a pulsed laser beam in toluene as a result of two-photon absorption of pump radiation and time-dependent interaction of a laser pulse with a medium is considered. (paper)
Fluid amplifier dynamic characteristics and Fortran program for numerical solution of time- dependent two-dimensional Navier-Stokes equation for viscous jet in arbitrary flow field
Studies and experiments to determine time dependency of plasma density, high frequency electron-phonon interactions in magnetic fields, and operation of stuffed-cusp plasma facility
The effective properties of composites are influenced by the time-dependent behavior of polymer matrices very sensitive to changes in temperature. Improved plate and beam models are required to efficiently design, and simulate composite structures when the long-term performance of large anisotropic composite structures is the matter of interest. In this work, mechanics of structure genome (MSG) is used to con-struct linear thermoviscoelastic plate and beam models that can homogenize three-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integral formulation for thermorheologically simple materials subject to finite temperature changes with the restriction that the strain is small. The reduced time parameter is introduced to relate the time-temperature dependency of the anisotropic material by means of master curves at reference conditions. The new formulation has been implemented in SwiftCompTM, a general-purpose multiscale constitutive modeling code based on MSG. Experimental data and three-dimensional direct numerical simulations of thin-ply high-strain composites (TP-HSC) using a commercial finite element analysis (FEA) package are conducted to verify the accuracy of SwiftCompTM results. The paper also analyzes the relationship between the shift factor of the polymer matrix and the temperature dependencies of the effective beam properties.
Evaluation of the electron-nuclear dynamics and relaxation mechanisms of gold and silver nanoclusters and their alloys is important for future photocatalytic, light harvesting, and photoluminescence applications of these systems. In this work, the effect of silver doping on the nonradiative excited state relaxation dynamics of the atomically precise thiolate-protected gold nanocluster [Au 25- n Ag n (SH) 18 ] –1 ( n = 1, 12, 25) is studied theoretically. Time-dependent density functional theory is used to study excited states lying in the energy range 0.0–2.5 eV. The fewest switches surface hopping method with decoherence correction was used to investigate the dynamics of these states. The HOMO–LUMO gap increases significantly upon doping of 12 silver atoms but decreases for the pure silver nanocluster. Doped clusters show a different response for ground state population increase lifetimes and excited state population decay times in comparison to the undoped system. The ground state recovery times of the S 1 –S 6 states in the first excited peak were found to be longer for [Au 13 Ag 12 (SH) 18 ] –1 than the corresponding recovery times of other studied nanoclusters, suggesting that this partially doped nanocluster is best for preserving electrons in an excited state. The decay time constants were in the range of 2.0–20 ps for the six lowest energy excited states. Among the higher excited states, S 7 has the slowest decay time constant although it occurs more quickly than S 1 decay. Overall, these clusters follow common decay time constant trends and relaxation mechanisms due to the similarities in their electronic structures.
Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.