Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SPHERICAL CONFIGURATION”

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.

441 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

On the character of turbulent-like flows in self-consistent models of core-collapse supernovae

Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities, especially development of convection, will determine the morphology of the explosion. We have performed 3D simulations using spherical-polar coordinates covering a reduced angular extent (90° computational domain), and with angular resolutions of 2°, 1°, 1/2°, and 1/4°, to study the development of turbulent-like flows in core-collapse supernova explosions on a time scale of order 100 ms. We have employed the multi-physics Chimera code that includes detailed nuclear physics and spectral neutrino transport. Coarse resolution models do not develop an inertial range, presumably due to the bottleneck effect, such that the energy is prevented from cascading down to small scales and tends to accumulate at large scales. High-resolution models instead, start to recover the k −5/3 scaling of Kolmogorov's theory. Stochasticity and few simulation samples limit our ability to predict the development of explosions. Over the simulated time period, our models show no clear trend in improving (or diminishing) conditions for explosion as the angular resolution is increased. However, we find that disordered flow provides an effective pressure, as characterized by the Reynolds stress, behind the shock. This contribution, ~40%–50% of the thermal pressure, aids shock revival and thus the development of the explosion. Finally, we show that the kinetic energy power spectrum of reduced angular extent and full 4π models are consistent, thus indicating that a 90° computational domain is an adequate configuration to study the character of turbulence in CCSNe.

Casanova Bustamante, Jordi↗
NASA NTRS2015

Optimization of a Wcl6 CVD System to Coat UO2 Powder with Tungsten

In order to achieve deep space exploration via Nuclear Thermal Propulsion (NTP), Marshall Space Flight Center (MSFC) is developing W-UO2 CERMET fuel elements, with focus on fabrication, testing, and process optimization. A risk of fuel loss is present due to the CTE mismatch between tungsten and UO2 in the W-60vol%UO2 fuel element, leading to high thermal stresses. This fuel loss can be reduced by coating the spherical UO2 particles with tungsten via H2/WCl6 reduction in a fluidized bed CVD system. Since the latest incarnation of the inverted reactor was completed, various minor modifications to the system design were completed, including an inverted frit sublimer. In order to optimize the parameters to achieve the desired tungsten coating thickness, a number of trials using surrogate HfO2 powder were performed. The furnace temperature was varied between 930 C and 1000degC, and the sublimer temperature was varied between 140 C and 200 C. Each trial lasted 73-82 minutes, with one lasting 205 minutes. A total of 13 trials were performed over the course of three months, two of which were re-coatings of previous trials. The powder samples were weighed before and after coating to roughly determine mass gain, and Scanning Electron Microscope (SEM) data was also obtained. Initial mass results indicated that the rate of layer deposition was lower than desired in all of the trials. SEM confirmed that while a uniform coating was obtained, the average coating thickness was 9.1% of the goal. The two re-coating trials did increase the thickness of the tungsten layer, but only to an average 14.3% of the goal. Therefore, the number of CVD runs required to fully coat one batch of material with the current configuration is not feasible for high production rates. Therefore, the system will be modified to operate with a negative pressure environment. This will allow for better gas mixing and more efficient heating of the substrate material, yielding greater tungsten coating per trial.

Belancik, Grace A.↗
NASA NTRS2013

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗
NASA NTRS2004

Surface Figure Measurement at 20K: Silicon Carbide and Cesic Mirrors

This report presents the facility, methods, and results of testing cryogenically-cooled spherical mirrors, using standard phase-shifting interferometry, at the Goddard Space Flight Center (GSFC). Two mirrors were supplied to GSFC by the European Space Technology Center, so that GSFC could render a second, independent cryo-measurement of their surface figures at 20K. These mirrors, produced by Galileo Avionica and its partners, demonstrate the technology of silicon carbide and SiC-composite lightweighted mirrors, designed for high accuracy (10 nm rms surface figure error) at both room temperature and 20K. The two mirrors provided for testing at GSFC include one made of sintered silicon carbide (mirror blank by Bettini), and a second made of the C/Si/SiC composite known as Cesic (supplied by ECM). Both mirrors are f/2 spheres with 150 mm clear aperture, and with integral mounts. At GSFC, the mirrors will be measured first at at room temperature, making use of standard techniques, with a predicted uncertainty of 2 nm. Then each mirror will be tested in the cryostat, down to 20K: first freely supported and, second, mounted to a Cesic plate that models a Cesic optical bench. The uncertainty of the resultant surface figure error is predicted to be about 5 nm rms. Details of the uncertainty budget are discussed in a related poster. GSFC's techniques and results can be usefully compared to the measurements performed on the same two mirrors by and for Galileo Avionica, who used a somewhat different test configuration and alignment approach.

Blake, Peter↗
NASA NTRS2010

Hydrogen Fire in a Storage Vessel

On October 23, 2007, the operations team began a procedure to sample the Liquid Hydrogen (LH2) storage vessels ("tanks"), and associated transfer system. This procedure was being performed to determine the conditions within the system, and if necessary, to purge the system of any excess Gaseous Hydrogen (GH2) in preparation for reactivation of the system. The system had not been used since 2003. The LH2 storage system contains two (2) spherical pressure vessels of 225,000 gallons in volume, with a maximum working pressure (MAWP) of 50 psig. Eight inch transfer piping connects them to the usage point. Operations began with activation of the burnstack for the LH2 storage area. Pneumatic (GN2) systems in the storage area were then activated and checked. Pressurization of storage tank number 1 with gaseous nitrogen (GN2) was initiated, with a target pressure of 10 psig, at which point samples were planned to be taken. At 5 psig, a loud noise was heard in the upper area of tank number 2. Smoke was seen exiting the burnstack and from the insulation on vent lines for both tanks. At this time tank number 1 was vented and the pressurization system was secured. The mishap resulted in physical damage to both storage tanks, as well as to some of the piping for both tanks. Corrective action included repair of the damaged hardware by a qualified contractor. Preventive action included documented organizational policy and procedures for establishing standby and mothball conditions for facilities and equipment, including provisions as detailed in the investigation report recommendations: Recommendation 1: The using organization should define necessary activities in order to place hydrogen systems in long term periods of inactivity. The defined activities should address requirements for rendering inert, isolation (i.e., physical disconnect, double block and bleed, etc.) and periodic monitoring. Recommendation 2: The using organization should develop a process to periodically monitor hazardous systems for proper configuration (i.e., a daily/weekly/monthly check sheet to verify critical purges are active).

Hester, Zena M.↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2024

Global Sensitivity Analysis of Simulated Remote Sensing Polarimetric Observations Over Snow

This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors. The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum. The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance. The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration. The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.

remote sensing↗
DOE OSTI2021

Design and Analysis of High Heat Flux Plasma-Facing Components for NSTX Upgrade

The National Spherical Tokamak Experiment (NSTX) upgrade (NSTX-U) requirements lead to enhanced heat loads on plasma-facing components (PFCs) especially in the divertor regions, where normal heat flux density can reach up to 8 MW/m 2 . For these high heat flux (HHF) regions, the PFCs were redesigned, to use castellations which relieved the thermal stresses associated with high incident heat flux. Another design feature of HHF PFCs is the absence of front face mounting holes which create localized areas of high temperature and thermal stress concentrations. Optimized shaping of the front face of the HHF PFCs eliminates regions where the front face is perpendicular to the heat flux direction during normal helicity operation and spreads the heat load. A unique mechanism was designed to connect PFC tiles to the NSTX-U center stack casing using locking features accessible from the low heat flux regions. Isotropic graphite was selected as the HHF PFC tile material. Initial tile design parameters were accessed using analytical expressions for pulsed heat flux loading. A working prototype of the locking mechanism was created during the initial stages of the design to prove the concept performance. This article presents an overview of the divertor HHF PFC design and the results of the thermal and structural analyses performed using ANSYS software. The results of the analyses cover normal operating conditions and disruptions which impose electromagnetic (EM) loads from eddy and halo currents. The 3-D, transient, nonlinear analyses took into account the temperature-dependent properties of the materials, friction interfaces between the parts, and variable electric properties of the parts and interfaces. The results confirm that the HHF PFC tiles remain within the allowable limits for the loads defined by the NSTX-U Recovery Project. In addition, tolerance stack up analyses were performed to ensure tile performance in the worst possible assembly configuration. Finally, the design process was completed successfully, and the NSTX HHF PFC tiles are currently in the production phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Hall Magnetohydrodynamic Power Generation and Drag Augmentation Using A Coaxial-Electrode Configuration During Hypersonic Entry

For interplanetary missions, a spacecraft must reduce its relative velocity from orbital speeds of multiple kilometers per second to zero in order to safely land on the surface of a planetary body. During this phase known as planetary reentry, the spacecraft undergoes high speed above Mach 5 up to Mach 30. This hyper-sonic flight regime imposes a strong bow shock in front of the vehicle that imparts extreme aerodynamic and thermal effects. In addition, during portions of this regime, a plasma is formed in the post-shock flow-field around the vehicle due to these high temperatures. Interaction between an applied magnetic field with this electrically conductive fluid, or magnetohydrodynamic (MHD) interaction, can be utilized to convert kinetic energy from the flow into storable electrical energy. Onboard MHD energy generation can benefit the spacecraft by increasing the power of active thermal control components and vehicle control systems during reentry. MHD interaction also exerts a Lorentz body force on the plasma flow that can be manipulated to act as an additional drag force through the vehicle body. MHD drag augmentation can benefit the spacecraft with an additional control mechanism that operates without moving parts. This can be utilized at higher altitudes and subsequently reduce convective heat transfer to the vehicle. MHD interaction with plasma flow suggests that both power generation and drag augmentation can be integrated together in one design. Modern numerical studies have only demonstrated the viability and potential of MHD energy generation and drag modulation separately for implementation into future spacecraft designs [1] [2] [3]. Steeves et. al. presented one possible generator design specifically for reentry vehicles that utilized a modular panel with two embedded electromagnets and two extruding electrodes for air plasma to flow between. Fujino et. al. presented results that indicated increased total drag with a permanent magnet embedded in the blunt body of the Orbital Reentry Experiment (OREX) trajectory. Furthermore, experimental investigations of these two MHD applications have been limited and separate due to the difficult technical nature of experimental de-sign and analysis. For MHD energy generation, there have been two designs which used a cylindrical fore-body with two embedded permanent magnets and two outward facing curved electrodes for tangency to an artificially ionized plasma that flowed along the body length [4] [5]. The first design was studied using an air microwave (MW) ionized supersonic plasma wind tunnel and the second design used radio frequency (RF) ionization of Argon, but both physically demonstrated the feasibility of MHD energy generation in reentry plasma conditions. For MHD drag modulation, one study involved a small-scale cylindrical blunt body embedded with permanent magnets immersed in plasma flow in an arcjet tunnel [6], and another utilized a spherical permanent magnet enclosed with a spherical forebody that was evaluated in a shock expansion tunnel [7]. Both experimental studies provided physical evidence of drag augmentation due to MHD interactions. Therefore, because of the multi-faceted effects of MHD along with both applications having been deemed feasible and beneficial during reentry, this indicates that a coupled power generation and drag modulation design could be achievable. The axi-symmetric nature of blunt bodies for reentry suggests that a spherical set of two ring electrodes near the nose along the body axis. Based on first principles, flow along the walls of the spherical body can then produce MHD interactions between the electrodes. As a result, two currents are generated: a Hall current that produces power and a primary inductive current that augments the drag force. The goal of this work is to provide a simulated design of a coaxial-electrode con-figuration for power generation and drag modulation for experimental testing.

E. Leong↗