Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient estimation”

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 577 records · Page 32

Thermocapillary Flow and Aggregation of Bubbles on a Solid Wall

During the electrolytic evolution of oxygen bubbles forming on a vertically oriented transparent tin oxide electrode, bubbles were found to be mutually attractive. The mechanism of the aggregation had never been explained satisfactorily until Guelcher et al. attributed it to thermocapillary flow. The gradient of surface tension of the liquid at the bubble's surface, which was established because of reaction heat and ohmic heat loss at the electrode wall, drives flow of the liquid adjacent to each bubble; the bubble "pumps" fluid along its surface away from the wall. Fluid flows toward the bubble to conserve mass and entrains nearby bubbles in the flow pattern. The same logic would apply when two bubbles of equal size are adjacent to each other on a warm wall. Each bubble drives thermocapillary flow and hence entrains the other in its flow pattern, which drives the aggregation. Our objective here is to perform experiments where the temperature gradient at the wall is well known and controlled. The theory can be quantitatively tested by studying aggregation of bubble pairs of equal size, and by varying system parameters such as temperature gradient, bubble size and fluid viscosity. The results are then compared with the theory in a quantitatively rigorous manner. We demonstrate that the theory without adjustable parameters is capable of quantitatively modeling the rate of aggregation of two bubbles. The equations governing the thermocapillary flow around a single stationary bubble on a heated or cooled wall in a semi-infinite domain were solved. Both Reynolds number and Marangoni number were much less than unity. The critical result is that liquid in the vicinity of a warm wall flows toward a stationary collector bubble. Consequently the thermocapillary flow around the stationary bubble entrains another bubble toward itself. The bubbles undergo hindered translation parallel to the wall with velocity U while the fluid flow field is described with u. Two velocities were equated by using a wall hindrance parameter q: U = qu [1] which shows the velocity of bubble is proportional to the entraining velocity. The hindrance parameter q can experimentally be measured independently. q can also be calculated by solving the equations of motion for a bubble translating parallel to a solid wall. The experimental cell is cylindrical with an ID of 10 cm and consists of a 1 cm deep main cell filled with silicone oil and flanked by two thermal reservoirs. The upper thermal reservoir was heated and the lower thermal reservoir was cooled so that the bubbles aggregate. Two types of silicone oil (eta = 0.02 and 0.50 Pa s) were used. Two equal sized air bubbles were injected into the cell with a syringe. The center-to-center distance of bubbles was observed through a microscope. Bubble radius ranged from 0.40 mm to 0.65 mm and the temperature gradients along with the cell ranged from 1400 to 5000 K/m. The bubbles aggregated when heat flows from the wall to the fluid. The velocities of bubbles were in the range of 1 - 10 microns/s. The separation r decreased more quickly when the temperature gradient was higher, bubble size was larger, and the oil viscosity was lower. r decreased more rapidly as the bubbles approached each other. Dimensionless time was arbitrarily set to be zero when the dimensionless center-to-center distance between the bubbles was 4. All the bubble trajectories fall onto one line, especially in the range of dimensionless distance from 4 to 3. This means the relative movement of the bubble pair is proportional to the temperature gradient and bubble size and it is inversely proportional to the viscosity of the oil. This result strongly suggests that the thermocapillary flow-based aggregation mechanism is correct. A value of q can be estimated by fitting the scaled data to Eq. [1]. A best fit value of q was obtained as q = 0.26 with a standard deviation of 0.03. Independent experimental results for q for a 0.5 mm radius bubble, give values of q in the range 0.11 to 0.23. The value of q obtained from solving the equations of motion reveals q has values in the range 0.23 - 0.30. Since the full scale of possible values of q is zero to one, the maximum deviation of independently determined values of q from the best fit value was 15% of this full scale. Thus reasonable quantitative agreement between theory and experiment has been obtained.

Kasumi, Hiroki↗

NASA Tech Briefs, August 2005

Topics include: Hidden Identification on Parts: Magnetic Machine-Readable Matrix Symbols; System for Processing Coded OFDM Under Doppler and Fading; Multipurpose Hyperspectral Imaging System; Magnetic-Flux-Compensated Voltage Divider; High-Performance Satellite/Terrestrial-Network Gateway; Internet-Based System for Voice Communication With the ISS; Stripline/Microstrip Transition in Multilayer Circuit Board; Dual-Band Feed for a Microwave Reflector Antenna; Quadratic Programming for Allocating Control Effort; Range Process Simulation Tool; Simulator of Space Communication Networks; Computing Q-D Relationships for Storage of Rocket Fuels; Contour Error Map Algorithm; Portfolio Analysis Tool; Glass Frit Filters for Collecting Metal Oxide Nanoparticles; Anhydrous Proton-Conducting Membranes for Fuel Cells; Portable Electron-Beam Free-Form Fabrication System; Miniature Laboratory for Detecting Sparse Biomolecules; Multicompartment Liquid-Cooling/Warming Protective Garments; Laser Metrology for an Optical-Path-Length Modulator; PCM Passive Cooling System Containing Active Subsystems; Automated Electrostatics Environmental Chamber; Estimating Aeroheating of a 3D Body Using a 2D Flow Solver; Artificial Immune System for Recognizing Patterns; Computing the Thermodynamic State of a Cryogenic Fluid; Safety and Mission Assurance Performance Metric; Magnetic Control of Concentration Gradient in Microgravity; Avionics for a Small Robotic Inspection Spacecraft; and Simulation of Dynamics of a Flexible Miniature Airplane.

Source record↗

Efficient On-Orbit Singularity-Free Geopotential Estimation

For vehicles in orbital cruise, either in orbit around a body or in a transfer orbit, gravity is the primary external force; therefore, the complexity of the geopotential model heavily influences the accuracy of the navigation state. As this type of spacecraft requires a highly accurate geopotential model and performs the computation numerous times a second, this calculation needs to be efficient. This calculation needs to include the first and second derivatives of the geopotential because navigation systems can require both the gravitational acceleration and gradient for inertial integration and state filtering. The most efficient method for calculating the geopotential, the forward column method, is discussed in detail, as well as a method to avoid the singularities that exist when using this method. This is shown to decrease the computation time of the geopotential compared to other popular methods. In addition, methods for first and second order propagation are discussed, which decrease the rate at which the full geopotential model needs to be calculated while maintaining the accuracy required of navigation systems. These estimation methods are shown to decrease the necessary computation of the gravity model by multiple orders of magnitude. This method decreases the computation time of the Exploration Upper Stage navigation system geopotential model by potentially an order of magnitude compared with the previous model without affecting the navigation error. A properly implemented geopotential model can have the accuracy of an 8 by 8 model while approaching the computational requirement as using only the J234 coefficients.

Amert, Joel↗

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗

Reproductive and leaf litterfall fluxes in forest ecosystem sites globally (1950-2022)

Forest allocation of net primary productivity (NPP) to reproduction is poorly quantified globally, despite its critical role in forest regeneration and a well-supported trade-off with allocation to growth. Although field measurements of total NPP are rare, our work finds that a proxy for reproductive carbon allocation constructed from leaf (L) and reproductive (R) litterfall fluxes, R/(R+L), is strongly correlated with R/NPP, facilitating analysis across a wide range of sites where biometric estimates of NPP are not available (R² = 0.85; Hanbury-Brown et al., 2022, Ward et al., in prep). To investigate relationships between ecosystem-scale reproductive allocation (RA) and climate, soil fertility, and stand age gradients, we conducted a literature search and synthesized 824 observations of annual average leaf and reproductive litterfall fluxes across forest sites globally. The zip file includes 1) a folder Data/ containing the litterfall data ("GlobalForestRA_data.csv") and metadata ("GlobalForestRA_metadata.doc") files. The data file includes geographic coordinates, long-term mean annual temperature and precipitation (1970-2000, extracted from WorldClim2.1), leaf and reproductive litterfall fluxes, sampling interval and protocols, forest characteristics (dominant leaf morphology, information pertaining to forest age and successional stage, and disturbance history) and soil properties (% sand, %silt, %clay, total phosphorus (P), nitrogen (N), cation exchange capacity (CEC) and pH) extracted from SoilGrids250 and from on-site measurements, where available. The metadata file contains information about each variable reported in the data file, including data sources, processing methods, and all references. The Data folder contains two additional files used to create Figure 1; these are described in greater detail in the README.2) R scripts GloalForestRA_analysis.r and GlobalForestRA_SI.r and a folder /Functions used to produce results, figures, and tables in the manuscript Ward et al. (in press)3) a README file describing how the data and R scripts can be used to reproduce statistical results, figures, and tables found in the manuscript. Ward et al. (in press)This repository can also be found at: https://github.com/r-ward/Global_Analysis_ForestRA.Ward, R.E., Zhang-Zheng, H. Aernethy, K., Adu-Bredu, S., Arroyo, L., Bailey, A. et al. (in press). Forest age rivals climate to explain reproductive allocation patterns in forest ecosystems globally. Ecology Letters. Hanbury-Brown, A.R., Ward, R.E. & Kueppers, L.M. (2022). Forest regeneration within Earth system models: current process representations and ways forward. New Phytol., 235, 20–40.Ward et al. (2025), Forest age rivals climate to explain reproductive allocation patterns in forest ecosystems globally, in prep.

54 ENVIRONMENTAL SCIENCES↗

Waterspout, Gust Fronts and Associated Cloud Systems

Nine waterspouts observed on five experimental days during the GATE period of observations are discussed. Primary data used are from 2 aircraft flying in different patterns, one above the other between 30 and 300 m. There is strong evidence associating whirl initiation with cumulus outflow. Computations prepared from estimates of convergence with the region suggest the possibility of vortex generation within 4 minutes. This analysis supports (1) the importance cumulus outflows may have in waterspout initiation and (2) the possibility that sea surface temperature gradients may be important in enabling waterspout development from modest size cumuli.

Simpson, J.↗

Hypersonic stability and transition

An evaluation is conducted of the development status of growing normal modes in the computation of flow fields over hypersonic vehicles, giving attention to the effects of surface cooling, pressure gradients, and suction on the growth rates of the normal modes. Also noted are the effects of bluntness on the stability properties of hypersonic vehicles, and the efficacy of the e exp N method of transition estimation at supersonic and hypersonic speeds. A survey is made of the precautions that must be taken to ensure the effectiveness of experimental and computational studies of hypersonic flow phenomena.

Reshotko, Eli↗

Mapping Pine Plantations in the Southeastern U.S. Using Structural, Spectral, and Temporal Remote Sensing Data

The southeastern U.S. produces the most industrial roundwood in the U.S. each year, largely from commercial pine plantations. The extent of plantation forests and management dynamics can be difficult to ascertain from periodic forest inventories, yet short-rotation tree plantations also present challenges for remote sensing. Here, we integrated spectral, temporal, and structural information from airborne and satellite platforms to distinguish pine plantations from natural forests and evaluate the contribution from planted forests to regional forest coverin the southeastern U.S. Within flight lines from NASA Goddard's Lidar, Hyperspectral, and Thermal (G-LiHT) Airborne Imager, lidar metrics of forest structure had the highest overall accuracy for pine plantations among single-source classifications (90%), but the combination of spectral and temporal metrics from Landsat generated comparable accuracy (91%). Combined structural, temporal, and spectral information from G-LiHT and Landsat had the highest accuracy for plantations (92%) and natural forests (88%). At a regional scale, classifications using Landsat spectral and temporal metrics had between 74 and 82% mean class accuracy for plantations.Regionally, plantations accounted for 28% of forest cover in the southeastern U.S., a result similar to plot-based estimates, albeit with greater spatial detail. Regional maps of plantation forests differed from existing map products, including the National Land Cover Database. Combining plantation extent in 2011 with Landsat based forest change data identified strong regional gradients in plantation dynamics since 1985, with distinct spatial patterns of rotation age (east-west) and plantation expansion (interior). Our analysis demonstrates the potential to improve the characterization of dynamic land cover classes, including economically important timber plantations, by integrating diverse remote sensing datasets. Critically, multi-source remote sensing provides an approach to leverage periodic forest inventory data for annual monitoring of managed forest landscapes.

Fagan, M. E.↗

Observational Constraint on Cloud Feedbacks Suggests Moderate Climate Sensitivity

Global climate models (GCMs) predict warming in response to increasing greenhouse gases, partly due to decreased tropical low-level cloud cover and reflectance. We use satellite observations that discriminate stratocumulus (Sc) from shallow cumulus (Cu) clouds to separately evaluate their sensitivity to warming and constrain the tropical contribution to low-cloud feedback. We find an observationally inferred low-level feedback two times smaller than a previous estimate. Cu are insensitive to warming whereas GCMs exhibit a large positive cloud feedback in Cu regions. In contrast, Sc show sensitivity to warming and the tropical inversion layer strength, controlled by the tropical Pacific SST gradient. Models fail to reproduce the historical SST gradient trends and therefore changes in inversion strength, generating an overestimate of the positive Sc cloud feedback. Continued weak east-Pacific warming would therefore produce a weaker low-cloud feedback and imply a more moderate climate sensitivity (3.47  0.33 K) than many models predict.

Cloud feedback↗

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗

Global particle buildup simulations with gas puff scan: application to WEST discharge

This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheath-limited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Operational Data Reduction Procedure for Determining Density and Vertical Structure of the Martian Upper Atmosphere from Mars Global Surveyor Accelerometer Measurements

The success of aerobraking by the Mars Global Surveyor (MGS) spacecraft was partly due to the analysis of MGS accelerometer data. Accelerometer data was used to determine the effect of the atmosphere on each orbit, to characterize the nature of the atmosphere, and to predict the atmosphere for future orbits. To interpret the accelerometer data, a data reduction procedure was developed to produce density estimations utilizing inputs from the spacecraft, the Navigation Team, and pre-mission aerothermodynamic studies. This data reduction procedure was based on the calculation of aerodynamic forces from the accelerometer data by considering acceleration due to gravity gradient, solar pressure, angular motion of the MGS, instrument bias, thruster activity, and a vibration component due to the motion of the damaged solar array. Methods were developed to calculate all of the acceleration components including a 4 degree of freedom dynamics model used to gain a greater understanding of the damaged solar array. The total error inherent to the data reduction procedure was calculated as a function of altitude and density considering contributions from ephemeris errors, errors in force coefficient, and instrument errors due to bias and digitization. Comparing the results from this procedure to the data of other MGS Teams has demonstrated that this procedure can quickly and accurately describe the density and vertical structure of the Martian upper atmosphere.

Cancro, George J.↗

A Validated All-Pressure Fluid Drop Model and Lewis Number Effects for a Binary Mixture

The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Harstad, K.↗

Data Fusion for Earth Science Remote Sensing

Beginning in 2004, NASA has supported the development of an international network of ground-based remote sensing installations for the measurement of greenhouse gas columns. This collaboration has been successful and is currently used in both carbon cycle investigations and in the efforts to validate the GOSAT space-based column observations of CO2 and CH4. With the support of a grant, this research group has established a network of ground-based column observations that provide an essential link between the satellite observations of CO2, CO, and CH4 and the extensive global in situ surface network. The Total Carbon Column Observing Network (TCCON) was established in 2004. At the time of this report seven sites, employing modern instrumentation, were operational or were expected to be shortly. TCCON is expected to expand. In addition to providing the most direct means of tying the in situ and remote sensing data sets together, TCCON provides a means of testing the retrieval algorithms of SCIAMACHY and GOSAT over the broadest variation in atmospheric state. TCCON provides a critically maintained and long timescale record for identification of temporal drift and spatial bias in the calibration of the space-based sensors. Finally, the global observations from TCCON are improving our understanding of how to use column observations to provide robust estimates of surface exchange of C02 and CH4 in advance of the launch of OCO and GOSAT. TCCON data are being used to better understand the impact of both regional fluxes and long-range transport on gradients in the C02 column. Such knowledge is essential for identifying the tools required to best use the space-based observations. The technical approach and methodology of retrieving greenhouse gas columns from near-IR solar spectra, data quality and process control are described. Additionally, the impact of and relevance to NASA of TCCON and satellite validation and carbon science are addressed.

data fusion↗

Shape Servoing of Deformable Objects using Adaptive Deformation Model Estimation

In this paper, we propose an adaptive shape servoing method to deform a soft object into a desired 3-D shape. The high dimensional representation and the unknown deformation properties of the soft object pose a challenge to actively manipulate its shape. To address this issue, we develop a method to compute the deformation Jacobian matrix in real-time. The Jacobian is estimated using a set of basis functions and its corresponding parameters to capture the dynamics of the system and relate the applied input motion to changes in the soft object's shape. An integral concurrent learning (ICL) based adaptive update law is derived using Lyapunov analysis to estimate the deformation parameters and prove its convergence. A physics-based simulation is used to validate the proposed method and controller by performing manipulation tasks with different desired configurations. The performance is compared with a standard gradient update law to demonstrate the accuracy and robustness of our approach.

Vrithik Raj Guthikonda↗

Study of Fluid Flow Control in Protein Crystallization using Strong Magnetic Fields

An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in "microgravity", researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined. Specific goals for the investigation are: 1. To develop an understanding of convection control in diamagnetic fluids with concentration gradients through experimentation and numerical modeling. Specifically solutal buoyancy driven convection due to crystal growth will be considered. 2. To develop predictive measures for successful crystallization in a magnetic field using analyses and numerical modeling for use in future protein crystal growth experiments. This will establish criteria that can be used to estimate the efficacy of magnetic field flow damping on crystallization of candidate proteins. 3. To demonstrate the understanding of convection damping by high magnetic fields to a class of proteins that is of interest and whose structure is as yet not determined. 4. To compare quantitatively, the quality of the grown crystals with and without a magnetic field. X-ray diffraction techniques will be used for the comparative studies. In a preliminary set of experiments, we studied crystal dissolution effects in a 5 Tesla magnet available at NASA Marshall Space Flight Center (MSFC). Using a Schlieren setup, a 1mm crystal of Alum (Aluminum-Potassium Sulfate) was introduced in a 75% saturated solution and the resulting dissolution plume was observed. The experiment was conducted both in the presence and absence of a magnetic field gradient. The magnet produces a gradient field of approx. 1 Tesla2/cm. Image analysis of the recorded images indicated an enhanced plume velocity that was of the order of the measurement limit. For this experiment, both the gradient and gravity fields are in the same direction resulting in an enhanced effective gravity that tends to accelerate the observed plume velocity. While the results are not conclusive, pending further tests, it clearly points out the inadequacy of the MSFC magnet for conducting protein crystallization experiments and the need for a stronger magnet. In spacebased experiments, however, where the gravitational effects are small, only a weak magnetic field will be required to control or mitigate the effects of convective contamination.

Ramachandran, Narayanan↗

Estimation of particle size based on LDV measurements in a de-accelerating flow field

The accuracy of velocity measurements made with a laser velocimeter is strongly dependent upon the response of the seeding particles to the dynamics of the flow field. The smaller the particle the better the response to flow fluctuations and gradients and therefore the more accurate velocity measurement. In direct conflict is the requirement of light scattering efficiency to obtain signals with the laser velocimeter which, in general, is better as the particle size is increased. In low speed flow fields these two requirements on particle size overlap and accurate measurements may be obtained. However in high speed flows, where the velocity gradients may be severe, very small particles are required to maintain sufficient dynamic response characteristics to follow the flow. Therefore if velocity measurements are to be made in these flows, the laser velocimeter must be designed with sufficient sensitivity to obtain signals from these small particles. An insitu determination of the size distribution of kaolin particles (Al2O3, .2 + or - SiO2 . 2H2O) in the 16-foot Transonic Tunnel and the sensitivity characteristics of the laser velocimeter system is described.

Meyers, J. F.↗

Characterizing and improving the performance of molten-salt-steam heat exchangers in concentrating solar power plants

Shell-and-tube heat exchangers (HXs) for steam generation from molten salts in concentrating solar power (CSP) plants experience thermal fatigue due to significant temperature gradients and inherent transient operation. Molten salt-steam HX design lifespans exceed actual lifespans, and, as a consequence, designers overpredict plant profitability and operators neglect appropriate prescriptions to optimize these lifetimes. Here, this study refines HX lifespan estimates with data benchmarked against thermal-fluid mechanical modeling of stress and accumulated fatigue. Reduced-order thermal models of the molten salt-steam, shell-and-tube evaporator and superheater predict transient temperature profiles along the two HXs salt-steam flow paths. The modeled evaporator and superheater temperature profiles enable assessment of cyclic stresses within the HX tubesheets, where molten-salt HX failures are most common. Evaporator and superheater performance data from a current 110 MW elec commercial CSP plant provide a basis for validating the reduced-order HX models. HX life predictions derived from stochastic failure distributions serve as inputs for simulating and optimizing existing plant operations. The impact of the updated lifespans on overall plant revenue depends on operating scenarios. This study suggests that typical ramping rates for a CSP plant with a high-temperature Rankine cycle result in an evaporator and superheater life of approximately 10 and 25 years, respectively, compared to the design target of 30 years. Reduced HX lifespans decrease operational plant revenue on average by 4.6-5.1%. Furthermore, there may be as many as four HX replacements over the 30-year lifetime of the plant; and, purchase agreement loss due to failure to meet contractual production requirements can have ramifications that include the risk of bankruptcy.

14 SOLAR ENERGY↗