Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global shape 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 235 records · Page 13

Following the Dust: Satellite Perspectives of Dust Sources, Transport, Deposition, and Impacts

Mineral dust is an integral component in the Earth system that interacts with the system's many other components involving the energy, water, and carbon cycles. Dust also degrades air quality and adversely affects human health. These interactions and impacts are not contained in regions nearby dust sources, but can reach very far because of the long-range transport on intercontinental and global scales. Satellite's routine sampling and extensive coverage in time and space makes it an ideal platform to follow the dust from sources to sinks and assess its impacts along the long journey. Dust particles are unique in their coarse size and irregular shape, which makes it feasible to distinguish them from other aerosol particles using remote sensing techniques. This talk will provide an overview of what we have learned from analyzing advanced satellite remote sensing measurements during the EOS-era supplemented by in situ observations and model simulations, including dust source characterization, seasonal and interannual variability, trans-Pacific and trans-Atlantic transport and deposition, and dust influences on the radiation budget, air quality, and ecosystems. The talk will also discuss challenges and opportunities to further improve the dust characterization and assessment of the impacts via remote sensing techniques.

Yu, Hongbin↗

Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models

Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.

Geometric imperfections↗

Applying transfer matrix method to the estimation of the modal characteristics of the NASA Mini-Mass Truss

It is beneficial to use a distributed parameter model for large space structures because the approach minimizes the number of model parameters. Holzer's transfer matrix method provides a useful means to simplify and standardize the procedure for solving the system of partial differential equations. Any large space structures can be broken down into sub-structures with simple elastic and dynamical properties. For each single element, such as beam, tether, or rigid body, we can derive the corresponding transfer matrix. Combining these elements' matrices enables the solution of the global system equations. The characteristics equation can then be formed by satisfying the appropriate boundary conditions. Then natural frequencies and mode shapes can be determined by searching the roots of the characteristic equation at frequencies within the range of interest. This paper applies this methodology, and the maximum likelihood estimation method, to refine the modal characteristics of the NASA Mini-Mast Truss by successively matching the theoretical response to the test data of the truss. The method is being applied to more complex configurations.

Shen, Ji-Yao↗

Effect of Wind Velocity on Flame Spread in Microgravity

A three-dimensional, time-dependent model is developed describing ignition and subsequent transition to flame spread over a thermally thin cellulosic sheet heated by external radiation in a microgravity environment. A low Mach number approximation to the Navier Stokes equations with global reaction rate equations describing combustion in the gas phase and the condensed phase is numerically solved. The effects of a slow external wind (1-20 cm/s) on flame transition are studied in an atmosphere of 35% oxygen concentration. The ignition is initiated at the center part of the sample by generating a line-shape flame along the width of the sample. The calculated results are compared with data obtained in the 10s drop tower. Numerical results exhibit flame quenching at a wind speed of 1.0 cm/s, two localized flames propagating upstream along the sample edges at 1.5 cm/s, a single line-shape flame front at 5.0 cm/s, three flames structure observed at 10.0 cm/s (consisting of a single line-shape flame propagating upstream and two localized flames propagating downstream along sample edges) and followed by two line-shape flames (one propagating upstream and another propagating downstream) at 20.0 cm/s. These observations qualitatively compare with experimental data. Three-dimensional visualization of the observed flame complex, fuel concentration contours, oxygen and reaction rate isosurfaces, convective and diffusive mass flux are used to obtain a detailed understanding of the controlling mechanism, Physical arguments based on lateral diffusive flux of oxygen, fuel depletion, oxygen shadow of the flame and heat release rate are constructed to explain the various observed flame shapes.

Prasad, Kuldeep↗

BrO and Inferred Bry Profiles over the Western Pacific: Relevance of Inorganic Bromine Sources and a Bry Minimum in the Aged Tropical Tropopause Layer

We report measurements of bromine monoxide (BrO) and use an observationally constrained chemical box model to infer total gas-phase inorganic bromine (Br(sub y)) over the tropical western Pacific Ocean (tWPO) during the CONTRAST field campaign (January-February 2014). The observed BrO and inferred Bry profiles peak in the marine boundary layer (MBL), suggesting the need for a bromine source from sea-salt aerosol (SSA), in addition to organic bromine (CBry ). Both profiles are found to be C-shaped with local maxima in the upper free troposphere (FT). The median tropospheric BrO vertical column density (VCD) was measured as 1.6 x 10(exp 13) molec cm(exp -2), compared to model predictions of 0.9 x 10(exp 13) molec cm(exp -2) in GEOS-Chem (CBr(sub y) but no SSA source), 0.4 x 10(exp 13) molec cm(exp -2) in CAM-Chem (CBr(sub y) and SSA), and 2.1 x 10(exp 13) molec cm(exp -2) in GEOS-Chem (CBry and SSA). Neither global model fully captures the Cshape of the Br(sun y) profile. A local Br(sub y) maximum of 3.6 ppt (2.9-4.4 ppt; 95% confidence interval, CI) is inferred between 9.5 and 13.5 km in air masses influenced by recent convective outflow. Unlike BrO, which increases from the convective tropical tropopause layer (TTL) to the aged TTL, gas-phase Br(sub y) decreases from the convective TTL to the aged TTL. Analysis of gas-phase Br(sub y) against multiple tracers (CFC-11, H2O/O3 ratio, and potential temperature) reveals a Br(sub y) minimum of 2.7 ppt (2.3-3.1 ppt; 95% CI) in the aged TTL, which agrees closely with a stratospheric injection of 2.6 +/- 0.6 ppt of inorganic Br(sub y) (estimated from CFC-11 correlations), and is remarkably insensitive to assumptions about heterogeneous chemistry. Bry increases to 6.3 ppt (5.6-7.0 ppt; 95% CI) in the stratospheric "middleworld" and 6.9 ppt (6.5-7.3 ppt; 95% CI) in the stratospheric "overworld". The local Br(sub y) minimum in the aged TTL is qualitatively (but not quantitatively) captured by CAM-Chem, and suggests a more complex partitioning of gas-phase and aerosol Br(sub y) species than previously recognized. Our data provide corroborating evidence that inorganic bromine sources (e.g., SSA-derived gas-phase Br(sub y) ) are needed to explain the gas-phase Br(sub y) budget in the upper free troposphere and TTL. They are also consistent with observations of significant bromide in Upper Troposphere-Lower Stratosphere aerosols. The total Br(sub y) budget in the TTL is currently not closed, because of the lack of concurrent quantitative measurements of gas-phase Br(sub y) species (i.e., BrO, HOBr, HBr, etc.) and aerosol bromide. Such simultaneous measurements are needed to (1) quantify SSA-derived Br(sub y) in the upper FT, (2) test Br(sub y) partitioning, and possibly explain the gas-phase Br(sub y) minimum in the aged TTL, (3) constrain heterogeneous reaction rates of bromine, and (4) account for all of the sources of Br(sub y) to the lower stratosphere.

Atmospheric Chemistry; bromine; troposphere; CONTR↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗

The MODIS Aerosol Algorithm, Products, Validation and Applications

The MODerate resolution Imaging Spectroradiometer (MODIS) currently aboard both the Terra and Aqua satellites produces a suite of products designed to characterize global aerosol distribution, optical thickness and particle size. Never before has a space-borne instrument been able to provide such detailed information, complementing field and modeling efforts to produce a comprehensive picture of aerosol characteristics. The three years of Terra-MODIS data have been validated by comparing with co-located AERONET observations of aerosol optical thickness and derivations of aerosol size parameters. Some 8000 comparison points located at 133 AERONET sites around the globe show that the MODIS aerosol optical thickness retrievals are accurate to within the pre-launch expectations. MODIS-derived size parameters are also compared with AERONET retrievals and found to agree well for fine-mode dominated aerosol regimes. Aerosol regimes dominated by dust aerosol are less accurate, attributed to what is thought to be nonsphericity. Errors due to nonsphericity will be reduced by introducing a new set of empirical phase functions, derived without any assumptions of particle shape. The major innovation that MODIS bring to the field of remote sensing of aerosol is the measure of particle size and the separation of finemode and coarsemode dominated aerosol regimes. Particle size can separate finemode man-made aerosols created during combustion, from larger natural aerosols originating from salt spray or wind erosion. This separation allows for the calculation of aerosol radiative effect and the estimation of the man-made aerosol radiative forcing. MODIS can also be used in regional studies of aerosol-cloud interaction that affect the global radiative and hydrological cycles.

Remer, L. A.↗

Transmantle flux tectonics

Venus, Earth, and Mars have surfaces that display topographic domes and depressions with quasi-circular planimetric shapes, relief of 0 to several km, and large spatial scales (10(exp 2) to 10(exp 4) km). Our morphostructural mapping reveals hierarchical arrangements of these features. They are explained by a model of long-acting mantle convection, as a particular case of convection in a stratified and random inhomogeneous medium, which develops the form of a hierarchy of different convective pattern scales, each arising from different levels in the mantle. The hypothesis of transmantle flux tectonics parsimoniously explains a diversity of seemingly unrelated terrestrial planetary phenomena, including Earth megaplumes, global resurfacing epochs on Venus, and cyclic ocean formation and global climate change for Mars. All these phenomenon are hypothesized to be parsimoniously explained by a process of transmantle flux tectonics in which long-acting mantle convection generates stresses in blocks of planetary lithosphere to produce distinctive quasi-circular global-hierarchical morphostructure (QGM) patterns. Transmantle flux tectonics differs from plume tectonics in that individual plumes are not considered in isolation. Rather, a wholly interactive process is envisioned in which various spatial and temporal scales of convection operate contemporaneously and hierarchically within other scales. This process of continual change by hierarchical convective cells affects the surface at varying temporal and spatial scales, and its effects are discernable through their relic geological manifestations, the QGM patterns.

Finn, V. J.↗

Global models: Robot sensing, control, and sensory-motor skills

Robotics research has begun to address the modeling and implementation of a wide variety of unstructured tasks. Examples include automated navigation, platform servicing, custom fabrication and repair, deployment and recovery, and science exploration. Such tasks are poorly described at onset; the workspace layout is partially unfamiliar, and the task control sequence is only qualitatively characterized. The robot must model the workspace, plan detailed physical actions from qualitative goals, and adapt its instantaneous control regimes to unpredicted events. Developing robust representations and computational approaches for these sensing, planning, and control functions is a major challenge. The underlying domain constraints are very general, and seem to offer little guidance for well-bounded approximation of object shape and motion, manipulation postures and trajectories, and the like. This generalized modeling problem is discussed, with an emphasis on the role of sensing. It is also discussed that unstructured tasks often have, in fact, a high degree of underlying physical symmetry, and such implicit knowledge should be drawn on to model task performance strategies in a methodological fashion. A group-theoretic decomposition of the workspace organization, task goals, and their admissible interactions are proposed. This group-mechanical approach to task representation helps to clarify the functional interplay of perception and control, in essence, describing what perception is specifically for, versus how it is generically modeled. One also gains insight how perception might logically evolve in response to needs of more complex motor skills. It is discussed why, of the many solutions that are often mathematically admissible to a given sensory motor-coordination problem, one may be preferred over others.

Schenker, Paul S.↗

Effect of Protuberance Shape and Orientation on Space Shuttle Orbiter Boundary-Layer Transition

This document describes an experimental study conducted to examine the effects of protuberances on hypersonic boundary-layer transition. The experiment was conducted in the Langley 20-Inch Mach 6 Tunnel on a series of 0.9%-scale Shuttle Orbiter models. The data were acquired to complement the existing ground-based boundary-layer transition database that was used to develop Version 1.0 of the boundary-layer transition RTF (return-to-flight) tool. The existing ground-based data were all acquired on 0.75%-scale Orbiter models using diamond-shaped ( pizza-box ) trips. The larger model scale facilitated in manufacturing higher fidelity protuberances. The end use of this experimental database will be to develop a technical basis (in the form of a boundary-layer transition correlation) to assess representative protrusion shapes, e.g., gap fillers and protrusions resulting from possible tile repair concepts. The primary objective of this study is to investigate the effects of protuberance-trip location and geometry on Shuttle Orbiter boundary-layer transition. Secondary goals are to assess the effects of gap-filler orientation and other protrusion shapes on boundary-layer transition. Global heat-transfer images using phosphor thermography of the Orbiter windward surface and the corresponding streamwise and spanwise heating distributions were used to infer the state of the boundary layer, i.e., laminar, transitional, or turbulent.

King, RUdolph A.↗

Optimized Non-Obstructive Particle Damping (NOPD) Treatment for Composite Honeycomb Structures

Non-Obstructive Particle Damping (NOPD) technology is a passive vibration damping approach whereby metallic or non-metallic particles in spherical or irregular shapes, of heavy or light consistency, and even liquid particles are placed inside cavities or attached to structures by an appropriate means at strategic locations, to absorb vibration energy. The objective of the work described herein is the development of a design optimization procedure and discussion of test results for such a NOPD treatment on honeycomb (HC) composite structures, based on finite element modeling (FEM) analyses, optimization and tests. Modeling and predictions were performed and tests were carried out to correlate the test data with the FEM. The optimization procedure consisted of defining a global objective function, using finite difference methods, to determine the optimal values of the design variables through quadratic linear programming. The optimization process was carried out by targeting the highest dynamic displacements of several vibration modes of the structure and finding an optimal treatment configuration that will minimize them. An optimal design was thus derived and laboratory tests were conducted to evaluate its performance under different vibration environments. Three honeycomb composite beams, with Nomex core and aluminum face sheets, empty (untreated), uniformly treated with NOPD, and optimally treated with NOPD, according to the analytically predicted optimal design configuration, were tested in the laboratory. It is shown that the beam with optimal treatment has the lowest response amplitude. Described below are results of modal vibration tests and FEM analyses from predictions of the modal characteristics of honeycomb beams under zero, 50% uniform treatment and an optimal NOPD treatment design configuration and verification with test data.

Panossian, H.↗

Canopy height and climate dryness parsimoniously explain spatial variation of unstressed stomatal conductance

The spatio-temporal variation of stomatal conductance directly regulates photosynthesis, water partitioning, and biosphere-atmosphere interactions. While many studies have focused on stomatal response to stresses, the spatial variation of unstressed stomatal conductance remains poorly determined, and is usually characterized in land surface models (LSMs) simply based on plant functional type (PFT). Here, we derived unstressed stomatal conductance at the ecosystem-scale using observations from 115 global FLUXNET sites. When aggregated by PFTs, the across-PFT pattern was highly consistent with the parameterizations of LSMs. However, PFTs alone captured only 17% of the variation in unstressed stomatal conductance across sites. Within the same PFT, unstressed stomatal conductance was negatively related to climate dryness and canopy height, which explained 45% of the total spatial variation. Our results highlight the importance of plant environment interactions in shaping stomatal traits. The trait-environment relationship established here provides an empirical approach for improved parameterizations of stomatal conductance in LSMs.

Yanlan Liu↗

Particle sizes and composition of Mars atmospheric dust based upon Viking and Mariner 9 observations

Mars atmospheric dust can play an important role in the thermal structure of the Mars atmosphere during periods of high dust loading. However, the radiative properties of Mars atmospheric dust remain uncertain due to uncertain definitions of the dust composition and size distribution. The analysis by Toon et al., of Mariner 9 IRIS spectra during the 1971-1972 global dust storm indicated a reasonable match between the modeled 9-micron absorption of montmorillinite and the observed 9-micron absorption. Toon et al. also determined that an effective (cross-section weighted) mean radius of 2.5 microns (R(sub mode) = 0.4 microns) provided a consistent fit of montmorillinite to the IRIS dust spectra at 9 microns. Pollack et al. analyzed Viking lander observations of atmospheric extinction and scattering at visible-near IR wavelengths (0.5-1.0 microns), and obtained consistency with the Toon et al. dust size distribution when the effects of nonspherical particle shapes were included. An additional, minor (1 percent) component of visible-ultraviolet absorbing material was required to model the derived visible (0.86) and ultraviolet (0.4-0.6) single-scattering albedos of the dust, since montmorillinite does not absorb sufficiently in this wavelength region. A combined analysis of the Viking IRTM and Mariner 9 observations was conducted to reassess the model of Mars atmospheric ultraviolet-to-infrared measurements of dust absorption and scattering. The optical constants for palagonite are incorporated in a doubling-adding radiative transfer model of the Mars atmosphere to simulate Mariner 9 IRIS spectra as well as the Viking IRTM IR band observations. Visible and ultraviolet single-scattering albedos based on the Hansen and Travis Mie scattering code were also derived. A tentative conclusion is that smaller dust particles (R(sub mode) = 0.15 microns, cross-section weighted mean R = 1.2 microns) composed of palagonite provide a much improved fit to the Mariner 9 IRIS spectra; agreement with the observed ratio of visible-to-infrared extinction opacities; and ultraviolet and visible single-scattering albedos comparable to their observed values.

Clancy, R. T.↗

Theoretical and Observational Determination of Global and Regional Radiation Budget, Forcing and Feedbacks at the Top-of-Atmosphere and Surface

Report consists of: 1. List of accomplishments 2. List of publications 3. Abstracts of published or submitted papers and 4. Subject invention disclosure. The accomplishments of the grant listed are: 1. Improved the third-order turbulence closure in cloud resolving models to remove the liquid water oscillation. 2. Used the University of California-Los Angeles (UCLA) large-eddy simulation (LES) model to provide data for radiation transfer testing. 3. Revised shortwave k-distribution models based on HITRAN 2000. 4. Developed a gamma-weighted two-stream radiative transfer model for radiation budget estimate applications. 5. Estimated the effect of spherical geometry to the earth radiation budget. 6. Estimated top-of-atmosphere irradiance over snow and sea ice surfaces. 7. Estimated the aerosol direct radiative effect at the top of the atmosphere. 8. Estimated the top-of-atmosphere reflectance of the clear-sky molecular atmosphere over ocean. 9. Developed and validated new set of Angular Distribution Models for the CERES TRMM satellite instrument (tropical) 10. Developed and validated new set of Angular Distribution Models for the CERES Terra satellite instrument (global) 11. Quantified the top-of-atmosphere direct radiative effect of aerosols over global oceans from merged CERES and MODIS observations 12 Clarified the definition of TOA flux reference level for radiation budget studies 13. Developed new algorithm for unfaltering CERES measured radiances 14. Used multiangle POLDER measurements to produce narrowband angular distribution models and examine the effect of scene identification errors on TOA albedo estimates 15. Developed and validated a novel algorithm called the Multidirectional Reflectance Matching (MRM) model for inferring TOA albedos from ice clouds using multi-angle satellite measurements. 16. Developed and validated a novel algorithm called the Multidirectional Polarized Reflectance Matching (MPRM) model for inferring particle shapes from ice clouds using multi-angle polarized satellite measurements. 17. Developed 4 advanced light scattering models including the three-dimensional (3D) uniaxial perfectly matched layer (UPML) finite-difference time-domain (FDTD) model. 18. Develop sunglint in situ measurement and study reflectance distribution in the sunglint area. 19. Lead a balloon-borne radiometer TOA albedo validation effort. 20. Developed a CERES surface UVB, UVA, and UV index product.

Loeb, Norman G.↗

Evaluation of GEOS Total Cloud Fraction with GLOBE Citizen Science Observations and Co-Located Satellite Data

Here we evaluate the total cloud fraction in cycled forecast experiments with the NASA Global Earth Observing System (GEOS) model. Forecasts were run for summer and winter periods of 2017/2018 and compared with ground-based and satellite observations. Citizen science observations from the Global Learning and Observations to Benefit the Environment (GLOBE) Program were matched with MODIS satellite overpasses and geostationary cloud fraction, yielding a dataset of coincident observations for comparison with hourly model output. The observations indicate a model bias toward overcast and clear conditions, with underestimation of intermediate cloud fractions. We investigate underestimation of variance in the sub-grid total water probability density function (PDF) as a possible cause. The PDF determines large-scale cloud fraction as the sub-grid fraction with total water exceeding the saturation specific humidity. We conduct model experiments in which the shape, width, and spatial dependence of the PDF are varied, and consider the impacts on cloud fraction relative to the combined observations.

Starke, Matthew↗

Photometric Correction of Hayabusa2's NIRS3 Spectra of Asteroid 162173 Ryugu Using Empirical Photometric Models

In 2018, Japanese Aerospace Exploration Agency’s spacecraft Hayabusa2 began a near infrared spectroscopic imaging survey of near-Earth asteroid 162173 Ryugu. Hayabusa2 is a successful sample-return mission with an overarching goal to provide a better understanding of the origin and evolution of our solar system. The target, Ryugu, is a low-albedo carbonaceous asteroid (Cb-type) that is linked to carbonaceous chondrite meteorites. It is thought to have originated in the main asteroid belt and migrated inward to become a near-Earth asteroid, and therefore reachable by spacecraft. Preliminary findings from processed NIRS3data have shown that hydroxyl-bearing minerals are present on the surface of Ryugu, and it is likely the result of impact fragments from an aqueously altered parent body.In this project, we have used newly calibrated and processed NIRS3 data using updated shape models of Ryugu. NIRS3 spectra have both a thermal and reflectance component. The thermal component (beyond 2.5 microns) was modeled and removed from all NIRS3 spectra. Ryugu spectra were taken at different viewing geometries, and a photometric model needed to be developed to normalize all of the spectra at a common geometry for each mission phase. We have used three empirical models: Minnaert, Lommel-Seeliger, and ROLO (RObotic LunarOrbiter). These models were chosen for their ability to relate the surface reflectance to the viewing geometry, as well as their compatibility with the asteroid’s albedo range. The three models provide the global light scattering properties of Ryugu’s surface and subsequently enable us to calculate the geometric albedo, phase integral, spherical bond albedo, and the average surface normal albedo for Ryugu.

Lucille Grace Williamson↗

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors.

Pifko, A. B.↗

Current plate motions

A global plate motion model, named NUVEL-1, which describes current plate motions between 12 rigid plates is described, with special attention given to the method, data, and assumptions used. Tectonic implications of the patterns that emerged from the results are discussed. It is shown that wide plate boundary zones can form not only within the continental lithosphere but also within the oceanic lithosphere; e.g., between the Indian and Australian plates and between the North American and South American plates. Results of the model also suggest small but significant diffuse deformation of the oceanic lithosphere, which may be confined to small awkwardly shaped salients of major plates.

Demets, C.↗