Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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 631 records · Page 35

Suborbital Lunar Gravity Experiment of an Electrodynamic Regolith Conveyor

NASA’s Swamp Works Electrostatics and Surface Physics Laboratory (ESPL) is planning to test a 4-phase electrodynamic regolith conveyor in simulated lunar gravity on a suborbital flight in 2023. A 4-phase electrodynamic regolith conveyor could convey regolith without the risk of rotating or vibratory actuation, which could jam or require regular maintenance due to the abrasive nature of lunar regolith. Another advantage of electrodynamic conveying is the reduction of conveying power, which is important considering the limited scale of available early-stage lunar power systems. The current state of the art(SOA) for lunar regolith conveying is based on recent NASA system studies for oxygen production plants. These plant designs require conveying rates around 100 kg/hr, to produce 10 mT/yr of oxygen from the regolith. To accomplish this, conventional augers and vibratory spiral conveyors have been identified as the SOA or the leading candidates due to their extensive use in the terrestrial material handling industry. At NASA KSC, the use of dynamic electric fields, generated by alternating high voltage on electrodes, has been developed as a dust mitigation solution known as the Electrodynamic Dust Shield (EDS). The EDS is being developed for lenses, solar panels, radiators, fabric and seals. ESPL researchers have shown the ability to move thin layers (~ a few mm) of dust with mW of power. In Academia, researchers have shown the ability to electrodynamically convey regolith up to ~1 kg/hr with a4-phase EDS and electromagnetically convey particles that would suggest a lunar regolith conveying rate of ~25 kg/hr. This paper will describe the design and expected operation of a 4-phase electrodynamic regolith conveyor in lunar gravity. This paper will also describe how an electrodynamic conveyor could be scaled for transporting ISRU relevant flow rates.

Aaron D S Olson↗

Acoustic and turbulence measurements of a tone-excited jet with and without flight simulation

Acoustic as well as turbulence measurements were made of tone-excited jets to obtain an understanding of the broadband noise augmentation mechanism. Results for both heated and unheated jets with and without flight simulation are presented for a range of excitation frequencies and levels, as well as for the zero order and first order spinning modes. It is argued that although the phase-locked large-scale turbulence structure is at the root of the noise amplification process, the actual noise generation mechanism lies in the small-scale turbulence. Results show that (1) the phase velocity of the excited large scale instability waves is subsonic relative to the ambient fluid and (2) broadband jet noise is almost uniform at all frequencies. Moreover, when the large-scale and small-scale structures are plotted against forward velocity (for a given upstream excitation level of 141 dB), the change in large-scale turbulence is negligible with forward velocity, whereas the small-scale turbulence decreases for both excited and unexcited jets. Also, the corresponding difference in far-field noise does not appear to change significantly, indicating that the changes in small-scale turbulence are responsible for jet noise amplification.

Ahuja, K. K.↗

Pan and Tilt Platform for a Thermal Vacuum Video Camera

The purpose of the Space Simulation Test Engineering (SSTE) Section at NASA Goddard Space Flight Center is "to provide full service test support, engineering expertise, and state-of-the-art facilities to perform simulated space and thermal cycle testing of Earth orbit and deep space flight hardware"(ref. 1). During the Integration & Testing phase of space flight projects, flight hardware must successfully demonstrate its functionality under conditions similar to the actual environment experienced in space. When testing hardware with mechanical components, such as antenna extensions, door openings, or solar panel deployments, visual observation would allow test personnel to verify mechanical operations in real-time. A video camera can supply data that the usual sensors inside a TV chamber cannot provide. Since mechanical functional testing could affect multiple regions of a test item, a camera with a fixed position may not be able to monitor and record the performance of different components of a mechanism, especially for large test items. By mounting the camera to a Pan and Tilt Platform, the camera can cover a larger area and its field-of-view (FOV) will increase significantly. This paper details the design process for the Pan and Tilt Platform.

Lui, Yan↗

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

compressor↗

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

unsteady↗

Numerical Simulation of Boundary Layer Ingesting (BLI) Inlet-Fan Interaction

Future civil transport designs may incorporate engine inlets integrated into the body of the aircraft to take advantage of efficiency increases due to weight and drag reduction. Additional increases in engine efficiency are predicted if the inlet ingests the lower momentum boundary layer flow. Previous studies have shown, however, that efficiency benefits of Boundary Layer Ingesting (BLI) ingestion are very sensitive to the magnitude of fan and duct losses, and blade structural response to the non-uniform flow field that results from a BLI inlet has not been studied in-depth. This paper presents an effort to extend the modeling capabilities of an existing rotating turbomachinery unsteady analysis code to include the ability to solve the external and internal flow fields of a BLI inlet. The TURBO code has been a successful tool in evaluating fan response to flow distortions for traditional engine/inlet integrations, such as the development of rotating stall and inlet distortion through compressor stages. This paper describes the first phase of an effort to extend the TURBO model to calculate the external and inlet flowfield upstream of fan so that accurate pressure distortions that result from BLI configurations can be computed and used to analyze fan aerodynamics and structural response. To validate the TURBO program modifications for the BLI flowfield, experimental test data obtained by NASA for a flushmounted S-duct with large amounts of boundary layer ingestion was modeled. Results for the flow upstream and in the inlet are presented and compared to experimental data for several high Reynolds number flows to validate the modifications to the solver. Quantitative data is presented that indicates good predictive capability of the model in the upstream flow. A representative fan is attached to the inlet and results are presented for the coupled inlet/fan model. The impact on the total pressure distortion at the AIP after the fan is attached is examined.

Inlet Fan Interactions↗

The interaction of collisionless shocks in astrophysical plasmas

The interaction between collisionless shocks, such as may arise in a number of situations in space plasmas, is examined by means of hybrid numerical simulations. The production of energetic particles by shock collisions is investigated. When the shocks are quasi-perpendicular, ions with energies greater than 10E(0) are produced, where E(0) is the kinetic energy of the unshocked ions. The effect is optimized when the two colliding shocks have roughly equal strengths. As the shocks become closer to quasi-parallel, the number and energy of the accelerated ions increases. Energies in excess of 30E(0) are obtained in the quasi-parallel regime, with a small fall-off when the shocks are exactly parallel. For quasi-parallel collisions, the ion energization can be attributed to the interaction of hot ions with both strong electromagnetic waves and the motional electric field located at the shocks. These results are applied to shock collisions in the distant solar wind, at planetary bow shocks, as well as to the impulsive phase of solar flares.

Cargill, P. J.↗

Regional Application of an Ecosystem Production Model for Studies of Biogeochemistry in the...

The degree to which primary production, soil carbon, and trace gas fluxes in tropical forests of the Amazon are limited by moisture availability and other environmental factors was examined using an ecosystem modeling application for the country of Brazil. A regional geographic information system (GIS) serves as the data source of climate drivers, satellite images, land cover, and soil properties for input to the NASA Ames-CASA (Carnegie-Ames-Stanford Approach) model over a 8-km grid resolution. Simulation results supports the hypothesis that net primary production (NPP) is limited by cloud interception of solar radiation over the humid northwestern portion of the region. Peak annual rates for NPP of nearly 1.4 kg C m-2yr -1are localized in the seasonally dry eastern Amazon in areas that we assume are primarily deep-rooted evergreen forest cover. Regional effects of forest conversion on NPP and soil carbon content are indicated in the model results, especially in seasonally dry areas. Comparison of model flux predictions along selected eco-climatic transects reveal moisture, soil, and land use controls on gradients of ecosystem production and soil trace gas emissions (CO2, N2O, and NO). These results are used to formulate a series of research hypotheses for testing in the next phase of regional modeling, which includes recalibration of the light-use efficiency term in CASA using field measurements of NPP, and refinements of vegetation index and soil property (texture and potential rooting depth) maps for the region.

Potter, C. S.↗

Modeling the Ionosphere-Thermosphere Response to a Geomagnetic Storm Using Physics-based Magnetospheric Energy Input: OpenGGCM-CTIM Results

The magnetosphere is a major source of energy for the Earth's ionosphere and thermosphere (IT) system. Current IT models drive the upper atmosphere using empirically calculated magnetospheric energy input. Thus, they do not sufficiently capture the storm-time dynamics, particularly at high latitudes. To improve the prediction capability of IT models, a physics-based magnetospheric input is necessary. Here, we use the Open Global General Circulation Model (OpenGGCM) coupled with the Coupled Thermosphere Ionosphere Model (CTIM). OpenGGCM calculates a three-dimensional global magnetosphere and a two-dimensional high-latitude ionosphere by solving resistive magnetohydrodynamic (MHD) equations with solar wind input. CTIM calculates a global thermosphere and a high-latitude ionosphere in three dimensions using realistic magnetospheric inputs from the OpenGGCM. We investigate whether the coupled model improves the storm-time IT responses by simulating a geomagnetic storm that is preceded by a strong solar wind pressure front on August 24, 2005. We compare the OpenGGCM-CTIM results with low-earth-orbit satellite observations and with the model results of Coupled Thermosphere-Ionosphere-Plasmasphere electrodynamics (CTIPe). CTIPe is an up-to-date version of CTIM that incorporates more IT dynamics such as a low-latitude ionosphere and a plasmasphere, but uses empirical magnetospheric input. OpenGGCMCTIM reproduces localized neutral density peaks at approx. 400 km altitude in the high-latitude dayside regions in agreement with in situ observations during the pressure shock and the early phase of the storm. Although CTIPe is in some sense a much superior model than CTIM, it misses these localized enhancements. Unlike the CTIPe empirical input models, OpenGGCM-CTIM more faithfully produces localized increases of both auroral precipitation and ionospheric electric fields near the high-latitude dayside region after the pressure shock and after the storm onset, which in turn effectively heats the thermosphere and causes the neutral density increase at 400 km altitude.

Connor, Hyunju K.↗

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Lin, Yuh-Lang↗

Computational Analyses of Pressurization in Cryogenic Tanks

A) Advanced Gas/Liquid Framework with Real Fluids Property Routines: I. A multi-fluid formulation in the preconditioned CRUNCH CFD(Registered TradeMark) code developed where a mixture of liquid and gases can be specified: a) Various options for Equation of state specification available (from simplified ideal fluid mixtures, to real fluid EOS such as SRK or BWR models). b) Vaporization of liquids driven by pressure value relative to vapor pressure and combustion of vapors allowed. c) Extensive validation has been undertaken. II. Currently working on developing primary break-up models and surface tension effects for more rigorous phase-change modeling and interfacial dynamics B) Framework Applied to Run-time Tanks at Ground Test Facilities C) Framework Used For J-2 Upper Stage Tank Modeling: 1) NASA MSFC tank pressurization: a) Hydrogen and oxygen tank pre-press, repress and draining being modeled at NASA MSFC. 2) NASA AMES tank safety effort a) liquid hydrogen and oxygen are separated by a baffle in the J-2 tank. We are modeling pressure rise and possible combustion if a hole develops in the baffle and liquid hydrogen leaks into the oxygen tank. Tank pressure rise rates simulated and risk of combustion evaluated.

Ahuja, Vineet↗

Analytic Spacecraft Attitude and Rate Estimation Performance During Attitude Sensor Outages

Analytic expressions for spacecraft attitude and rate estimation performance of an attitude estimation filter in terms of sensor specifications are useful tools for spacecraft design. Farrenkopf (1978) famously found analytic expressions for steady-state pre-update and post-update attitude and gyro bias estimate error variances for an attitude estimation filter for a single-axis spacecraft with a Rate Output Gyro (ROG). Markley and Reynolds (2000) extended the analysis for a Rate-Integrating Gyro (RIG) with angle white noise. These expressions allow for the rapid evaluation of system performance during preliminary mission design phases. One contribution of this paper is the analytic calculation of the steady-state pre-update and post-update angular rate estimate uncertainty for both the ROG and RIG cases. The primary contribution of this paper is the extension of the results for both the ROG and the RIG cases to the situation of an attitude sensor outage. This situation arises frequently in practice; for example when a star sensor’s field of view is occluded, when a star sensor’s readings are unreliable during a thruster burn that vibrates the spacecraft, or during star sensor outages due to radiation upsets. Analytic expressions for the attitude estimate uncertainty, gyro bias estimate uncertainty, and angular rate estimate uncertainty are given in terms of the attitude sensor outage interval, the star tracker measurement noise, and gyro noise parameters. Validity of the analytic results is demonstrated via Monte Carlo simulation.

Galante, Joseph M.↗

Analytic Spacecraft Attitude and Rate Estimation Performance During Attitude Sensor Outages

Analytic expressions for spacecraft attitude and rate estimation performance of an attitude estimation filter in terms of sensor specifications are useful tools for spacecraft design. Farrenkopf (1978) famously found analytic expressions for steady-state pre-update and post-update attitude and gyro bias estimate error variances for an attitude estimation filter for a single-axis spacecraft with a Rate Output Gyro (ROG). Markley and Reynolds (2000) extended the analysis for a Rate-Integrating Gyro (RIG) with angle white noise. These expressions allow for the rapid evaluation of system performance during preliminary mission design phases. One contribution of this paper is the analytic calculation of the steady-state pre-update and post-update angular rate estimate uncertainty for both the ROG and RIG cases. The primary contribution of this paper is the extension of the results for both the ROG and the RIG cases to the situation of an attitude sensor outage. This situation arises frequently in practice; for example when a star sensor’s field of view is occluded, when a star sensor’s readings are unreliable during a thruster burn that vibrates the spacecraft, or during star sensor outages due to radiation upsets. Analytic expressions for the attitude estimate uncertainty, gyro bias estimate uncertainty, and angular rate estimate uncertainty are given in terms of the attitude sensor outage interval, the star tracker measurement noise, and gyro noise parameters. Validity of the analytic results is demonstrated via Monte Carlo simulation.

Galante, Joseph M.↗

Finite Element Simulation of Solid Rocket Booster Separation Motors During Motor Firing

One of the toughest challenges facing Solid Rocket Booster (SRB) engineers is to ensure that any design changes made to the Shuttle-Derived Booster Separation Motors (BSM) for future space exploration vehicles is able to withstand the increasingly hostile motor firing environment without cracking its critical component - the graphite throat. This paper presents a critical analysis methodology and techniques for assessing effects of BSM design changes with great accuracy and precision. For current Space Shuttle operation, the motor firing occurs at SRB separation - approximately 125 seconds after Shuttle launch at an altitude of about 28 miles. The motor operation event lasts about two seconds, however, the surface temperature of the graphite throat increases approximately 3400 F in less than one second with a corresponding increase in surface pressure of approximately 2200 pounds per square inch (psi) in less than one-tenth of a second. To capture this process fully and accurately, a two-phase sequentially coupled thermal-mechanical finite element approach was developed. This method allows the time- and location-dependent pressure fields to interact with the spatial-temporal thermal fields throughout the operation. The material properties of graphite throat are orthotropic and temperature-dependent. The analysis involves preload and multiple body contacts.

Yu. Weiping↗

A Historical Overview of the NASA Orbital Debris Program Office’s Laboratory Optical Measurements

The NASA Orbital Debris Program Office (ODPO) has used laboratory measurements to help bring ground-based measurements together with models to ascertain Earth-orbiting target parameters of interest to support various orbital debris models. In 2005, the Optical Measurement Center (OMC) was established to simulate space-based illumination conditions using equipment and techniques that recreate telescopic observations, particularly source-target-sensor orientations. The intent was to recreate light curves using known aspect angles of known targets and phase angles (angle is defined by the vertex between illumination source-object-detector) to complement telescopic observations that could be used to update the current optical size estimation model (OSEM) – a model that converts object brightness into size for orbital debris models. To support the above goals, the laboratory has undergone several equipment upgrades to increase capabilities over almost 20 years of operation. The primary instrumentation acquires reflectance measurements and includes a solar-like light source, CCD camera with astrometric filters, and robotic arm. A rotary arm was added approximately five years after full operation to allow acquisition through a full 360° range of phase angles. Another part of the OMC instrumentation is a field spectrometer, predominately used for field operations to acquire pre- and post-flight spacecraft material spectral measurements. Additionally, reflectance spectroscopy of various materials is also of interest resulting from hypervelocity impact tests, pristine spacecraft materials, or samples of materials that are used in spacecraft design. These measurements are stored in NASA’s Spectral Material Database, a resource that is still being populated today. The study of spectral measurements also enabled the development of spectral unmixing routines to support the identification of spacecraft materials from spectral data gathered by ground based telescopes. Preliminary OMC investigations focused on feasibility studies to acquire 360° rotation light curves of simple shapes at a single-phase angle and extended to measurements of representative fragments from ground-based explosion tests. To correlate the light curves with ground-based optical measurements, a focused study on high area to mass materials was conducted in support of a newly identified population (at the time) in geosynchronous orbit (GEO) consisting of multi-layered insulation. To further characterize orbital debris, a larger selection of materials was analyzed using laboratory photometric measurements that included representative targets from pristine spacecraft materials and ground-based impact tests. Around 2012, an initiative was requested to understand the feasibility of active debris removal (ADR) of larger targets using grappling methods for spent rocket bodies. Using a priori information on selected targets, scaled-down versions of rocket bodies were generated thanks to improvements in 3D printing technology and machining. These targets were studied in the OMC to understand rotation characteristics. These were compared with telescopic data to determine if the tumble and rotation angles would allow ADR. In 2013, the OMC focused on combining spectral measurements with photometric data to characterize GEO orbital debris. Several years later, NASA acquired a Titan III Transtage test article from “The Boneyard” with a high-resemblance to on-orbit Titan III Transtage rocket bodies, allowing physical access to a representative rocket body that suffered fragmentations in GEO. This prompted the creation of 3D models using lidar technology and spectral measurements of the materials. Focused research also transitioned to specific materials (i.e., solar cells) when telescopic surveys requested characterization of specific GEO targets. In the different research products presented, the focus has been to understand the various parameters that influence optical size estimation, including albedo, phase functions, and brightness variations. Work in this area continues with newer sources of data, including DebriSat, a high-fidelity 56-kg spacecraft replica representative of a modern low Earth orbit (LEO) satellite subjected to a laboratory hypervelocity impact test to understand fragmentation events and to support updates to satellite breakup models and size estimation models. Utilizing the vast population of fragments from DebriSat and prior laboratory impact experiments, the ODPO has focused on acquiring bidirectional reflectance distribution function (BRDF) data to characterize targets in the laboratory, thus removing aspect angle dependencies. Additionally, the DebriSat project has provided improved processes for measuring size via image acquisition, such that a true fragment size can be directly compared to the derived size using the OSEM. The team continues to assess BRDFs and use spectral measurement data to investigate the parameters used in the OSEM, specifically magnitudes, albedo, and phase functions.

Heather Cowardin↗

Characterization of Meteorites by Focused Ion Beam Sectioning: Recent Applications to CAIs and Primitive Meteorite Matrices

Focused ion beam (FIB) sectioning has revolutionized preparation of meteorite samples for characterization by analytical transmission electron microscopy (TEM) and other techniques. Although FIB is not "non-destructive" in the purest sense, each extracted section amounts to no more than nanograms (approximately 500 cubic microns) removed intact from locations precisely controlled by SEM imaging and analysis. Physical alteration of surrounding material by ion damage, fracture or sputter contamination effects is localized to within a few micrometers around the lift-out point. This leaves adjacent material intact for coordinate geochemical analysis by SIMS, microdrill extraction/TIMS and other techniques. After lift out, FIB sections can be quantitatively analyzed by electron microprobe prior to final thinning, synchrotron x-ray techniques, and by the full range of state-of-the-art analytical field-emission scanning transmission electron microscope (FE-STEM) techniques once thinning is complete. Multiple meteorite studies supported by FIB/FE-STEM are currently underway at NASA-JSC, including coordinated analysis of refractory phase assemblages in CAIs and fine-grained matrices in carbonaceous chondrites. FIB sectioning of CAIs has uncovered epitaxial and other overgrowth relations between corundum-hibonite-spinel consistent with hibonite preceding corundum and/or spinel in non-equilibrium condensation sequences at combinations of higher gas pressures, dust-gas enrichments or significant nebular transport. For all of these cases, the ability of FIB to allow for coordination with spatially-associated isotopic data by SIMS provides immense value for constraining the formation scenarios of the particular CAI assemblage. For carbonaceous chondrites matrix material, FIB has allowed us to obtain intact continuous sections of the immediate outer surface of Murchison (CM2) after it has been experimentally ion processed to simulate solar wind space weathering. The surface amorphization and loss of OH produced by the irradiation provides important clues regarding space weathering on primitive asteroids such as the OSIRIS-Rex target 101955 Bennu.

Christoffersen, Roy↗

Global ocean tide mapping using TOPEX/Poseidon altimetry

The investigation's main goals are to produce accurate tidal maps of the main diurnal, semidiurnal, and long-period tidal components in the world's deep oceans. This will be done by the application of statistical estimation techniques to long time series of altimeter data provided by the TOPEX/POSEIDON mission, with additional information provided by satellite tracking data. In the prelaunch phase, we will use in our simulations and preliminary work data supplied by previous oceanographic missions, such as Seasat and Geosat. These results will be of scientific interest in themselves. The investigation will also be concerned with the estimation of new values, and their uncertainties, for tidal currents and for the physical parameters appearing in the Laplace tidal equations, such as bottom friction coefficients and eddy viscosity coefficients. This will be done by incorporating the altimetry-derived charts of vertical tides as boundary conditions in the integration of those equations. The methodology of the tidal representation will include the use of appropriate series expansions such as ocean-basin normal modes and spherical harmonics. The results of the investigation will be space-determined tidal models of coverage and accuracy superior to that of the present numerical models of the ocean tides, with the concomitant benefits to oceanography and associated disciplinary fields.

Sanchez, Braulio V.↗

Symmetry, Statistics and Structure in MHD Turbulence

Here, we examine homogeneous MHD turbulence in terms of truncated Fourier series. The ideal MHD equations and the associated statistical theory of absolute equilibrium ensembles are symmetric under P, C and T. However, the presence of invariant helicities, which are pseudoscalars under P and C, dynamically breaks this symmetry. This occurs because the surface of constant energy in phase space has disjoint parts, called components: while ensemble averages are taken over all components, a dynamical phase trajectory is confined to only one component. As the Birkhoff-Khinchin theorem tells us, ideal MHD turbulence is thus non-ergodic. This non-ergodicity manifests itself in low-wave number Fourier modes that have large mean values (while absolute ensemble theory predicts mean values of zero). Therefore, we have coherent structure in ideal MHD turbulence. The level of non-ergodicity and amount of energy contained in the associated coherent structure depends on the values of the helicities, as well as on the presence, or not, of a mean magnetic field and/or overall rotation. In addition to the well known cross and magnetic helicities, we also present a new invariant, which we call the parallel helicity, since it occurs when mean field and rotation axis are aligned. The question of applicability of these results to real (i.e., dissipative) MHD turbulence is also examined. Several long-time numerical simulations on a 64(exp 3) grid are given as examples. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection of these results with inverse spectral cascades, selective decay, and magnetic dynamos is also discussed.

Shebalin, John V.↗