Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Momentum Imaging”

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 361 records · Page 20

Hybridized Exciton-Photon-Phonon States in a Transition Metal Dichalcogenide van der Waals Heterostructure Microcavity

Excitons in atomically thin transition-metal dichalcogenides (TMDs) have been established as an attractive platform to explore polaritonic physics, owing to their enormous binding energies and giant oscillator strength. Basic spectral features of exciton polaritons in TMD microcavities, thus far, were conventionally explained via two-coupled-oscillator models. This ignores, however, the impact of phonons on the polariton energy structure. Here we establish and quantify the threefold coupling between excitons, cavity photons, and phonons. For this purpose, we employ energy-momentum-resolved photoluminescence and spatially resolved coherent two-dimensional spectroscopy to investigate the spectral properties of a high-quality-factor microcavity with an embedded WSe 2 van der Waals heterostructure at room temperature. Our approach reveals a rich multibranch structure which thus far has not been captured in previous experiments. Simulation of the data reveals hybridized exciton-photon-phonon states, providing new physical insight into the exciton polariton system based on layered TMDs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Plasma Actuators for Turbomachinery Flow Control

This report is Part I of the final report of NASA Cooperative Agreement contract no. NNX07AC02A. The period of performance was January 1, 2007 to December 31, 2010. This report includes the project summary, a list of publications and reprints of the publications that appeared in archival journals. Part II of the final report includes a Ph.D. dissertation and is published separately as NASA/CR-2012-2172655. The research performed under this project was focused on the operation of surface dielectric barrier discharge (DBD) devices driven by high voltage, nanosecond scale pulses plus constant or time varying bias voltages. The main interest was in momentum production and the range of voltages applied eliminated significant heating effects. The approach was experimental supplemented by computational modeling. All the experiments were conducted at Princeton University. The project provided comprehensive understanding of the associated physical phenomena. Limitations on the performance of the devices for the generation of high velocity surface jets were established and various means for overcoming those limitations were proposed and tested. The major limitations included the maximum velocity limit of the jet due to electrical breakdown in air and across the dielectric, the occurrence of backward breakdown during the short pulse causing reverse thrust, the buildup of surface charge in the dielectric offsetting the forward driving potential of the bias voltage, and the interaction of the surface jet with the surface through viscous losses. It was also noted that the best performance occurred when the nanosecond pulse and the bias voltage were of opposite sign. Solutions include the development of partially conducting surface coatings, the development of a semiconductor diode inlaid surface material to suppress the backward breakdown. Extension to long discharge channels was studied and a new ozone imaging method developed for more quantitative determination of surface jet properties.

Miles, Richard, B↗

Lateral Shock of the R Aquarii Jet

The R Aqr jet was observed with the VLA B-configuration at two epochs separated by approximately 13.2 yr. Comparison of the resulting 6 cm continuum images show that the radio jet has undergone a lateral counterclockwise rotation of approximately 6 deg-12 deg on the plane of the sky. The model of jet parcels on independent trajectories is difficult to reconcile with these observations and leads us to consider a path-oriented jet (i.e., younger parcels follow the same path as older parcels). Comparison of the most recent radio image with a nearly contemporaneous HST/FOC ultraviolet image at approximately 2330 Angstroms suggests that the ultraviolet emission lies along the leading side of the rotating radio jet. In conjunction with a proper motion analysis of the jet material that yields empirical space-velocity and resulting acceleration-magnitude relationships as a function of distance from the central source, we evaluate the observational results in terms of a schematic model in which the jet emission consists of plane-parallel isothermal shocks along the leading edge of rotation. In such a radiating shock, the ultraviolet-emitting region is consistent with the adiabatic region in the form of a high-temperature, low-density sheath that surrounds the cooled postshock radio-emitting region. Within the context of the schematic model, we obtain the temperatures, densities, and pressures within the preshock, adiabatic, and postshock regions as a function of distance from the central source; the physical parameters so derived compare favorably to previously published estimates. We obtain a total jet mass of 3.1 x 10(exp -5) solar mass and an age of approximately 115 yr. We evaluate the model in the context of its density-boundary condition, its applicability to an episodic or quasi-continuous jet, and angular momentum considerations.

Hollis, J. M.↗

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗

Double Asteroid Redirection Test (DART) Mission

NASA’s Double Asteroid Redirection Test (DART) mission was humanity’s first attempt to move a celestial body, demonstrating the capability to perform a kinetic impact on a planetary defense–relevant sized asteroid. DART was part of the international collaboration known as the Asteroid Impact & Deflection Assessment (AIDA), involving NASA, the European Space Agency (ESA), the Agenzia Spaziale Italiana (ASI), and scientists around the world. DART was a key step to demonstrating preparedness to respond to planetary defense scenarios, and it provides a crucial data point for likely outcomes. Near-Earth objects (NEOs) greater than 140 m in size are of particular interest to planetary defense because they have the potential to cause significant damage if they were to impact Earth, and also because they are difficult to detect, with less than 50% of the predicted population discovered as of 2022 (National Academies Press, 2022). With an appropriately sized spacecraft and enough warning (typically many years to decades), a kinetic impact can slightly alter the orbit of an asteroid in a way that, over time, prevents the asteroid from colliding with Earth in the future. DART’s target was Dimorphos, the smaller (~150-m-diameter) member of the binary asteroid system (65803) Didymos, which is a near-Earth, potentially hazardous, and well-characterized asteroid system. By simply observing changes to the system after impact and comparing them with a pre-impact reference, it was possible to use ground-based telescopes to observe the deflection in the orbit of Dimorphos after impact. Developed and operated by the Johns Hopkins University Applied Physics Laboratory (APL), the mission entered formulation in 2015 after multiple years of concept development. The project was administered according to NPR 7120.5, with technical oversight and funding through the Planetary Missions Program Office (PMPO) at Marshall Space Flight Center (MSFC) and overall support as a directed mission from NASA’s Planetary Defense Coordination Office (PDCO). The DART spacecraft hosted a singular payload, the Didymos Reconnaissance and Asteroid Camera for Optical navigation (DRACO), and a deployable CubeSat contributed by ASI named the Light Italian CubeSat for Imaging of Asteroids (LICIACube). On 11 September 2022, DART deployed LICIACube, which subsequently followed the DART spacecraft at a safe distance and observed the immediate aftermath of the DART impact. DART was designed to autonomously detect, navigate to, and impact Dimorphos. This autonomous design was chosen to maximize the probability of impact, since commanding from the ground could result in course corrections arriving too late. On the day of impact, 26 September 2022, the spacecraft’s autonomous systems successfully detected and locked on to Dimorphos, impacting its surface within 2 m of the center of the illuminated figure (Jensenius et al., 2023). No human intervention was required for a successful impact, demonstrating that humanity possesses the technology to perform a kinetic impact. Within 2 weeks of impact, it was clear that the orbit of Dimorphos had been significantly altered. On 11 October 2022, NASA Administrator Bill Nelson announced that the new orbital period of Dimorphos was shortened by approximately 32 ± 2 min, from 11 h and 55 min before impact to 11 h and 23 min after impact. With additional observations over the following months, the accuracy of this measurement improved to a –33.24 min ± 1.4 s orbital period change (Naidu et al., 2023; Scheirich et al., 2023), and Beta (β), the momentum transfer enhancement parameter, was reported to be 3.6 (Cheng et al., 2023). Subsequent studies examined the details of DART’s impact site, modeled the impact event, investigated the ejecta produced, and analyzed the dynamics of the Didymos system. These combined results clearly demonstrate that the project met all Level 1 mission requirements.

asteroid↗

Imaging X-Ray Polarimetry Explorer Mission Attitude Determination and Control Concept

The goal of the Imaging X-Ray Polarimetry Explorer (IXPE) Mission is to expand understanding of high-energy astrophysical processes and sources, in support of NASA's first science objective in Astrophysics: "Discover how the universe works." X-ray polarimetry is the focus of the IXPE science mission. Polarimetry uniquely probes physical anisotropies-ordered magnetic fields, aspheric matter distributions, or general relativistic coupling to black-hole spin-that are not otherwise measurable. The IXPE Observatory consists of Spacecraft and Payload modules. The Payload includes three polarization sensitive, X-ray detector units (DU), each paired with its corresponding grazing incidence mirror module assemblies (MMA). A deployable boom provides the correct separation (focal length) between the DUs and MMAs. These Payload elements are supported by the IXPE Spacecraft. A star tracker is mounted directly with the deployed Payload to minimize alignment errors between the star tracker line of sight (LoS) and Payload LoS. Stringent pointing requirements coupled with a flexible structure and a non-collocated attitude sensor-actuator configuration requires a thorough analysis of control-structure interactions. A non-minimum phase notch filter supports robust control loop stability margins. This paper summarizes the IXPE mission science objectives and Observatory concepts, and then it describes IXPE attitude determination and control implementation. IXPE LoS pointing accuracy, control loop stability, and angular momentum management are discussed.

Attitude Determination and Contro↗

Investigation of Unsteady Flow Behavior in Transonic Compressor Rotors with LES and PIV Measurements

In the present study, unsteady flow behavior in a modern transonic axial compressor rotor is studied in detail with large eddy simulation (LES) and particle image velocimetry (PIV). The main purpose of the study is to advance the current understanding of the flow field near the blade tip in an axial transonic compressor rotor near the stall and peak-efficiency conditions. Flow interaction between the tip leakage vortex and the passage shock is inherently unsteady in a transonic compressor. Casing-mounted unsteady pressure transducers have been widely applied to investigate steady and unsteady flow behavior near the casing. Although many aspects of flow have been revealed, flow structures below the casing cannot be studied with casing-mounted pressure transducers. In the present study, unsteady velocity fields are measured with a PIV system and the measured unsteady flow fields are compared with LES simulations. The currently applied PIV measurements indicate that the flow near the tip region is not steady even at the design condition. This self-induced unsteadiness increases significantly as the compressor rotor operates near the stall condition. Measured data from PIV show that the tip clearance vortex oscillates substantially near stall. The calculated unsteady characteristics of the flow from LES agree well with the PIV measurements. Calculated unsteady flow fields show that the formation of the tip clearance vortex is intermittent and the concept of vortex breakdown from steady flow analysis does not seem to apply in the current flow field. Fluid with low momentum near the pressure side of the blade close to the leading edge periodically spills over into the adjacent blade passage. The present study indicates that stall inception is heavily dependent on unsteady behavior of the flow field near the leading edge of the blade tip section for the present transonic compressor rotor.

Hah, Chunill↗

Concentration Measurements in Self-Excited, Momentum-Dominated Helium Jets

Flow structure of momentum-dominated pure helium jets discharged vertically into ambient air was investigated using high-speed rainbow schlieren deflectometry (RSD) technique. Effects of the operating parameters, i.e., Reynolds number (Re) and Richardson number (Ri), on the oscillatory behavior of the flow were examined over a range of experimental conditions. To seek the individual effect of these parameters, one of them was fixed and the other was varied with certain constraints. Measurements revealed highly periodic oscillations in the laminar region as well as high regularity in transition and turbulent regions. Maximum spectral power profiles at different axial locations indicated the oscillation amplitude increasing until the breakdown of the jet in the turbulent regime. The transition from the laminar to turbulent flow was also investigated. Fast Fourier transform analysis performed in the transition regime showed that the flow oscillates at a unique frequency, which was the same in the upstream laminar flow region. Measured deflection angle data were used in Abel inversion algorithm to construct the helium concentration fields. Instantaneous helium concentration contours revealed changes in the flow structure and evolution of vortical structures during an oscillation cycle. Temporal evolution plots of helium concentration at different axial location showed repeatable oscillations at all axial and radial locations up to the turbulent regime. A cross-correlation technique, applied to find the spatial displacements of the vortical structures, provided correlation coefficient peaks between consecutive schlieren images. Results show that the vortical structure convected and accelerated only in the axial direction.

Yildirim, Bekir Sedat↗

Attenuation of Gas Turbulence by a Nearly Stationary Dispersion of Fine Particles

Turbulence attenuation by greater than a factor of two has been observed in many practical gas flows carrying volume fractions as small as 0.01% of dispersed particles. Particles which cause such attenuation usually are smaller than the smallest scales of the turbulence and have time constants 5 to 10 times greater than the time scale of a typical turbulent eddy. That is, strongly attenuating particles usually have Stokes numbers in the range of 5 to 10, indicating that they do not respond to the turbulent fluctuations, but instead just fall through the flow responding only to the mean flow. There are two mechanisms by which free falling particles may attenuate turbulence. First, the unresponsive particles act as a drag on the turbulent eddies, passing energy from the turbulent eddies to the small scale wakes of the particles where it is quickly dissipated by viscosity. The second mechanism is more complicated. Particles falling under gravity convert gravitational potential energy to turbulent velocity fluctuations. If the particles are large, this mechanism increases the overall turbulence level. However, with moderate size particles, the small scale turbulence generated apparently distorts the turbulent eddies leading to more rapid dissipation. Unfortunately, this conclusion is supported only by circumstantial evidence to date. The objectives of the experiment are to use microgravity to separate the two mechanisms. A region of nearly-isotropic decaying turbulence with zero mean flow will be formed in a box in the microgravity environment. Different sets of particles with Stokes numbers in the range of 2 to 20 will be dispersed in the flow. With zero gravity and no mean fluid velocity the particles will have zero mean velocity. With the large Stokes numbers, the fluctuating velocities will also be small. Therefore, the only attenuation mechanism will be the direct action of the particles on the turbulence. Control experiments will also be done in which the particles fall through the measurement volume. Measurements will be acquired using a high resolution image velocimetry (PIV) system being developed specifically for work in particle-laden flows. The measurements will include the decay of the turbulence kinetic energy under various particle loadings. The spatial spectra of the turbulence will also be measured. In a second set of experiments, the interaction of a single eddy with a collection of nearly stationary particles will be examined. The eddy will be a vortex ring emitted by a jet pulse through an orifice. The distortion of the vortex under the influence of the particles will be examined to gain a better understanding of how fine particles can cause such large reductions in turbulence levels. This experiment could not be conducted in terrestrial gravity because the high particle velocities would overwhelm the relatively low speed motion of the vortex ring. This experimental program is just getting underway. The initial challenge is to build a closed facility containing reasonably homogeneous and isotropic turbulence with zero mean velocity. Our approach is to use a set of synthetic jets mounted on the periphery of a transparent plexiglass box to create the turbulence. A synthetic jet is a plenum chamber with an orifice open to the volume of interest. The volume of the chamber fluctuates periodically so alternately a jet is ejected from the volume or flow is drawn back in as a sink. The asymmetry of this situation results in a net transport of momentum and kinetic energy into the volume of interest. The present apparatus includes eight synthetic jets each powered independently by a six inch loudspeaker. The synthetic jets discharge through ejector tubes to increase the scale of the turbulence. Construction of the apparatus is now complete and preliminary flow visualization studies have been conducted. The PIV system is also under development. A compact dual-pulse YAG laser has been acquired as the light source and special software is under development to allow simultaneous measurements of both the particle phase and the fluid phase (marked by fine tracers).

Paris, A. D.↗

AOCS Performance and Stability Validation for a 160-m Solar Sail with Control-Structure Interactions

Future solar sail missions, such as NASA's Solar Polar Imager Vision, will require sails with dimensions on the order of 50-500 m. We are examining a square sail design with moving mass (trim control mass, TCM) and quadrant rotation primary actuators plus pulsed plasma thrusters (PPTs) at the mast tips for backup attitude control. Quadrant rotation is achieved via roll stabilizer bars (RSB) at the mast tips. At these sizes, given the gossamer nature of the sail supporting structures, flexible modes may be low enough to interact with the control system, especially as these actuators are located on the flexible structure itself and not on the rigid core. This paper develops a practical analysis of the flexible interactions using state-space systems and modal data from finite element models of the system. Torsion and bending of the masts during maneuvers could significantly affect the function of the actuators while activation of the membrane modes could adversely affect the thrust vector direction and magnitude. Analysis of the RSB and TCM dynamics for developing high-fidelity simulations is included. For control analysis of the flexible system, standard finite-element models of the flexible sail body are loaded and the modal data is used to create a modal coordinate state-space system. Key parameters include which modes to include, which nodes are of interest for force inputs and displacement outputs, connecting nodes through which external forces and torques are applied from the flex body to the core, any nominal momentum in the system, and any steady rates. The system is linearized about the nominal attitude and rate. The state-space plant can then be analyzed with a state-space controller, and Bode, Nyquist, step and impulse responses generated. The approach is general for any rigid core with a flexible appendage. This paper develops a compensator for a simple two-mass flex system and extrapolates the results to the solar sail. A finite element model of the 20 m solar sail by ATK Space Systems, recently validated in ground tests, is used to demonstrate the sail analysis approach.

Wie, Bong↗

Lunar Reconnaissance Orbiter (LRO) Guidance, Navigation and Control (GN&C) Overview

The National Aeronautics and Space Administration s (NASA) Lunar Reconnaissance Orbiter (LRO) launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle and into a direct insertion trajectory to the oon. LRO, which was designed, built, and operated by the NASA Goddard Space Flight Center in Greenbelt, MD, is gathering crucial data on the lunar environment that will help astronauts prepare for long-duration lunar expeditions. The mission has a nominal life of 1 year as its seven instruments find safe landing sites, locate potential resources, characterize the radiation environment, and test new technology. To date, LRO has been operating well within the bounds of its requirements and has been collecting excellent science data images taken from the LRO Camera Narrow Angle Camera of the Apollo landing sites appeared on cable news networks. A significant amount of information on LRO s science instruments is provided at the LRO mission webpage. LRO s Guidance, Navigation and Control (GN&C) subsystem is made up of an onboard attitude control system (ACS) and a hardware suite of sensors and actuators. The LRO onboard ACS is a collection of algorithms based on high level and derived requirements, and reflect the science and operational events throughout the mission lifetime. The primary control mode is the Observing mode, which maintains the lunar pointing orientation and any offset pointing from this baseline. It is within this mode that all science instrument calibrations, slews and science data is collected. Because of a high accuracy requirement for knowledge and pointing, the Observing mode makes use of star tracker (ST) measurement data to determine an instantaneous attitude pointing. But even the star trackers alone do not meet the tight requirements, so a six-state Kalman Filter is employed to improve the noisy measurement data. The Observing mode obtains its rate information from an inertial reference unit (IRU) and in the event of an IRU failure, the rate data is be derived from the star tracker, but with degraded pointing performance. The Delta-V control mode responsibility is to maintain attitude pointing during the cruise trajectory, insertion burns and lunar orbit maintenance by adjustments made to the spacecraft s velocity magnitude and vector direction. The ACS also provides for a thruster based system momentum management algorithm (known as Delta-H) to maintain the system and wheel momentum to within acceptable levels. In the event an anomaly causes the LRO spacecraft to lose the ability to maintain its current attitude pointing, a Sun Safe mode is included in the ACS for the purpose of providing a known power and thermally safe coarse inertial sun attitude for an indefinite period of time, within the manageable limits of the reaction wheels. The Sun Safe mode is also the initial spacecraft control mode off of the launch vehicle and provides for a means to null tip-off rates immediately after separation. The nominal configuration is to use the IRU for rate information in the controller. In the event of a gyro failure a gyroless control mode was developed that computes rate information from the CSS data.

Garrick, Joseph↗

In-Flight Guidance, Navigation, and Control Performance Results for the GOES-16 Spacecraft

The Geostationary Operational Environmental Satellite-R Series (GOES-R), which launched in November 2016, is the first of the next generation geostationary weather satellites. GOES-R provides 4 times the resolution, 5 times the observation rate, and 3 times the number of spectral bands for Earth observations compared with its predecessor spacecraft. Additionally, Earth relative and Sun-relative pointing and pointing stability requirements are maintained throughout reaction wheel desaturation events and station keeping activities, allowing GOES-R to provide continuous Earth and sun observations. This paper reviews the pointing control, pointing stability, attitude knowledge, and orbit knowledge requirements necessary to realize the ambitious Image Navigation and Registration (INR) objectives of GOES-R. This paper presents a comparison between low-frequency on-orbit pointing results and simulation predictions for both the Earth Pointed Platform (EPP) and Sun Pointed Platform (SPP). Results indicate excellent agreement between simulation predictions and observed on-orbit performance, and compliance with pointing performance requirements. The EPP instrument suite includes 6 seismic accelerometers sampled at 2 KHz, allowing in-flight verification of jitter responses and comparison back to simulation predictions. This paper presents flight results of acceleration, shock response spectrum (SRS), and instrument line of sight responses for various operational scenarios and instrument observation modes. The results demonstrate the effectiveness of the dual-isolation approach employed on GOES-R. The spacecraft provides attitude and rate data to the primary Earth-observing instrument at 100 Hz, which are used to adjust instrument scanning. The data must meet accuracy and latency numbers defined by the Integrated Rate Error (IRE) requirements. This paper discusses the on-orbit IRE results, showing compliance to these requirements with margin. During the spacecraft checkout period, IRE disturbances were observed and subsequently attributed to thermal control of the Inertial Measurement Unit (IMU) mounting interface. Adjustments of IMU thermal control and the resulting improvements in IRE are presented. Orbit knowledge represents the final element of INR performance. Extremely accurate orbital position is achieved by GPS navigation at Geosynchronous Earth Orbit (GEO). On-orbit performance results are shown demonstrating compliance with the stringent orbit position accuracy requirements of GOES-R, including during station keeping activities and momentum desaturation events. As we show in this paper, the on-orbit performance of the GNC design provides the necessary capabilities to achieve GOES-R mission objectives.

Chapel, Jim↗

Impact-Seismic Investigations of the InSight Mission

Impact investigations will be an important aspect of the InSight mission. One of the scientific goals of the mission is a measurement of the current impact rate at Mars. Impacts will additionally inform the major goal of investigating the interior structure of Mars. In this paper, we review the current state of knowledge about seismic signals from impacts on the Earth, Moon, and laboratory experiments. We describe the generalized physical models that can be used to explain these signals. A discussion of the appropriate source time function for impacts is presented, along with spectral characteristics including the cutoff frequency and its dependence on impact momentum. Estimates of the seismic efficiency (ratio between seismic and impact energies) vary widely. Our preferred value for the seismic efficiency at Mars is 5 × 10−4, which we recommend using until we can measure it during the InSight mission, when seismic moments are not used directly. Effects of the material properties at the impact point and at the seismometer location are considered. We also discuss the processes by which airbursts and acoustic waves emanate from bolides, and the feasibility of detecting such signals. We then consider the case of impacts on Mars. A review is given of the current knowledge of present-day cratering on Mars: the current impact rate, characteristics of those impactors such as velocity and directions, and the morphologies of the craters those impactors create. Several methods of scaling crater size to impact energy are presented. The Martian atmosphere, although thin, will cause fragmentation of impactors, with implications for the resulting seismic signals. We also benchmark several different seismic modeling codes to be used in analysis of impact detections, and those codes are used to explore the seismic amplitude of impactinduced signals as a function of distance from the impact site. We predict a measurement of the current impact flux will be possible within the timeframe of the prime mission (one Mars year) with the detection of ∼ a few to several tens of impacts. However, the error bars on these predictions are large. Specific to the InSight mission, we list discriminators of seismic signals from impacts that will be used to distinguish them from marsquakes. We describe the role of the InSight Impacts Science Theme Group during mission operations, including a plan for possible night-time meteor imaging. The impacts detected by these methods during the InSight mission will be used to improve interior structure models, measure the seismic efficiency, and calculate the size frequency distribution of current impacts.

Ingrid Daubar↗

Extreme Longitudinal Compression of Optimized Beams for MEV Ultrafast Electron Diffraction (Final Technical Report)

We worked out the design of a high repetition rate MeV energy ultrafast electron diffraction instrument based on the existing Cornell photoinjector, which can readily be applied to the presented findings. This example is a blueprint of other similarly arranged UED setups. Using particle tracking simulations in conjunction with multiobjective genetic algorithm optimization, we explored the smallest bunch lengths, emittance, and probe spot sizes achievable. As two limits, we defined stroboscopic conditions (with single electrons per pulse) and operation with 10 5 electrons per bunch which may be suitable for single-shot diffraction images. In the stroboscopic case, the flexibility provided by the many cavity bunching and acceleration allows for longitudinal phase space linearization without a higher harmonic field, providing sub-fs bunch lengths at the sample. Given low emittance photoemission conditions, these small bunch lengths can be maintained with probe transverse sizes at the single micron (1 μm) scale and below. In the case of 10 5 electrons per pulse, we simulated state-of-the-art 5D brightness conditions: rms bunch lengths of 10 fs with 3-nm normalized emittances, while permitting repetition rates as high as 1.3 GHz. We showed that in conjunction with collimating apertures, a novel focusing scheme achieves very high-quality emittance compensation for the central core of the beam composing 40% of particles, for a resulting beam size of 5 μm (rms). Finally, to aid in the design of new SRF-based ultrafast electron diffraction machines, we simulated the trade-off between the number of cavities used and achievable bunch length and emittance. In the longitudinal dimension, we made use of the fact that MeV UED requires much lower energy than the 15-MeV maxi mum energy of Cornell’s CBETA injector, and we may therefore use several of the SRF cavities for bunch length compression. In practice, we used a genetic optimization algorithm to choose the phases and amplitudes of the cavities appropriately for optimal bunching. In the zero space charge case, we found that bunching and acceleration are distributed across the six cavities in a way that produces a linearizing effect. And we showed that the ultimate bunch length can be limited by time-of-flight differences arising from transverse size and transverse momentum spread. The space charge code developed and used for this development is now permanent part of the Bmad accelerator simulation code and has already contributed to other developments, e.g., for the EIC electron cooler design.

43 PARTICLE ACCELERATORS↗

Three-dimensional ballistocardiography and respiratory motion in sustained microgravity

BACKGROUND: We measured the three-dimensional ballistocardiogram (BCG) in a free-floating subject in sustained microgravity during spaceflight to test the usefulness of such measurements for future non-invasive monitoring of cardiac function, and to examine the effects of respiratory movement on the BCG in three axes. METHODS: Acceleration was measured using a three-axis accelerometer fastened to the lumbar region of the subject while simultaneous recordings of ECG, and respiratory motion via impedance plethysmography were also made. Data were recorded during a 146-s period of inactivity on the part of the subject during which time there was no contact with the spacecraft. RESULTS: Total body motion due to respiratory activity was consistent with that calculated from the known action of the diaphragm and conservation of momentum. The accelerations due to cardiac activity, ensemble averaged over the R-R interval, were greatest along the head-to-foot axis. Maximum amplitude of the HIJK complex of the BCG generated by ventricular ejection was greatest in the head to foot axis (approximately 70 x 10(-3) m x s(-2)), but there were also substantial accelerations along the dorsoventral axis of up to 43 10(-3) m x s(-2), that are not measured interrestrial two-dimensional studies. The amplitude of the BCG was strongly affected by lung volume, with accelerations being reduced 50 to 70% between end-inspiration and end-expiration. CONCLUSIONS: These data suggest a greatly reduced transmission of the cardiac motion to the body at end-expiration (FRC) than at higher lung volumes. The BCG might be further developed as a non-invasive means of monitoring parameters such as stroke volume in microgravity.

Flight Experiment↗

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗

Combined Experimental/Numerical Study of the Soot Formation Process in a Gasoline Direct-Injection Spray in the Presence of Laser-Induced Plasma Ignition

Here, combustion issued from an eight-hole, direct-injection spray was experimentally studied in a constant-volume pre-burn combustion vessel using simultaneous high-speed diffused back-illumination extinction imaging (DBIEI) and OH* chemiluminescence. DBIEI has been employed to observe the liquid-phase of the spray and to quantitatively investigate the soot formation and oxidation taking place during combustion. The fuel-air mixture was ignited with a plasma induced by a single-shot Nd:YAG laser, permitting precise control of the ignition location in space and time. OH* chemiluminescence was used to track the high-temperature ignition and flame. The study showed that increasing the delay between the end of injection and ignition drastically reduces soot formation without necessarily compromising combustion efficiency. For long delays between the end of injection and ignition (1.9 ms) soot formation was eliminated in the main downstream charge of the fuel spray. However, poorly atomized and large droplets formed at the end of injection (dribble) eventually do form soot near the injector even when none is formed in the main charge. The quantitative soot measurements for these spray and ignition scenarios, resolved in time and space, represents a significant new achievement. Reynolds-averaged Navier-Stokes (RANS) simulations were performed to assess spray mixing and combustion. An analysis of the predicted fuel-air mixture in key regions, defined based upon experimental observations, was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing. The mixture analysis indicates that soot production can be avoided if the flame propagates into regions where the equivalence ratio (Φ) is already below 2. Reactive RANS simulations have also been performed, but with a poor match against the experiment, as the flame speed and heat-release rate are largely over estimated. This modeling weakness appears related to a very high level of turbulent viscosity predicted for the high-momentum spray in the RANS simulations, which is an important consideration for modeling ignition and flame propagation in mixtures immediately created by the spray.

42 ENGINEERING↗

Efficiency Studies of Fast Neutron Tracking Using MCNP

Fast neutron identification and spectroscopy is of great interest to nuclear physics experiments. Using the neutron elastic scattering, the fast neutron momentum can be measured. Wang and Morris introduced the theoretical concept that the initial fast neutron momentum can be derived from up to three consecutive elastic collisions between the neutron and the target, including the information of two consecutive recoil ion tracks and the vertex position of the third collision or two consecutive elastic collisions with the timing information. Here, we also include the additional possibility of measuring the deposited energies from the recoil ions. In this paper, we simulate the neutron elastic scattering using the Monte Carlo N-Particle Transport Code (MCNP) and study the corresponding neutron detection and tracking efficiency. The corresponding efficiency and the scattering distances are simulated with different target materials, especially natural silicon (92.23%28Si, 4.67%29Si, and 3.1%30Si) and helium-4 (4He). The timing of collision and the recoil ion energy are also investigated, which are important characters for the detector design. We also calculate the ion traveling range for different energies using the software, “The Stopping and Range of Ions in Matter (SRIM)”, showing that the ion track can be most conveniently observed in 4He unless sub-micron spatial resolution can be obtained in silicon.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗