Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 811 records · Page 45

Fluorescence Studies of Protein Crystal Nucleation

We have postulated that, in the case of tetragonal chicken egg white lysozyme, crystal growth occurs by the addition of pre-critical nuclei sized n-mers that form in the bulk solution, and that the n-mer growth units were multiples of the tetrameric 4(sub 3) helical structure. These have the strongest intermolecular bonds in the crystal and are therefore likely to be the first species formed. High resolution AFM studies provide strong supporting evidence for this model, but the data also suggest that the actual species in solution may not be identical in structure to that found in the crystal. We are using fluorescence resonance energy transfer (FRET) to study the initial solution phase self-assembly process, using covalent fluorescent derivatives which crystallize in the characteristic P4(sub 3)2(sub 1)2(sub 1) space group. FRET studies are being carried out between the cascade blue (CB-lys, donor, Ex(sub max) 366 nm, Em 420 nm) and lucifer yellow (LY-lys, acceptor, Ex(sub max) 430 nm, Em 528 nm) asp101 derivatives. The estimated R(sub 0) for this probe pair, the distance where 50% of the donor energy is transferred to the acceptor, is approx. 1.2 nm, compared to 2.2 nm between the side chain carboxyls of adjacent asp101's in the crystalline 4(sub 3) helix. The short donor lifetime of 2.80 ns (chi(sup 2) = 0.644), coupled with the large average distances between the molecules (greater than or equal to 50 nm) in solution, ensure that any energy transfer observed is not due to random diffusive interactions. Lifetime data show that CB-lys has a single lifetime when it is the only species in solution. Similarly, LY-lys also exhibits a single lifetime of 4.63 ns (chi(sup 2) = 0.42) when alone in solution. Addition of LY-lys to CB-lys results in the appearance of a third lifetime component of 0.348ns for the CB-lys. The fractional intensities of the different species present can be used to estimate the distribution of monomer and n-mers in solution. The self-association process is a function of the protein concentration relative to the saturation concentration, and observing it in dilute solution (conc. less than or equal to 10(exp -5)M) requires that the experiments be performed under low solubility conditions, i.e., low temperatures and high salt concentrations. Data from preliminary steady state FRET studies with N-terminal bound pyrene acetic acid (PAA-lys, donor, Ex 340 nm, Em 376 nm) and asp101 LY-lys as an acceptor showed a consistent trend of decreasing donor fluorescence intensity with increasing total protein concentration. The FRET data have been obtained at pH 4.6, 0.1M NaAc buffer, at 5 and 7% NaCl, 4 C. The corresponding C(sub sat) values are 0.471 and 0.362 mg/ml (approx. 3.3 and approx. 2.5 x 10(exp -5)M respectively). The donor fluorescence decrease is more pronounced at7% NaCl, consistent with the expected increased intermolecular interactions at higher salt concentrations as reflected in the lower solubility. Results from these and other ongoing studies will be discussed in conjunction with an emerging model for how tetragonal lysozyme crystals nucleate and the relevance of that model to other proteins.

Pusey, Marc↗

Thermal Edge-Effects Model for Automated Tape Placement of Thermoplastic Composites

Two-dimensional thermal models for automated tape placement (ATP) of thermoplastic composites neglect the diffusive heat transport that occurs between the newly placed tape and the cool substrate beside it. Such lateral transport can cool the tape edges prematurely and weaken the bond. The three-dimensional, steady state, thermal transport equation is solved by the Green's function method for a tape of finite width being placed on an infinitely wide substrate. The isotherm for the glass transition temperature on the weld interface is used to determine the distance inward from the tape edge that is prematurely cooled, called the cooling incursion Delta a. For the Langley ATP robot, Delta a = 0.4 mm for a unidirectional lay-up of PEEK/carbon fiber composite, and Delta a = 1.2 mm for an isotropic lay-up. A formula for Delta a is developed and applied to a wide range of operating conditions. A surprise finding is that Delta a need not decrease as the Peclet number Pe becomes very large, where Pe is the dimensionless ratio of inertial to diffusive heat transport. Conformable rollers that increase the consolidation length would also increase Delta a, unless other changes are made, such as proportionally increasing the material speed. To compensate for premature edge cooling, the thermal input could be extended past the tape edges by the amount Delta a. This method should help achieve uniform weld strength and crystallinity across the width of the tape.

Costen, Robert C.↗

A Coupled Evaluation of Operational MODIS and Model Aerosol Products for Maritime Environments Using Sun Photometry: Evaluation of the Fine and Coarse Mode

Although satellite retrievals and data assimilation have progressed to where there is a good skill for monitoring maritime Aerosol Optical Depth (AOD), there remains uncertainty in achieving further degrees of freedom, such as distinguishing fine and coarse mode dominated species in maritime environments (e.g., coarse mode sea salt and dust versus fine mode terrestrial anthropogenic emissions, biomass burning, and maritime secondary production). For the years 2016 through 2019, we performed an analysis of 550 nm total AOD550, fine mode AOD (FAOD550; also known as FM AOD in the literature), coarse mode AOD (CAOD550), and fine mode fraction (η550) between Moderate Resolution Spectral Imaging Radiometer (MODIS) V6.1 MOD/MYD04 dark target aerosol retrievals and the International Cooperative for Aerosol Prediction (ICAP) core four multi-model consensus (C4C) of analyses/short term forecasts that assimilate total MODIS AOD550. Differences were adjudicated by the global shipboard Maritime Aerosol Network (MAN) and selected island AERONET sun photometer observations with the application of the spectral deconvolution algorithm (SDA). Through a series of conditional and regional analyses, we found divergence included regions of terrestrial influence and latitudinal dependencies in the remote oceans. Notably, MODIS and the C4C and its members, while having good correlations overall, have a persistent +0.04 to +0.02 biases relative to MAN and AERONET for typical AOD550 values (84th% < 0.28), with the C4C underestimating significant events thereafter. Second, high biases in AOD550 are largely associated with the attribution of the fine mode in satellites and models alike. Thus, both MODIS and C4C members are systematically overestimating AOD550 and FAOD550 but perform better in characterizing the CAOD550. Third, for MODIS, findings are consistent with previous reports of a high bias in the retrieved Ångström Exponent, and we diagnosed both the optical model and cloud masking as likely causal factors for the AOD550 and FAOD550 high bias, whereas for the C4C, it is likely from secondary overproduction and perhaps numerical diffusion. Fourth, while there is no wind-speed-dependent bias for surface winds <12 m s−1, the C4C and MODIS AOD550s also overestimate CAOD550 and FAOD550, respectively, for wind speeds above 12 m/s. Finally, sampling bias inherent in MAN, as well as other circumstantial evidence, suggests biases in MODIS are likely MODIS and the C4C products have their own strengths and challenges for a given climate application and discuss needed research. even larger than what was diagnosed here. We conclude with a discussion on how

Aerosols↗

SAO and Kelvin Waves in the EuroGRIPS GCMS and the UK Meteorological Offices Analyses

This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.

Amodei, M.↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made substantial progress toward the goal of producing a unified model of the basic physical processes responsible for solar wind acceleration. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a 1D model of plasma heating and acceleration. The accomplishments in Year 2 are divided into these two categories: 1a. Focused Study of Kinetic Magnetohydrodynamic (MHD) Turbulence. lb. Focused Study of Non - WKB Alfven Wave Rejection. and 2. The Unified Model Code. We have continued the development of the computational model of a time-study open flux tube in the extended corona. The proton-electron Monte Carlo model is being tested, and collisionless wave-particle interactions are being included. In order to better understand how to easily incorporate various kinds of wave-particle processes into the code, the PI performed a detailed study of the so-called "Ito Calculus", i.e., the mathematical theory of how to update the positions of particles in a probabilistic manner when their motions are governed by diffusion in velocity space.

Cranmer, Steven R.↗

The 3D Navier-Stokes analysis of a Mach 2.68 bifurcated rectangular mixed-compression inlet

The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a three-dimensional (3D) Navier-Stokes flow solver. A two-equation turbulence model, and a porous bleed model based on unchoked bleed hole discharge coefficients were used. Comparisons were made with experimental data, inviscid theory, and two-dimensional Navier-Stokes analyses. The main objective was to gain insight into the inlet fluid dynamics. Examination of the computational results along with the experimental data suggest that the cowl shock-sidewall boundary layer interaction near the leading edge caused a substantial separation in the wind tunnel inlet model. As a result, the inlet performance may have been compromised by increased spillage and higher bleed mass flow requirements. The internal flow contained substantial waves that were not in the original inviscid design. 3D effects were fairly minor for this inlet at on-design conditions. Navier-Stokes analysis appears to be an useful tool for gaining insight into the inlet fluid dynamics. It provides a higher fidelity simulation of the flowfield than the original inviscid design, by taking into account boundary layers, porous bleed, and their interactions with shock waves.

Mizukami, M.↗

Boundary-Layer-Ingesting Inlet Flow Control

This paper gives an overview of a research study conducted in support of the small-scale demonstration of an active flow control system for a boundary-layer-ingesting (BLI) inlet. The effectiveness of active flow control in reducing engine inlet circumferential distortion was assessed using a 2.5% scale model of a 35% boundary-layer-ingesting flush-mounted, offset, diffusing inlet. This experiment was conducted in the NASA Langley 0.3-meter Transonic Cryogenic Tunnel at flight Mach numbers with a model inlet specifically designed for this type of testing. High mass flow actuators controlled the flow through distributed control jets providing the active flow control. A vortex generator point design configuration was also tested for comparison purposes and to provide a means to examine a hybrid vortex generator and control jets configuration. Measurements were made of the onset boundary layer, the duct surface static pressures, and the mass flow through the duct and the actuators. The distortion and pressure recovery were determined by 40 total pressure measurements on 8 rake arms each separated by 45 degrees and were located at the aerodynamic interface plane. The test matrix was limited to a maximum free-stream Mach number of 0.85 with scaled mass flows through the inlet for that condition. The data show that the flow control jets alone can reduce circumferential distortion (DPCPavg) from 0.055 to about 0.015 using about 2.5% of inlet mass flow. The vortex generators also reduced the circumferential distortion from 0.055 to 0.010 near the inlet mass flow design point. Lower inlet mass flow settings with the vortex generator configuration produced higher distortion levels that were reduced to acceptable levels using a hybrid vortex generator/control jets configuration that required less than 1% of the inlet mass flow.

Owens, Lewis R.↗

Boundary-Layer-Ingesting Inlet Flow Control

This paper gives an overview of a research study conducted in support of the small-scale demonstration of an active flow control system for a boundary-layer-ingesting (BLI) inlet. The effectiveness of active flow control in reducing engine inlet circumferential distortion was assessed using a 2.5% scale model of a 35% boundary-layer-ingesting flush-mounted, offset, diffusing inlet. This experiment was conducted in the NASA Langley 0.3-meter Transonic Cryogenic Tunnel at flight Mach numbers with a model inlet specifically designed for this type of testing. High mass flow actuators controlled the flow through distributed control jets providing the active flow control. A vortex generator point design configuration was also tested for comparison purposes and to provide a means to examine a hybrid vortex generator and control jets configuration. Measurements were made of the onset boundary layer, the duct surface static pressures, and the mass flow through the duct and the actuators. The distortion and pressure recovery were determined by 40 total pressure measurements on 8 rake arms each separated by 45 degrees and were located at the aerodynamic interface plane. The test matrix was limited to a maximum free-stream Mach number of 0.85 with scaled mass flows through the inlet for that condition. The data show that the flow control jets alone can reduce circumferential distortion (DPCP(sub avg)) from 0.055 to about 0.015 using about 2.5% of inlet mass flow. The vortex generators also reduced the circumferential distortion from 0.055 to 0.010 near the inlet mass flow design point. Lower inlet mass flow settings with the vortex generator configuration produced higher distortion levels that were reduced to acceptable levels using a hybrid vortex generator/control jets configuration that required less than 1% of the inlet mass flow.

Owens, Lewis R.↗

Coupled effects of conduction in the crystal and thermo-solutal convection in a rectangular inclined enclosure

To date modeling of crystal growth of optoelectronic materials using Physical Vapor Transport has been limited to the study of the fluid phase. To achieve it, the equations of coupled heat, mass and momentum transfer in the gas have to be solved. The first objective of this study is to examine the effect of heat conduction in the crystal on the fluid flow in the neighborhood of the interface. Heat transfer boundary conditions on both interfaces were modified to take into account the additional heat flux between gas and solid. It is proved that heat conduction does not affect the fluid flow. In the presence of gravity, density gradients in the fluid phase generate convection responsible for the problem of a nonplanar growth of the interface. The second objective is to study systematically under one-g the different possible flows in order to solve this problem. Depending on the parameters, a diffusive mode and three convective modes (thermal, solutal and thermo-solutal) are observed. The competition between thermal and solutal convections leads to a mathematical condition which can be used to achieve a planar growth. It is proven that, under the physical conditions chosen, this mathematical condition cannot be thermodynamically satisfied.

Mennetrier, Christophe↗

Simulation of fluid flows during growth of organic crystals in microgravity

Several counter diffusion type crystal growth experiments were conducted in space. Improvements in crystal size and quality are attributed to reduced natural convection in the microgravity environment. One series of experiments called DMOS (Diffusive Mixing of Organic Solutions) was designed and conducted by researchers at the 3M Corporation and flown by NASA on the space shuttle. Since only limited information about the mixing process is available from the space experiments, a series of ground based experiments was conducted to further investigate the fluid dynamics within the DMOS crystal growth cell. Solutions with density differences in the range of 10 to the -7 to 10 to the -4 power g/cc were used to simulate microgravity conditions. The small density differences were obtained by mixing D2O and H2O. Methylene blue dye was used to enhance flow visualization. The extent of mixing was measured photometrically using the 662 nm absorbance peak of the dye. Results indicate that extensive mixing by natural convection can occur even under microgravity conditions. This is qualitatively consistent with results of a simple scaling analysis. Quantitave results are in close agreement with ongoing computational modeling analysis.

Roberts, Gary D.↗

Generic hypersonic vehicle performance model

An integrated computational model of a generic hypersonic vehicle was developed for the purpose of determining the vehicle's performance characteristics, which include the lift, drag, thrust, and moment acting on the vehicle at specified altitude, flight condition, and vehicular configuration. The lift, drag, thrust, and moment are developed for the body fixed coordinate system. These forces and moments arise from both aerodynamic and propulsive sources. SCRAMjet engine performance characteristics, such as fuel flow rate, can also be determined. The vehicle is assumed to be a lifting body with a single aerodynamic control surface. The body shape and control surface location are arbitrary and must be defined. The aerodynamics are calculated using either 2-dimensional Newtonian or modified Newtonian theory and approximate high-Mach-number Prandtl-Meyer expansion theory. Skin-friction drag was also accounted for. The skin-friction drag coefficient is a function of the freestream Mach number. The data for the skin-friction drag coefficient values were taken from NASA Technical Memorandum 102610. The modeling of the vehicle's SCRAMjet engine is based on quasi 1-dimensional gas dynamics for the engine diffuser, nozzle, and the combustor with heat addition. The engine has three variable inputs for control: the engine inlet diffuser area ratio, the total temperature rise through the combustor due to combustion of the fuel, and the engine internal expansion nozzle area ratio. The pressure distribution over the vehicle's lower aft body surface, which acts as an external nozzle, is calculated using a combination of quasi 1-dimensional gas dynamic theory and Newtonian or modified Newtonian theory. The exhaust plume shape is determined by matching the pressure inside the plume, calculated from the gas dynamic equations, with the freestream pressure, calculated from Newtonian or Modified Newtonian theory. In this manner, the pressure distribution along the vehicle after body expansion surface is then determined. The aerodynamic modeling, the engine modeling, and the exhaust plume analysis are described in more detail. A description of the computer code used to perform the above calculations is given and an input/output example is then given. The computer code is available on a Macintosh floppy disk.

Chavez, Frank R.↗

Bridgman crystal growth

The objective of this theoretical research effort was to improve the understanding of the growth of Pb(x)Sn(1-x)Te and especially how crystal quality could be improved utilizing the microgravity environment of space. All theoretical growths are done using the vertical Bridgman method. It is believed that improved single crystal yields can be achieved by systematically identifying and studying system parameters both theoretically and experimentally. A computational model was developed to study and eventually optimize the growth process. The model is primarily concerned with the prediction of the thermal field, although mass transfer in the melt and the state of stress in the crystal were of considerable interest. The evolution is presented of the computer simulation and some of the important results obtained. Diffusion controlled growth was first studied since it represented a relatively simple, but nontheless realistic situation. In fact, results from this analysis prompted a study of the triple junction region where the melt, crystal, and ampoule wall meet. Since microgravity applications were sought because of the low level of fluid movement, the effect of gravitational field strength on the thermal and concentration field was also of interest. A study of the strength of coriolis acceleration on the growth process during space flight was deemed necessary since it would surely produce asymmetries in the flow field if strong enough. Finally, thermosolutal convection in a steady microgravity field for thermally stable conditions and both stable and unstable solutal conditions was simulated.

Carlson, Frederick↗

Trioxane-Air Counterflow Diffusion Flames in Normal and Microgravity

Trioxane, a weakly bound polymer of formaldehyde (C3H6O3, m.p. 61 C, b.p. 115 C), is a uniquely suited compound for studying material flammability. Like many of the more commonly used materials for such tests (e.g., delrin, polyethylene, acrylic sheet, wood, and paper), it displays relevant phenomena (internal heat conduction, melting, vaporization, thermal decomposition, and gas phase reaction of the decomposition products). Unlike the other materials, however, it is non-sooting and has simple and well-known chemical kinetic pathways for its combustion. Hence it should prove to be much more useful for numerical modeling of surface combustion than the complex fuels typically used. We have performed the first exploratory tests of trioxane combustion in the counterflow configuration to determine its potential as a surrogate solid fuel which allows detailed modeling. The experiments were performed in the spring and summer of 1998 at the National Institute of Standards and Technology in Gaithersburg, MD, and at NASA-GRC in Cleveland. Using counterflow flames at 1-g, we measured the fuel consumption rate and the extinction conditions with added N2 in the air; at mg conditions, we observed the ignition characteristics and flame shape from video images. We have performed numerical calculations of the flame structure, but these are not described here due to space limitations. This paper summarizes some burning characteristics of trioxane relevant to its use for studying flame spread and fire suppression.

Linteris, Gregory T.↗

Effects of a Forward-swept Front Rotor on the Flowfield of a Counterrotation Propeller

The effects of a forward-swept front rotor on the flowfield of a counterrotation model propeller at takeoff conditions at zero degree angle of attack are studied by solving the unsteady three-dimensional Euler equations. The configuration considered is an uneven blade count counterrotation model with twelve forward-swept blades on the fore rotor and ten aft-swept blades on the aft rotor. The flowfield is compared with that of a reference aft-swept counterrotation geometry and Laser Doppler Velocimeter (LDV) measurements. At the operating conditions considered, the forward-swept blade experiences a higher tip loading and produces a stronger tip vortex compared to the aft-swept blade, consistent with the LDV and acoustic measurements. Neither the solution nor the LDV data indicated the formation of a leading edge vortex. The predicted radial distribution of the circumferentially averaged axial velocity at the measurement station agreed very closely with LDV data, while crossflow velocities showed poor agreement. The discrepancy between prediction and LDV data of tangential and radial velocities is due in part to the insufficient mesh resolution in the region between the rotors and in the tip region to track the tip vortex. The vortex is diffused by the time it arrives at the measurement station. The uneven blade count configuration requires the solution to be carried out for six blade passages of the fore rotor and five passages of the aft rotor, thus making grid refinement prohibitive.

Nallasamy, M.↗

Geologic studies of outer solar system satellites: Iapetus and Europa

A model of regolith evolution on bodies subjected to an asymmetric impact flux is presented. The general effect of an asymmetric impact flux is ballistic diffusion of regolith from areas of high impact flux to areas of lower flux. The effect is most pronounced on low gravity bodies where ejecta travel distances are large, and on bodies with strong flux asymmetries. It may be especially important on Iapetus, where retrograde debris from Phoebe striking the leading hemisphere may produce a flux asymmetry as large as 100 to 1 from apex to antapex. We find that the net amount of material transported by ballistic diffusion is inversely proportional to satellite density. Theoretical calculations and Voyager observations support the view that tidal heating is sufficient to maintain a liquid layer throughout the history of Europa. Photosynthetic, thermal, electrical, and chemical energy sources are considered. The calculations suggest that there may be regions on Europa, very limited in space and time, with physical conditions that are within the range of adaptation of life on Earth. Calculations for the electrical currents within Europa for a thick ice crust are very small.

Squyres, S. W.↗

Direct simulation of isothermal-wall supersonic channel flow

The motivation for this work is the fact that in turbulent flows where compressibility effects are important, they are often poorly understood. A few examples of such flows are those associated with astrophysical phenomena and those found in combustion chambers, supersonic diffusers and nozzles, and over high-speed airfoils. For this project, we are primarily interested in compressibility effects near solid surfaces. Our main objective is an improved understanding of the fundamentals of compressible wall-bounded turbulence, which can in turn be used to cast light upon modeling concepts such as the Morkovin hypothesis and the Van Driest transformation. To this end, we have performed a direct numerical simulation (DNS) study of supersonic turbulent flow in a plane channel with constant-temperature walls. All of the relevant spatial and temporal scales are resolved so that no sub grid scale or turbulence model is necessary. The channel geometry was chosen so that finite Mach number effects can be isolated by comparing the present results to well established incompressible channel data. Here the fluid is assumed to be an ideal gas with constant specific heats, constant Prandtl number, and power-law temperature-dependent viscosity. Isothermal-wall boundary conditions are imposed so that a statistically stationary state may be obtained. The flow is driven by a uniform (in space) body force (rather than a mean pressure gradient) to preserve stream wise homogeneity, with the body force defined so that the total mass flux is constant.

Coleman, Gary N.↗

Application of a two-stream radiative transfer model for leaf lignin and cellulose concentrations from spectral reflectance measurements, part 2

We used the Kubelka-Munk theory of diffuse spectral reflectance in layers to analyze influences of multiple chemical components in leaves. As opposed to empirical approaches to estimation of plant chemistry, the full spectral resolution of laboratory reflectance data was retained in an attempt to estimate lignin or other constituent concentrations from spectral band positions. A leaf water reflectance spectrum was derived from theoretical mixing rules, reflectance observations, and calculations from theory of intrinsic k- and s-functions. Residual reflectance bands were then isolated from spectra of fresh green leaves. These proved hard to interpret for composition in terms of simple two component mixtures such as lignin and cellulose. We next investigated spectral and dilution influences of other possible components (starch, protein). These components, among others, added to cellulose in hypothetical mixtures, produce band displacements similar to lignin, but will disguise by dilution the actual abundance of lignin present in a multicomponent system. This renders interpretation of band positions problematical. Knowledge of end-members and their spectra, and a more elaborate mixture analysis procedure may be called for. Good observational atmospheric and instrumental conditions and knowledge thereof are required for retrieval of expected subtle reflectance variations present in spectra of green vegetation.

Conel, James E.↗

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗