Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape Matching”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

A Technique to Measure Coronal Electron Density, Temperature, and Velocity Above 2.5 R from Sun Center Using Polarized Brightness Spectrum

The current model for the polarized brightness (pB) spectrum has a decades-long history of progressively incorporating its dependence on electron density N(e), temperature T(e), and flow velocity in the radial direction V(e). The pBNe spectrum follows the exact shape of the photosphere spectrum, which is not smooth, which is expected from the thermal Doppler broadening of the photosphere spectrum due to the high coronal T(e). The pBN(e) spectrum is smooth, but the free coronal electrons remain static and unaffected by solar wind, and the pBN(e)T(e)V(e) spectrum is red-shifted by electrons seeing a red-shifted photosphere spectrum as they flow away from the Sun as solar wind, which takes a radial direction above 2.5 R from Sun center. In this article, we review the progress of the above three model pB spectra in describing the observations and highlight the differences, first by comparing the three model pB spectra against wavelength using a model for Ne and constant values for T(e) and V (e), and second by generating three model 2D pB maps by integrating over a selected wavelength region in the three model pB spectra along lines of sight passing through the 14 July 2000 (“Bastille Day”) coronal mass ejection (CME) model, which contains 3D information on N(e), T(e), and V e. In this regard, the COronal Diagnostic EXperiment (CODEX) on the International Space Station (ISS) in 2024 will measure N(e), T(e), and V(e) by matching the measured pB with modeled pBN(e)T(e)V(e) in selected wavelength regions using multiple filters.

Electron density↗
NASA NTRSDate not supplied

An Approach for Practical Grid-Resolved Roughness Aerothermodynamic Simulations for Woven Heatshield Surfaces

An approach for simulating the aerothermal environment over a patterned roughness surface within an otherwise axisymmetric flow is developed. The patterned roughness is simplified into a sinusoidal shape that captures the basic features of a charred 3MDCP ablator. Two orientations of this pattern, which are equally flight relevant due to the forming process, enable periodic boundary conditions along the edges of a grid spanning a single roughness element. This simplified three-dimensional grid represents the minimum problem size for a patterned roughness simulation, which enables computationally efficient grid-resolved roughness simulations. This approach is validated by simulating the Langley Mach 6 measurements made on a similar sinusoidal surface, which results in convective heating within 5% of the experimental data. Applying this approach to Mars Sample Return (MSR) Earth Entry System (EES) flight cases, considering simulations with both 11 species air and 30 species air with ablation products, results in heating augmentation that is slightly lower than the widely-used Dahm correlation approach, with the difference being dependent on pattern orientation. This result provides the only insight into the behavior of the roughness augmentation at the high temperature reacting flow conditions present for EES flight cases, which are not captured by available ground test measurements. For shear augmentation, the simulations revealed that the surface-parallel pressure component is dominant. For arc-jet tests targeted to flight values for shear, the higher arc-jet pressures required to match flight values for shear result in up to a 400% increase in the rough wall shear value, which is not captured by the Dahm correlation. This indicates that, to avoid significant over-testing, grid-resolved roughness simulations are required to determine arc-jet conditions that produce flight shear levels in rough-wall scenarios.

Christopher O Johnston↗
NASA NTRS2023

A Technique to Measure Coronal Electron Density, Temperature, and Velocity Above 2.5 R⨀ from Sun Center using Polarized Brightness Spectrum

The current model for the polarized brightness (pB) spectrum has a decades-long history of progressively incorporating its dependence on electron density Ne, temperature Te, and flow velocity in the radial direction V e . The pB N e spectrum follows the exact shape of the photosphere spectrum, which is not smooth, which is expected from the thermal Doppler broadening of the photosphere spectrum due to the high coronal T e ; the pB N e T e spectrum is smooth, but the free coronal electrons remain static and unaffected by solar wind, and the pB N e T e V e spectrum is red-shifted by electrons seeing a red-shifted photosphere spectrum as they flow away from the Sun as solar wind, which takes a radial direction above 2.5 R ⨀ from Sun center. In this article, we review the progress of the above three model pB spectra in describing the observations and highlight the differences, first by comparing the three model pB spectra against wavelength using a model for Ne and constant values for T e and V e , and second by generating three model 2D pB maps by integrating over a selected wavelengthregion in the three model pB spectra along lines of sight passing through the 14 July 2000 (“Bastille Day”) coronal mass ejection (CME) model, which contains 3D information on N e , T e , and V e . In this regard, the COronal Diagnostic EXperiment (CODEX) on the International Space Station (ISS) in 2024 will measure N e , T e , and V e by matching the measured pB with modeled pB N e T e V e in selected wavelength regions using multiple filters.

Electron density↗
NASA NTRS1989

Monolithic Solid Oxide Fuel Cell development

The Monolithic Solid Oxide Fuel Cell (MSOFC) is an oxide-ceramic structure in which appropriate electronic and ionic conductors are fabricated in a honeycomb shape similar to a block of corrugated paperboard. These electronic and ionic conductors are arranged to provide short conduction paths to minimize resistive losses. The power density achievable with the MSOFC is expected to be about 8 kW/kg or 4 kW/L, at fuel efficienceis over 50 percent, because of small cell size and low resistive losses in the materials. The MSOFC operates in the range of 700 to 1000 C, at which temperatures rapid reform of hydrocarbon fuels is expected within the nickel-YSZ fuel channels. Tape casting and hot roll calendering are used to fabricate the MSOFC structure. The performance of the MSOFC has improved significantly during the course of development. The limitation of this system, based on materials resistance alone without interfacial resistances, is 0.093 ohm-sq cm area-specific resistance (ASR). The current typical performance of MSOFC single cells is characterized by ASRs of about 0.4 to 0.5 ohm-sq cm. With further development the ASR is expected to be reduced below 0.2 ohm-sq cm, which will result in power levels greater than 1.4 W/sq cm. The feasibility of the MSOFC concept was proven, and the performance was dramatically improved. The differences in thermal expansion coefficients and firing shrinkages among the fuel cell materials were minimized. As a result of good matching of these properties, the MSOFC structure was successfully fabricated with few defects, and the system shows excellent promise for development into a practical power source.

Myles, K. M.↗
DOE OSTI2023

Understanding Biases in Sample Preparation Techniques for Coupled Scanning Electron Microscopy and MAMA PuO 2 Morphological Analysis

In this project, the scanning electron microscopy (SEM) sampling method used during the statistical design study (SDS) was investigated to determine if any sampling biases were present in the analyzed data. Using standard particle size distribution powders from the National Institute of Standards and Technology (NIST 1984 standard reference material) with the origin wet dispersion method, it was determined that a bias to smaller particles was present. This was supported by theoretical calculations using Stokes’ law to determine the settling rate of spherical particles of roughly the same size and mass as those found in the SDS. Based on the theoretical calculations, it was determined that the settling rate for each of the 76 powder sets in the SDS could be unique based on specific particle shape and mass distributions, making a universal correction factor/formula not applicable. Therefore, priority shifted to developing an improved wet dispersion method that significantly reduced the particle settling rate for all particle size and shapes. This was achieved by replacing the original solvent (isopropyl alcohol) with a heavy liquid (lithium heteropolytungstates), which dramatically slowed the settling rate and allowed for the capture of a suitable homogeneous aliquot. SEM imaging and Morphological Analysis for Material Attribution (MAMA) software analysis were conducted on the NIST standard, and the SEM/MAMA data were compared to data captured by a dynamic image analysis particle size analyzer. The resulting data confirmed that the new wet dispersion method does indeed deliver an improved representative aliquot to the SEM stub. For instance, in the NIST certificate, the average particle size is ~17.1 µm ± 2.2 µm with a normal distribution. The initial wet dispersion method resulted in a drastically reduced average particle size of 6.1 µm in addition to a non-representative heavy bi-modal distribution whereas the improved LST wet dispersion method resulting in an average particle size that was much closer to the NIST certificate (12.7 µm) with a similar normal distribution. Although the improved method was still short of the NIST certificate average, atomic force microscopy analysis determined that the resulting ~20-25% reduction in size was due to particles sinking into the carbon sticky tape used for SEM imaging. It is believed that that this bias can be calibrated in a much more predicable manner than the original settling rate bias. In addition, the matching normal distribution curves between the NIST certificate and the heavy liquid method indicate a much-improved representative aliquot has been sampled and imaged. A surrogate CeO 2 powder was used to reflect PuO 2 more accurately and to aid in implementing radiological controls and shielding. The resulting data sets from the SEM/MAMA method and the particle size analyzer give almost identical average particle sizes and particle distribution statistics. Future work will re-analyze several select runs from the SDS to determine if morphological signatures can be found with the improved sampling method.

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

Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption in the ASCENDS 2011 Airborne Campaign

We have previously demonstrated an efficient pulsed, wavelength-resolved IPDA lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission. Our team participated in the 2010 ASCENDS airborne campaigns we flew airborne version of the CO2 and O2 lidar on the NASA DC-8. The CO2 lidar measures the atmospheric backscatter profiles and shape of the 1572.33 nm absorption line using 250 mW average laser power, 30 wavelength samples per scan and 300 scans per second. Most flights had 5-6 altitude steps to greater than 12 km, and clear CO2 line shapes were observed at all altitudes. Our post-flight analysis estimated the Iidar range and pulse energies at each wavelength every second. We then solved for the best-fit CO2 absorption line shape, and calculated the Differential Optical Depth (DOD) at the line peak. We compared these to CO2 DODs calculated from spectroscopy based on HITRAN 2008 and the conditions from airborne in-situ readings. Analysis of the 2010 measurements over the Pacific Ocean and Lamont OK shows the expected -linear change of the peak DOD with altitude. For measurements at altitudes greater than 6 km the random errors were approximately 0.3 ppm for 80 sec averaging times. After the 2010 flights we improved the airborne lidar's scan uniformity, calibration and receiver sensitivity. Our team participated in the seven ASCENDS science flights during late July and August 2011. These flights were made over a wide variety of surface and cloud conditions near the US, including over the central valley of California, over several mountain ranges, over both broken and solid stratus cloud deck over the Pacific Ocean, snow patches on mountain tops, over thin and broken clouds above the US Southwest and Iowa, and over forests near the WLEF tower in Wisconsin. Analyses show the retrievals of lidar range and CO2 column absorption, as wen as estimates of CO2 mixing ratio worked well when measuring over topography with rapidly changing height and reflectivity, through thin clouds and to stratus cloud tops. For regions where the CO2 concentration was relatively constant, the measured CO2 absorption profile (averaged for 50 sec) matched the predicted profile to better than 1% RMS error for all flight altitudes. For 1 & 10 second averaging, the scatter in the retrievals was limited by signal shot noise (i.e. the sigual photon count). Analysis to date shows the decrease in CO2 due to vegetation when flying easterward over the Great Plains as well as the increase in CO2 concentration in the vicinity ofthe coal-fired power plant in New Mexico. Examples of these and other results will be presented.

Abshire, James B.↗
DOE OSTI2025

Backlighting extended gas halos around luminous red galaxies: Kinematic Sunyaev-Zel’dovich effect from DESI Y1 and ACT data

The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel’dovich effect (kSZ). We precisely measure this effect (𝑆/𝑁 ≈ 10) via velocity stacking with 825,283 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over ≥ 4,000 deg 2 . We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift but clear trends with stellar mass and absolute magnitude in 𝑔, 𝑟, and 𝑧 bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5% confidence level, i.e., probability to exceed (PTE)⁢ = 0.005). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris 𝑧 = 0.5, PTE⁢ = 0.37) over weaker-feedback cases (IllustrisTNG 𝑧 = 0.8, PTE⁢ = 0.045), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around the DESI Y1 bright galaxy sample (BGS) and the emission-line galaxies (ELG) whose features match qualitative expectations. Finally, we forecast 𝑆/𝑁 ∼ 50 for future stacked kSZ measurements using data from the ACT, the DESI Y3, and the Rubin Observatory. Furthermore, these measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Characterization and controllability of radiated power via extrinsic impurity seeding in strongly negative triangularity plasmas in DIII-D

Experiments with extrinsic impurity seeding in strongly negative triangularity shapes in DIII-D achieved radiated power fractions (relative to input power) of up to ≈85% total radiation and ≈55% core radiation in steady-operating conditions. The relationship between core and total radiation was sensitive to impurity species and input power. Attempts to reach higher radiation levels via higher impurity flows resulted in radiative collapse disruptions. Nitrogen, neon, argon, and krypton were tested. Injection was by gas puffing, usually controlled by feeding back real-time estimates for core or total radiating fractions ($P_\textrm{rad} / P_\textrm{input}$). Argon and krypton were controlled by feeding back total radiating fraction, whereas neon was controlled by feeding back the core radiating fraction, due to the lower efficiency of neon as a divertor radiator. Nitrogen flows were pre-programmed. Reasonable total $P_\textrm{rad}$ control target following was achieved with argon or krypton. Control was more challenging with the neon/core radiation configuration, which was more prone to slow response and overshooting of the control target. Poor particle removal contributed to the control challenge: neon particle inventory within the last closed flux surface was roughly constant for up to 1 s (the longest duration tested) after neon injection was halted. However, separate experiments with laser-blow off of non-recycling impurities measured a short impurity confinement time, on the scale of the energy confinement time of ∼100 ms. Modeling with the Aurora impurity transport simulation matched experimental neon density profiles with full recycling (R = 1.0) and weak pumping but predicted rapid decreases in neon inventory if pumping were increased or recycling decreased. This indicates that changes outside the confined plasma (adding a well-placed pump) would improve controllability for all highly recycling species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

Mapping From an Instrumented Glove to a Robot Hand

An algorithm has been developed to solve the problem of mapping from (1) a glove instrumented with joint-angle sensors to (2) an anthropomorphic robot hand. Such a mapping is needed to generate control signals to make the robot hand mimic the configuration of the hand of a human attempting to control the robot. The mapping problem is complicated by uncertainties in sensor locations caused by variations in sizes and shapes of hands and variations in the fit of the glove. The present mapping algorithm is robust in the face of these uncertainties, largely because it includes a calibration sub-algorithm that inherently adapts the mapping to the specific hand and glove, without need for measuring the hand and without regard for goodness of fit. The algorithm utilizes a forward-kinematics model of the glove derived from documentation provided by the manufacturer of the glove. In this case, forward-kinematics model signifies a mathematical model of the glove fingertip positions as functions of the sensor readings. More specifically, given the sensor readings, the forward-kinematics model calculates the glove fingertip positions in a Cartesian reference frame nominally attached to the palm. The algorithm also utilizes an inverse-kinematics model of the robot hand. In this case, inverse-kinematics model signifies a mathematical model of the robot finger-joint angles as functions of the robot fingertip positions. Again, more specifically, the inverse-kinematics model calculates the finger-joint commands needed to place the fingertips at specified positions in a Cartesian reference frame that is attached to the palm of the robot hand and that nominally corresponds to the Cartesian reference frame attached to the palm of the glove. Initially, because of the aforementioned uncertainties, the glove fingertip positions calculated by the forwardkinematics model in the glove Cartesian reference frame cannot be expected to match the robot fingertip positions in the robot-hand Cartesian reference frame. A calibration must be performed to make the glove and robot-hand fingertip positions correspond more precisely. The calibration procedure involves a few simple hand poses designed to provide well-defined fingertip positions. One of the poses is a fist. In each of the other poses, a finger touches the thumb. The calibration subalgorithm uses the sensor readings from these poses to modify the kinematical models to make the two sets of fingertip positions agree more closely.

Goza, Michael↗
NASA NTRS2020

SWIRP (Submm-Wave and Long Wave InfraRed Polarimeter); Development and Characterization of a Sub-Mm Polarimeter for Ice Cloud Investigations

A major source of uncertainty in climate models is the presence, shape and distribution of ice particles in the uppermost layers of the clouds. The effects of this component are poorly constrained, turning ice particles into an almost-free variable in many climate models.NASA-GSFC is developing a new instrument aimed at measuring the size and shape of ice particles. The instrument consists of two sub-mm polarimeters (at 220 and 670 GHz) coupled with a long-wave infrared polarimeter at 10 micron. Each polarimeter has identical V-pol and H-pol channels; the axes of polarization are defined geometrically by the orientation of the waveguide elements, and the purity has been measured in the lab. The instrument is configured as a conical scanner, suitable for deployment as a payload on a small satellite or on a high-altitude sub-orbital platform. From a 400 km orbit, the instrument has a 3dB spatial resolution of 20 (10) km at 220 (670) GHz and a swath of 600 km over 180 degrees of view.The BAPTA (Bearing And Power Transfer Assembly) carries heritage from the SSMIS design, now in its 22nd year of on-orbit operation, but with a much reduced SWaP (Size Weight and Power) footprint, suitable for a small satellite.The main components of the instrument have been fabricated and are undergoing final testing prior to their integration as a single unit. The sub-mm channels have dedicated secondary reflectors which illuminate a shared primary reflector. The receiving units are placed behind the focal point of the optical arrangement, so that all beams equally illuminate the primary reflector and are almost co-located on the ground (within a single 220 GHz footprint). Primary and secondary beam patterns have been measured and verified to match the as-designed expectations. A Zytex (TM) window is deployed to protect the secondary reflectors and the feed horns from debris and other contaminants, and to reduce the heat load from the active (hot) IR calibration unit. The insertion loss of Zytex has been measured and is accounted in the calibration equation of the sub-mm channels.The radiometric performance of the sub-mm receivers has been characterized in the lab and under operational conditions of temperature and pressure.This paper discusses the design constraints on the sub-mm components, details of the scientific goals and their flowdown, and describes the characterization of the polarimeters. Options to optimize the layout and distribution of the masses within the assembly, with the goal of making the instrument even more compact and fully-compatible with cubesat-class satellites will be presented.

De Amici, G.↗
NASA NTRS1993

Meteorological analysis of the November 25 FIRE cirrus-2 case: A well-defined ridge-crest cirrus system over Oklahoma

On the morning of 25 Nov. 1991, a cirrus cloud system formed in a region extending eastward from the continental divide into central Oklahoma and southwestern Kansas. The system moved slowly eastward during the day and maintained a relatively constant scale and shape. From the satellite perspective, this cloud system had a number of morphological similarities to the 27-28 Oct. FIRE (First International Satellite Cloud Climatology Project Regional Experiment) Cirrus-1 case analyzed by Starr and Wylie. By comparison, the present case had much better rawinsonde coverage. The rawinsonde network was much more extensive and the cloud system was relatively slow moving and well-located with respect to the enhanced sounding network and the National Weather Service (NWS) Wind Profiler Network. In addition to the rawinsonde soundings (1200 and 1800 UTC on 25 Nov. and 0000 and 0600 UTC on 26 Nov.) and satellite observations, four aircraft sorties were flown in the northern portion of this cirrus cloud system from mid-morning to mid afternoon. These included microphysical-radiation profiling missions by the NCAR Sabreliner near Tulsa (about midway between COF and OUN) and the NCAR King Air near Ponca City (about 80 km upwind of Tulsa), a microphysical profiling mission by the UND Citation near Tulsa, and a NASA ER-2 overflight. Thus, very good data are available for this case although, unfortunately, the cloud system remained south of the FIRE Hub site at Coffeyville (COF) where extensive ground-based remote sensor was located. The northern edge of the cloud system was visible from COF during most of the day. Given the wealth of available in situ and remote sensing observations, the similarity to the FIRE Cirrus-1 focal case study, and the particularly good match in scale and location of the cloud system with the rawinsonde and profiler networks, analysis of this case may yield very useful results on the life-cycle of this cirrus cloud system. The quality of the meteorological description for this case was assessed based on our preliminary analysis of the rawinsonde observations.

Starr, David↗
DOE OSTI2020

Determining Compression Characteristics of Honeycomb Material - 19674

The objective is to determine compression test characteristics of stainless steel honeycomb material to be able to represent honeycomb structures accurately in analytical models used to simulate hypothetical accident scenarios of shipping packages. Honeycomb is a material used primary in the aerospace industry due to its high strength-to-weight ratio. Because it doesn't have a shelf life and can withstand high heats, it is an excellent candidate for a structural material in package designs. Honeycomb has orthotropic material properties. The honeycomb currently being evaluated is made of metal ribbons spot welded together to form hexagon pattern between metal plates. The hexagon ribbons are brazed to the metal plates. The complexity of the honeycomb significantly lengthens the simulation time, which is further complicated when the brazing, welding, and imperfections of the material is considered. This means to affectively represent the honeycomb in simulation programs, such as Abaqus, the material needs to be approximated as a uniform orthotropic material. This requires material properties in each of the three directions, T, W, and L. The T direction is defined as the 'strong' direction, perpendicular to the honeycomb sheet. The L direction is parallel the ribbon and the W direction is perpendicular to the ribbon. Honeycomb has 3 different compression stages. Stage 1 is the initial compression where the honeycomb maintains its structural integrity and does not permanently deform from loads from normal operations. Stage 2 is when the initial buckling causes irreversible damage to the honeycomb. Stage 2 is the one we are most interested in because it absorbs the most energy from a hypothetical accident scenario. During Stage 2, the honeycomb fails layer-by-layer, indicating that the more layers, the longer material crushing is sustained, and thus the more energy that is absorbed. Stage 3 is final compression, similar to compressing solid metal. During stage 3 the stress increases with a diminishing rate of elongation and the honeycomb is completely failed where the plates between the honeycomb core sandwich the crushed honeycomb ribbon. The stress-strain graph bellow shows compression test of 3 different samples in the T direction. The red lines divide the three stages. The far left is Stage 1, the middle stage 2, and the right stage 3. The data collected is displayed in the table below. The expressions in the table represent the regression formula that represents each stage of the compression on a stress-strain graph. Stage 1 is assumed to intersect with the origin; data for stage 3 for the crush in W and L directions where unable to be gathered due to the nature of failure for those directions. Stage 1 and 2 are linear regressions while stage 3 is a degree 2, to best match the curve. Linear regression is used on stage 2 because ultimately that represents the energy absorbed. With more data, a degree n x 2 regression would be a more appropriate regression for stage 2, where n is the number of layers. The table below is the approximation of each scenario and stage. Each stage starts and ends at the intersection of the next stage. There where two main objectives for these tests: firstly to understand how honeycomb performs under extreme compression, and secondly to be able to numerically represent the honeycomb structure. Both of these where completed to varying degrees. It is important understood how the honeycomb fails. If it is crushed in the T direction it retains its integrity even after being crushed, however when crushed in the L and W crush direction, if it fails, it disintegrates, and loses all integrity. Also, the more layers the more time the material spends in stage 2. Failure is started by buckling, thus if there is any imperfection, the stress will not spike but transition straight into stage 2. The data collected and aggregated can be used for initial simulation of honeycomb used in packages. The initial testing has set the ground work for more data to be collected in order to verify results and to allow more confidence in the simulation results. This is only the initial data collected. The next steps is to continue to collect more data to verify results and to test more variants of honeycomb, with different brazing, layers, and shape.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Hybrids of Solar Sail, Solar Electric, and Solar Thermal Propulsion for Solar-System Exploration

Solar sails have long been known to be an attractive method of propulsion in the inner solar system if the areal density of the overall spacecraft (S/C) could be reduced to approx.10 g/sq m. It has also long been recognized that the figure (precise shape) of useful solar sails needs to be reasonably good, so that the reflected light goes mostly in the desired direction. If one could make large reflective surfaces with reasonable figure at an areal density of approx.10 g/sq m, then several other attractive options emerge. One is to use such sails as solar concentrators for solar-electric propulsion. Current flight solar arrays have a specific output of approx. 100W/kg at 1 Astronomical Unit (AU) from the sun, and near-term advances promise to significantly increase this figure. A S/C with an areal density of 10 g/sq m could accelerate up to 29 km/s per year as a solar sail at 1 AU. Using the same sail as a concentrator at 30 AU, the same spacecraft could have up to approx. 45 W of electric power per kg of total S/C mass available for electric propulsion (EP). With an EP system that is 50% power-efficient, exhausting 10% of the initial S/C mass per year as propellant, the exhaust velocity is approx. 119 km/s and the acceleration is approx. 12 km/s per year. This hybrid thus opens attractive options for missions to the outer solar system, including sample-return missions. If solar-thermal propulsion were perfected, it would offer an attractive intermediate between solar sailing in the inner solar system and solar electric propulsion for the outer solar system. In the example above, both the solar sail and solar electric systems don't have a specific impulse that is near-optimal for the mission. Solar thermal propulsion, with an exhaust velocity of the order of 10 km/s, is better matched to many solar system exploration missions. This paper derives the basic relationships between these three propulsion options and gives examples of missions that might be enabled by such hybrids.

In-space Propulsion↗
DOE OSTI2022

OPALCO - Decatur Island Solar and Energy Storage Project: An Assessment of Battery Technical Performance

Orcas Power & Light Cooperative (OPALCO) is a member-owned, nonprofit cooperative utility that provides energy services to approximately 11,200 customers across 20 islands in San Juan County, Washington. OPALCO’s mostly hydroelectric power is generated by Bonneville Power Administration and delivered to the islands by submarine cables. In 2016, as part of the second round of funding from the Washington state Clean Energy Fund, OPALCO received a $1 million matching grant to support a project that deployed a 504-kW LG community photovoltaic (PV) system in combination with a 1 MW/2 MWh lithium-iron-phosphate battery energy storage system (BESS) on Decatur Island, Washington. The Decatur Island Substation is essential to ensuring reliable energy for the residents of the San Juan Islands as it is the point of interconnection with the mainland transmission system. The BESS, in combination with the community solar array, will deliver an innovative method to both defer the costly upgrade of the transmission system and allow for other high-value applications intended to benefit the utility and its customers. In 2018, Pacific Northwest National Laboratory (PNNL) completed a preliminary economic assessment for several identified use cases in collaboration with OPALCO. Between August 2021 and May 2022, extensive testing was conducted, and the results were used to assess the technical performance of the BESS subjected to actual field operations. Both reference performance and use case tests were performed: (a) Reference performance tests assess the general technical capabilities of the BESS, such as energy capacity, round-trip efficiency (RTE), ramp rate, and signal tracking capability. These are the first tests performed (baseline) and are repeated after use case tests (post cycle). A standardized U.S. Department of Energy (DOE) energy storage performance protocol was used to characterize the BESS, including representative duty cycle profiles, test procedure guidance, and calculation guidance for determining key characteristics. (b) Use case tests examine the performance of the BESS for specific use cases using duty cycles developed by PNNL in collaboration with OPALCO. Four use cases were selected for testing: 1) demand charge reduction, 2) load shaping, 3) outage mitigation, and 4) transmission deferral. The use case duty cycles were developed based on utility and site-specific characteristics in addition to the technical characteristics and physical capabilities of the BESS. Use case tests were performed between the baseline and post cycle tests. This report describes the BESS and its components, presents testing and performance analysis results, and shares key insights and lessons learned from this project. Outcomes of the tests and analyses will help OPALCO understand the performance of the Decatur Island BESS in its current state and design appropriate operational strategies for this and other BESSs over the long term.

14 SOLAR ENERGY↗
DOE OSTI2021

Low-Cost Aero Technology Demonstrations

The main focus of this work was to demonstrate the use of polymeric additive manufacturing (AM) to create tooling for both preforming and consolidation. Polymeric tooling was utilized where both modest and higher pressures are used for part consolidation. The key focus for the AM tooling development was for fabrication of complex structures such as ducting, C-channel stiffened skins, and airfoils where conventional male tooling would typically be trapped in the cured part. The AM tooling was evaluated for use as a tool master used to fabricate and re-shape deformable/re-formable mandrels based on SpinTech’s shape memory composite technology known as Smart Tooling. The AM tooling was also evaluated for use as a mold for composite infusion and consolidation. Key performance parameters were tracked for project schedule completion with each step comprising of “art to part” cycle time, cost, and model fidelity for dimensions, performance, and cost. The primary focus of this demonstration was to determine if a 50% cost reduction was achievable, for each AM tooling-set, as compared to conventional processes. UDRI leveraged project partner SpinTech, who manufactures tools and parts in these categories and thus provided a baseline regarding current best practices and provided valuable feedback during the entirety of this demonstration. This demonstration primarily focused on the use of AM tooling for fabrication of three composite component structures which are typically utilized in aircraft and comprise salient geometric features of broad interest. These components are often tooling intensive and have features requiring extraction of male tools which are usually trapped by the geometry. The three structures selected by the team included: 1) A one-piece airfoil shell comprised of compound contours where male tooling would be trapped unless the part were manufactured in two halves as is typically the case. 2) A one-piece duct used for air handling, comprised of compound contours where male tooling would be trapped unless the part were manufactured in two halves, or a washout mandrel were to be used. 3) A co-cured C-channel stiffened skin where typically C-channels would be individually manufactured and then bonded to a cured skin. The demonstration was comprised of three main tasks: • Task 1: AM Tool Feasibility Study – ensure the AM tooling meets the performance requirements as specified by SpinTech to match baseline performance. • Task 2: Complex Tool Demonstration – Fabricate tooling, preforms, and parts representative of an airfoil and duct. • Task 3: Large Aerostructure Fabrication Demonstration – Fabricate tooling, preforms, and part representative of a C-channel stiffened skin. With the conclusion of this project, a decision tree was developed to determine the key considerations necessary to determine if use of AM tooling for the three selected structures was able to attain the same quality as historically achieved on metallic tooling, while providing a significant cost reduction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗
DOE OSTI2022

Energy Northwest - Horn Rapids Solar and Storage: An Assessment of Battery Technical Performance

Chartered in 1957 as a joint action agency of the state, Energy Northwest (ENW) is a consortium of 27 public utility districts and municipalities across Washington state. ENW takes advantage of economies of scale and shared services to help utilities run their operations more efficiently and at lower cost, to the benefit of more than 1.5 million customers. ENW develops, owns, and operates a diverse mix of electricity generating resources, including hydro, solar, and wind projects – and the Northwest’s only active nuclear energy facility. These projects provide enough reliable, affordable, and environmentally responsible energy to power more than a million homes each year, and that carbon-free electricity is provided at the cost of generation. The agency continually explores new generation projects to meet its members’ needs. In 2017, as part of the second round of funding from the Washington state Clean Energy Fund, the Washington State Department of Commerce granted up to $3 million in matching funds to develop an estimated $6.5 million project that deployed a 4-MW, 20-acre solar generating array of photovoltaic (PV) panels coupled with a 1 MW/5.5 MWh lithium-iron-phosphate battery energy storage system (BESS) in Richland, Washington. The combination of PV and BESS will provide a predictable, renewable generating source and will also serve as a training ground for solar and battery technicians throughout the nation. The City of Richland will purchase the power from the project and utilize the benefits of the energy storage. The project provides Washington state with its first opportunity to integrate a large-scale solar and storage facility into its clean mix of hydro, nuclear, and wind resources. This first-of-its-kind facility combines solar generation with battery storage and technician training. In 2019, Pacific Northwest National Laboratory (PNNL) worked with ENW to assess the integrated PV and BESS in representative use cases that could benefit the City of Richland. Between March and May 2022, extensive testing was conducted, and the results were used to assess the technical performance of the BESS subjected to actual field operations. Both reference performance and use case tests were performed: (A) Reference performance tests assess the general technical capabilities of the BESS, such as energy capacity, round-trip efficiency (RTE), ramp rate, and signal tracking capability. These are the first tests performed (baseline), and they are repeated after use case tests (post cycle). A standardized U.S. Department of Energy (DOE) energy storage performance protocol was used to characterize the BESS, including representative duty cycle profiles, test procedure guidance, and calculation guidance for determining key characteristics. (B) Use case tests examine the performance of the BESS for specific use cases using duty cycles developed by PNNL in collaboration with ENW. Five use cases were selected for testing: 1) demand charge reduction, 2) load shaping, 3) transmission charge reduction, 4) Volt-VAR service, and 5) outage mitigation. The use case duty cycles were developed based on utility and site-specific characteristics in addition to the technical characteristics and physical capabilities of the BESS. Use case tests were performed between the baseline and post cycle tests. This report describes the BESS and its components, presents testing and performance analysis results, and shares key insights and lessons learned from this project. Outcomes of the tests and analyses will help ENW understand the performance of the Horn Rapids BESS in its current state and design appropriate operational strategies for this and other BESSs over the long term.

14 SOLAR ENERGY↗