SEARCH · Engineering Papers

Results for “Numerical approximation & analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Post-Test Inspection of NASA's Evolutionary Xenon Thruster Long-Duration Test Hardware: Discharge Chamber

The NEXT Long-Duration Test is part of a comprehensive thruster service life assessment intended to demonstrate overall throughput capability, validate service life models, quantify wear rates as a function of time and operating condition, and identify any unknown life-limiting mechanisms. The test was voluntarily terminated in February 2014 after demonstrating 51,184 hours of high-voltage operation, 918 kg of propellant throughput, and 35.5 MN-s of total impulse. The post-test inspection of the thruster hardware began shortly afterwards with a combination of non-destructive and destructive analysis techniques, and is presently nearing completion. This paper presents relevant results of the post-test inspection for the discharge chamber as well as other miscellaneous components such as the high-voltage propellant isolators and electrical cabling. Comparison of magnetic field measurements taken during pretest and post-test inspections indicate that the field strength did not degrade, consistent with performance data obtained during the test. Inspection of discharge chamber mesh samples show a deposition coating primarily composed of grid material that is approximately 15 micrometers in thickness. This thickness is well within the retention capability of the mesh and is therefore not expected to present any issues. Approximately 3.1 grams of deposition flakes were found at the bottom of the discharge chamber, composed primarily of grid material and carbon. Calculated size histograms of these flakes indicate that 99% have a maximum dimension of 200 micrometers or smaller, which is significantly less than the ion optics grid gap. Larger flakes that are capable of causing a grid-to-grid short will be analyzed to determine if their formation will occur in flight or is a facility effect. The high-voltage propellant isolators as well as numerous other electrical insulators were inspected and no evidence of arcing or any other issues were found.

Electric Propulsion↗

NASA NTRS2021

Post-test Inspection of NASA’s Evolutionary Xenon Thruster Long-Duration Test Hardware: Discharge Chamber

The NEXT Long-Duration Test is part of a comprehensive thruster service life assessment intended to demonstrate overall throughput capability, validate service life models, quantify wear rates as a function of time and operating condition, and identify any unknown life-limiting mechanisms. The test was voluntarily terminated in February 2014after demonstrating 51,184 hours of high-voltage operation, 918 kg of propellant throughput, and 35.5 MN-s of total impulse. The post-test inspection of the thruster hardware began shortly afterwards with a combination of non-destructive and destructive analysis techniques, and is presently nearing completion. This paper presents relevant results of the post-test inspection for the discharge chamber as well as other miscellaneous components such as the high-voltage propellant isolators and electrical cabling. Comparison of magnetic field measurements taken during pretest and post-test inspections indicate that the field strength did not degrade, consistent with performance data obtained during the test. Inspection of discharge chamber mesh samples show a deposition coating primarily composed of grid material that is approximately 15 μm in thickness. This thickness is well within the retention capability of the mesh and is therefore not expected to present any issues. Approximately 3.1 grams of deposition flakes were found at the bottom of the discharge chamber, composed primarily of grid material and carbon. Calculated size histograms of these flakes indicate that 99% have a maximum dimension of 200 μm or smaller, which is significantly less than the ion optics grid gap. Larger flakes that are capable of causing a grid-to-grid short will be analyzed to determine if their formation will occur in flight or is a facility effect. The high-voltage propellant isolators as well as numerous other electrical insulators were inspected and no evidence of arcing or any other issues were found.

Rohit Shastry↗

DOE OSTI2024

Evaluation of the Radioactive Material Release in the Harborview Research and Training Building and Implications for Emergency Response

On May 2, 2019, during the 137 Cs source recovery operation, a source capsule in a research irradiator containing approximately 77.1 TBq was breached. Based on a geometric reconstruction analysis of the damage to the capsule, approximately 46.3 GBq (0.04%) was impacted by the chop saw (grinder) inside a Mobile Hot Cell (MHC) on the loading dock at the University of Washington Harborview Research and Training (HRT) Building. A very small fraction of the material impacted, less than 1%) was released from the Mobile Hot Cell and then to the rest of the HRT Building. The objectives of this project were to assess the accidental release of 137 CsCl and its implications related to emergency response methods and the ramifications of 137 CsCl transport. The phenomenology of this event was also compared with past alkali halide dispersal events. The vast number of measurements and samples collected by the remediation contractors, the Department of Energy's Nuclear Emergency Support Team, and the small number of retrospective samples collected by the authors informed the analysis. The techniques included (1) autoradiography and electron microscopy of samples collected from the HRT Building and the irradiator, (2) 3D visualization of deposition on surfaces and within the ventilation system, and (3) a study of the damage to the source capsule to evaluate the Cs particle size and particle composition due to the grinding accident. Subsequently, the cesium contaminant transport through the numerous pathways in the building was reconstructed to assess the deposition on surfaces as a function of particle size. The implications for emergency response are relevant to data quality and management. A Data Quality Objective (DQO) guides data collection methods so that they have appropriate accuracy and precision for the intended application. Recommendations were made with respect to the sample collection protocols and sample archival.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2021

Sensitivity Analysis of Tunable Equation of State Material Model In Pulsed Mercury Target Simulation

A pulsed neutron spallation target is subjected to very short but intense loads from repeated proton pulses. Approximately 60% of the energy from each proton pulse is deposited into the mercury target material and the stainless-steel target structure, leading to a high-pressure region in both the stationary target structure and the flowing mercury. The high-pressure region propagates and leads to fluid-structure interaction. The resultant loading on the target structure containing liquid mercury is difficult to predict, although various simulation approaches and material models for the mercury have been tried. To date, the best match of simulation to experimental data is obtained by using an equation of state (EOS) material model with a specified tensile cutoff pressure, which simulates the cavitation threshold. The inclusion of a threshold to represent cavitation is key to the successful predictions of stress waves triggered by the high-energy pulse striking the mercury and vessel. However, recent measurements of target structure strain show that significant discrepancies remain between the measured and simulated strain values in the EOS mercury model. These differences grow when noncondensable helium gas is intentionally injected into the flowing mercury to reduce the loading on the structure. An EOS-based proportional–integral–derivative (PID) mercury model has been proposed to reduce the gap between the measured and simulated vessel strain responses for targets with gas injection. The conceptual and numerical description and initial investigation of the PID model are presented in previous work. Further studies of this PID model — including the sensitivity of the structure’s strain response to model parameters (the tensile cutoff, PID parameters Kp, Ki, and Kd) — are reported in this article. Results show the strain response is more sensitive to changes in the tensile cutoff value than to changes in the model parameters Kp, Ki, and Kd. These results will aid in future work where the model parameters will be optimized to match simulation data to strain measurements.

Lin, Lianshan↗

NASA NTRS1993

A study of the large-scale infrared emission from a selected dark cloud

An investigation of the infrared emission energetics and embedded population in the rho Ophiuchi dark cloud is summarized. With a distance of approximately 140 pc, the rho Ophiuchi cloud is one of the closest regions of recent star formation. It is also one of the best studied such regions with numerous observations at all wavelengths. The Infrared Astronomy Satellite (IRAS) data of the cloud provided a new glimpse of the overall structure of the cloud. In particular, the interaction of radiation from the Sco-Oph OB Association on the external heating of the cloud was very evident on Skyflux and Survey CO-Add images produced by IRAS. The infrared survey also revealed a number of new embedded sources in the cloud which have subsequently been observed from the ground. In earlier study, the overall energies of the cloud using the IRAS data was explored. The main conclusions of that work were: (1) the overall luminosity of the cloud is well explained by the emission of the known B-stars, HD 147889, SR-3, and S1, along with a 15 percent contribution from the external radiation field; (2) the dust physical temperatures were significantly lower than the observed CO gas temperatures; and (3) dust grains are heated to only 10 percent to 20 percent of the total depth into the cloud. This analysis was extended by drawing on data from large-scale CO maps of Loren (1989) and from near-infrared surveys of the embedded population.

Young, Erick T.↗

NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗

NASA NTRSDate not supplied

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low earth orbit for scientific, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with 8 Teledyne dosimeters over varying shielding effectiveness. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 –AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat 3 unit (U) models. Using NOVICE model radiation analysis coding, the shielding effectiveness's, based on a generic CubeSat 3U structure with 4 electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cmto4.44 cm (0.550 g/cm2–12.0 g/cm2) thick aluminum. For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID)increases by nearly a magnitude for thin-walled aluminum 0.550g/cm2-0.686 g/cm2(0.204 cm –0.254 cm) typical CubeSat structures. The shielding effectiveness by NOVICE Sigma estimates, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from Earth's magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast, to a thin walled generic 3U CubeSat, Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm2from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the spacecraft’s mission value versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗

DOE OSTI2025

Kolmogorov-Arnold wavefunctions

Here, this work investigates Kolmogorov-Arnold network-based (KAN) wave-function Ansätz as viable representations for quantum Monte Carlo simulations. Through systematic analysis of one-dimensional model systems, we evaluate their computational efficiency and representational power against established methods. Our numerical experiments suggest some efficient training methods and we explore how the computational cost scales with desired precision, particle number, and system parameters. Roughly speaking, KANs seem to be 10 times cheaper computationally than other neural-network-based Ansätz . We also introduce a novel approach for handling strong short-range potentials—a persistent challenge for many numerical techniques—which generalizes efficiently to higher-dimensional, physically relevant systems with short-ranged strong potentials common in atomic and nuclear physics.

1-dimensional systems↗

DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned paper is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low Earth orbit for science, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low Earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with eight Teledyne µDosimeters over varying shielding effectivenesses. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 – AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat three unit (3U) models. Using NOVICE model radiation analysis coding, the shielding effectivenesses, based on a generic CubeSat 3U structure with four electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cm to 4.44 cm (0.550 g/cm 2 to 12.0 g/cm 2 ). For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID) increases by nearly a magnitude for thin-walled aluminum 0.550-g/cm 2 - 0.686-g/cm 2 (0.204-cm – 0.254-cm) typical CubeSat structures. The shielding effectiveness estimated by NOVICE Sigma, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one-year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from the Earth’s magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast to a thin-walled generic 3U CubeSat, the Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm 2 from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the mission value of the spacecraft versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗

NASA NTRS1997

Neural Network Prediction of Failure of Damaged Composite Pressure Vessels from Strain Field Data Acquired by a Computer Vision Method

This effort used a new and novel method of acquiring strains called Sub-pixel Digital Video Image Correlation (SDVIC) on impact damaged Kevlar/epoxy filament wound pressure vessels during a proof test. To predict the burst pressure, the hoop strain field distribution around the impact location from three vessels was used to train a neural network. The network was then tested on additional pressure vessels. Several variations on the network were tried. The best results were obtained using a single hidden layer. SDVIC is a fill-field non-contact computer vision technique which provides in-plane deformation and strain data over a load differential. This method was used to determine hoop and axial displacements, hoop and axial linear strains, the in-plane shear strains and rotations in the regions surrounding impact sites in filament wound pressure vessels (FWPV) during proof loading by internal pressurization. The relationship between these deformation measurement values and the remaining life of the pressure vessels, however, requires a complex theoretical model or numerical simulation. Both of these techniques are time consuming and complicated. Previous results using neural network methods had been successful in predicting the burst pressure for graphite/epoxy pressure vessels based upon acoustic emission (AE) measurements in similar tests. The neural network associates the character of the AE amplitude distribution, which depends upon the extent of impact damage, with the burst pressure. Similarly, higher amounts of impact damage are theorized to cause a higher amount of strain concentration in the damage effected zone at a given pressure and result in lower burst pressures. This relationship suggests that a neural network might be able to find an empirical relationship between the SDVIC strain field data and the burst pressure, analogous to the AE method, with greater speed and simplicity than theoretical or finite element modeling. The process of testing SDVIC neural network analysis and some encouraging preliminary results are presented in this paper. Details are given concerning the processing of SDVIC output data such that it may be used as back propagation neural network (BPNN) input data. The software written to perform this processing and the BPNN algorithm are also discussed. It will be shown that, with limited training, test results indicate an average error in burst pressure prediction of approximately six percent,

Russell, Samuel S.↗

DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2008

Simulation of the Impact of New Ocean Surface Wind Measurements on H*Wind Analyses

The H*Wind analysis, a product of the Hurricane Research Division of NOAA's Atlantic Oceanographic and Meteorological Laboratory, brings together wind measurements from a variety of observation platforms into an objective analysis of the distribution of surface wind speeds in a tropical cyclone. This product is designed to improve understanding of the extent and strength of the wind field, and to improve the assessment of hurricane intensity. See http://www.aoml.noaa.gov/hrd/data sub/wind.html. The Hurricane Imaging Radiometer (HIRAD) is a new passive microwave remote sensor for hurricane observations that is currently under development by NASA Marshall Space Flight Center, NOAA Hurricane Research Division, the University of Central Florida and the University of Michigan. HIRAD is being designed to enhance the current real-time airborne ocean surface winds observation capabilities of NOAA and USAF Weather Squadron hurricane hunter aircraft using the operational airbome Stepped Frequency Microwave Radiometer (SFMR). Unlike SFMR, which measures wind speed and rain rate along the ground track directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (approximately 3 x the aircraft altitude, or approximately 2 km from space). The instrument is described in a separate paper presented at this conference. The present paper describes a set of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing instruments (air, surface, and space-based) are simulated from the output of a numerical model from the University of Miami, and those results are used to construct H*Wind analyses. Evaluations will be presented on the relative impact of HIRAD and other instruments on H*Wind analyses, including the use of HIRAD from 2 aircraft altitudes and from a space-based platform.

Miller, Timothy↗

DOE OSTI2023

Aerodynamic-whistles-based ultrasonic tone generators for bat deterrence

Novel ultrasonic bat deterrents based on aerodynamic whistles are proposed and investigated experimentally and numerically for their ability to generate ultrasound. The baseline deterrent, a single-whistle design inspired by Beeken [“Fluid ultrasonic generator,” U.S. patent 3,432,804 (1969)], is examined first. It consists of two resonating cavities/chambers. The whistle is “powered” by a regulated high-pressure air supply, and the performance of the whistle is examined for a range of supply pressures. Far-field acoustic measurements in the 20 Hz–50 kHz frequency range are made in an anechoic chamber. The noise measurements are supplemented with two- and three-dimensional unsteady Reynolds-averaged Navier–Stokes (uRANS) simulations to investigate the mechanisms of ultrasound generation. The Ffowcs Williams–Hawkings acoustic analogy is used with the three-dimensional uRANS results to predict the far-field radiation. The far-field acoustic predictions are in good agreement with the measurements in the anechoic chamber. The peak frequency (fundamental) of the radiated ultrasound of the baseline deterrent is approximately 23 kHz. Harmonics and a sub-harmonic of the fundamental tone are also observed. The numerical simulations show that the two resonating chambers of the baseline deterrent operate out-of-phase and Helmholtz resonance determines the whistling frequency over the range of supply air pressure considered. A six-whistle deterrent targeting a broad spectral coverage in the 20–50 kHz frequency range is designed, fabricated, and tested in the anechoic chamber. Each whistle in the deterrent is obtained by geometrically scaling the baseline whistle. Measurements show that the six-whistle ultrasound deterrent generates ultrasound with six dominant peaks in the designed frequency range when the supply air pressure exceeds 5 psig. The proposed ultrasound devices can be used for a variety of purposes including bat deterrence at wind turbines.

42 ENGINEERING↗

NASA NTRS2000

A Mechanistic Study of Nucleate Boiling Heat Transfer Under Microgravity Conditions

Experimental studies of growth and detachment processes of a single bubble and multiple bubbles formed on a heated surface have been conducted in the parabola flights of KC-135 aircraft. Distilled water and PF5060 were used as the test liquids. A micro-fabricated test surface was designed and built. Artificial cavities of diameters 10 microns, 7 microns and 4 microns were made on a thin polished Silicon wafer that was electrically heated by a number of small heating elements on the back side in order to control the surface superheat. Bubble growth period, bubble size and shape from nucleation to departure were measured under subcooled and saturation conditions. Significantly larger bubble departure diameters and bubble growth periods than those at earth normal gravity were observed. Bubble departure diameters as large as 20 mm for water and 6 mm for PF5060 were observed as opposed to about 3 mm for water and less than 1 mm for PF5060 at earth normal gravity respectively. It is found that the bubble departure diameter can be approximately related to the gravity level through the relation D(sub d) proportional 1/g(exp 1/2). For water,the effect of wall superheat and liquid subcooling on bubble departure diameter is found to be small.The growth periods are found to be very sensitive to liquid subcooling at a given wall superheat. However,the preliminary results of single bubble dynamics using PF5060 showed that the departure diameter increases when wall superheat is elevated at the same gravity and subcooling. Growth period of single bubbles in water has been found to vary as t(sub g) proportional g(exp -.93). For water, when the magnitude of horizontal gravitational components was comparable to that of gravity normal to the surface, single bubbles slid along the heater surface and departed with smaller diameter at the same gravity level in the direction normal to the surface. For PF5060, even a very small horizontal gravitational component caused the sliding of bubble along the surface. The numerical simulation has been carried out by solving under the condition of axisymmetry, the mass, momentum, and energy equations for the vapor and the liquid phases. In the model the contribution of micro-layer has been included and instantaneous shape of the evolving vapor-liquid interface is determined from the analysis. Consistent with the experimental results, it is found that effect of reduced gravity is to stretch the growth period and bubble diameter It is found that effect of reduced gravity is to stretch the growth period and bubble diameter at departure. The numerical simulations are in good agreement with the experimental data for both the departure diameters and the growth periods. In the study on dynamics of multiple bubbles, horizontal merger of 2,3 4,and 5 bubbles was observed. It is found that after merger of 2 and 3 bubbles the equivalent diameter of the detached bubble is smaller than that of a single bubble departing at the same gravity level. During and after bubble merger, liquid still fills the space between the vapor stems so as to form mushroom type bubbles. The experimental and numerical studies conducted so far have brought us a step closer to prediction of nucleate boiling heat fluxes under low gravity conditions. Preparations for a space flight are continuing.

Dhir, V. K.↗

NASA NTRS2007

An Analysis of Peak Wind Speed Data from Collocated Mechanical and Ultrasonic Anemometers

This study compared peak wind speeds reported by mechanical and ultrasonic anemometers at Cape Canaveral Air Force Station and Kennedy Space Center (CCAFS/KSC) on the east central coast of Florida and Vandenberg Air Force Base (VAFB) on the central coast of California. Launch Weather Officers, forecasters, and Range Safety analysts need to understand the performance of wind sensors at CCAFS/KSC and VAFB for weather warnings, watches, advisories, special ground processing operations, launch pad exposure forecasts, user Launch Commit Criteria (LCC) forecasts and evaluations, and toxic dispersion support. The legacy CCAFS/KSC and VAFB weather tower wind instruments are being changed from propeller-and-vane (CCAFS/KSC) and cup-and-vane (VAFB) sensors to ultrasonic sensors under the Range Standardization and Automation (RSA) program. Mechanical and ultrasonic wind measuring techniques are known to cause differences in the statistics of peak wind speed as shown in previous studies. The 45th Weather Squadron (45 WS) and the 30th Weather Squadron (30 WS) requested the Applied Meteorology Unit (AMU) to compare data between the RSA ultrasonic and legacy mechanical sensors to determine if there are significant differences. Note that the instruments were sited outdoors under naturally varying conditions and that this comparison was not designed to verify either technology. Approximately 3 weeks of mechanical and ultrasonic wind data from each range from May and June 2005 were used in this study. The CCAFS/KSC data spanned the full diurnal cycle, while the VAFB data were confined to 1000-1600 local time. The sample of 1-minute data from numerous levels on five different towers on each range totaled more than 500,000 minutes of data (482,979 minutes of data after quality control). The ten towers were instrumented at several levels, ranging from 12 ft to 492 ft above ground level. The ultrasonic sensors were collocated at the same vertical levels as the mechanical sensors and typically within 15 ft horizontally of each another. Data from a total of 53 RSA ultrasonic sensors, collocated with mechanical sensors were compared. The 1- minute average wind speed/direction and the 1-second peak wind speed/direction were compared.

Short, David A.↗

NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

NASA NTRS2015

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗