Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “object oriented modeling”

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 595 records · Page 33

Acquisition of Long-Duration, Low-Gravity Slosh Data Utilizing Existing ISS Equipment (SPHERES) for Calibration of CFD Models of Coupled Fluid-Vehicle Behavior

Accurate prediction of coupled fluid slosh and launch vehicle or spacecraft dynamics (e.g., nutation/precessional movement about various axes, attitude changes, ect.) requires Computational Fluid Dynamics (CFD) models calibrated with low-gravity, long duration slosh data. Recently completed investigations of reduced gravity slosh behavior have demonstrated the limitations of utilizing parabolic flights on specialized aircraft with respect to the specific objectives of the experiments. Although valuable data was collected, the benefits of longer duration low-gravity environments were clearly established. The proposed research provides the first data set from long duration tests in zero gravity that can be directly used to benchmark CFD models, including the interaction between the sloshing fluid and the tank/vehicle dynamics. To explore the coupling of liquid slosh with the motion of an unconstrained tank in microgravity, NASA's Kennedy Space Center, Launch Services Program has teamed up with the Florida Institute of Technology (FIT), Massachusetts Institute of Technology (MIT) and the Office of the Chief Technologist (OCT) to perform a series of slosh dynamics experiments on the International Space Station using the SPHERES platform. The Synchronized Position Hold Engage Reorient Experimental Satellites (SPHERES) testbed provides a unique, free-floating instrumented platform on ISS that can be utilized in a manner that would solve many of the limitations of the current knowledge related to propellant slosh dynamics on launch vehicle and spacecraft fuel tanks. The six degree of freedom (6-DOF) motion of the SPHERES free-flyer is controlled by an array of cold-flow C02 thrusters, supplied from a built-in liquid C02 tank. These SPHERES can independently navigate and re-orient themselves within the ISS. The intent of this project is to design an externally mounted tank to be driven inside the ISS by a set of two SPHERES devices (Figure 1 ). The tank geometry simulates a launch vehicle upper stage propellant tank and the maneuvers replicate those of real vehicles. The design includes inertial sensors, data acquisition, image capture and data storage interfaces to the SPHERES VERTIGO computer system on board the flight article assembly. The design also includes mechanical and electronic interfaces to the existing SPHERES hardware, which include self-contained packages that can operate in conjunction with the existing SPHERES electronics.

Schallhorn, Paul↗

Acquisition of Long-Duration, Low-Gravity Slosh Data Utilizing Existing ISS Equipment (SPHERES) for Calibration of CFD Models of Coupled Fluid-Vehicle Behavior

Accurate prediction of coupled fluid slosh and launch vehicle or spacecraft dynamics (e.g., nutation/precessional movement about various axes, attitude changes, ect.) requires Computational Fluid Dynamics (CFD) models calibrated with low-gravity, long duration slosh data. Recently completed investigations of reduced gravity slosh behavior have demonstrated the limitations of utilizing parabolic flights on specialized aircraft with respect to the specific objectives of the experiments. Although valuable data was collected, the benefits of longer duration low-gravity environments were clearly established. The proposed research provides the first data set from long duration tests in zero gravity that can be directly used to benchmark CFD models, including the interaction between the sloshing fluid and the tank/vehicle dynamics. To explore the coupling of liquid slosh with the motion of an unconstrained tank in microgravity, NASA's Kennedy Space Center, Launch Services Program has teamed up with the Florida Institute of Technology (FIT), Massachusetts Institute of Technology (MIT) and the NASA Game Changing Development Program (GCD) to perform a series of slosh dynamics experiments on the International Space Station using the SPHERES platform. The Synchronized Position Hold Engage Reorient Experimental Satellites (SPHERES) testbed provides a unique, free-floating instrumented platform on ISS that can be utilized in a manner that would solve many of the limitations of the current knowledge related to propellant slosh dynamics on launch vehicle and spacecraft fuel tanks. The six degree of freedom (6-DOF) motion of the SPHERES free-flyer is controlled by an array of cold-flow C02 thrusters, supplied from a built-in liquid C02 tank. These SPHERES can independently navigate and re-orient themselves within the ISS. The intent of this project is to design an externally mounted tank to be driven inside the ISS by a set of two SPHERES devices (Figure 1). The tank geometry simulates a launch vehicle upper stage propellant tank and the maneuvers replicate those of real vehicles. The design includes inertial sensors, data acquisition, image capture and data storage interfaces to the SPHERES VERTIGO computer system on board the flight article assembly. The design also includes mechanical and electronic interfaces to the existing SPHERES hardware, which include self-contained packages that can operate in conjunction with the existing SPHERES electronics

Schallhorn, Paul↗

Perceived visual speed constrained by image segmentation

Little is known about how or where the visual system parses the visual scene into objects or surfaces. However, it is generally assumed that the segmentation and grouping of pieces of the image into discrete entities is due to 'later' processing stages, after the 'early' processing of the visual image by local mechanisms selective for attributes such as colour, orientation, depth, and motion. Speed perception is also thought to be mediated by early mechanisms tuned for speed. Here we show that manipulating the way in which an image is parsed changes the way in which local speed information is processed. Manipulations that cause multiple stimuli to appear as parts of a single patch degrade speed discrimination, whereas manipulations that perceptually divide a single large stimulus into parts improve discrimination. These results indicate that processes as early as speed perception may be constrained by the parsing of the visual image into discrete entities.

NASA Discipline Neuroscience↗

Structural Health Management of Damaged Aircraft Structures Using the Digital Twin Concept

The development of multidisciplinary integrated Structural Health Management (SHM) tools will enable accurate detection, and prognosis of damaged aircraft under normal and adverse conditions during flight. As part of the digital twin concept, methodologies are developed by using integrated multiphysics models, sensor information and input data from an in-service vehicle to mirror and predict the life of its corresponding physical twin. SHM tools are necessary for both damage diagnostics and prognostics for continued safe operation of damaged aircraft structures. The adverse conditions include loss of control caused by environmental factors, actuator and sensor faults or failures, and structural damage conditions. A major concern in these structures is the growth of undetected damage/cracks due to fatigue and low velocity foreign object impact that can reach a critical size during flight, resulting in loss of control of the aircraft. To avoid unstable, catastrophic propagation of damage during a flight, load levels must be maintained that are below a reduced load-carrying capacity for continued safe operation of an aircraft. Hence, a capability is needed for accurate real-time predictions of damage size and safe load carrying capacity for structures with complex damage configurations. In the present work, a procedure is developed that uses guided wave responses to interrogate damage. As the guided wave interacts with damage, the signal attenuates in some directions and reflects in others. This results in a difference in signal magnitude as well as phase shifts between signal responses for damaged and undamaged structures. Accurate estimation of damage size, location, and orientation is made by evaluating the cumulative signal responses at various pre-selected sensor locations using a genetic algorithm (GA) based optimization procedure. The damage size, location, and orientation is obtained by minimizing the difference between the reference responses and the responses obtained by wave propagation finite element analysis of different representative cracks, geometries, and sizes.

Seshadri, Banavara R.↗

Comparing viewer and array mental rotations in different planes

Participants imagined rotating either themselves or an array of objects that surrounded them. Their task was to report on the egocentric position of an item in the array following the imagined rotation. The dependent measures were response latency and number of errors committed. Past research has shown that self-rotation is easier than array rotation. However, we found that imagined egocentric rotations were as difficult to imagine as rotations of the environment when people performed imagined rotations in the midsagittal or coronal plane. The advantages of imagined self-rotations are specific to mental rotations performed in the transverse plane.

NASA Center ARC↗

Development and application of virtual reality for man/systems integration

While the graphical presentation of computer models signified a quantum leap over presentations limited to text and numbers, it still has the problem of presenting an interface barrier between the human user and the computer model. The user must learn a command language in order to orient themselves in the model. For example, to move left from the current viewpoint of the model, they might be required to type 'LEFT' at a keyboard. This command is fairly intuitive, but if the viewpoint moves far enough that there are no visual cues overlapping with the first view, the user does not know if the viewpoint has moved inches, feet, or miles to the left, or perhaps remained in the same position, but rotated to the left. Until the user becomes quite familiar with the interface language of the computer model presentation, they will be proned to lossing their bearings frequently. Even a highly skilled user will occasionally get lost in the model. A new approach to presenting type type of information is to directly interpret the user's body motions as the input language for determining what view to present. When the user's head turns 45 degrees to the left, the viewpoint should be rotated 45 degrees to the left. Since the head moves through several intermediate angles between the original view and the final one, several intermediate views should be presented, providing the user with a sense of continuity between the original view and the final one. Since the primary way a human physically interacts with their environment should monitor the movements of the user's hands and alter objects in the virtual model in a way consistent with the way an actual object would move when manipulated using the same hand movements. Since this approach to the man-computer interface closely models the same type of interface that humans have with the physical world, this type of interface is often called virtual reality, and the model is referred to as a virtual world. The task of this summer fellowship was to set up a virtual reality system at MSFC and begin applying it to some of the questions which concern scientists and engineers involved in space flight. A brief discussion of this work is presented.

Brown, Marcus↗

Development and testing of tip devices for horizontal axis wind turbines

A theoretical and field experimental program has been carried out to investigate the use of tip devices on horizontal axis wind turbine rotors. The objective was to improve performance by the reduction of tip losses. While power output can always be increased by a simple radial tip extension, such a modification also results in an increased gale load both because of the extra projected area and longer moment arm. Tip devices have the potential to increase power output without such a structural penalty. A vortex lattice computer model was used to optimize three basic tip configuration types for a 25 kW stall limited commercial wind turbine. The types were a change in tip planform, and a single-element and double-element nonplanar tip extension (winglets). A complete data acquisition system was developed which recorded three wind speed components, ambient pressure, temperature, and turbine output. The system operated unattended and could perform real-time processing of the data, displaying the measured power curve as data accumulated in either a bin sort mode or polynomial curve fit. Approximately 270 hr of perormance data were collected over a three-month period. The sampling interval was 2.4 sec; thrus over 400,000 raw data points were logged. Results for each of the three new tip devices, compared with the original tip, showed a small decrease (of the order of 1 kW) in power output over the measured range of wind speeds from cut-in at about 4 m/s to over 20 m/s, well into the stall limiting region. Changes in orientation and angle-of-attack of the winglets were not made. For aircraft wing tip devices, favorable tip shapes have been reported and it is likely that the tip devices tested in this program did not improve rotor performance because they were not optimally adjusted.

Gyatt, G. W.↗

Optomechanical Alignment of the Grating Wheel Mechanism for a Ground-based, Cryogenic, Near-Infrared Astronomy Instrument

We describe the population, optomechanical alignment, and alignment verification of near-infrared gratings on the grating wheel mechanism (GWM) for the Infrared Multi- Object Spectrometer (IRMOS). IRMOS is a cryogenic (80 K) facility instrument for the Mayall Telescope (3.8 m) at Kitt Peak National Observatory and a MEMS spectrometer concept demonstrator for NASA's James Webb Space Telescope. The IRMOS optics, bench, and mechanisms are predominantly made of Al 6061 -T651. The GWM consists of 13 planar diffraction gratings and one flat imaging mirror (58 x 57 mm), each mounted at a unique compound angle on a 31.8 cm diameter gear. The Al 6061 grating substrates are stress relieved for enhanced cryogenic performance and the optical surface is replicated from an off-the-shelf master. The imaging mirror is diamond turned and post-polished. The grating mechanism spans a projected diameter of approximately 48cm when fully assembled, utilizes several flexure designs throughout the system to accommodate thermal gradient situations, and is controlled using custom software with an off-the-shelf controller. Each optic is aligned in six degrees of freedom relative to the GWM coordinate system, which is defined relative to an optical alignment cube mounted at the center of the gear. The tip/tilt (Rx, Ry) orientation of a given grating is measured using the zero-order return from an autocollimating theodolite. Each optic's mount includes a one-piece shim located between the optic and the gear. The shim is machined to fine align each optic. We also describe alignment verification, where grating diffractive properties are compared to model predictions.

Gutkowski, Sharon M.↗

Alignment of the Grating Wheel Mechanism for a Ground-Based, Cryogenic, Near-Infrared Astronomy Instrument

We describe the population, optomechanical alignment, and alignment verification of near-infrared gratings on the grating wheel mechanism (GWM) for the Infrared Multi-Object Spectrometer (IRMOS). IRMOS is a cryogenic (80 K), principle investigator-class instrument for the 2.1 m and Mayall 3.8 m telescopes at Kitt Peak National Observatory, and a MEMS spectrometer concept demonstrator for the James Webb Space Telescope. The GWM consists of 13 planar diffraction gratings and one flat imaging mirror (58 x 57 mm), each mounted at a unique compound angle on a 32 cm diameter gear. The mechanism is predominantly made of Al 6061. The grating substrates are stress relieved for enhanced cryogenic performance. The optical surfaces are replicated from off-the-shelf masters. The imaging mirror is diamond turned. The GWM spans a projected diameter of approx. 48 cm when fully assembled, utilizes several flexure designs to accommodate potential thermal gradients, and is controlled using custom software with an off-the-shelf controller. Under ambient conditions, each grating is aligned in six degrees of freedom relative to a coordinate system that is referenced to an optical alignment cube mounted at the center of the gear. The local tip/tilt (Rx/Ry) orientation of a given grating is measured using the zero-order return from an autocollimating theodolite. The other degrees of freedom are measured using a two-axis cathetometer and rotary table. Each grating's mount includes a one-piece shim located between the optic and the gear. The shim is machined to fine align each grating. We verify ambient alignment by comparing grating difractive properties to model predictions.

Gutkowski, Sharon M.↗

Deformation and fracture of thin sheet aluminum-lithium alloys: The effect of cryogenic temperatures

The objective is to characterize the fracture behavior and to define the fracture mechanisms for new Al-Li-Cu alloys, with emphasis on the role of indium additions and cryogenic temperatures. Three alloys were investigated in rolled product form: 2090 baseline and 2090 + indium produced by Reynolds Metals, and commercial AA 2090-T81 produced by Alcoa. The experimental 2090 + In alloy exhibited increases in hardness and ultimate strength, but no change in tensile yield strength, compared to the baseline 2090 composition in the unstretched T6 condition. The reason for this behavior is not understood. Based on hardness and preliminary Kahn Tear fracture experiments, a nominally peak-aged condition was employed for detailed fracture studies. Crack initiation and growth fracture toughness were examined as a function of stress state and microstructure using J(delta a) methods applied to precracked compact tension specimens in the LT orientation. To date, J(delta a) experiments have been limited to 23 C. Alcoa 2090-T81 exhibited the highest toughness regardless of stress state. Fracture was accompanied by extensive delamination associated with high angle grain boundaries normal to the fatigue precrack surface and progressed microscopically by a transgranular shear mechanism. In contrast the two peak-aged Reynolds alloys had lower toughness and fracture was intersubgranular without substantial delamination. The influences of cryogenic temperature, microstructure, boundary precipitate structure, and deformation mode in governing the competing fracture mechanisms will be determined in future experiments. Results contribute to the development of predictive micromechanical models for fracture modes in Al-Li alloys, and to fracture resistant materials.

Wagner, John A.↗

HOTCFGM-2D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Bi-Directionally Components with Bi-Directionally Graded Microstructures

The objective of this two-year project was to develop and deliver to the NASA-Glenn Research Center a two-dimensional higher-order theory, and related computer codes, for the analysis and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, brisk blades). To satisfy this objective, two-dimensional version of the higher-order theory, HOTCFGM-2D, and four computer codes based on this theory, for the analysis and design of structural components functionally graded in the radial and circumferential directions were developed in the cylindrical coordinate system r-Theta-z. This version of the higher-order theory is a significant generalization of the one-dimensional theory, HOTCFGM-1D, developed during the FY97 for the analysis and design of cylindrical structural components with radially graded microstructures. The generalized theory is applicable to thin multi-phased composite shells/cylinders subjected to steady-state thermomechanical, transient thermal and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial direction, and arbitrarily distributed in the radial and circumferential direction, thereby allowing functional grading of the internal reinforcement in the r-Theta plane. The four computer codes fgmc3dq.cylindrical.f, fgmp3dq.cylindrical.f, fgmgvips3dq.cylindrical.f, and fgmc3dq.cylindrical.transient.f are research-oriented codes for investigating the effect of functionally graded architectures, as well as the properties of the multi-phase reinforcement, in thin shells subjected to thermomechanical and inertial loading, on the internal temperature, stress and (inelastic) strain fields. The reinforcement distribution in the radial and circumferential directions is specified by the user. The thermal and inelastic properties of the individual phases can vary with temperature. The inelastic phases are presently modeled by the power-law creep model generalized to multi-directional loading (within fgmc3dq.cylindrical.f and fgmc3dq.cylindrical.transient.f for steady-state and transient thermal loading, respectively), and incremental plasticity and GVIPS unified viscoplasticity theories (within the steady-state loading versions fgmp3dq.cylindrical.f and fgmgvips3dq.cylindrical.f).

Pindera, Marek-Jerzy↗

Thermal Considerations for 2039 Opposition Class Nuclear Electric Propulsion/Chemical Propulsion Crewed Mars Mission

The high specific impulse (Isp) of Nuclear Electric Propulsion (NEP) technology offers the potential for advanced space mission capabilities. However, the five critical technology elements of NEP vehicles have yet to prove technical maturity levels for consideration into mission design. In response to the critical reviews by the NASA Engineering and Safety Center (NESC) and the National Academies of Sciences, Engineering, and Medicine (NASEM), NASA’s Space Nuclear Propulsion (SNP) project created an NEP Technology Maturation Plan (TMP) for focused development of NEP technology. The TMP called for a coordinated set of technology development efforts to meet this objective. The Modular Assembled Radiators for NEP VehicLes (MARVL) Early Career Initiative (ECI) project was initiated to develop a portion of the fifth Critical Technology Element (CTE) of the NEP vehicle: the Primary Heat Rejection Subsystem (PHRS). A target application of a 2039 human-rated Mars mission was outlined in the TMP. For the outlined mission, a NEP vehicle will experience several thermal environments which will impact the design and operation of the PHRS. To maintain radiator temperatures within the required effective temperature range, the effect of natural, induced, and NEP internally generated heat loads on the radiator panel must be well understood. Furthermore, this analysis is critical for analyzing the influence of various orientations and positions of the NEP vehicle relative to nearby celestial bodies throughout the mission. This study conducted a complete enveloping analysis of the thermal environments influencing the NEP vehicle throughout the mission. Thermal analysis was conducted for the radiator panels based on the defined mission environments. This thermal analysis concludes with the selection of ideal radiator orientations for the NEP vehicle, and the identification of worst case hot and cold environmental sink temperatures throughout the mission. For the target application, the environmental sink temperature while the reactor is powered OFF or powered ON ranges from 30 K to 353 K and 2.7 K to 243 K respectively. When considering interplanetary space, the minimum environmental sink temperature when the reactor is powered OFF and the radiators are oriented “edge to Sun” is 2.7 K. The environmental thermal models generated in this study will be used for future studies with the full vehicle system model. The environmental sink temperature curves generated will be used for future radiator and component analysis to predict transient performance in the space environment. The environmental sink temperature and heat rejection capability curves will inform the trade between commissioning orbits that are in consideration. The model may also serve as a useful tool as reference for future crewed space missions, missions involving radiators or temperature sensitive equipment, or other missions requiring analysis of natural orbital thermal environments.

Nuclear Electric Propulsion↗

Maneuver Design and Calibration for the Genesis Spacecraft

Genesis is the fifth mission selected as part of NASA's Discovery Program. The objective of Genesis is to collect solar wind samples for a period of approximately two years while in a halo orbit about the Earth-Sun L I point. At the end of this period, the samples are to be returned to a specific recovery point on the Earth for subsequent analysis. This goal has never been attempted before and presents a formidable challenge in terms of mission design and operations, particularly planning and execution of propulsive maneuvers. To achieve a level of cost-effectiveness consistent with a Discovery-class mission, the Genesis spacecraft design was adapted to the maximum extent possible from designs used on earlier missions, such as Mars Global Surveyor (MGS) and Stardust, another sample collection mission. The spacecraft design for Genesis is shown. Spin stabilization was chosen for attitude control, in lieu of three-axis stabilization, with neither reaction wheels nor accelerometers included. This precludes closed-loop control of propulsive maneuvers and implies that any attitude changes, including spin changes and precessions, will behave like translational propulsive maneuvers and affect the spacecraft trajectory. Moreover, to minimize contamination risk to the samples collected, all thrusters were placed on the side opposite the sample collection canister. The orientation and characteristics of thrusters are indicated. For large maneuvers (>2.5 m/s), two 5 lbf thrusters will be used for delta v, with precession to the burn attitude, followed by spin-up from 1.6 to 10 rpm before the burn and spin down to 1.6 rpm afterwards, then precession back to the original spin attitude. For small maneuvers (<2.5 m/s), no spin change is needed and four 0.2 lbf thrusters are used for Av. Single or double 360 deg. precession changes are required whenever the desired delta v falls inside the two-way turn circle (about 0.4 m/s) based on the mass properties, spin rate and lever arm lengths based on thruster locations. In such instances, delta v resulting from spacecraft precession cannot be used effectively as a component of the desired delta v, and must therefore be removed by precessing at least one complete revolution around the turn circle. To eliminate cross-track execution errors, a second revolution in the opposite direction would also be performed. This paper will address the design of propulsive maneuvers in light of the aforementioned challenges and other constraints. Maneuver design will be performed jointly by the Navigation Team at JPL and the Spacecraft Team at LMA, based on the process indicated . Typical maneuver timelines will be presented which address considerations introduced by attitude changes. These include nutation, which is introduced by precessing or spinning down and must be given sufficient time to damp out prior to execution of subsequent events, as well as sun and earth pointing constraints, which must be considered to ensure sufficient spacecraft power and to minimize telecommunications interruptions, respectively. The paper will include a description of how individual propulsive maneuvers are resolved into components to account for delta v from translational burns and spacecraft attitude changes required to carry out such maneuvers. Contributions to maneuver delta v arising from attitude changes, based on mass properties for the period just after launch, are indicated. Similar curves will be presented spanning all mission phases from launch through return. A set of closed-form equations for resolving maneuver components, base on a specific delta v required for correction or deterministic changes to the spacecraft trajectory will be presented, as well. In addition to nominal maneuvers, special calibration maneuvers are planned to improve open-loop modeling of maneuvers and to reduce execution errors. Uncalibrated execution errors are indicated. Such errors could be reduced by 50% or more over the course of the mission. Special calibrations are of particular importance for the return leg of the mission, since the sample canister must be returned to a specific location within the Utah Test and Training Range (UTTR) for mid-air retrieval. An entry angle tolerance of no less than +/- 0.08 deg. is required to achieve this objective. Biasing of the final return maneuvers coupled with a specific maneuver mode to use a series of well-characterized spin changes to effect these maneuvers is part of the current Genesis baseline mission plan. Another important objective of calibrations is to better characterize precession maneuvers. Such maneuvers are part of most propulsive maneuvers, but are also required periodically to maintain sun-pointing for power or daily during solar-wind pointing during collection periods. Although relatively small, such maneuvers will have a significant cumulative impact on orbit determination, particularly in the halo portion of the mission. The current mission design also calls for three stationkeeping maneuvers during each halo orbit of approximately six months duration. These stationkeeping maneuvers may be sufficiently small that single or double 360 deg. precession changes may be required. Because there are no accelerometers on board the spacecraft, calibration can only be performed with the aid of ground-based radiometric tracking. To establish a high degree of accuracy in characterizing the magnitude of burns, the spacecraft spin axis should be along the line of sight to the Earth, providing Doppler measurements with <1 mm/sec accuracy in S-Band. Emission constraints allow such alignment only during certain portions of the mission when the Earth-spacecraft-sun geometry is favorable. The impact of precessions, or burns at times when geometry is not favorable, can be assessed by reconstruction of the spacecraft trajectory using tracking arcs of several days before and after the event.

Williams, Kenneth E.↗

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis↗

Single-Frame Terrain Mapping Software for Robotic Vehicles

This software is a component in an unmanned ground vehicle (UGV) perception system that builds compact, single-frame terrain maps for distribution to other systems, such as a world model or an operator control unit, over a local area network (LAN). Each cell in the map encodes an elevation value, terrain classification, object classification, terrain traversability, terrain roughness, and a confidence value into four bytes of memory. The input to this software component is a range image (from a lidar or stereo vision system), and optionally a terrain classification image and an object classification image, both registered to the range image. The single-frame terrain map generates estimates of the support surface elevation, ground cover elevation, and minimum canopy elevation; generates terrain traversability cost; detects low overhangs and high-density obstacles; and can perform geometry-based terrain classification (ground, ground cover, unknown). A new origin is automatically selected for each single-frame terrain map in global coordinates such that it coincides with the corner of a world map cell. That way, single-frame terrain maps correctly line up with the world map, facilitating the merging of map data into the world map. Instead of using 32 bits to store the floating-point elevation for a map cell, the vehicle elevation is assigned to the map origin elevation and reports the change in elevation (from the origin elevation) in terms of the number of discrete steps. The single-frame terrain map elevation resolution is 2 cm. At that resolution, terrain elevation from 20.5 to 20.5 m (with respect to the vehicle's elevation) is encoded into 11 bits. For each four-byte map cell, bits are assigned to encode elevation, terrain roughness, terrain classification, object classification, terrain traversability cost, and a confidence value. The vehicle s current position and orientation, the map origin, and the map cell resolution are all included in a header for each map. The map is compressed into a vector prior to delivery to another system.

Rankin, Arturo L.↗

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

Processing Motion Signals in Complex Environments

Motion information is critical for human locomotion and scene segmentation. Currently we have excellent neurophysiological models that are able to predict human detection and discrimination of local signals. Local motion signals are insufficient by themselves to guide human locomotion and to provide information about depth, object boundaries and surface structure. My research is aimed at understanding the mechanisms underlying the combination of motion signals across space and time. A target moving on an extended trajectory amidst noise dots in Brownian motion is much more detectable than the sum of signals generated by independent motion energy units responding to the trajectory segments. This result suggests that facilitation occurs between motion units tuned to similar directions, lying along the trajectory path. We investigated whether the interaction between local motion units along the motion direction is mediated by contrast. One possibility is that contrast-driven signals from motion units early in the trajectory sequence are added to signals in subsequent units. If this were the case, then units later in the sequence would have a larger signal than those earlier in the sequence. To test this possibility, we compared contrast discrimination thresholds for the first and third patches of a triplet of sequentially presented Gabor patches, aligned along the motion direction. According to this simple additive model, contrast increment thresholds for the third patch should be higher than thresholds for the first patch.The lack of a measurable effect on contrast thresholds for these various manipulations suggests that the pooling of signals along a trajectory is not mediated by contrast-driven signals. Instead, these results are consistent with models that propose that the facilitation of trajectory signals is achieved by a second-level network that chooses the strongest local motion signals and combines them if they occur in a spatio-temporal sequence consistent with a trajectory. These results parallel the lack of increased apparent contrast along a static contour made up of similarly oriented elements.

Verghese, Preeti↗

CAPRI: A Geometric Foundation for Computational Analysis and Design

CAPRI is a software building tool-kit that refers to two ideas; (1) A simplified, object-oriented, hierarchical view of a solid part integrating both geometry and topology definitions, and (2) programming access to this part or assembly and any attached data. A complete definition of the geometry and application programming interface can be found in the document CAPRI: Computational Analysis PRogramming Interface appended to this report. In summary the interface is subdivided into the following functional components: 1. Utility routines -- These routines include the initialization of CAPRI, loading CAD parts and querying the operational status as well as closing the system down. 2. Geometry data-base queries -- This group of functions allow all top level applications to figure out and get detailed information on any geometric component in the Volume definition. 3. Point queries -- These calls allow grid generators, or solvers doing node adaptation, to snap points directly onto geometric entities. 4. Calculated or geometrically derived queries -- These entry points calculate data from the geometry to aid in grid generation. 5. Boundary data routines -- This part of CAPRI allows general data to be attached to Boundaries so that the boundary conditions can be specified and stored within CAPRI s data-base. 6. Tag based routines -- This part of the API allows the specification of properties associated with either the Volume (material properties) or Boundary (surface properties) entities. 7. Geometry based interpolation routines -- This part of the API facilitates Multi-disciplinary coupling and allows zooming through Boundary Attachments. 8. Geometric creation and manipulation -- These calls facilitate constructing simple solid entities and perform the Boolean solid operations. Geometry constructed in this manner has the advantage that if the data is kept consistent with the CAD package, therefore a new design can be incorporated directly and is manufacturable. 9. Master Model access This addition to the API allows for the querying of the parameters and dimensions of the model. The feature tree is also exposed so it is easy to see where the parameters are applied. Calls exist to allow for the modification of the parameters and the suppression/unsuppression of nodes in the tree. Part regeneration is performed by a single API call and a new part becomes available within CAPRI (if the regeneration was successful). This is described in a separate document. Components 1-7 are considered the CAPRI base level reader.

Haimes, Robert↗