Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear potential flow”

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 397 records · Page 22

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT 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. 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 (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. 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 two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗

Aerodynamic control of NASP-type vehicles through vortex manipulation, volume 4

Forebody Vortex Control (FVC) is an emerging technology that has received widespread and concentrated attention by many researchers for application on fighter aircraft to enhance aerodynamic controllability at high angles of attack. This research explores potential application of FVC to a NASP-type configuration. The configuration investigated is characterized by a slender, circular cross-section forebody and a 78 deg swept delta wing. A man-in-the-loop, six-degress-of-freedom, high-fidelity simulation was developed that demonstrates the implementation and advantages of pneumatic forebody vortex control. Static wind tunnel tests were used as the basis for the aerodynamic characteristics modeled in the simulation. Dynamic free-to-roll wind tunnel tests were analyzed and the wing rock motion investigated. A non-linear model of the dynamic effects of the bare airframe and the forebody vortex control system were developed that closely represented the observed behavior. Multiple state-of-the-art digital flight control systems were developed that included different utilizations of pneumatic vortex control. These were evaluated through manned simulation. Design parameters for a pneumatic forebody vortex control system were based on data collected regarding the use of blowing and the mass flow required during realistic flight maneuvers.

Smith, Brooke C.↗

Exploring Airfoil Table Generation using XFOIL and OVERFLOW

The rotorcraft design process is a continuously evolving field of research that incorporates a number of software programs. An accurate airfoil table is critical in the design and testing process for rotorcraft. With multiple flow solvers available and flow conditions of multirotor UAM vehicles potentially covering a wide range of Reynolds and Mach numbers, a documented approach for developing airfoil tables is needed. Using benchmark data from legacy airfoil tables and wind tunnel tests for comparison, simulations for a comprehensive test matrix could guide rotorcraft design engineers in generating their own airfoil tables using the XFOIL and OVERFLOW solvers. The motivation for this study is to investigate flow solver features to develop a best practices document for airfoil table generation. The study uses the OVERFLOW and XFOIL flow solvers, coupled with the airfoil table generator AFTGen, to analyze three airfoils for a specific Reynolds numbers flow regime and provide details on how well each flow solver performs within a specific angle of attack range, Mach number range, Reynolds number range, and in different flow conditions, such as turbulent and transitional flow. OVERFLOW analyses in AFTGen for fully turbulent and transition flow are compared with XFOIL results and experimental test data for the section lift, section drag, and pressure coefficients. XFOIL ultimately yields results that are accurate within the linear angle of attack range and below a Mach number of 0.4 but tends to overpredict lift and underpredict drag unless the flow is in the compressible regime. XFOIL cannot accurately model stall and post-stall conditions due to the nature of the solver. This is evident in nearly every case run with XFOIL, where the linear range is usually predicted acceptably and the lift coefficient is overpredicted as the stall angle of attack is approached (with the exception being the generally poor correlation with most of the SSC-A09 cases). OVERFLOW is limited at low Mach numbers, and appears to perform best at Mach numbers of 0.4 and above. The exploration of airfoil table generation using XFOIL and OVERFLOW yielded moderately successful results for the NACA 0012 airfoil, reasonably good results for the RC(4)-10 airfoil, and less accurate results for the SSC-A09 airfoil.

Airfoil Table Generation↗

Stability investigations of airfoil flow by global analysis

As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.

Morzynski, Marek↗

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗

Evaluating GCM land surface hydrology parameterizations by computing river discharges using a runoff routing model: Application to the Mississippi basin

To relate general circulation model (GCM) hydrologic output to readily available river hydrographic data, a runoff routing scheme that routes gridded runoffs through regional- or continental-scale river drainage basins is developed. By following the basin overland flow paths, the routing model generates river discharge hydrographs that can be compared to observed river discharges, thus allowing an analysis of the GCM representation of monthly, seasonal, and annual water balances over large regions. The runoff routing model consists of two linear reservoirs, a surface reservoir and a groundwater reservoir, which store and transport water. The water transport mechanisms operating within these two reservoirs are differentiated by their time scales; the groundwater reservoir transports water much more slowly than the surface reservior. The groundwater reservior feeds the corresponding surface store, and the surface stores are connected via the river network. The routing model is implemented over the Global Energy and Water Cycle Experiment (GEWEX) Continental-Scale International Project Mississippi River basin on a rectangular grid of 2 deg X 2.5 deg. Two land surface hydrology parameterizations provide the gridded runoff data required to run the runoff routing scheme: the variable infiltration capacity model, and the soil moisture component of the simple biosphere model. These parameterizations are driven with 4 deg X 5 deg gridded climatological potential evapotranspiration and 1979 First Global Atmospheric Research Program (GARP) Global Experiment precipitation. These investigations have quantified the importance of physically realistic soil moisture holding capacities, evaporation parameters, and runoff mechanisms in land surface hydrology formulations.

Liston, G. E.↗

Spaceborne radar observations: A guide for Magellan radar-image analysis

Geologic analyses of spaceborne radar images of Earth are reviewed and summarized with respect to detecting, mapping, and interpreting impact craters, volcanic landforms, eolian and subsurface features, and tectonic landforms. Interpretations are illustrated mostly with Seasat synthetic aperture radar and shuttle-imaging-radar images. Analogies are drawn for the potential interpretation of radar images of Venus, with emphasis on the effects of variation in Magellan look angle with Venusian latitude. In each landform category, differences in feature perception and interpretive capability are related to variations in imaging geometry, spatial resolution, and wavelength of the imaging radar systems. Impact craters and other radially symmetrical features may show apparent bilateral symmetry parallel to the illumination vector at low look angles. The styles of eruption and the emplacement of major and minor volcanic constructs can be interpreted from morphological features observed in images. Radar responses that are governed by small-scale surface roughness may serve to distinguish flow types, but do not provide unambiguous information. Imaging of sand dunes is rigorously constrained by specific angular relations between the illumination vector and the orientation and angle of repose of the dune faces, but is independent of radar wavelength. With a single look angle, conditions that enable shallow subsurface imaging to occur do not provide the information necessary to determine whether the radar has recorded surface or subsurface features. The topographic linearity of many tectonic landforms is enhanced on images at regional and local scales, but the detection of structural detail is a strong function of illumination direction. Nontopographic tectonic lineaments may appear in response to contrasts in small-surface roughness or dielectric constant. The breakpoint for rough surfaces will vary by about 25 percent through the Magellan viewing geometries from low to high Venusian latitudes. Examples of anomalies and system artifacts that can affect image interpretation are described.

Ford, J. P.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Cohen, S. C.↗

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Cohen, S. C.↗

Manipulating Neutral Atoms in Chip-Based Magnetic Traps

Several techniques for manipulating neutral atoms (more precisely, ultracold clouds of neutral atoms) in chip-based magnetic traps and atomic waveguides have been demonstrated. Such traps and waveguides are promising components of future quantum sensors that would offer sensitivities much greater than those of conventional sensors. Potential applications include gyroscopy and basic research in physical phenomena that involve gravitational and/or electromagnetic fields. The developed techniques make it possible to control atoms with greater versatility and dexterity than were previously possible and, hence, can be expected to contribute to the value of chip-based magnetic traps and atomic waveguides. The basic principle of these techniques is to control gradient magnetic fields with suitable timing so as to alter a trap to exert position-, velocity-, and/or time-dependent forces on atoms in the trap to obtain desired effects. The trap magnetic fields are generated by controlled electric currents flowing in both macroscopic off-chip electromagnet coils and microscopic wires on the surface of the chip. The methods are best explained in terms of examples. Rather than simply allowing atoms to expand freely into an atomic waveguide, one can give them a controllable push by switching on an externally generated or a chip-based gradient magnetic field. This push can increase the speed of the atoms, typically from about 5 to about 20 cm/s. Applying a non-linear magnetic-field gradient exerts different forces on atoms in different positions a phenomenon that one can exploit by introducing a delay between releasing atoms into the waveguide and turning on the magnetic field.

Aveline, David↗

On the stability of the boundary layer on a transonic swept wing

Both incompressible and compressible linear stability theory are applied to the three-dimensional compressible boundary layer on a particular transonic sweptback wing of infinite span. A spatial stability theory is used which identifies the growth direction with the real part of the complex angle of the group velocity. It is found that in the forward, but not the rear, crossflow instability region, the maximum amplification rates of the steady disturbances may be calculated to within about 10% by the incompressible stability theory. There is little difference between the sixth and eighth-order compressible theories. The maximum amplification rate of the steady disturbances at any chordwise station is closely related to the maximum crossflow at that station independent of the Reynolds number. For other than crossflow instability, there can be large differences between the incompressible and compressible theories, both as to the amplification rate and the angle of the wavenumber vector for maximum instability. Amplitude ratios of individual wave components are obtained by integrating the spatial amplification rate along the growth direction subject to the constraint that the wavenumber vector is irrotational. This procedure yields steady disturbances aligned with the local potential flow direction whose wavelengths are nearly independent of downstream distance.

Mack, L. M.↗

Individual Variability and Predictive Biomarkers of Human Sensitivity to Ionizing Radiation

High-linear energy transfer (LET) ionizing radiation is a major health hazard for astronauts who will be exposed to galactic cosmic rays during upcoming lunar and Mars missions. Predicting and mitigating personalized health risks from space radiation exposure requires understanding the factors underlying individual radiation sensitivity, which currently remain unknown. We started to address this challenge by identifying the genomic and demographic associations with sensitivity to low and high-LET ionizing radiation in the largest human sample cohort to date and by evaluating baseline DNA damage as a potential predictor of radiation responses. We have collected and irradiated peripheral blood mononuclear cells (PBMCs) from 768 healthy donors of matched ethnicity, 50/50 male/female, 18-70 years old. We have analyzed the responses of 576 donor PBMCs to irradiation with simulated galactic cosmic ray components: 350MeV/n 28Si, 350MeV/n 40Ar and 600MeV/n 56Fe, at 1.1 and 3 particles/100μm2 fluences, as well as 0.1Gy and 1Gy doses of gamma rays, at 4 and 24 hours post irradiation. We quantified post-irradiation DNA repair based on radiation-induced 53BP1+ foci formation, together with oxidative stress and cell death using a flow cytometry-based assay. We also assessed baseline levels of DNA repair in PBMCs at the time of sample collection prior to irradiation experiments. We observed a wide variability of subject- and LET-dependent radiation responses, with radiation-induced DNA repair foci increasing with LET, though oxidative stress being notably reduced by high-LET irradiation, potentially due to a switch between hydrogen peroxide and oxygen radical-based mechanisms. We also identified a relationship between few DNA repair foci at baseline and increased DNA repair after irradiation, accompanied by an alteration in immunoregulatory cytokine secretion, which might be adapted as biomarkers to predict ionizing radiation sensitivity. Among demographic variables, baseline foci levels were positively associated with age and latent cytomegalovirus infection. Finally, we have performed low-throughput whole genome sequencing of all samples and are currently in process of identifying the genes and pathways associated with low and high-LET ionizing radiation sensitivity in humans.

Space radiation, human spaceflight, DNA damage↗

Lee waves: Benign and malignant

The flow of an incompressible fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as dynamically incompressible, even though there exists a mean static upward density gradient. Even in the linear approximation - i.e., for small disturbances - this model explains a great many of the flow phenomena observed in the lee of mountains. However, nonlinearities do arise importantly, in three ways: (1) through amplification due to the decrease of mean density with height; (2) through the large (scaled) size of the obstacle, such as a mountain range; and (3) from dynamically singular levels in the fluid field. These effects produce a complicated array of phenomena - large departure of the streamlines from their equilibrium levels, high winds, generation of small scales, turbulence, etc. - that present hazards to aircraft and to lee surface areas. The nonlinear disturbances also interact with the larger-scale flow in such a manner as to impact global weather forecasts and the climatological momentum balance. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km), where recent observations show them to be of a length scale that must involve the coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the phenomena associated with these wave dynamics are being studied with a view to their potential impact on high performance aircraft, including the projected National Aerospace Plane (NASP). The presentation shows the results of analysis and of state-of-the-art numerical simulations, validated where possible by observational data, and illustrated with photographs from nature.

Wurtele, M. G.↗

Land-Atmosphere Interactions: Successes, Problems and Prospects

After two decades of active research, a much better understanding of the broader role of biospheric processes on the local climate has emerged. A surface-albedo increase, particularly in desert border regions of the subtropics (as well as the deforested tropical regions), leads to a net surface energy deficit, which in turn leads to a relative sinking and reduced rainfall. On the other hand, studies of the influence of altered ratios of evapotranspiration and sensible fluxes, in situations where the net solar income is unchanged, show that evapotranspiration is a more desirable flux for increased precipitation and vitality of the biosphere. Besides providing water vapor and convective available potential energy (CAPE) to the lower troposphere, evapotranspiration helps in building larger CAPE before "turning on" the moist-convection. Larger CAPE in the lower troposphere enables convection to reach into the deeper atmosphere thereby heating the upper troposphere; indeed, moist-convection is also accompanied by the evaporation of falling precipitation that cools and moistens the lower atmosphere. While convective, as opposed to stratiform, precipitation reduces the fractional cloud cover; it also allows more solar radiation to reach the surface thereby invigorating surface fluxes. These, together with moist convection and associated downdrafts help to maintain the characteristic upper temperature limit(s) of the moist-land as well as oceanic regions. Regardless of the above understanding, several important problems continue to hinder the accurate simulation of a realistic land atmosphere interaction in a numerical model (both GCM and/or Meso-scale models). Among the unsolved problems are parameterization of sub-grid scale land processes that include small-scale variability of soil moisture, snow-cover and snow-physics, the biodiversity of the biosphere, orography, local drainage characteristics under natural conditions, and surface flow over the natural terrain. A well-known non-linear response of surface fluxes to these variations makes the problem of parameterizing land-atmosphere interaction processes hard-to-address and simulate, particularly in a GCM. In our presentation, we will discuss how orographic, snow-cover, and water table interactions can be included into a Simple Biosphere Model such as SiB/SSiB. Figure I shows how, in the Russian region, spring snowmelt affects the soil moisture profile. Corresponding figure 2 shows how interaction with the water table decreases the natural evapotranspiration in the Sahel region simulation. While these simulations need better validation with data, the simulations reveal that surface processes are sensitive to these parameterizations. With these developments, we continue to advance our understanding of the interaction of land with the atmosphere aloft, but the intrinsic variability of the newer parameters, e. g., hydraulic properties of the soil, diminish the positive influences of these advances on the improved climate simulation with GCMs.

Sud, Y. C.↗

Getting granite dikes out of the source region

Whether a dike can propagate far from a magma reservoir depends upon the competition between the rate at which propagation widens the dike and the rate at which freezing constricts the aperture available for magma flow. Various formulations are developed for a viscous fluid at temperature T(sub m) intruding a growing crack in an elastic solid. The initial solid temperature equals T(sub m) at the source and decreases linearly with distance from the source. If T(sub m) is the unique freezing temperature of the fluid, dike growth is initially self-similar and an essentially exact solution is obtained; if T(sub m) is above the solidus temperature, the solution is approximate but is designed to overestimate the distance the dike may propagate. The ability of a dike to survive thermally depends primarily upon a single parameter that is a measure of the ratio of the dike frozen margin thickness to elastic thickness. Perhaps more intuitively, one may define a minimum distance from the essentially solid reservoir wall to the point at which the host rock temperature drops below the solidus, necessary for dikes to propagate far into subsolidus rock. It is concluded that for reasonable material properties and source conditions, most basalt dikes will have little difficulty leaving the source region, but most rhyolite dikes will be halted by freezing soon after the magma encounters rock at temperatures below the magma solidus. While these results can explain why granitic dikes are common near granitic plutons but rare elsewhere, the potentially large variation in magmatic systems makes it premature to rule out the possibility that most granites are transported through the crust in dikes. Nonetheless, these results highlight difficulties with such proposals and suggest that it may also be premature to rule out the possibility that most granite plutons ascend as more equidimensional bodies.

Rubin, Allan M.↗

Considerations on Electrolytic Conductivity Measurement for Monitoring of Ionic Silver Biocide Dosing

NASA interest in ionic silver (Ag+) as a biocide for spacecraft potable water systems motivates the development of Ag+ concentration sensors to ensure nominal dosing. The electrolytic conductivity change of highly-purified potable water is linearly related to the concentration of chemically-dosed Ag+ and could serve as a useful proxy measurement, while the conductivity change during electrolytic dosing may be less so, depending on the influent water chemistry and electrolytic efficiency. Understanding and mitigating the potentially deleterious effects of Ag+ interaction with conductivity measurement systems requires investigation. Issues associated with traditional conductivity cells, which rely on wetted(typically metal or graphite) electrodes, and capacitively-coupled contactless conductivity detection (C4D) for this application are considered. Traditional conductivity cells may potentially be subject to significant excitation-induced or auto-galvanic reduction of Ag+. These may result in Ag+ depletion or electrode fouling and associated measurement error. Proper selection of electrode material, excitation parameters, and cell geometry may limit such effects. C4D uses electrodes placed outside an inert, insulating material, with dielectric polarization enabling the production of an electric field and resultant current across the analyte solution. This approach could potentially mitigate problems with Ag+ depletion and fouling by eliminating the possibility of auto-galvanic deposition and reducing the Faradaic current density. However, it is necessary to confirm C4D performance at very low conductivity levels and to determine if long-term operation produces conductive deposits, which could result in measurement error. A commercial C4D system with 1/16” (1.59 mm)outer diameter flow tubing and claimed performance in the conductivity range appropriate for this application was identified, and its sensitivity in the upper parts of this range was confirmed in a preliminary experiment. A concept for a C4D detector that could allow for full rate flow-through using planar electrodes and thin-film dielectric layers is discussed.

silver↗

Particle Loading Tests on HEPA Flat Sheet Media at Sub-Ambient Pressures Using a Lunar Dust Simulant

When humans return to the moon under the NASA Artemis program, their activities on the lunar surface will inevitably lead to the intrusion of some level of lunar dust into the lander cabin. Therein the crew would be exposed to the potential hazards of lunar dust and the possibility of subsequent transfer into orbiting segments after docking. The spacecraft’s cabin filtration system will need to be effective at removing the airborne lunar dust to properly purify the breathable cabin air and minimize the spreading of the dust throughout the vehicle and orbital segments through the mission duration. The fine nature of the lunar dust will require high efficiency filtration, such as HEPA. A series of tests were performed in a specially designed recirculating sealed flow loop, for testing filter media and filter elements at the NASA Glenn Research Center. The flow loop was used to assess the performance and capacity of flat sheet filter media at representative cabin pressures using JSC-1AF lunar dust simulant and at high rates of particle loading. The pressure drop across the filter media was measured as a function of accumulated particle mass load at ambient pressure and at two sub-ambient pressures, 0.0703MPa and 0.0565MPa, and at a media velocity that was scaled relative to its pleated configuration. The challenge particle flows were generated by a custom designed particle generator that introduces dispersed particles of the lunar simulant at high concentrations. An optical particle counting instrument provided filter efficiency measurements within the sealed environment. The pressure drop was found to increase linearly with the amount of dust load on the media, for all test conditions, while the starting pressure drop was found to be lower at the lowest sub-ambient pressure case. High filter efficiency was maintained after high particle loads on the media.

Sub-ambient↗