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.
Reduction in the complexity of gas turbine fuel controls without loss of control accuracy, reliability, or effectiveness as a method for reducing engine costs is discussed. A description and analysis of hydromechanical approach are presented. A computer simulation of the control mechanism is given and performance of a physical model in engine test is reported.
The Hg/Cd/Te solidification in a Bridgman cell was thermally analyzed by Continuum's VAST code. The energy equation is solved in an axisymmetric, quasi steady domain for both the molten and solid alloy regions. Alloys composition is calculated by a simplified one dimensional model to estimate its effect on melt thermal conductivity and, consequently, on the temperature field within the cell. Solidification is assumed to occur at a fixed temperature of 979 K. Simplified boundary conditions are included to model both the radiant and conductive heat exchange between the furnace walls and the alloy. Calculations are performed to show how the steady state isotherms are affected by: (1) the hot and cold furnace temperatures; (2) boundary condition parameters; and (3) the growth rate which affects the calculated alloy's composition.
Three dimensional incompressible vortex rings are used to model downburst flow field, and modifications are made to express the induced velocities analytically. As a result, this three dimensional downburst model can be used in trajectory optimization, as well as in flight simulations. Two such vortex rings are superposed to match the JAWS AB Corridor wind profiles and the DFW Downburst wind profiles. Model parameters are determined through nonlinear least square fit.
As the number of inverter-interfaced distributed generators (IIDGs) increases, the decrease in the inertia of the microgrid affects the stability of microgrid voltage and current. The virtual synchronous generator (VSG) control method simulates the output characteristics of the synchronous machine to improve the inertia and damping of the microgrid. However, low-frequency oscillations can occur in systems with multi-VSGs connected in parallel, primarily when operating in island mode. Therefore, this article develops the small-signal state-space model for multi-inverters based on VSG and analyses the robust stability region (RSR) of the control parameters. Firstly, an extensible small-signal state-space model is established and simplified in the dq axis, which is a necessary preprocessing for the subsequent stability analysis. Secondly, the guardian map method based on bialternate product operation is constructed, and the RSR is established through translation mapping and rotation mapping. Therefore, the solution of RSR is transformed into the problem of identifying the Hurwitz matrix. Finally, the simulations and experiments verify that the calculated region can ensure stable operation during startup and meet the requirement to suppress low-frequency oscillations. In conclusion, the method of solving RSR proposed in this paper provides a basis for selecting controller parameters in practical engineering.
The modeling of porous bleed regions as boundary conditions in computational fluid dynamics (CFD) simulations of supersonic inlet flows has been improved through a scaling of sonic flow coefficient data for 90deg bleed holes. The scaling removed the Mach number as a factor in computing the sonic flow coefficient and allowed the data to be fitted with a quadratic equation, with the only factor being the ratio of the plenum static pressure to the surface static pressure. The implementation of the bleed model into the Wind-US CFD flow solver was simplified by no longer requiring the evaluation of the flow properties at the boundary-layer edge. The quadratic equation can be extrapolated to allow the modeling of small amounts of blowing, which can exist when recirculation of the bleed flow occurs within the bleed region. The improved accuracy of the bleed model was demonstrated through CFD simulations of bleed regions on a flat plate in supersonic flow with and without an impinging oblique shock. The bleed model demonstrated good agreement with experimental data and three-dimensional CFD simulations of bleed holes.
The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.
In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.
Here we present a wave generalization of the classic Schwarzschild method for constructing self-consistent halos—such a halo consists of a suitable superposition of waves instead of particle orbits, chosen to yield a desired mean density profile. As an illustration, the method is applied to spherically symmetric halos. We derive an analytic relation between the particle distribution function and the wave superposition amplitudes and show how it simplifies in the high-energy (WKB) limit. We verify the stability of such constructed halos by numerically evolving the Schrödinger-Poisson system. The algorithm provides an efficient and accurate way to simulate the time-dependent halo substructures from wave interference. We use this method to construct halos with a variety of density profiles, all of which have a core from the ground-state wave function, though the core-halo relation need not be the standard one.
This article proposes a high-level engineering guide to develop an integrated system model for power quality analysis in electrical submersible pump (ESP) applications with variable frequency drives (VFDs). Such analyses are troublesome in the industry due to the complexity of these systems. In this article, simple steps to perform system integration analyses of such arrangements, including torsional analysis, are developed. A simplified VFD-ESP model suitable for coupled electrical and mechanical analysis in steady state is proposed in all their configurations of practical interest. Such a model can be easily implemented in common simulation software, significantly reducing engineering efforts for implementation and analysis. The focus of the proposed model is the prediction of ESP failures that might result from a poor power quality caused by VFDs. Analytical expressions of different types of harmonics in these systems, as well as their accurate locations in the frequency domain, including their interharmonics and common-mode harmonics, are derived for this purpose. Further, the effectiveness of the proposed model is verified through offline and real-time hybrid simulation results. Finally, a comparison between simulation results obtained using the proposed model and measurements collected on a down-scale laboratory prototype is carried out to demonstrate the accuracy of the suggested modeling approach.
Quantifying changes in soil organic carbon (SOC) stocks and other soil properties is essential for understanding how soils will respond to land management practices and global change. Although they are widely used, comparisons of SOC stocks at fixed depth (FD) intervals are subject to errors when changes in bulk density or soil organic matter occur. The equivalent soil mass (ESM) method has been recommended in lieu of FD for assessing changes in SOC stocks in mineral soils, but ESM remains underutilized for SOC stocks and has rarely been used for other soil properties. In this paper, we draw attention to the limitations of the FD method and demonstrate the advantages of the ESM approach. We provide illustrations to show that the FD approach is susceptible to errors not only for quantifying SOC stocks but also for soil mass-based properties such as SOC mass percent, C:N mass ratio, and δ 13 C. Here we describe the ESM approach and show how it mitigates the FD method limitations. Using bulk density change simulations applied to an empirical dataset from bioenergy cropping systems, we show that the ESM method provides consistently lower errors than FD when quantifying changes in SOC stocks and other soil properties. To simplify the use of ESM, we detail how the method can be integrated into sampling schemes, and we provide an example R computer script that can perform ESM calculations on large datasets. We encourage future studies, whether temporal or comparative, to utilize sampling methods that are amenable to the ESM approach. Overall, we agree with previous recommendations that ESM should be the standard method for evaluating SOC stock changes in mineral soils, but we further suggest that ESM may also be preferred for comparisons of other soil properties including mass percentages, elemental mass ratios, and stable isotope composition.
Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in turn provides us with information about the nature of dark energy. Conventional methods generate shear maps from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser–Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations—an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly scaling shear map construction alternative using a scalable Gaussian process prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest neighbors and fits hyperparameters using batched leave-one-out cross-validation. This work is the first step toward a full, scalable mass mapping method. We work in a simplified regime where we validate our method by interpolating and analyzing maps given noisy point-estimate data from all three shear fields, taken from a suite of N -body ray-tracing simulations. We also show that we can perform these operations at the scale of billions of galaxies on high-performance computing platforms.
The direct simulation Monte Carlo (DSMC) method has been modified to deal with the very small degree of ionization that occurs in the early stages of the 'communications blackout' phase of re-entry. The simplified model of the process assumes that the ionization is entirely a consequence of the reaction N + O yields NO(+) + e-. The measurements from the RAM C-II test vehicle at an altitude of 81 km have been chosen as the test case. It was found that the form of the electron distribution along the shock layer was in agreement with the measurements. That normal to the surface was also in agreement with the measurements if the surface is regarded as fully catalytic. However, quantitative agreement for the fully catalytic case requires a reaction rate about an order of magnitude greater than those generally associated with this reaction.
The Urban Air Mobility concept envisions a new future for transportation in urban environments, primarily focused around vertiport operations. A vertiport refers to an identifiable ground or elevated area used for the vertical takeoff and landing of an aircraft. One critical role within vertiport operations that has been proposed but is under-researched is the vertiport manager. The vertiport manager is expected to manage the ground-to-air operations at the vertiport, but their precise responsibilities, tasks, and challenges are currently undefined. The High Density Vertiplex Subproject at NASA is responsible for developing and testing technologies, concepts, and architectures that will support the infrastructure needed for terminal environments around vertiports. Our previous work investigated ground control station operators performing simulated and live flight operations within a remote operations environment. In the current study, we extended from our previous work to explore the role of a vertiport manager. Specifically, we investigated five participants serving as a vertiport manager across both simulated and live flight operations. We performed multiple knowledge elicitation techniques, a thematic content analysis of qualitative data, and naturalistic observations to create a list of insights, design and training recommendations, and a simplified cognitive task diagram of the vertiport manager’s primary task. Findings from this study are discussed within the context of the current debate on the degree to which this role should be automated. A third alternative is suggested to allow a human vertiport manager to be successful by teaming with increasingly autonomous systems, designing the role to improve safety and efficiency of Urban Air Mobility operations.
The Urban Air Mobility concept envisions a new future for transportation in urban environments, primarily focused around vertiport operations. A vertiport refers to an identifiable ground or elevated area used for the vertical takeoff and landing of an aircraft. One critical role within vertiport operations that has been proposed but is under-researched is the vertiport manager. The vertiport manager is expected to manage the ground-to-air operations at the vertiport, but their precise responsibilities, tasks, and challenges are currently undefined. The High Density Vertiplex Subproject at NASA is responsible for developing and testing technologies, concepts, and architectures that will support the infrastructure needed for terminal environments around vertiports. Our previous work investigated ground control station operators performing simulated and live flight operations within a remote operations environment. In the current study, we extended from our previous work to explore the role of a vertiport manager. Specifically, we investigated five participants serving as a vertiport manager across both simulated and live flight operations. We performed multiple knowledge elicitation techniques, a thematic content analysis of qualitative data, and naturalistic observations to create a list of insights, design and training recommendations, and a simplified cognitive task diagram of the vertiport manager’s primary task. Findings from this study are discussed within the context of the current debate on the degree to which this role should be automated. A third alternative is suggested to allow a human vertiport manager to be successful by teaming with increasingly autonomous systems, designing the role to improve safety and efficiency of Urban Air Mobility operations.
The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.
The yield of chromosomal aberrations has been shown to increase in the lymphocytes of astronauts after long-duration missions of several months in space. Chromosome exchanges, especially translocations, are positively correlated with many cancers and are therefore a potential biomarker of cancer risk associated with radiation exposure. Although extensive studies have been carried out on the induction of chromosomal aberrations by low- and high-LET radiation in human lymphocytes, fibroblasts, and epithelial cells exposed in vitro, there is a lack of data on chromosome aberrations induced by low dose-rate chronic exposure and mixed field beams such as those expected in space. Chromosome aberration studies at NSRL will provide the biological validation needed to extend the computational models over a broader range of experimental conditions (more complicated mixed fields leading up to the galactic cosmic rays (GCR) simulator), helping to reduce uncertainties in radiation quality effects and dose-rate dependence in cancer risk models. These models can then be used to answer some of the open questions regarding requirements for a full GCR reference field, including particle type and number, energy, dose rate, and delivery order. In this study, we designed a simplified mixed field beam with a combination of proton, helium, oxygen, and iron ions with shielding or proton, helium, oxygen, and titanium without shielding. Human fibroblasts cells were irradiated with these mixed field beam as well as each single beam with acute and chronic dose rate, and chromosome aberrations (CA) were measured with 3-color fluorescent in situ hybridization (FISH) chromosome painting methods. Frequency and type of CA induced with acute dose rate and chronic dose rates with single and mixed field beam will be discussed. A computational chromosome and radiation-induced DNA damage model, BDSTRACKS (Biological Damage by Stochastic Tracks), was updated to simulate various types of CA induced by acute exposures of the mixed field beams used for the experiments. The chromosomes were simulated by a polymer random walk algorithm with restrictions to their respective domains in the nucleus [1]. The stochastic dose to the nucleus was calculated with the code RITRACKS [2]. Irradiation of a target volume by a mixed field of ions was implemented within RITRACKs, and the fields of ions can be delivered over specific periods of time, allowing the simulation of dose-rate effects. Similarly, particles of various types and energies extracted from a pre-calculated spectra of galactic cosmic rays (GCR) can be used in RITRACKS. The number and spatial location of DSBs (DNA double-strand breaks) were calculated in BDSTRACKS using the simulated chromosomes and local (voxel) dose. Assuming that DSBs led to chromosome breaks, and simulating the rejoining of damaged chromosomes occurring during repair, BDSTRACKS produces the yield of various types of chromosome aberrations as a function of time (only final yields are presented). A comparison between experimental and simulation results will be shown.
A 14-degree-of-freedom (14-DOF) flexible multibody dynamics model is developed and verified for a semi-submersible floating offshore wind turbine (FOWT). The model considers the coupled dynamics of the platform, tower, nacelle, blades, and mooring subjected to external wind and wave loads. The platform is simplified as a rigid body connected to the seabed by the mooring system. The tower and blade are modeled as flexible cantilever beams. The energy method is used to derive the governing equations of motion, where the kinetic energy, potential energy, and work done by external wind-wave forces are all deduced in a global coordinate system. The 5-MW baseline semi-submersible FOWT is used to verify the derived model against the results simulated from FAST developed by the US National Renewable Energy Laboratory (NREL) at two scenarios: free decay state and different wind-wave load cases. Results show that the established FOWT model can well reflect the vibration characteristics of FAST model. The application of the model to the control of platform pitch with a tuned mass damper is studied. Furthermore, the simplified model could provide a low-order method for the structural dynamics analysis and advanced vibration control design for the multi-body components of the semi-submersible FOWT in the future.
Understanding air propagation and condensation following a catastrophic vacuum break in particle accelerator beamlines cooled by liquid helium is essential for ensuring operational safety. This review summarizes experimental and theoretical work conducted in our cryogenics lab to address this issue. Systematic measurements were performed to study nitrogen gas propagation in uniform copper tubes cooled by both normal liquid helium (He I) and superfluid helium (He II). These experiments revealed a nearly exponential deceleration of the gas front, with stronger deceleration observed in He II-cooled tubes. To interpret these results, a one-dimensional (1D) theoretical model was developed, incorporating gas dynamics, heat transfer, and condensation mechanisms. The model successfully reproduced key experimental observations in the uniform tube system. However, recent experiments involving a bulky copper cavity designed to mimic the geometry of a superconducting radiofrequency (SRF) cavity revealed strong anisotropic flow patterns of nitrogen gas within the cavity, highlighting limitations in extrapolating results from simplified tube geometries to real accelerator beamlines. To address these complexities, we outline plans for systematic studies using tubes with multiple bulky cavities and the development of a two-dimensional (2D) model to simulate gas dynamics in these more intricate configurations. As a result, these efforts aim to provide a comprehensive understanding of vacuum breaks in particle accelerators and improve predictive capabilities for their operational safety.