Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster computing”

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 901 records · Page 50

Scattering and Radiative Properties of Morphologically Complex Carbonaceous Aerosols: A Systematic Modeling Study

This paper provides a thorough modeling-based overview of the scattering and radiative properties of a wide variety of morphologically complex carbonaceous aerosols. Using the numerically-exact superposition T-matrix method, we examine the absorption enhancement, absorption Angstroem exponent (AAE), backscattering linear depolarization ratio (LDR), and scattering matrix elements of black-carbon aerosols with 11 different model morphologies ranging from bare soot to completely embedded soot-sulfate and soot-brown carbon mixtures. Our size-averaged results show that fluffy soot particles absorb more light than compact bare-soot clusters. For the same amount of absorbing material, the absorption cross section of internally mixed soot can be more than twice that of bare soot. Absorption increases as soot accumulates more coating material and can become saturated. The absorption enhancement is affected by particle size, morphology, wavelength, and the amount of coating. We refute the conventional belief that all carbonaceous aerosols have AAEs close to 1.0. Although LDRs caused by bare soot and certain carbonaceous particles are rather weak, LDRs generated by other soot-containing aerosols can reproduce strong depolarization measured by Burton et al. for aged smoke. We demonstrate that multi-wavelength LDR measurements can be used to identify the presence of morphologically complex carbonaceous particles, although additional observations can be needed for full characterization. Our results show that optical constants of the host/coating material can significantly influence the scattering and absorption properties of soot-containing aerosols to the extent of changing the sign of linear polarization. We conclude that for an accurate estimate of black-carbon radiative forcing, one must take into account the complex morphologies of carbonaceous aerosols in remote sensing studies as well as in atmospheric radiation computations.

carbonaceous aerosols; scattering matrix; polariza↗

MPIRUN: A Portable Loader for Multidisciplinary and Multi-Zonal Applications

Multidisciplinary and multi-zonal applications are an important class of applications in the area of Computational Aerosciences. In these codes, two or more distinct parallel programs or copies of a single program are utilized to model a single problem. To support such applications, it is common to use a programming model where a program is divided into several single program multiple data stream (SPMD) applications, each of which solves the equations for a single physical discipline or grid zone. These SPMD applications are then bound together to form a single multidisciplinary or multi-zonal program in which the constituent parts communicate via point-to-point message passing routines. One method for implementing the message passing portion of these codes is with the new Message Passing Interface (MPI) standard. Unfortunately, this standard only specifies the message passing portion of an application, but does not specify any portable mechanisms for loading an application. MPIRUN was developed to provide a portable means for loading MPI programs, and was specifically targeted at multidisciplinary and multi-zonal applications. Programs using MPIRUN for loading and MPI for message passing are then portable between all machines supported by MPIRUN. MPIRUN is currently implemented for the Intel iPSC/860, TMC CM5, IBM SP-1 and SP-2, Intel Paragon, and workstation clusters. Further, MPIRUN is designed to be simple enough to port easily to any system supporting MPI.

Fineberg, Samuel A.↗

System for Automated Calibration of Vector Modulators

Vector modulators are used to impose baseband modulation on RF signals, but non-ideal behavior limits the overall performance. The non-ideal behavior of the vector modulator is compensated using data collected with the use of an automated test system driven by a LabVIEW program that systematically applies thousands of control-signal values to the device under test and collects RF measurement data. The technology innovation automates several steps in the process. First, an automated test system, using computer controlled digital-to-analog converters (DACs) and a computer-controlled vector network analyzer (VNA) systematically can apply different I and Q signals (which represent the complex number by which the RF signal is multiplied) to the vector modulator under test (VMUT), while measuring the RF performance specifically, gain and phase. The automated test system uses the LabVIEW software to control the test equipment, collect the data, and write it to a file. The input to the Lab - VIEW program is either user-input for systematic variation, or is provided in a file containing specific test values that should be fed to the VMUT. The output file contains both the control signals and the measured data. The second step is to post-process the file to determine the correction functions as needed. The result of the entire process is a tabular representation, which allows translation of a desired I/Q value to the required analog control signals to produce a particular RF behavior. In some applications, corrected performance is needed only for a limited range. If the vector modulator is being used as a phase shifter, there is only a need to correct I and Q values that represent points on a circle, not the entire plane. This innovation has been used to calibrate 2-GHz MMIC (monolithic microwave integrated circuit) vector modulators in the High EIRP Cluster Array project (EIRP is high effective isotropic radiated power). These calibrations were then used to create correction tables to allow the commanding of the phase shift in each of four channels used as a phased array for beam steering of a Ka-band (32-GHz) signal. The system also was the basis of a breadboard electronic beam steering system. In this breadboard, the goal was not to make systematic measurements of the properties of a vector modulator, but to drive the breadboard with a series of test patterns varying in phase and amplitude. This is essentially the same calibration process, but with the difference that the data collection process is oriented toward collecting breadboard performance, rather than the measurement of output from a network analyzer.

Lux, James↗

Improve Data Mining and Knowledge Discovery Through the Use of MatLab

Data mining is widely used to mine business, engineering, and scientific data. Data mining uses pattern based queries, searches, or other analyses of one or more electronic databases/datasets in order to discover or locate a predictive pattern or anomaly indicative of system failure, criminal or terrorist activity, etc. There are various algorithms, techniques and methods used to mine data; including neural networks, genetic algorithms, decision trees, nearest neighbor method, rule induction association analysis, slice and dice, segmentation, and clustering. These algorithms, techniques and methods used to detect patterns in a dataset, have been used in the development of numerous open source and commercially available products and technology for data mining. Data mining is best realized when latent information in a large quantity of data stored is discovered. No one technique solves all data mining problems; challenges are to select algorithms or methods appropriate to strengthen data/text mining and trending within given datasets. In recent years, throughout industry, academia and government agencies, thousands of data systems have been designed and tailored to serve specific engineering and business needs. Many of these systems use databases with relational algebra and structured query language to categorize and retrieve data. In these systems, data analyses are limited and require prior explicit knowledge of metadata and database relations; lacking exploratory data mining and discoveries of latent information. This presentation introduces MatLab(R) (MATrix LABoratory), an engineering and scientific data analyses tool to perform data mining. MatLab was originally intended to perform purely numerical calculations (a glorified calculator). Now, in addition to having hundreds of mathematical functions, it is a programming language with hundreds built in standard functions and numerous available toolboxes. MatLab's ease of data processing, visualization and its enormous availability of built in functionalities and toolboxes make it suitable to perform numerical computations and simulations as well as a data mining tool. Engineers and scientists can take advantage of the readily available functions/toolboxes to gain wider insight in their perspective data mining experiments.

Shaykhian, Gholam Ali↗

Improve Data Mining and Knowledge Discovery through the use of MatLab

Data mining is widely used to mine business, engineering, and scientific data. Data mining uses pattern based queries, searches, or other analyses of one or more electronic databases/datasets in order to discover or locate a predictive pattern or anomaly indicative of system failure, criminal or terrorist activity, etc. There are various algorithms, techniques and methods used to mine data; including neural networks, genetic algorithms, decision trees, nearest neighbor method, rule induction association analysis, slice and dice, segmentation, and clustering. These algorithms, techniques and methods used to detect patterns in a dataset, have been used in the development of numerous open source and commercially available products and technology for data mining. Data mining is best realized when latent information in a large quantity of data stored is discovered. No one technique solves all data mining problems; challenges are to select algorithms or methods appropriate to strengthen data/text mining and trending within given datasets. In recent years, throughout industry, academia and government agencies, thousands of data systems have been designed and tailored to serve specific engineering and business needs. Many of these systems use databases with relational algebra and structured query language to categorize and retrieve data. In these systems, data analyses are limited and require prior explicit knowledge of metadata and database relations; lacking exploratory data mining and discoveries of latent information. This presentation introduces MatLab(TradeMark)(MATrix LABoratory), an engineering and scientific data analyses tool to perform data mining. MatLab was originally intended to perform purely numerical calculations (a glorified calculator). Now, in addition to having hundreds of mathematical functions, it is a programming language with hundreds built in standard functions and numerous available toolboxes. MatLab's ease of data processing, visualization and its enormous availability of built in functionalities and toolboxes make it suitable to perform numerical computations and simulations as well as a data mining tool. Engineers and scientists can take advantage of the readily available functions/toolboxes to gain wider insight in their perspective data mining experiments.

Shaykahian, Gholan Ali↗

DeepSAT: A Deep Learning Approach to Tree-Cover Delineation in 1-m NAIP Imagery for the Continental United States

High resolution tree cover classification maps are needed to increase the accuracy of current land ecosystem and climate model outputs. Limited studies are in place that demonstrates the state-of-the-art in deriving very high resolution (VHR) tree cover products. In addition, most methods heavily rely on commercial softwares that are difficult to scale given the region of study (e.g. continents to globe). Complexities in present approaches relate to (a) scalability of the algorithm, (b) large image data processing (compute and memory intensive), (c) computational cost, (d) massively parallel architecture, and (e) machine learning automation. In addition, VHR satellite datasets are of the order of terabytes and features extracted from these datasets are of the order of petabytes. In our present study, we have acquired the National Agriculture Imagery Program (NAIP) dataset for the Continental United States at a spatial resolution of 1-m. This data comes as image tiles (a total of quarter million image scenes with ~60 million pixels) and has a total size of ~65 terabytes for a single acquisition. Features extracted from the entire dataset would amount to ~8-10 petabytes. In our proposed approach, we have implemented a novel semi-automated machine learning algorithm rooted on the principles of "deep learning" to delineate the percentage of tree cover. Using the NASA Earth Exchange (NEX) initiative, we have developed an end-to-end architecture by integrating a segmentation module based on Statistical Region Merging, a classification algorithm using Deep Belief Network and a structured prediction algorithm using Conditional Random Fields to integrate the results from the segmentation and classification modules to create per-pixel class labels. The training process is scaled up using the power of GPUs and the prediction is scaled to quarter million NAIP tiles spanning the whole of Continental United States using the NEX HPC supercomputing cluster. An initial pilot over the state of California spanning a total of 11,095 NAIP tiles covering a total geographical area of 163,696 sq. miles has produced true positive rates of around 88 percent for fragmented forests and 74 percent for urban tree cover areas, with false positive rates lower than 2 percent for both landscapes.

Imagery↗

MODIS NDVI Change Detection Techniques and Products Used in the Near Real Time ForWarn System for Detecting, Monitoring, and Analyzing Regional Forest Disturbances

This presentation discusses MODIS NDVI change detection methods and products used in the ForWarn Early Warning System (EWS) for near real time (NRT) recognition and tracking of regionally evident forest disturbances throughout the conterminous US (CONUS). The latter has provided NRT forest change products to the forest health protection community since 2010, using temporally processed MODIS Aqua and Terra NDVI time series data to currently compute and post 6 different forest change products for CONUS every 8 days. Multiple change products are required to improve detectability and to more fully assess the nature of apparent disturbances. Each type of forest change product reports per pixel percent change in NDVI for a given 24 day interval, comparing current versus a given historical baseline NDVI. EMODIS 7 day expedited MODIS MOD13 data are used to obtain current and historical NDVIs, respectively. Historical NDVI data is processed with Time Series Product Tool (TSPT); and 2) the Phenological Parameters Estimation Tool (PPET) software. While each change products employ maximum value compositing (MVC) of NDVI, the design of specific products primarily differs in terms of the historical baseline. The three main change products use either 1, 3, or all previous years of MVC NDVI as a baseline. Another product uses an Adaptive Length Compositing (ALC) version of MVC to derive an alternative current NDVI that is the freshest quality NDVI as opposed to merely the MVC NDVI across a 24 day time frame. The ALC approach can improve detection speed by 8 to 16 days. ForWarn also includes 2 change products that improve detectability of forest disturbances in lieu of climatic fluctuations, especially in the spring and fall. One compares current MVC NDVI to the zonal maximum under the curve NDVI per pheno-region cluster class, considering all previous years in the MODIS record. The other compares current maximum NDVI to the mean of maximum NDVI for all previous MODIS years.

Spruce, Joseph P.↗

Efficiency Improvements to the Displacement Based Multilevel Structural Optimization Algorithm

Multilevel Structural Optimization (MSO) continues to be an area of research interest in engineering optimization. In the present project, the weight optimization of beams and trusses using Displacement based Multilevel Structural Optimization (DMSO), a member of the MSO set of methodologies, is investigated. In the DMSO approach, the optimization task is subdivided into a single system and multiple subsystems level optimizations. The system level optimization minimizes the load unbalance resulting from the use of displacement functions to approximate the structural displacements. The function coefficients are then the design variables. Alternately, the system level optimization can be solved using the displacements themselves as design variables, as was shown in previous research. Both approaches ensure that the calculated loads match the applied loads. In the subsystems level, the weight of the structure is minimized using the element dimensions as design variables. The approach is expected to be very efficient for large structures, since parallel computing can be utilized in the different levels of the problem. In this paper, the method is applied to a one-dimensional beam and a large three-dimensional truss. The beam was tested to study possible simplifications to the system level optimization. In previous research, polynomials were used to approximate the global nodal displacements. The number of coefficients of the polynomials equally matched the number of degrees of freedom of the problem. Here it was desired to see if it is possible to only match a subset of the degrees of freedom in the system level. This would lead to a simplification of the system level, with a resulting increase in overall efficiency. However, the methods tested for this type of system level simplification did not yield positive results. The large truss was utilized to test further improvements in the efficiency of DMSO. In previous work, parallel processing was applied to the subsystems level, where the derivative verification feature of the optimizer NPSOL had been utilized in the optimizations. This resulted in large runtimes. In this paper, the optimizations were repeated without using the derivative verification, and the results are compared to those from the previous work. Also, the optimizations were run on both, a network of SUN workstations using the MPICH implementation of the Message Passing Interface (MPI) and on the faster Beowulf cluster at ICASE, NASA Langley Research Center, using the LAM implementation of UP]. The results on both systems were consistent and showed that it is not necessary to verify the derivatives and that this gives a large increase in efficiency of the DMSO algorithm.

Plunkett, C. L.↗

Gene Expression Dynamics Inspector (GEDI): for integrative analysis of expression profiles

Genome-wide expression profiles contain global patterns that evade visual detection in current gene clustering analysis. Here, a Gene Expression Dynamics Inspector (GEDI) is described that uses self-organizing maps to translate high-dimensional expression profiles of time courses or sample classes into animated, coherent and robust mosaics images. GEDI facilitates identification of interesting patterns of molecular activity simultaneously across gene, time and sample space without prior assumption of any structure in the data, and then permits the user to retrieve genes of interest. Important changes in genome-wide activities may be quickly identified based on 'Gestalt' recognition and hence, GEDI may be especially useful for non-specialist end users, such as physicians. AVAILABILITY: GEDI v1.0 is written in Matlab, and binary Matlab.dll files which require Matlab to run can be downloaded for free by academic institutions at http://www.chip.org/~ge/gedihome.html Supplementary information: http://www.chip.org/~ge/gedihome.html.

Database Management Systems↗

Reactive Collision Avoidance Algorithm

The reactive collision avoidance (RCA) algorithm allows a spacecraft to find a fuel-optimal trajectory for avoiding an arbitrary number of colliding spacecraft in real time while accounting for acceleration limits. In addition to spacecraft, the technology can be used for vehicles that can accelerate in any direction, such as helicopters and submersibles. In contrast to existing, passive algorithms that simultaneously design trajectories for a cluster of vehicles working to achieve a common goal, RCA is implemented onboard spacecraft only when an imminent collision is detected, and then plans a collision avoidance maneuver for only that host vehicle, thus preventing a collision in an off-nominal situation for which passive algorithms cannot. An example scenario for such a situation might be when a spacecraft in the cluster is approaching another one, but enters safe mode and begins to drift. Functionally, the RCA detects colliding spacecraft, plans an evasion trajectory by solving the Evasion Trajectory Problem (ETP), and then recovers after the collision is avoided. A direct optimization approach was used to develop the algorithm so it can run in real time. In this innovation, a parameterized class of avoidance trajectories is specified, and then the optimal trajectory is found by searching over the parameters. The class of trajectories is selected as bang-off-bang as motivated by optimal control theory. That is, an avoiding spacecraft first applies full acceleration in a constant direction, then coasts, and finally applies full acceleration to stop. The parameter optimization problem can be solved offline and stored as a look-up table of values. Using a look-up table allows the algorithm to run in real time. Given a colliding spacecraft, the properties of the collision geometry serve as indices of the look-up table that gives the optimal trajectory. For multiple colliding spacecraft, the set of trajectories that avoid all spacecraft is rapidly searched on-line. The optimal avoidance trajectory is implemented as a receding-horizon model predictive control law. Therefore, at each time step, the optimal avoidance trajectory is found and the first time step of its acceleration is applied. At the next time step of the control computer, the problem is re-solved and the new first time step is again applied. This continual updating allows the RCA algorithm to adapt to a colliding spacecraft that is making erratic course changes.

Scharf, Daniel↗

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.

NASA Discipline Exobiology↗

Numerical experiments with model monophyletic and paraphyletic taxa

The problem of how accurately paraphyletic taxa versus monophyletic (i.e., holophyletic) groups (clades) capture underlying species patterns of diversity and extinction is explored with Monte Carlo simulations. Phylogenies are modeled as stochastic trees. Paraphyletic taxa are defined in an arbitrary manner by randomly choosing progenitors and clustering all descendants not belonging to other taxa. These taxa are then examined to determine which are clades, and the remaining paraphyletic groups are dissected to discover monophyletic subgroups. Comparisons of diversity patterns and extinction rates between modeled taxa and lineages indicate that paraphyletic groups can adequately capture lineage information under a variety of conditions of diversification and mass extinction. This suggests that these groups constitute more than mere "taxonomic noise" in this context. But, strictly monophyletic groups perform somewhat better, especially with regard to mass extinctions. However, when low levels of paleontologic sampling are simulated, the veracity of clades deteriorates, especially with respect to diversity, and modeled paraphyletic taxa often capture more information about underlying lineages. Thus, for studies of diversity and taxic evolution in the fossil record, traditional paleontologic genera and families need not be rejected in favor of cladistically-defined taxa.

Non-NASA Center↗

Aeroacoustic Computations of a Transonic Truss-Braced Wing Aircraft: Part 2 – Acoustic Signature and Noise Source Identification

High-fidelity, time-dependent simulations of a Boeing-designed, transonic, truss-braced-wing aircraft in cruise (clean) and landing configurations are leveraged to generate synthetic microphone-phased-array data for airframe noise prediction and assessment. These data sets are used to compute source localization (beamform) maps to determine the location and strength of primary and secondary airframe noise sources associated with this unique configuration. The synthetic phased-array implementation mimics the setup of a flight test. As this study is ongoing, preliminary integrated far-field spectra for the cruise configuration obtained at multiple spatial resolutions revealed significant tonal content that lacked convergence with increased resolution. The origin of several of these tones and their unusual convergence behavior was traced to the larger-than-normal trailing-edge thickness of the “as-tested” cruise model being simulated. Reducing the trailing-edge thickness to more realistic values eliminated most of the tones at low to moderate frequencies and improved spectrum convergence significantly. Applying lessons learned from the cruise simulations, several modifications to the geometry of the landing configuration were made and are described in this work. Results from permeable and solid Ffowcs-Williams and Hawkings surfaces at two different spatial resolutions (coarse and medium) are used to illustrate the major noise sources and determine convergence of the CLEAN integrated noise levels for the entire aircraft as well as major subcomponents. We demonstrate that the low-frequency content of the far-field spectrum is dominated by noise generated from the main landing gear, while the medium- and high-frequency content is dominated by the wing-leading-edge Krueger flaps. Since analysis of the acoustic maps for the landing configuration revealed several clusters of multiple sources along the wing leading edge, “high resolution” processing of the array data was used to distinguish more accurately the locations of sources.

Transonic Truss-Braced Wing↗

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗

A Stochastic Model of Space-Time Variability of Tropical Rainfall: I. Statistics of Spatial Averages

Global maps of rainfall are of great importance in connection with modeling of the earth s climate. Comparison between the maps of rainfall predicted by computer-generated climate models with observation provides a sensitive test for these models. To make such a comparison, one typically needs the total precipitation amount over a large area, which could be hundreds of kilometers in size over extended periods of time of order days or months. This presents a difficult problem since rain varies greatly from place to place as well as in time. Remote sensing methods using ground radar or satellites detect rain over a large area by essentially taking a series of snapshots at infrequent intervals and indirectly deriving the average rain intensity within a collection of pixels , usually several kilometers in size. They measure area average of rain at a particular instant. Rain gauges, on the other hand, record rain accumulation continuously in time but only over a very small area tens of centimeters across, say, the size of a dinner plate. They measure only a time average at a single location. In making use of either method one needs to fill in the gaps in the observation - either the gaps in the area covered or the gaps in time of observation. This involves using statistical models to obtain information about the rain that is missed from what is actually detected. This paper investigates such a statistical model and validates it with rain data collected over the tropical Western Pacific from ship borne radars during TOGA COARE (Tropical Oceans Global Atmosphere Coupled Ocean-Atmosphere Response Experiment). The model incorporates a number of commonly observed features of rain. While rain varies rapidly with location and time, the variability diminishes when averaged over larger areas or longer periods of time. Moreover, rain is patchy in nature - at any instant on the average only a certain fraction of the observed pixels contain rain. The fraction of area covered by rain decreases, as the size of a pixel becomes smaller. This means that within what looks like a patch of rainy area in a coarse resolution view with larger pixel size, one finds clusters of rainy and dry patches when viewed on a finer scale. The model makes definite predictions about how these and other related statistics depend on the pixel size. These predictions were found to agree well with data. In a subsequent second part of the work we plan to test the model with rain gauge data collected during the TRMM (Tropical Rainfall Measuring Mission) ground validation campaign.

Kundu, Prasun K.↗

Development of a Ground Multi-Mission Low-Cost Optical Terminal (LCOT) for Free-Space Optical Communications

Once confined to the realm of laboratory experiments and theoretical papers, space-based laser communications (lasercomm) are on the verge of achieving mainstream status. Organizations from Facebook to NASA, and missions from cubesats to Orion are employing lasercomm to achieve gigabit communication speeds at mass and power requirements lower than that of traditional radio frequency (RF) methods. Since first demonstrating free-space optical communications services with Lunar Laser Communications Demonstration (LLCD) in 2013, NASA has invested in developing optical communications technologies and capabilities to enhancing its space communications networks. Along with evolving optical space terminals, NASA is also developing lasercomm grounds stations capable of meeting the rapidly increasing data volume demands of upcoming missions, from low-earth to lunar orbits and beyond and integrating these advanced capabilities into its Near-Space and Deep Space Networks. To meet this emerging need, the Low-Cost Optical Terminal (LCOT) project at NASA’s Goddard Space Flight Center (GSFC) is designing, building and validating a prototype for a flexible, multi-mission, and economical optical ground terminal that could be used as a blueprint for a global network of optical ground stations, capable of supporting a wide variety of missions. To date a major impediment to widespread adoption of laser communication has been the lack of an existing ground network infrastructure. A mission that wishes to take advantage of laser communication not only needs to invest in an optical space terminal, but it must also finance the creation of ground terminals to receive the downlink signal. This adds significant additional cost. Missions that do decide to incur the cost of financing a network of ground terminals end up building highly specialized optical receivers that are operable as receivers for that specific mission only. Significant Non-Recurring Engineering (NRE) cost is invested to build highly specialized one-of-a-kind ground terminals that go into storage after that particular mission is over. This is not an economical approach and does nothing to grow the number of optical ground stations available to future missions. In essence each mission that wants to take advantage of the benefits of lasercom has to start from scratch to provide a ground terminal network to support it. As long as this is the case, the cost for using laser communications will be too high for most missions to consider. LCOT intends to close this gap in technology by designing and developing a standard optical ground terminal design that is flexible enough to serve as a receiver for a wide range of future missions – a ground terminal that can be quickly reconfigured to receive downlinks at different wavelengths using different signal formats. Not only does LCOT have the industry-building objectives of utilizing commercial-off-the-shelf (COTS) components to the maximum extent possible, but also spurs the commercial development of other necessary lasercomm components not currently offered by industry. Finally, LCOT will give NASA scientists and engineers a facility where they can gain real-world experience with optical communications. It will give engineers a cost-effective way to try out new concepts and processes by providing the infrastructure for such testing. In this way it is hoped LCOT will serve as a stimulus for innovation in optical communications and speed its widespread adoption by future missions. The LCOT is comprised of five subsystems: Free-Space Optical, Transceiver, Amplifier, Monitor and Control, and Observatory Infrastructure. In August 2021, the LCOT team installed a 70 cm telescope, developed by Planewave Instruments that was optimized for optical communications. Free-Space Optical subsystem comprises of the telescope and its associated hardware, including a transmitter optical assembly, wide field cameras, two optical benches, and an adaptive optics subsystem. The transmit optical assembly, a unique concept design, is a cluster of four functionally independent transmit subassemblies located on the receive telescope. In addition to receiving optical signals and directing the expanded beam with high precision to the space terminal, it also performs tracking functions. The transmit optical assembly will support operations from Low Earth Orbit (LEO) through lunar and will be used as a template for industry manufacturing. The Optical Infrastructure subsystem is responsible for providing environmentally controlled shelters for LCOT equipment and various other systems. To maintain the safety and proper functionality of the telescope, a 16 ft Astrohaven clamshell dome procured which provides all-sky coverage without the need to rotate the dome. Additionally, Atmospheric Monitoring Assembly (AMA) will be part of optical infrastructure subsystem to ensure accurate performance of the LCOT. Like existing optical ground stations, LCOT will measure standard weather station parameters, infrared all sky image of cloud cover, and cloud height. LCOT, however, adds requirements for measuring night time seeing and, in the future, daytime seeing. Unlike other optical ground terminals, the LCOT is transceiver agnostic; user transceivers may be duplex transceivers, standalone receivers, or standalone transmitters with or without acquisition beacon functionality. As such, the LCOT project accommodates testing with external customer transceivers in a flexible manor, further complimenting its intended multi-mission goals. Another unique component of LCOT is the use of a new amplifier technology – the Very Large Mode Area (VLMA) amplifiers. This new technology allows LCOT to avoid the issues faced by previous laser communications ground terminals, gives users more flexibility and modular capability, and is capable of reaching an order of magnitude higher peak power than traditional High Power Optical Amplifiers (HPOA). One drawback of the VLMA HPOA approach is that the amplified light is output into free-space. The solution developed by LCOT is an optics train that couples the output of the VLMA amplifier into a short fiber for transport to the transmit telescopes with high efficiency. Like many of the LCOT components, a set of detailed manufacturing drawings have been created for the optics train to allow any machine shop with a multi-axis Computer Numerical Control (CNC) machine to fabricate the piece parts from commonly available materials. In line with the goals of LCOT, the monitor and control functions are developed as a modular and flexible system with the ability to support future hardware or algorithm changes, minimizing disruptions. A main priority of development in the Monitor and Control Subsystem (MCS) is the safety monitor system.

laser communication↗

Development of a Ground Multi-mission Low Cost Optical Terminal(LCOT) for Free-Space Optical Communication

Once confined to the realm of laboratory experiments and theoretical papers, space-based laser communications (lasercomm) are on the verge of achieving mainstream status. Organizations from Facebook to NASA, and missions from cubesats to Orion are employing lasercomm to achieve gigabit communication speeds at mass and power requirements lower than that of traditional radio frequency (RF) methods. Since first demonstrating free-space optical communications services with Lunar Laser Communications Demonstration (LLCD) in 2013, NASA has invested in developing optical communications technologies and capabilities to enhancing its space communications networks. Along with evolving optical space terminals, NASA is also developing lasercomm grounds stations capable of meeting the rapidly increasing data volume demands of upcoming missions, from low-earth to lunar orbits and beyond and integrating these advanced capabilities into its Near-Space and Deep Space Networks. To meet this emerging need, the Low-Cost Optical Terminal (LCOT) project at NASA’s Goddard Space Flight Center (GSFC) is designing, building and validating a prototype for a flexible, multi-mission, and economical optical ground terminal that could be used as a blueprint for a global network of optical ground stations, capable of supporting a wide variety of missions. To date a major impediment to widespread adoption of laser communication has been the lack of an existing ground network infrastructure. A mission that wishes to take advantage of laser communication not only needs to invest in an optical space terminal, but it must also finance the creation of ground terminals to receive the downlink signal. This adds significant additional cost. Missions that do decide to incur the cost of financing a network of ground terminals end up building highly specialized optical receivers that are operable as receivers for that specific mission only. Significant Non-Recurring Engineering (NRE) cost is invested to build highly specialized one-of-a-kind ground terminals that go into storage after that particular mission is over. This is not an economical approach and does nothing to grow the number of optical ground stations available to future missions. In essence each mission that wants to take advantage of the benefits of lasercom has to start from scratch to provide a ground terminal network to support it. As long as this is the case, the cost for using laser communications will be too high for most missions to consider. LCOT intends to close this gap in technology by designing and developing a standard optical ground terminal design that is flexible enough to serve as a receiver for a wide range of future missions – a ground terminal that can be quickly reconfigured to receive downlinks at different wavelengths using different signal formats. Not only does LCOT have the industry-building objectives of utilizing commercial-off-the-shelf (COTS) components to the maximum extent possible, but also spurs the commercial development of other necessary lasercomm components not currently offered by industry. Finally, LCOT will give NASA scientists and engineers a facility where they can gain real-world experience with optical communications. It will give engineers a cost-effective way to try out new concepts and processes by providing the infrastructure for such testing. In this way it is hoped LCOT will serve as a stimulus for innovation in optical communications and speed its widespread adoption by future missions. The LCOT is comprised of five subsystems: Free-Space Optical, Transceiver, Amplifier, Monitor and Control, and Observatory Infrastructure. In August 2021, the LCOT team installed a 70 cm telescope, developed by Planewave Instruments that was optimized for optical communications. Free-Space Optical subsystem comprises of the telescope and its associated hardware, including a transmitter optical assembly, wide field cameras, two optical benches, and an adaptive optics subsystem. The transmit optical assembly, a unique concept design, is a cluster of four functionally independent transmit subassemblies located on the receive telescope. In addition to receiving optical signals and directing the expanded beam with high precision to the space terminal, it also performs tracking functions. The transmit optical assembly will support operations from Low Earth Orbit (LEO) through lunar and will be used as a template for industry manufacturing. The Optical Infrastructure subsystem is responsible for providing environmentally controlled shelters for LCOT equipment and various other systems. To maintain the safety and proper functionality of the telescope, a 16 ft Astrohaven clamshell dome procured which provides all-sky coverage without the need to rotate the dome. Additionally, Atmospheric Monitoring Assembly (AMA) will be part of optical infrastructure subsystem to ensure accurate performance of the LCOT. Like existing optical ground stations, LCOT will measure standard weather station parameters, infrared all sky image of cloud cover, and cloud height. LCOT, however, adds requirements for measuring night time seeing and, in the future, daytime seeing. Unlike other optical ground terminals, the LCOT is transceiver agnostic; user transceivers may be duplex transceivers, standalone receivers, or standalone transmitters with or without acquisition beacon functionality. As such, the LCOT project accommodates testing with external customer transceivers in a flexible manor, further complimenting its intended multi-mission goals. Another unique component of LCOT is the use of a new amplifier technology – the Very Large Mode Area (VLMA) amplifiers. This new technology allows LCOT to avoid the issues faced by previous laser communications ground terminals, gives users more flexibility and modular capability, and is capable of reaching an order of magnitude higher peak power than traditional High Power Optical Amplifiers (HPOA). One drawback of the VLMA HPOA approach is that the amplified light is output into free-space. The solution developed by LCOT is an optics train that couples the output of the VLMA amplifier into a short fiber for transport to the transmit telescopes with high efficiency. Like many of the LCOT components, a set of detailed manufacturing drawings have been created for the optics train to allow any machine shop with a multi-axis Computer Numerical Control (CNC) machine to fabricate the piece parts from commonly available materials. In line with the goals of LCOT, the monitor and control functions are developed as a modular and flexible system with the ability to support future hardware or algorithm changes, minimizing disruptions. A main priority of development in the Monitor and Control Subsystem (MCS) is the safety monitor system.

Haleh Safavi↗

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗