SEARCH · Search NASA

Results for “First principles calculations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

196 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Interactions of Ku70/80 with Double-Strand DNA: Energetic, Dynamics, and Functional Implications

Space radiation is a proficient inducer of DNA damage leading to mutation, aberrant cell signaling, and cancer formation. Ku is among the first responding proteins in nucleus to recognize and bind the DNA double strand breaks (DSBs) whenever they are introduced. Once loaded Ku works as a scaffold to recruit other repair factors of non-homologous end joining and facilitates the following repair processes. The crystallographic study of the Ku70/80 heterodimer indicate the core structure of this protein shows virtually no conformational change after binding with DNA. To investigate the dynamical features as well as the energetic characteristics of Ku-DNA binding, we conduct multi-nanosecond molecular dynamics simulations of a modeled Ku70/80 structure and several complexes with two 24-bp DNA duplexes. Free energy calculations show significant energy differences between the complexes with Ku bound at DSBs and those with Ku associated at an internal site of a chromosome. The results also reveal detailed interactions between different nucleotides and the amino acids along the DNA-binding cradle of Ku, indicating subtle binding preference of Ku at specific DNA sequences. The covariance matrix analyses along the trajectories demonstrate the protein is stimulated to undergo correlated motions of different domains once bound to DNA ends. Additionally, principle component analyses identify these low frequency collective motions suitable for binding with and translocation along duplex DNA. It is proposed that the modification of dynamical properties of Ku upon binding with DSBs may provide a signal for the further recruitment of other repair factors such as DNA-PKcs, XLF, and XRCC4.

Hu, Shaowen↗

NASA NTRS2017

Design and Principles Enabling the Space Reference FOM

A first complete draft of the Simulation Interoperability Standards Organization (SISO) Space Reference Federation Object Model (FOM) has now been produced. This paper provides some insights into its capabilities and discusses the opportunity for reuse in other domains. The focus of this first version of the standard is execution control, time management and coordinate systems, well-known reference frames, as well as some basic support for physical entities. The biggest part of the execution control is the coordinated start-up process. This process contains a number of steps, including checking of required federates, handling of early versus late joiners, sharing of federation wide configuration data and multi-phase initialization. An additional part of Execution Control is the coordinated and synchronized transition between Run mode, Freeze mode and Shutdown. For time management, several time lines are defined, including real-time, scenario time, High Level Architecture (HLA) logical time and physical time. A strategy for mixing simulations that use different time steps is introduced, as well as an approach for finding common boundaries for fully synchronized freeze. For describing spatial information, a mechanism with a set of reference frames is specified. Each reference frame has a position and orientation related to a parent reference frame. This makes it possible for federates to perform calculations in reference frames that are convenient to them. An operation on the Moon can be performed using lunar coordinates whereas an operation on Earth can be performed using Earth coordinates. At the same time, coordinates in one reference frame have an unambiguous relationship to a coordinate in another reference frame. While the Space Reference FOM is originally being developed for Space operations, the authors believe that many parts of it can be reused for any simulation that has a focus on physical processes with one or more coordinate systems, and require high fidelity and repeatability.

Moeller, Bjoern↗

NASA NTRS2014

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing↗

NASA NTRS2015

Reliability and Probabilistic Risk Assessment - How They Play Together

PRA methodology is one of the probabilistic analysis methods that NASA brought from the nuclear industry to assess the risk of LOM, LOV and LOC for launch vehicles. PRA is a system scenario based risk assessment that uses a combination of fault trees, event trees, event sequence diagrams, and probability and statistical data to analyze the risk of a system, a process, or an activity. It is a process designed to answer three basic questions: What can go wrong? How likely is it? What is the severity of the degradation? Since 1986, NASA, along with industry partners, has conducted a number of PRA studies to predict the overall launch vehicles risks. Planning Research Corporation conducted the first of these studies in 1988. In 1995, Science Applications International Corporation (SAIC) conducted a comprehensive PRA study. In July 1996, NASA conducted a two-year study (October 1996 - September 1998) to develop a model that provided the overall Space Shuttle risk and estimates of risk changes due to proposed Space Shuttle upgrades. After the Columbia accident, NASA conducted a PRA on the Shuttle External Tank (ET) foam. This study was the most focused and extensive risk assessment that NASA has conducted in recent years. It used a dynamic, physics-based, integrated system analysis approach to understand the integrated system risk due to ET foam loss in flight. Most recently, a PRA for Ares I launch vehicle has been performed in support of the Constellation program. Reliability, on the other hand, addresses the loss of functions. In a broader sense, reliability engineering is a discipline that involves the application of engineering principles to the design and processing of products, both hardware and software, for meeting product reliability requirements or goals. It is a very broad design-support discipline. It has important interfaces with many other engineering disciplines. Reliability as a figure of merit (i.e. the metric) is the probability that an item will perform its intended function(s) for a specified mission profile. In general, the reliability metric can be calculated through the analyses using reliability demonstration and reliability prediction methodologies. Reliability analysis is very critical for understanding component failure mechanisms and in identifying reliability critical design and process drivers. The following sections discuss the PRA process and reliability engineering in detail and provide an application where reliability analysis and PRA were jointly used in a complementary manner to support a Space Shuttle flight risk assessment.

Safie, Fayssal M.↗

NASA NTRS2003

Bit-Serial Adder Based on Quantum Dots

A proposed integrated circuit based on quantum-dot cellular automata (QCA) would function as a bit-serial adder. This circuit would serve as a prototype building block for demonstrating the feasibility of quantum-dots computing and for the further development of increasingly complex and increasingly capable quantum-dots computing circuits. QCA-based bit-serial adders would be especially useful in that they would enable the development of highly parallel and systolic processors for implementing fast Fourier, cosine, Hartley, and wavelet transforms. The proposed circuit would complement the QCA-based circuits described in "Implementing Permutation Matrices by Use of Quantum Dots" (NPO-20801), NASA Tech Briefs, Vol. 25, No. 10 (October 2001), page 42 and "Compact Interconnection Networks Based on Quantum Dots" (NPO-20855), which appears elsewhere in this issue. Those articles described the limitations of very-large-scale-integrated (VLSI) circuitry and the major potential advantage afforded by QCA. To recapitulate: In a VLSI circuit, signal paths that are required not to interact with each other must not cross in the same plane. In contrast, for reasons too complex to describe in the limited space available for this article, suitably designed and operated QCA-based signal paths that are required not to interact with each other can nevertheless be allowed to cross each other in the same plane without adverse effect. In principle, this characteristic could be exploited to design compact, coplanar, simple (relative to VLSI) QCA-based networks to implement complex, advanced interconnection schemes. To enable a meaningful description of the proposed bit-serial adder, it is necessary to further recapitulate the description of a quantum-dot cellular automation from the first-mentioned prior article: A quantum-dot cellular automaton contains four quantum dots positioned at the corners of a square cell. The cell contains two extra mobile electrons that can tunnel (in the quantum-mechanical sense) between neighboring dots within the cell. The Coulomb repulsion between the two electrons tends to make them occupy antipodal dots in the cell. For an isolated cell, there are two energetically equivalent arrangements (denoted polarization states) of the extra electrons. The cell polarization is used to encode binary information. Because the polarization of a nonisolated cell depends on Coulomb-repulsion interactions with neighboring cells, universal logic gates and binary wires could be constructed, in principle, by arraying QCA of suitable design in suitable patterns. Again, for reasons too complex to describe here, in order to ensure accuracy and timeliness of the output of a QCA array, it is necessary to resort to an adiabatic switching scheme in which the QCA array is divided into subarrays, each controlled by a different phase of a multiphase clock signal. In this scheme, each subarray is given time to perform its computation, then its state is frozen by raising its inter-dot potential barriers and its output is fed as the input to the successor subarray. The successor subarray is kept in an unpolarized state so it does not influence the calculation of preceding subarray. Such a clocking scheme is consistent with pipeline computation in the sense that each different subarray can perform a different part of an overall computation. In other words, QCA arrays are inherently suitable for pipeline and, moreover, systolic computations. This sequential or pipeline aspect of QCA would be utilized in the proposed bit-serial adders.

Fijany, Amir↗

NASA NTRS2010

GOATS - Orbitology Component

The GOATS Orbitology Component software was developed to specifically address the concerns presented by orbit analysis tools that are often written as stand-alone applications. These applications do not easily interface with standard JPL first-principles analysis tools, and have a steep learning curve due to their complicated nature. This toolset is written as a series of MATLAB functions, allowing seamless integration into existing JPL optical systems engineering modeling and analysis modules. The functions are completely open, and allow for advanced users to delve into and modify the underlying physics being modeled. Additionally, this software module fills an analysis gap, allowing for quick, high-level mission analysis trades without the need for detailed and complicated orbit analysis using commercial stand-alone tools. This software consists of a series of MATLAB functions to provide for geometric orbit-related analysis. This includes propagation of orbits to varying levels of generalization. In the simplest case, geosynchronous orbits can be modeled by specifying a subset of three orbit elements. The next case is a circular orbit, which can be specified by a subset of four orbit elements. The most general case is an arbitrary elliptical orbit specified by all six orbit elements. These orbits are all solved geometrically, under the basic problem of an object in circular (or elliptical) orbit around a rotating spheroid. The orbit functions output time series ground tracks, which serve as the basis for more detailed orbit analysis. This software module also includes functions to track the positions of the Sun, Moon, and arbitrary celestial bodies specified by right ascension and declination. Also included are functions to calculate line-of-sight geometries to ground-based targets, angular rotations and decompositions, and other line-of-site calculations. The toolset allows for the rapid execution of orbit trade studies at the level of detail required for the early stage of mission concept development.

Haber, Benjamin M.↗

NASA NTRS1999

Studies of Radiation and Microphysics in Cirrus and Marine Stratocumulus Clouds

Two tasks were completed during this period. In the first, we examined the polarization of millimeter-wavelength radar beams scattered by ice crystals. Because of their non-spherical shape and size, ice crystals depolarize the incident polarized radar beam. In principle, this depolarization can be used to identify ice from liquid water, as well as provide some information on size. However, the amount of de-polarization is small, producing only a weak signal at the receiver. Our task was to determine the magnitude of such a signal and decide if our radar would be capable of measuring it under typical cirrus conditions. The theoretical study was carried out by Henrietta Lemke, a visiting graduate student from Germany. She had prior experience using a discrete dipole code to compute scattering depolarization. Dr. Kultegin Aydin of the Penn State Electrical Engineering Department, who is also expert in this area, consulted with us on this project at no cost to the project. Our conclusion was that the depolarization signal is too weak to be usefully measured by our system. Therefore we proceeded no further in this study. The second task involved the study of the effect of stratus microphysics on surface cloud forcing. Manajit Sengupta, a graduate student, and the project PI jointly carried out this task. The study used data culled from over a year of continuous radar and radiometer observations at the Atmospheric Radiation Measurement (ARM) site in Oklahoma. The study compared solar radiation calculations made using constant microphysics with calculations made using a retrieved mean particle size. The results showed that on average the constant microphysics produced the correct forcing when compared with the observed forcing. We conclude, therefore, that there is little impetus on radiation grounds alone to include explicit microphysics in climate models. The question of pollutant particle emission impacts on microphysics remains to be resolved. A manuscript is in preparation and will be submitted this year.

Source record↗

NASA NTRS1995

Development of Fast Algorithms Using Recursion, Nesting and Iterations for Computational Electromagnetics

In the first phase of our work, we have concentrated on laying the foundation to develop fast algorithms, including the use of recursive structure like the recursive aggregate interaction matrix algorithm (RAIMA), the nested equivalence principle algorithm (NEPAL), the ray-propagation fast multipole algorithm (RPFMA), and the multi-level fast multipole algorithm (MLFMA). We have also investigated the use of curvilinear patches to build a basic method of moments code where these acceleration techniques can be used later. In the second phase, which is mainly reported on here, we have concentrated on implementing three-dimensional NEPAL on a massively parallel machine, the Connection Machine CM-5, and have been able to obtain some 3D scattering results. In order to understand the parallelization of codes on the Connection Machine, we have also studied the parallelization of 3D finite-difference time-domain (FDTD) code with PML material absorbing boundary condition (ABC). We found that simple algorithms like the FDTD with material ABC can be parallelized very well allowing us to solve within a minute a problem of over a million nodes. In addition, we have studied the use of the fast multipole method and the ray-propagation fast multipole algorithm to expedite matrix-vector multiplication in a conjugate-gradient solution to integral equations of scattering. We find that these methods are faster than LU decomposition for one incident angle, but are slower than LU decomposition when many incident angles are needed as in the monostatic RCS calculations.

Chew, W. C.↗

NASA NTRS2022

Chapter 10 - Remote Sensing Measurements of Aerosol Properties

Satellite instruments have proven especially capable at monitoring the quantity of airborne particles in columns of atmosphere, globally. This chapter describes the principles of satellite measurements and retrieval algorithms, and surveys current instruments and their capabilities. We outline the issues associated with retrieval algorithms, such as surface characterization and aerosol proximity to clouds, and the challenges with interpretation of the results. The relationship between measured aerosol properties and climate-relevant aerosol properties simulated in models is outlined, as well as how measurements are used to evaluate models. Most space-based aerosol instruments are passive sensors that measure reflected sunlight at multiple wavelengths, some at multiple viewing angles. A few are active sensors that send out their own laser light and measure the returned signal. Except when clouds are present, the excess amount of light scattered back to space, beyond that expected from the surface and atmospheric gas, is attributed to aerosol. Satellite measurements are used in many ways in aerosol research. They often provide the only method for monitoring hazardous phenomena such as major wildfire and volcanic eruption plumes, especially in remote areas. Stable, long-term, near-global-scale satellite data records make it possible to identify regional and global aerosol trends. Aerosol radiative effects on climate can be quantified on a near-global scale and used to estimate the strength of aerosol–radiation and aerosol–cloud interactions as well as to evaluate climate model simulations of these interactions. Aerosol-type mapping from satellite imagery is helpful for source attribution, model validation, and to constrain particle light-absorption properties that are essential for radiative forcing calculations. The range of aerosol properties retrieved from satellite observations has grown considerably since the first global estimates of aerosol optical depth (τ a) over ocean were made in the late 1970s. Methods for retrieving particle size and light-absorption properties were explored in the 1990s using multispectral, multi-angle observations, and polarization in visible and near-infrared wavelengths. Sensitivity to particle light absorption, primarily from black or brown carbon content, improved with the inclusion of UV channels, and sensitivity to very thin aerosol layers in the upper troposphere and lower stratosphere was advanced with the use of limb-sounding instruments and active sensors. There are limitations to every measurement technique, including satellite aerosol remote sensing. For wide-swath, passive instruments, aerosol retrievals near clouds can present substantial challenges as far as 15 km away due to cloud-scattered light contaminating the signal. In nearly all cases, retrievals over bright snow and ice surfaces are precluded because surface reflectance uncertainties can overwhelm the aerosol signal. Similarly, meteorological cloud is identified and masked out where possible. Data from passive sensors also lack vertical resolution except those that view toward the limb or where multi-angle imagery is acquired over plumes from wildfires, erupting volcanoes, and wind-blown dust. Yet, passive sensors provide vastly more coverage than the active instruments that mitigate these issues. Particle microphysical information is qualitative from all remote sensing techniques, relying on proxies to infer particle composition, hygroscopicity, and the amount of light-absorbing material. Further, particles smaller than about 200 nm diameter cannot be distinguished from atmospheric gas molecules with remote sensing, which hinders studies of cloud condensation nuclei and their effects on clouds. Most satellite instruments dedicated to aerosol observations are in low-Earth, near-polar, sun-synchronous orbits, which means they cross the equator at the same local time each day. Most are set on cycles that repeat approximately every 16 days, which makes it difficult to monitor aerosol evolution locally. Geostationary satellites make it possible to observe changes occurring from minutes to hours over regions up to 8000 km in size, but lack coverage of high latitudes, and often provide more limited constraints on aerosol properties. Ground-truth data are vital for satellite aerosol-retrieval validation. The AErosol RObotic NETwork (AERONET) of sun photometers was created in 1993 and has become an established global network of over 350 instruments for validating satellite measurements. The network, as well as global networks of ground-based lidars, solar flux radiometers and other sun photometers, are widely used for evaluating global satellite retrievals and model simulations. NASA's Earth Observing System (EOS) program beginning in 1999 led to improvements in reliability, spatial resolution, and spectral resolution (and hence, to improved particle size discrimination and light absorption properties). Satellite payloads include advanced broad-swath and multi-angle imagers, along with the first space-based active sensor focused largely on long-term aerosol monitoring. Since about 2002, Europe's SENTINEL and operational meteorological satellite fleets are also providing sustained aerosol observations, with planned continuation until at least 2030. Satellite remote sensing instruments offer valuable data for evaluating aerosol representations in global climate models. They have been used to assess aerosol optical and physical properties, trends and distributions, and are applied increasingly as direct model constraints in data assimilation to create global aerosol reanalysis products. Aerosol optical depth is the most common quantity adopted for routine model evaluation, including multiwavelength data to loosely constrain particle-size distributions. These evaluations of multiple models have revealed general biases in their regional aerosol amounts and seasonal patterns of transport and removal. Although satellite measurements have near-global coverage, substantial errors can be introduced into the model observation comparison unless attention is paid to spatial and temporal collocation, cloud screening, subgrid-scale variability, and measurement uncertainties that vary with retrieval conditions.

aerosol properties↗

NASA NTRSDate not supplied

Synergy of Satellite Radiation, Precipitation, and Other Meteorological Variable Observations for Global Mean Sea Surface Turbulent Heat Flux Estimation

Sea surface turbulent heat flux is one of the key components in global freshwater and energy balances and plays an important role in atmospheric dynamics, thermodynamics and general circulation. This turbulent heat flux is dominantly decided by sea surface latent heat release with some contribution from sensible heat exchange. The flux and its anomaly could significantly affect ocean heat storage and ocean circulation. They are part of and have great impacts on climate variability. Though the turbulent heat flux is extremely important for the climate and weather systems, only very limited ship and buoy observations of the flux are available over open oceans. There are no global, operational, direct measurements of this crucial meteorological variable. Global estimates are, basically, indirectly calculated from bulk turbulent flux parameterization with a combination of satellite column water vapor, sea surface water temperature, and wind speed observations. Critical parameters such as sea surface air temperature and humidity are estimated from empirical relations of theses variables with the water vapor and sea surface water temperature, respectively. Lacking accurate knowledge on surface air temperature and humidity, uncertainties in the parameterization and potential changes in the non-linear turbulent processes with long-term climate variations could cause large errors in estimated long-term turbulent fluxes from the indirect method as shown in current global sea surface turbulent flux datasets. This study uses synergized data of satellite global radiation, precipitation, and other meteorological variable observations to estimate sea surface turbulent heat fluxes. The radiation observations are made by the satellite Clouds and the Earth’s Radiant Energy System (CERES) sensors, while the global precipitation data is from NASA’s satellite Global Precipitation Climatology Project (GPCP). Other data includes satellite sea surface water temperature and wind speed observations. These datasets are obtained from a wide range of space sensors from passive to active instruments and from visible and near infrared to thermal infrared and microwave spectral sounders. One significant feature of these datasets is that they have multi-decades long climate records. Top-of-atmosphere (TOA) radiation and its anomaly represents the net heat energy input to the climate system, and the oceanic precipitation and ocean-land moisture transport can be used to quantify sea surface latent heat release. Based on the synergized datasets and the principle of global water and energy balances, global mean turbulent heat fluxes are estimated. The turbulent heat anomalies for the first two decades of the 21st century are, then, obtained mainly from global CERES radiation and GPCP precipitation anomalies, along with CERES derived ocean-land heat transports. Analysis indicates that the uncertainties in the estimated turbulent flux anomalies may be reduced considerably. The results suggest a strong needs in synchronized and synergized observations of atmospheric radiation, precipitation, and oceanic meteorological variables for long-term climate studies.

Bing Lin↗

NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗

NASA NTRS2019

Experimental Demonstration and System Analysis for Plasmonic Force Propulsion

This project developed a higher-fidelity model of a recently envisioned small spacecraft propulsion system for precision pointing and proximity control. Plasmonic force propulsion harnesses solar light focused onto plasmon reactive subwavelength nanostructures to accelerate and expel nanoparticle propellant via strong optical forces. The goal of the project was to show that plasmonic space propulsion can provide the level of proximity and attitude control envisioned for future NASA nano/picosatellite missions, a level that is better than state-of-the-art approaches. We achieved this goal by showing that plasmonic force thrusters are feasible for a range of advanced mission concepts requiring swarm formations in a deep space environment. We performed three case studies that evaluated the performance of the plasmonic force propulsion thruster in a deep space, microsatellite swarm formation. These case studies assumed the propulsion system could generate thrust at the level predicted from our Phase 1 study (1.6 μN). Through these cases we were able to analyze the concept within a mission specific context through detailed orbital dynamics calculations. Results indicate that, with the Phase 1 estimated thrust level, the approach is promising for providing attitude control to swarm formation spacecraft. Further, we achieved goals related to technology development. Specifically, we experimentally demonstrated nanoparticle acceleration due to plasmonic forces with asymmetric nanostructures excited by focused laser light. Additionally, we investigated the thrust sensitivity and nanoparticle propellant injection dependencies upon thermal effects. As a result of our study, plasmonic force propulsion is at an early TRL 3. Active research and design has been conducted analytically and in the laboratory. Furthermore, practical applications such as the three case studies have been identified for the scientific basic principles that were observed. Future efforts related to fundamental understanding of these techniques should focus on 1) developing a standalone array of asymmetric nanostructures that can effectively interact with a stream or reservoir of particles or 2) experimentally evaluate a dielectrophoretic injector for nanoparticle propellant. The main limitation discovered about plasmonic propulsion regards performance estimates significantly below the Phase 1 estimations. Specifically, original assumptions in the Phase 1 project (notably, a linear array of asymmetric nanostructures) is not a viable approach to achieving significant acceleration, high exhaust velocity, of nanoparticles. More specifically, we assumed in Phase 1 that nanoparticles would be accelerated in series by a long linear array of asymmetric nanostructures. That is, the acceleration of the nanoparticle would build and increase with the kick received by each subsequent nanostructure. This is fundamentally flawed. The potential profile of a single nanostructure is such that it prohibits this phenomenon. The potential energy associated with the plasmon-generated dielectrophoretic force is a potential well, which is good for trapping nanoparticles, but cannot provide significant acceleration of particles to expel them out and away from the nanostructure. Further, a nanoparticle expelled from the first nanostructure would need to overcome the potential barrier for entry into the next nanostructure accelerating stage. Fundamentally, this effect means that a linear array of nanostructures is not a viable accelerating structure. Correspondingly then, acceleration can, or should, only be provided by one nanostructure, and the net acceleration and thrust force of a single nanostructure is small (~cm/s exhaust velocities, sub-nN level thrust vs. the 100’s m/s, μN originally envisioned). While our experiments demonstrated acceleration and manipulation of a nanoparticle using laser light in aqueous environment, the achievable energy and momentum addition to the nanoparticle from a single nanostructure stage is too low for useful propulsion. In terms of thrust prediction, the estimated thrust of 1.6 μN in Phase 1 is reduced to a few nN of thrust with this new insight and understanding of the concept. This thrust level is too small to achieve attitude control of swarms as originally envisioned.

Plasmonic Force↗

NASA NTRS1994

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗

NASA NTRS1962

PYROLYTIC GRAPHITE ROCKET THRUST CHAMBER INVESTIGATION PROJECT. MONTHLY LETTER

This is the Monthly Letter Progress Report for the eighth month of a twelve month technical study directed toward the following objectives: A. Determining analytically and experimentally the feasibility and attractiveness of using free-standing pyrolytic graphite for a radiation cooled liquid rocket thrust chamber and exit nozzle for space application. B. Establishing through analysis and experimental demonstration, a set of design principles to facilitate the design of a reliable pyrolytic graphite thrust chamber assembly for a specific range of operating conditions. C. Demonstrate by means of test firings of complete thrust chamber assemblies, the validity of the design data and principles developed during the program. Preparations have been completed for test firing the two 100 lb. thrust Boron Pyralloy chambers at Marquardt during the first week in March 1962. These two chambers will first be static pressure checked with water to a pressure of 135 psig. The chambers will then be test fired with N204/N2H4-MMH propellants . Test samples of eleven different pyrolytic carbide and nitride materials have been ordered for high temperature oxidation tests and should be delivered by the end of March for evaluation in the plasma torch facility. The specific materials ordered in free- standing condition in thicknesses from 0.030 inch to 0.065 inch were as follows: (a) Pyrolytic graphite (on hand) (b) Pyrolytic graphite-boron alloy (on hand) (c) Pyrolytic graphite-tungsten alloy (d)Boron nitride (e) Titanium nitride (f) Silicon carbide (g) Titanium carbide (h) Zirconium carbide (i) Hafnium carbide (j) Niobium carbide (k) Tantalum carbide (l) High density (1 1-90) graphite (m) High density graphite (heat treated) These materials were ordered from Raytheon, HTM, and American Meta 1 Products . A summary of the available data on the oxidation rates of pyrolytic materials under slow air and rocket motor conditions is shown in Figure 1. High Temperature Materials, Inc. has been contacted relative to providing Pyrographite thrust chambers with extra-thick walls . They believe they can deliver chambers with substantially thicker walls than were possible a short time ago. Such a configuration, with acceptable residual stresses, would provide a longer firing life since it is ultimately limited by the oxidation rate of the chamber wall material. Stress and structural studies have continued. The analytical stress studies point out some of the problems associated with the experimental determination of the elastic constants of this highly anisotropic material. Further analysis may shed light on improved techniques for making these experimental determinations. The large potential weight saving available with the use of pyrolytic graphite as a thrust chamber wall material in a radiation cooled rocket motor has been calculated and is shown graphically in Figure 2. Figure A presents the program progress schedule and the actual expenditures as of 7 February 1962.

ROCKET ENGINE↗

NASA NTRS2000

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗

NASA NTRS2013

Put X-Ray Polarimetry on the MAP!

With Prof. R. Novick and others at the Columbia Astrophysics Laboratory I help found the field of X-ray polarimetry in the early 1970s. Currently I have more experience with the design, construction, calibration, and space flight of such instruments than anyone on the planet. The early probing beginnings saw only one definitive measurement (that of the integrated low-energy X-ray emission from the Crab Nebula sans pulsar) and a number of upper limits. The limited success did nevertheless inspire a number of detailed theoretical calculations, concentrating at first on neutron stars and black holes showing how precise measurements (e.g. degree of polarization and position angle as a function of pulse phase) would provide definitive limitations on otherwise equally plausible theoretical models. Over time the theoretical foundation has grown (e.g. the proceedings of the X-Ray Polarimetry Workshop held at SLAC in 2004). I will outline these foundations. It is important to understand the history of X-ray polarimetry beyond the early excitement. A polarimeter was at the focus of the original Einstein mission but was dropped during the restructuring. A polarimeter was successfully proposed (R. Novick PI, I was a Co-I) and built for the original Spectrum-X mission. During the years before the cancellation of Spectrum-X, the potential flight of this device stood in the way of other space flights for polarimeters --- "let us wait and see what it finds". This was unfortunate as there were a number of reasons why that polarimeter should not have been flown on the mission. Perhaps the most significant (but not only) reason was that a shared focal plane provided very little observing time. This is an extremely important point in considering the Roadmap. It is doubtful that there many 100%-polarized sources and so the "signal" is more typically a small fraction of the source flux. Thus, the source itself provides a substantial background, making continuum polarimetry even more difficult than narrow-line spectroscopy. X-ray polarimetry thus requires a dedicated mission that can, without programmatic pressures from other instruments/users, devote the integration time to perform meaningful measurements. The recently cancelled GEMS might have been such a mission. At least it was dedicated to polarimetry. Performing meaningful measurements is not going to be easy. In part because of the long hiatus and lack of experience, there appears to be too much pressure to "sell" polarimetry missions by the number of sources for which one might answer the simple question is, or is not, the integrated and time averaged emission from the source polarized at some confidence level? This was a fine question for the 1970 s but, I maintain, it is not today. One simply doesn't want to measure the time averaged polarization of the Crab s pulsar, but one wants to know the polarization as a function of energy and pulse phase to compare, e.g. to optical and radio measurements which divide even the primary pulse into dozens of phase bins. Such observations can distinguish amongst competing theories for the pulsed emission. The Roadmap should define what meaningful experiments are. I will pose some suggestions. Note that , because X-rays are usually believed originate in either non-thermal or highly aspherical situations we expect X-ray polarimetry to be much more important and rich in astrophysical information as opposed to the visible, where starlight often dominates the emission. One has often dreamt about an instrument that does polarimetry whilst it does other things, and I will discuss this. Even in this case, one needs to realize that the observing time will be driven by the polarimetry, otherwise no useful polarization measurements will be made. Finally, I will discuss some misconceptions that have appeared in the literature and at conferences which indicate to me that certain fundamental principles of polarimeter design and performance are not clearly understoodt is important for the Roadmap that we understand these principles and not advocate instruments that will not (or even worse cannot) accomplish the promised scientific return. I will discuss these principles.

Weisskopf, Martin C.↗