Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 1,333 records · Page 74

Trimming an aircraft model for flight simulation

Real-time piloted aircraft simulations with digital computers have been performed at Ames Research Center (ARC) for over two decades. For the simulation of conventional aircraft models, the establishment of initial vehicle and control orientations at various operational flight regimes has been adequately handled by either analog techniques or simple inversion processes. However, exotic helicopter configurations have been introduced recently that require more sophisticated techniques because of their expanded degrees of freedom and environmental vibration levels. At ARC, these techniques are used for the backward solutions to real-time simulation models as required for the generation of trim points. These techniques are presented in this paper with examples from a blade-element helicopter simulation model.

Mcfarland, Richard E.↗

Comprehensive Understanding for Vegetated Scene Radiance Relationships

Directional reflectance distributions spanning the entire existent hemisphere were measured in two field studies; one using a Mark III 3-band radiometer and one using the rapid scanning bidirectional field instrument called PARABOLA. Surfaces measured included corn, soybeans, bare soils, grass lawn, orchard grass, alfalfa, cotton row crops, plowed field, annual grassland, stipa grass, hard wheat, salt plain shrubland, and irrigated wheat. Analysis of field data showed unique reflectance distributions ranging from bare soil to complete vegetation canopies. Physical mechanisms causing these trends were proposed. A 3-D model was developed and is unique in that it predicts: (1) the directional spectral reflectance factors as a function of the sensor's azimuth and zenith angles and the sensor's position above the canopy; (2) the spectral absorption as a function of location within the scene; and (3) the directional spectral radiance as a function of the sensor's location within the scene. Initial verification of the model as applied to a soybean row crop showed that the simulated directional data corresponded relatively well in gross trends to the measured data. The model was expanded to include the anisotropic scattering properties of leaves as a function of the leaf orientation distribution in both the zenith and azimuth angle modes.

Kimes, D. S.↗

Dynamical evolution of the Oort cloud

New studies of the dynamical evolution of cometary orbits in the Oort cloud are made using a revised version of Weissman's (1982) Monte Carlo simulation model, which more accurately mimics the perturbation of comets by the giant planets. It is shown that perturbations by Saturn and Jupiter provide a substantial barrier to the diffusion of cometary perihelia into the inner solar system. Perturbations by Uranus and Neptune are rarely great enough to remove comets from the Oort cloud, but do serve to scatter the comets in the cloud in initial energy. The new model gives a population of 1.8 to 2.1 x 10 to the 12th comets for the present-day Oort cloud, and a mass of 7 to 8 earth masses. Perturbation of the Oort cloud by giant molecular clouds in the galaxy is discussed, as is evidence for a massive 'inner Oort cloud' internal to the observed one. The possibility of an unseen solar companion orbiting in the Oort cloud and causing periodic comet showers is shown to be dynamically plausible but unlikely, based on the observed cratering rate on the earth and moon.

Weissman, P. R.↗

Thoracic Pressure Does Not Impact CSF Pressure via Compartment Compliance

Space acquired neuro-ocular syndrome (SANS) remains a difficult risk to characterize due to the complex multi-factorial etiology related to physiological responses to the spaceflight environment. Fluid shift and the resultant change on the Cardiovascular (CV) and cerebral spinal fluid systems (CSF) in the absence of gravity continue to be considered a contributing factor to the progression of SANS. In this study, we utilize a computational model of the CSF and CV interface to establish the sensitivity that intracranial pressure, and subsequently the optic nerve sheath pressure, exhibits due to variations in thoracic pressure, assuming the cranial perfusion pressure, i.e. mean arterial pressure (MAP) to central venous pressure (CVP), is known. Methods: The GRC Cross cutting computational modeling project created as model of the CSF and CV interaction within the cranial vault by extending the work of Stevens et al. [1] by modifying the representative anatomy to include a separate venous sinus, jugular veins, secondary veins and extra jugular pathways [2-3] to more adequately represent the vascular drainage pathways from the cranial vault (Figure 1). Assuming the MAP, CVP and thoracic pressure are known, we initiated this enhanced computational model assuming a supine positon and utilized a linear ramp to vary the thoracic pressure from the assumed supine state to the target pressure corresponding to set MAP and CVP values. The model generates the time based CSF pressure values (Figure2). Results and Conclusions: Following this analysis, CSF pressure shows significant independence from thoracic pressure changes (16 mmHg in thoracic pressure produces < 1mmHg change in CSF pressure), being mostly dependent on perfusion pressure. Similarly fluid redistribution is not predicted to be impacted over a level of 1mL. We note that this simulation represents an acute changes (order of 10's of minutes) and does not represent the long term effects.

Munster, Drayton W.↗

Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study

Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.

Convery, P. D.↗

Three-Dimensional, Transgenic Cell Models to Quantify Space Genotoxic Effects

The space environment contains radiation and chemical agents known to be mutagenic and carcinogenic to humans. Additionally, microgravity is a complicating factor that may modify or synergize induced genotoxic effects. Most in vitro models fail to use human cells (making risk extrapolation to humans more difficult), overlook the dynamic effect of tissue intercellular interactions on genotoxic damage, and lack the sensitivity required to measure low-dose effects. Currently a need exists for a model test system that simulates cellular interactions present in tissue, and can be used to quantify genotoxic damage induced by low levels of radiation and chemicals, and extrapolate assessed risk to humans. A state-of-the-art, three-dimensional, multicellular tissue equivalent cell culture model will be presented. It consists of mammalian cells genetically engineered to contain multiple copies of defined target genes for genotoxic assessment,. NASA-designed bioreactors were used to coculture mammalian cells into spheroids, The cells used were human mammary epithelial cells (H184135) and Stratagene's (Austin, Texas) Big Blue(TM) Rat 2 lambda fibroblasts. The fibroblasts were genetically engineered to contain -a high-density target gene for mutagenesis (60 copies of lacl/LacZ per cell). Tissue equivalent spheroids were routinely produced by inoculation of 2 to 7 X 10(exp 5) fibroblasts with Cytodex 3 beads (150 micrometers in diameter). at a 20:1 cell:bead ratio, into 50-ml HARV bioreactors (Synthecon, Inc.). Fibroblasts were cultured for 5 days, an equivalent number of epithelial cells added, and the fibroblast/epithelial cell coculture continued for 21 days. Three-dimensional spheroids with diameters ranging from 400 to 600 micrometers were obtained. Histological and immunohistochemical Characterization revealed i) both cell types present in the spheroids, with fibroblasts located primarily in the center, surrounded by epithelial cells; ii) synthesis of extracellular matrix; and iii,, mitotic cells located throughout the spheroids. Spheroidal integrity and cell viability were retained for the 30-day test period after removal of spheroids from the bioreactor. Potential utility of this three-dimensional, transgenic model for genotoxicity was initially assessed by exposure of spheroids to 0-2 Gy neon at dose rates of 0.3 to 1.5 Gy/min (National Institute of Radiological Sciences, Chiba, Japan). Quantification of mutation at the lacl gene revealed a linear dose response for mutation induction. Limited sequencing analysis of mutant clones revealed higher frequencies of deletions and multiple base sequence changes with increasing dose. These results suggest that our three-dimensional, transgenic model is applicable to a wide variety of studies involving the quantification, identification, and characterization of genotoxicity incurred in space and on Earth. This model uniquely allows investigation of the interaction of relevant factors, namely cell-to-cell interactions and the mechanistic interaction of microgravity with radiation insults and DNA repair. Using this three-dimensional model will allow us to obtain dual genotoxic information (i.e., mutation rate plus chromosome aberration data) from the same system so that one endpoint can be used to reference the other, thereby increasing the fidelity of the data set. Moreover, the tissue-equivalent nature of the three-dimensional model provides high confidence for relevance of risk assessment, i.e., the establishment of quality factors directly applicable to the microgravity environment.

Gonda, S. R.↗

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗

LAMPS software and mesoscale prediction studies

The full-physics version of the LAMPS model has been implemented on the Perkin-Elmer computer. In Addition, LAMPS graphics processors have been rewritten to the run on the Perkin-Elmer and they are currently undergoing final testing. Numerical experiments investigating the impact of convective parameterized latent heat release on the evolution of a precipitating storm have been performed and the results are currently being evaluated. Curent efforts include the continued evaluation of the impact of initial conditions on LAMPS model results. This work will help define measurement requirements for future research field projects as well as for observations in support of operational forecasts. Also, the impact of parameterized latent heat on the evolution of precipitating systems is continuing. This research is in support of NASA's proposed Earth Observation Mission (EOM).

Perkey, D. J.↗

The Atmosphere Observing System (AOS): Synergistic Aerosol, Cloud, Convection and Precipitation Measurement and Modeling Systems

The 2017 Decadal Survey (DS) highlighted Earth System Science themes, science and application questions, and several high priority objectives that have led to the inclusion of Aerosols (A) and Clouds-Convection-Precipitation (CCP) as Designated Observables (DOs). On June 1, 2018, several NASA centers (GSFC, LaRC, JPL, MSFC, GRC and ARC) submitted a joint Study Plan to the NASA Earth Science Division for the Aerosol (A) and Cloud, Convection, and Precipitation (CCP) Pre-formulation Study (ACCP), with the ACCP Study concluding in early 2021. The new mission now in pre-phase A is being referred to as the Atmosphere Observing System (AOS), an integral part of NASA’s Earth System Observatory (ESO) strategy. The DS and the ACCP team recognized the science merit in combining the A and CCP DOs for both enhancing the ability to address a number of science objectives and also to provide an expanded capability to address additional objectives beyond those of the individual DOs. A critical element of the ACCP observing strategy is to make extensive use of new passive and active sensors as well as of the so-called Program-of-Record (PoR), complemented by a fully integrated sub-orbital component. In order to achieve maximum benefit, all these observations need to be integrated into comprehensive observing and modeling/data assimilation systems. Such an approach requires comprehensive model-data synthesis capabilities that needs to be conceived in conjunction with the space-based and suborbital components of AOS. In this presentation we will summarize the major science goals of AOS including cloud feedbacks, atmospheric convection, emphasizing aerosol processes and aerosol radiative effects, and the synergistic aspects of clouds-precipitation-aerosol interactions. We will describe examples of how AOS data will be used across space and time to better initialize forecasts and train modeling systems, and to infuse models and data assimilation systems with AOS data for advancing operational predictions and to generate expanded hindcasts and reconstruction of the climate record.

Arlindo da Silva↗

Prediction of Calcium Oxalate Stone Risk with Interventions

Introduction Mitigations intended to reduce bone mineral loss in-flight also influence the level of renal stone risk to astronauts. Extending a recent study completed by our group [1], we investigate the likely influence of combined fluid intake and Ca excretion mitigation strategies. Methods A computational biochemistry model representing CaOx crystal precipitation, growth, and agglomeration [2] is combined with a probabilistic analysis to predict the in-flight CaOx renal stone incidence risk ratio (IRR) relative to pre-flight values using 1517 astronaut 24-hour urine chemistries. The output of this simulation is used to assess the combined impact of controlling fluid intake and Ca excretion on maintaining IRR < 1.2 for > 95% of the population, a threshold determined from initial analysis of the model output compared to referent populations [1]. Results and Conclusions Figure 1 illustrates the simulated in-flight astronaut population (n=40000) with CaOx IRR > 1.2. The red lines indicate the known clinical risk boundary of each factor, and the white box is considered the area of high combined clinical risk. The simulation and identified risk threshold conform well to the expected clinical risk boundaries. The results, assuming insensible water losses equal or exceed 0.7 L/day, imply mitigations achieve risks similar to pre-flight levels when either Ca < 150 mg/day or fluid volume intake > 3.2 L/day individually, or combined Ca < 200 and fluid volume intake > 2.5 L/day. We continue to investigate similar compounded effects and look to apprise engineering trades by informing relative risk changes during mission design. References [1] D. A. Goodenow-Messman, S. A. Gokoglu, M. Kassemi, and J. Myers, Jerry G., “Characterization of Astronaut CaOx Renal Stone Incidence Rates to Quantify In-flight and Post-flight Relative Risk,” npj Microgravity, in review, 2021. [2] M. Kassemi and D. Thompson, “Prediction of renal crystalline size distributions in space using a PBE analytic model. 1. Effect of microgravity-induced biochemical alterations,” Am. J. Physiol. Physiol., vol. 311, no. 3, pp. F520–F530, Sep. 2016, doi: 10.1152/ajprenal.00401.2015.

Debra A Goodenow↗

Observations of Land Surface Variability Using Passive Microwave Sensing

Understanding the global variability of land surface wetness (soil moisture), skin temperature, and related surface fluxes of heat and moisture is key to assessing the importance of the land surface in influencing climate. The feasibility of producing model estimates of these quantities is being studied as part of the International Satellite Land Surface Climatology Project (ISLSCP) Global Soil Wetness Project (GSWP). In the GSWP approach, meteorological observations and analyses are used to drive global circulation models. Satellite measurements can provide independent estimates of key land surface parameters that are needed for initializing and validating the climate models and for monitoring long-term change. Satellite observations of the land surface can also be assimilated into soil models to estimate moisture in the root zone. In our research, passive microwave satellite data recorded during 1978-1987 from the Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR) are being used to examine spatial and temporal trends in surface soil moisture, vegetation, and temperature. These data include observations at C and X bands (6.6 and 10.7 GHz), which are not available on the current Special Sensor Microwave/Imager (SSM/I) and are precursors to data that will become available from the Advanced Microwave Scanning Radiometer (AMSR) on Advanced Earth Observing Satellite (ADEOS-II) and Earth Observing System (EOS) PM1 in the year 2000. A chart shows a time-series of SMMR-derived surface temperature, T-e and surface soil moisture M, retrieved on a 0.5 deg x 0.5 deg grid and further averaged over a 4 deg x 10 deg study region in the African Sahel. Also shown are National Center for Environmental Prediction (NCEP) model outputs of surface temperature, T-sfc, and soil wetness, Soil-w. The variables have been scaled to have similar dynamic ranges on the plots. The NCEP data from the NCEP Reanalysis Project are monthly averages on a 2.5 deg x 2.5 deg grid averaged over the 4 deg x 10 deg study area. Comparisons of SMMR retrievals with forecast model output show the potential of the satellite data for validating model output and monitoring long-term trends. Continuing work will extend these results to other regions to validate the retrievals more quantitatively. In preparation for the launch of AMSR, field experiments are planned in collaboration with the Global Energy and Water Cycle Experiment (GEWEX) Coordinated Enhanced Observing Period (CEOP) experiments to evaluate the satellite-derived soil moisture measurements and to demonstrate their usefulness for land surface hydrology and climate. Additional information is contained in the original.

Njoku, Eni G.↗

SMC: SCENIC Model Control

NASAs Space Communications and Navigation (SCaN) program manages three active networks: the Near Earth Network, the Space Network, and the Deep Space Network. These networks simultaneously support NASA missions and provide communications services to customers worldwide. To efficiently manage these resources and their capabilities, a team of student interns at the NASA Glenn Research Center is developing a distributed system to model the SCaN networks. Once complete, the system shall provide a platform that enables users to perform capacity modeling of current and prospective missions with finer-grained control of information between several simulation and modeling tools. This will enable the SCaN program to access a holistic view of its networks and simulate the effects of modifications in order to provide NASA with decisional information. The development of this capacity modeling system is managed by NASAs Strategic Center for Education, Networking, Integration, and Communication (SCENIC). Three primary third-party software tools offer their unique abilities in different stages of the simulation process. MagicDraw provides UMLSysML modeling, AGIs Systems Tool Kit simulates the physical transmission parameters and de-conflicts scheduled communication, and Riverbed Modeler (formerly OPNET) simulates communication protocols and packet-based networking. SCENIC developers are building custom software extensions to integrate these components in an end-to-end space communications modeling platform. A central control module acts as the hub for report-based messaging between client wrappers. Backend databases provide information related to mission parameters and ground station configurations, while the end user defines scenario-specific attributes for the model. The eight SCENIC interns are working under the direction of their mentors to complete an initial version of this capacity modeling system during the summer of 2015. The intern team is composed of four students in Computer Science, two in Computer Engineering, one in Electrical Engineering, and one studying Space Systems Engineering.

Simulation↗

The Staring OBservations of the Atmosphere (SOBA) Mission Concept

The Staring OBservations of the Atmosphere (SOBA) Mission is a concept that was developed and matured under the guidance of the NASA Ames Project EXellence (APEX) program. If funded, it will provide an unprecedented opportunity to improve ash transport forecasts and climate model simulations associated with volcanic eruptions. NASA and National science objectives require a better understanding of volcanic aerosol and trace gas emissions, transport, chemical transformation, and deposition, since they impact Earth's climate and atmospheric composition, human health, and aviation safety. Natural hazards such as the 2010 eruption of the Eyjafjallajökull volcano in Iceland demonstrated how existing remote-sensing assets were inadequate for individual volcanic event monitoring. During this eruption, available instruments were unable to provide data necessary to initialize volcanic plume transport models so that they could accurately predict the quantity and location of volcanic ash. As a result, thousands of flights around the world were grounded unnecessarily, at great expense. Volcanoes can also play a large role in regulation of the Earth's climate, so SOBA observations will also be used to evaluate and improve volcanic aerosol and trace gas simulation in chemical transport models (CTMs) and global climate models (GCMs). We propose the development of an airborne remote sensing concept and field campaign that will respond to an eruption and provide near real time observations of a volcanic plume, specifically ash injection height, transport, aerosol microphysical physical properties, and the location and concentration of sulfur dioxide (SO2) (sulfate (SO42-) aerosol precursor). This airborne system will utilize a depolarization sensitive, downward looking Light Detection And Ranging (lidar) instrument and an ultraviolet (UV) imaging spectrometer, and will provide data to be ingested by volcanic ash advisory models. Furthermore, the lessons learned in the development of this system could eventually guide regular deployment of similar systems by NASA or other government agencies.

Aerosols↗

Verifying Implementation of the Dryden Turbulence Model and MIL-F-8785 Gust Gradient

Turbulence modeling in human-in-the-loop simulation is important to assessing aircraft handling qualities and pilot performance and to provide additional realism for pilot training. In the simulation community, the Dryden turbulence spectra is a popular choice for modeling the linear turbulent gusts because its rational form is efficiently reproduced by passing white noise through linear filters. The MIL-F-8785 gust gradients similarly use additional linear filters to model the gradient of the turbulent gust over the wing, and it represents the gust gradients as perturbations to the air-relative rotational rates. The Cockpit Motion Facility at NASA Langley Research Center (LaRC) models continuous random turbulence using the Dryden one-dimensional spectra and MIL-F-8785 gust gradient. The facility recently reviewed and updated its verification of these models as part of an initiative to improve motion cueing under turbulence. This exercise introduced improved methods for verifying the turbulence models and led to rediscovery of model assumptions that informed improvements to implementation.

Madden, Michael M.↗

PolarMERRA: A collaborative effort for NASA’s polar reanalysis

The polar atmosphere is often poorly represented in global models, adversely impacting analyses in these regions and forecasts on much larger scales. PolarMERRA is an initiative that brings together the modeling expertise from NASA’s global modeling and assimilation office (GMAO) and observational expertise from NASA’s Cryospheric sciences lab to assess and improve the representation of polar processes in the NASA’s modelling system in preparation for upcoming and next generation reanalyses. Model improvements in resolution, physics, and data assimilation have been guided by assessment of model skill in the polar regions. Preliminary results from a novel validation framework indicate a reduction in near surface air temperature, turbulent fluxes, and precipitation biases.

Chelsea Parker↗

Modeling of Ultrasonically Generated Liquid-Liquid Dispersions During Controlled Directional Solidification

There are innumerable two-component systems in which two very different liquid phases co-exist in equilibrium over a range of temperature and composition, e.g., oil and water, salt fluxes and solders, aluminum and lead. Often it is of practical concern to fabricate a solid component consisting of a uniform dispersion of one phase in the other. Unfortunately, uniform microstructural development during solidification of two immiscible liquids is hampered by inherent, often large, density differences between the phases that lead to severe segregation. Uniformity is also compromised by preferential wetting and coalescence phenomena. It is, however, well known that ultrasonic energy can initiate and maintain a fine liquid-liquid dispersion. The work presented here extends that observation by application of ultrasonic energy to promote uniform phase incorporation during controlled directional solidification. To this end experiments with the transparent organic, immiscible, succinonitrile-glycerol system were conducted and the numerous processing parameters associated with this technique were evaluated in view of optimizing dispersion uniformity. In view of the initial experimental results a model that predicts the dispersed liquid droplet size as a function of material properties, sample geometry, and applied energy has been developed. In the mathematical model we consider the ultrasonic field in an experimental ampoule of length L and diameter D induced by a probe having a vibration frequency of f=2OKhz (circular frequency omega = 2 pi f). The amplitude is adjustable from A=65 to 13Omicrons. The probe tip diameter is d, the liquid has a density of p, in which the speed of sound and surface tension are, respectively, c and sigma. The mathematical model and numerical investigation for the experiments [1] is done using the following assumptions: (i) The droplet size is small in comparison to the sound wave length; (ii) The forces between droplets are neglected (relative concentration is small); (iii) The droplet is stable if the kinetic energy, E(sub K), of the liquid motion due to ultrasonic field influence is less then the binding energy, E(sub S), due to the surface tension (it is easy to show that the surface energy of two droplets resulting from one is larger by about a factor of two.); (iv) The stability limit is characterized by E(sub S) to approx. E(sub K).

Grugel, R. N.↗

Partial Overhaul and Initial Parallel Optimization of KINETICS, a Coupled Dynamics and Chemistry Atmosphere Model

KINETICS is a coupled dynamics and chemistry atmosphere model that is data intensive and computationally demanding. The potential performance gain from using a supercomputer motivates the adaptation from a serial version to a parallelized one. Although the initial parallelization had been done, bottlenecks caused by an abundance of communication calls between processors led to an unfavorable drop in performance. Before starting on the parallel optimization process, a partial overhaul was required because a large emphasis was placed on streamlining the code for user convenience and revising the program to accommodate the new supercomputers at Caltech and JPL. After the first round of optimizations, the partial runtime was reduced by a factor of 23; however, performance gains are dependent on the size of the data, the number of processors requested, and the computer used.

KINETICS↗

Monte Carlo Modeling of Lunar Megaregolith Development

A Monte Carlo model of lunar megaregolith development is proposed. Minimum megaregolith depths are obtained in the preliminary model. The model impacts a gridded, initially flat surface area randomly, with points not constrained to grid intersections. Information is recorded at grid points on evolving topography, deepest depth impacted to, and number of hits at that location. Two types of lunar surfaces ar simulated, maria and highlands. Crater production functions are after run to densities that are multiples of observed lunar crater densities for craters of 8 km in diameter and larger. all craters are considered to have formed in a gravity scaled regime. Edge effect problems are corrected in this model by allowing craters forming outside the grid area or along its edges to effect the gridded area.

Cashore, J.↗