Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional random fields”

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 379 records · Page 21

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗

{sup 90}Y -PET imaging: Exploring limitations and accuracy under conditions of low counts and high random fraction

Purpose: {sup 90}Y -positron emission tomography (PET) imaging is becoming a recognized modality for postinfusion quantitative assessment following radioembolization therapy. However, the extremely low counts and high random fraction associated with {sup 90}Y -PET may significantly impair both qualitative and quantitative results. The aim of this work was to study image quality and noise level in relation to the quantification and bias performance of two types of Siemens PET scanners when imaging {sup 90}Y and to compare experimental results with clinical data from two types of commercially available {sup 90}Y microspheres. Methods: Data were acquired on both Siemens Biograph TruePoint [non-time-of-flight (TOF)] and Biograph microcomputed tomography (mCT) (TOF) PET/CT scanners. The study was conducted in three phases. The first aimed to assess quantification and bias for different reconstruction methods according to random fraction and number of true counts in the scan. The NEMA 1994 PET phantom was filled with water with one cylindrical insert left empty (air) and the other filled with a solution of {sup 90}Y . The phantom was scanned for 60 min in the PET/CT scanner every one or two days. The second phase used the NEMA 2001 PET phantom to derive noise and image quality metrics. The spheres and the background were filled with a {sup 90}Y solution in an 8:1 contrast ratio and four 30 min acquisitions were performed over a one week period. Finally, 32 patient data (8 treated with Therasphere{sup ®} and 24 with SIR-Spheres{sup ®}) were retrospectively reconstructed and activity in the whole field of view and the liver was compared to theoretical injected activity. Results: The contribution of both bremsstrahlung and LSO trues was found to be negligible, allowing data to be decay corrected to obtain correct quantification. In general, the recovered activity for all reconstruction methods was stable over the range studied, with a small bias appearing at extremely high random fraction and low counts for iterative algorithms. Point spread function (PSF) correction and TOF reconstruction in general reduce background variability and noise and increase recovered concentration. Results for patient data indicated a good correlation between the expected and PET reconstructed activities. A linear relationship between the expected and the measured activities in the organ of interest was observed for all reconstruction method used: a linearity coefficient of 0.89 ± 0.05 for the Biograph mCT and 0.81 ± 0.05 for the Biograph TruePoint. Conclusions: Due to the low counts and high random fraction, accurate image quantification of {sup 90}Y during selective internal radionuclide therapy is affected by random coincidence estimation, scatter correction, and any positivity constraint of the algorithm. Nevertheless, phantom and patient studies showed that the impact of number of true and random coincidences on quantitative results was found to be limited as long as ordinary Poisson ordered subsets expectation maximization reconstruction algorithms with random smoothing are used. Adding PSF correction and TOF information to the reconstruction greatly improves the image quality in terms of bias, variability, noise reduction, and detectability. On the patient studies, the total activity in the field of view is in general accurately measured by Biograph mCT and slightly overestimated by the Biograph TruePoint.

60 APPLIED LIFE SCIENCES↗

Power demand risk models on milling machines

The measurement of power demand risks in manufacturing power systems will benefit manufacturers and the wider society. If the risks can be characterized using manufacturing parameters, manufacturers can better control risks originating in those parameters, select less power-risky production plans, and reduce utility costs and resource consumption. The measurement of risk can also help manufacturers and power suppliers to protect their power systems from unexpected disturbances. Existing measures of risk, however, do not consider time duration, and thus cannot accurately quantify the risks in manufacturing power systems; the risks of a period of high power demand must be evaluated with the duration of the surge. Therefore, new methods of measuring power demand risks are proposed, adapting measures drawn from the field of finance. With a focus on milling operations, processing power is shown to be a function of processing amount and processing time, and a power demand distribution is directly derived as a joint distribution of and A and T. A bivariate random variable model with copulas is applied to examine the correlation in the joint distribution. Then, based on evaluation of a probability distribution of power demand from and, new risk measures are introduced. Herein illustrative examples are provided to show how the proposed measures can quantify the power demand risks from milling machines, based on manufacturing parameters. Certain manufacturing parameters are found to affect overall power demand risks, including i) raw material type, ii) variability in processing time, and iii) correlation between and. In the examples, these three factors increase power demand risks by up to 108%, 67%, and 1% respectively.

42 ENGINEERING↗

Stochastic Galerkin methods for the steady-state Navier–Stokes equations

We study the steady-state Navier–Stokes equations in the context of stochastic finite element discretizations. Specifically, we assume that the viscosity is a random field given in the form of a generalized polynomial chaos expansion. For the resulting stochastic problem, we formulate the model and linearization schemes using Picard and Newton iterations in the framework of the stochastic Galerkin method, and we explore properties of the resulting stochastic solutions. We also propose a preconditioner for solving the linear systems of equations arising at each step of the stochastic (Galerkin) nonlinear iteration and demonstrate its effectiveness for solving a set of benchmark problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗

A Critical Examination of Wind-Wave Spectral Functional Form

Traditionally, data from random ocean waves are presented in spectral functions. The spectra are the result of Fourier analysis. Fourier spectral analysis has dominated data analysis for, at least, the last hundred years. It has been the standard method for is examining the global amplitude-frequency distributions. Although Fourier transform valid under extremely general conditions, there are some crucial restrictions for the Fourier spectral analysis. The system must be linear, and the data must be stationary- otherwise, the resulting spectrum will make little physical sense. The stationarity requirement is also a common required criterion for most of other available data analysis methods. Nevertheless, few, if any, natural phenomena are linear and stationary. To compound these complications is the imperfection of our probes or numerical schemes the interactions of the imperfect probes even with a perfect linear system can make the final data nonlinear. Furthermore, all the available data are usually of finite duration. Under these conditions, Fourier analysis is of limited use, For lack of alternatives, however, Fourier analysis is still used to process such data. The loose application of Fourier analysis and the insouciant adoption of the stationary and linear assumptions may lead to misleading conclusions. Ocean waves are know to be nonlinear, and the wind system generating the wave field are seldom stationary- As a result, the traditional examination of the spectral form hardly made physical sense. A new method for analyzing nonlinear and nonstationary data has been developed. The key part is the Empirical Mode Decomposition (EMD) method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF) that serve as the basis of the representation of the data, This decomposition method is adaptive, and, therefore, highly efficient. The IMFs admit well-behaved Hilbert transforms, and yield instantaneous energy and frequency as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Among the main conceptual innovations is the introduction of the instantaneous frequencies for complicated data sets, which eliminate the need of spurious harmonics to represent nonlinear and nonstationary signals. The spectral forms of the ocean waves are very different. This new method offers the first physical meaningful view of the wave spectrum. Data from laboratory and field will be presented to illustrate the differences.

Huang, Norden E.↗

X-Ray Stereo Digital Image Correlation

Digital Image Correlation (DIC) is a well-established, non-contact diagnostic technique used to measure shape, displacement and strain of a solid specimen subjected to loading or deformation. Yet, measurements using standard DIC can have significant errors or be completely infeasible in challenging experiments, such as explosive, combustion, or fluid-structure interaction applications, where beam-steering due to index of refraction variation biases measurements or where the sample is engulfed in flames or soot. To address these challenges, we propose using X-ray imaging instead of visible light imaging for stereo-DIC, since refraction of X-rays is negligible in many situations, and X-rays can penetrate occluding material. Two techniques of creating an appropriate pattern for X-ray DIC are presented, both based on adding a dense material in a random speckle pattern on top of a less-dense specimen. A standard dot-calibration target is adapted for X-ray imaging, allowing the common bundle-adjustment calibration process in commercial stereo-DIC software to be used. High-quality X-ray images with sufficient signal-to-noise ratios for DIC are obtained for aluminum specimens with thickness up to 22.2 mm, with a speckle pattern thickness of only 80 μm of tantalum. The accuracy and precision of X-ray DIC measurements are verified through simultaneous optical and X-ray stereo-DIC measurements during rigid in-plane and out-of-plane translations, where errors in the X-ray DIC displacements were approximately 2–10 μm for applied displacements up to 20 mm. Finally, a vast reduction in measurement error—5–20 times reduction of displacement error and 2–3 times reduction of strain error—is demonstrated, by comparing X-ray and optical DIC when a hot plate induced a heterogeneous index of refraction field in the air between the specimen and the imaging systems. Subsequently, these findings show the feasibility of using X-ray-based stereo-DIC for non-contact measurements in exacting experimental conditions, where optical DIC cannot be used.

47 OTHER INSTRUMENTATION↗

Magnetism in Mixed Valence, Defect, Cubic Perovskites: BaIn 1– x Fe x O 2.5+δ , x = 0.25, 0.50, and 0.75. Local and Average Structures

The series BaIn 1–x Fe x O 2.5+δ , x = 0.25, 0.50, and 0.75, has been prepared under air-fired and argon-fired conditions and studied using X-ray diffraction, d.c. and a.c. susceptibility, Mössbauer spectroscopy, neutron diffraction, X-ray near edge absorption spectroscopy (XANES), and X-ray pair distribution (PDF) methods. While Ba 2 In 2 O 5 (BaInO 2.5 ) crystallizes in an ordered brownmillerite structure, Ibm2, and Ba 2 Fe 2 O 5 (BaFeO 2.5 ) crystallizes in a complex monoclinic structure, P2 1 /c, showing seven Fe 3+ sites with tetrahedral, square planar, and octahedral environments, all phases studied here crystallize in the cubic perovskite structure, Pm$\bar{}3}$m, with long-range disorder on the small cation and oxygen sites. 57 Fe Mössbauer studies indicate a mixed valency, Fe 4+ /Fe 3+ , for both the air-fired and argon-fired samples. The increased Fe 3+ content for the argon-fired samples is reflected in increased cubic cell constants and in the increased Mössbauer fraction. It appears that the Pm$\bar{}3}$m phases are only metastable when fired in argon. From a slightly modified percolation theory for a primitive cubic lattice (taking into account the presence of random O atom vacancies), long-range spin order is permitted for the x = 0.50 and 0.75 phases. Instead, the d.c. susceptibility shows only zero-field-cooled (ZFC) and field-cooled (FC) divergences at ~6 K [5 K] for x = 0.50 and at ~22 K [21 K] for x = 0.75, with values for the argon-fired samples in [ ]. Neutron diffraction data for the air-fired samples confirm the absence of long-range magnetic order at any studied temperature. For the air-fired x = 0.50, a.c. susceptibility data show a frequency-dependent χ'(max) and spin glass behavior, while for x = 0.75, χ'(max) is invariant with frequency, ruling out either a spin glass or a superparamagnetic ground state. These behaviors are discussed in terms of competing Fe 3+ –Fe 3+ antiferromagnetic exchange and ferromagnetic Fe 3+ –Fe 4+ exchange. The PDF and 57 Fe Mössbauer data indicate a local structure at short interatomic distances, which deviates strongly from the average Pm$\bar{}3}$m model. Fe Mössbauer, PDF, and XANES data show a systematic dependence on x and indicate that the Fe 3+ sites are largely fourfold-coordinated and Fe 4+ sites are fivefold- or sixfold-coordinated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility. Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson’s R and Spearman’s ρ coefficients as well as Kendall’s τ and Top–Down correlation. We also test the behavior of association measures, including the coefficient of determination R 2 , Kendall’s W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spear man’s ρ, Kendall’s τ, and Kendall’s W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R 2 coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships. Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

59 BASIC BIOLOGICAL SCIENCES↗

Application of a Planar Doppler Velocimetry System to a High Reynolds Number Compressible Jet

A Planar Doppler Velocimetry (PDV) system has been constructed and used to investigate the instantaneous turbulent velocity structure of a round high-speed compressible air jet with a low-speed co-flow. The exit condition was Mach=0.85 at ambient pressure, yielding a Reynolds number of about 650,000 on diameter. The PDV system was installed at NASA Langley Research Center in the Small Anechoic Jet Facility (SAJF), a chamber in which both the acoustic and aerodynamic properties of jets can be studied. For this test, the goal was to gather data which can be used to relate the turbulence structure of the jet to the levels and character of the acoustic noise produced by the jet. The current PDV system can acquire single-velocity-component, single-shot, planar images (15ns exposures) at 30 Hz. For this paper, the primary data set consists of 240 frames of velocity data acquired with both the jet and the low-speed co-flow seeded with light-scattering articles. Thus, velocities could be measured everywhere in the jet shear layer, both in the jet fluid and in the entrained co-flow. Some data were also taken with only the jet flow seeded. These provided mixing concentration images along with the reduced velocity fields. Other images were taken with only the co-flow seeded. These produced unique quantitative images of high speed entrainment. Optical "laser speckle" noise is the largest source of random noise in pulsed PDV systems. Components for the PDV imaging system were specifically selected to minimize speckle noise. To reduce systematic velocity errors due to laser drift, a frequency monitoring reference leg with a temperature-tuned reference iodine cell, was employed. In the course of this study, a novel flow seeder was developed. It enabled continuously variable seeding of the flow with particles of Sheared Pyrogenic Amorphous Hydrophobic Silica (SPAHS). The seeder comprised a dry fluidized bed hopper and a supersonic nozzle "pickup." Shearing action in the pickup dispersed the seed material in an exceptionally fine cloud (approximately 0.3 micron). These particles followed the flow well, did not clump or cake on screens or model surfaces, and were not susceptible to evaporation. Because of the refractory nature of the particles, SPAHS seeding should also be applicable to anticipated future testing at high temperatures.

Smith, Michael W.↗

Investigation of Flow Instabilities in the Inlet Ducts of DP-1C VTOL Aircraft

An investigation of flow instabilities in the inlet ducts of a two-engine vertical takeoff and landing aircraft DP-1C is described in this report. Recent tests revealed that the engines stall during run ups while the aircraft is operating on the ground. These pop stalls occurred at relatively low power levels, sometimes as low as 60 percent of the engine full speed. Inability to run the engines up to the full speed level is attributed to in-ground effects associated with hot gas ingestion. Such pop stalls were never experienced when the aircraft was tested on a elevated grid platform, which ensured that the aircraft was operating in out-of-the-ground-effect conditions. Based on available information on problems experienced with other vertical takeoff and landing aircraft designs, it was assumed that the engine stalls were caused by partial ingestion of hot gases streaming forward from the main exit nozzle under the aircraft inlets, which are very close to the ground. It was also suggested that the nose wheel undercarriage, located between the inlets, may generate vortices or an unstable wake causing intense mixing of hot exit gases with incoming inlet flow, which would enhance the hot gas ingestion. After running a short three-day series of tests with fully instrumented engine inlets, it is now believed the most probable reason for engine pop stalls are random ingestions of a vortex generated between the two streams moving in opposite directions: outbound hot gas stream from the main nozzle close to the ground and inbound inlet flow above. Originally, the vortex is in a horizontal plane. However, at a certain velocity ratio of these two streams, the vortex attaches either to the ground or the aircraft surface at one end and the other end is swallowed by one of the aircraft inlets. Once the vortex enters the inlet duct, a puff of hot air can be sucked through the vortex core into the engine, which causes a serious inlet flow field distortion followed by an engine stall. Once the engine stalls, the outflow from the inlet pushes the vortex away and the engine resumes normal operation. This hypothesis needs to be verified experimentally; e.g., by extensive smoke flow visualization ahead of the aircraft inlets.

Lepicovsky, Jan↗

Fuzzy self-learning control for magnetic servo system

It is known that an effective control system is the key condition for successful implementation of high-performance magnetic servo systems. Major issues to design such control systems are nonlinearity; unmodeled dynamics, such as secondary effects for copper resistance, stray fields, and saturation; and that disturbance rejection for the load effect reacts directly on the servo system without transmission elements. One typical approach to design control systems under these conditions is a special type of nonlinear feedback called gain scheduling. It accommodates linear regulators whose parameters are changed as a function of operating conditions in a preprogrammed way. In this paper, an on-line learning fuzzy control strategy is proposed. To inherit the wealth of linear control design, the relations between linear feedback and fuzzy logic controllers have been established. The exercise of engineering axioms of linear control design is thus transformed into tuning of appropriate fuzzy parameters. Furthermore, fuzzy logic control brings the domain of candidate control laws from linear into nonlinear, and brings new prospects into design of the local controllers. On the other hand, a self-learning scheme is utilized to automatically tune the fuzzy rule base. It is based on network learning infrastructure; statistical approximation to assign credit; animal learning method to update the reinforcement map with a fast learning rate; and temporal difference predictive scheme to optimize the control laws. Different from supervised and statistical unsupervised learning schemes, the proposed method learns on-line from past experience and information from the process and forms a rule base of an FLC system from randomly assigned initial control rules.

Tarn, J. H.↗

Cell-wall architecture and lignin composition of wheat developed in a microgravity environment

The microgravity environment encountered during space-flight has long been considered to affect plant growth and developmental processes, including cell wall biopolymer composition and content. As a prelude to studying how microgravity is perceived - and acted upon - by plants, it was first instructive to investigate what gross effects on plant growth and development occurred in microgravity. Thus, wheat seedlings were exposed to microgravity on board the space shuttle Discovery (STS-51) for a 10 day duration, and these specimens were compared with their counterparts grown on Earth under the same conditions (e.g. controls). First, the primary roots of the wheat that developed under both microgravity and 1 g on Earth were examined to assess the role of gravity on cellulose microfibril (CMF) organization and secondary wall thickening patterns. Using a quick freeze/deep etch technique, this revealed that the cell wall CMFs of the space-grown wheat maintained the same organization as their 1 g-grown counterparts. That is, in all instances, CMFs were randomly interwoven with each other in the outermost layers (farthest removed from the plasma membrane), and parallel to each other within the individual strata immediately adjacent to the plasma membranes. The CMF angle in the innermost stratum relative to the immediately adjacent stratum was ca 80 degrees in both the space and Earth-grown plants. Second, all plants grown in microgravity had roots that grew downwards into the agar; they did not display "wandering" and upward growth as previously reported by others. Third, the space-grown wheat also developed normal protoxylem and metaxylem vessel elements with secondary thickening patterns ranging from spiral to regular pit to reticulate thickenings. Fourthly, both the space- and Earth-grown plants were essentially of the same size and height, and their lignin analyses revealed no substantial differences in their amounts and composition regardless of the gravitational field experienced, i.e. for the purposes of this study, all plants were essentially identical. These results suggest that the microgravity environment itself at best only slightly affected either cell wall biopolymer synthesis or the deposition of CMFs, in contrast to previous assertions.

manned↗

Lagrangian statistics in weakly forced two-dimensional turbulence

Measurements of Lagrangian single-point and multiple-point statistics in a quasi-two-dimensional stratified layer system are reported. The system consists of a layer of salt water over an immiscible layer of Fluorinert and is forced electromagnetically so that mean-squared vorticity is injected at a well-defined spatial scale r i . Simultaneous cascades develop in which enstrophy flows predominately to small scales whereas energy cascades, on average, to larger scales. Lagrangian correlations and one- and two-point displacements are measured for random initial conditions and for initial positions within topological centers and saddles. Some of the behavior of these quantities can be understood in terms of the trapping characteristics of long-lived centers, the slow motion near strong saddles, and the rapid fluctuations outside of either centers or saddles. We also present statistics of Lagrangian velocity fluctuations using energy spectra in frequency space and structure functions in real space. We compare with complementary Eulerian velocity statistics. We find that simultaneous inverse energy and enstrophy ranges present in spectra are not directly echoed in real-space moments of velocity difference. Nevertheless, the spectral ranges line up well with features of moment ratios, indicating that although the moments are not exhibiting unambiguous scaling, the behavior of the probability distribution functions is changing over short ranges of length scales. Furthermore, implications for understanding weakly forced 2D turbulence with simultaneous inverse and direct cascades are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

What can solve the strong CP problem?

Three possible strategies have been advocated to solve the strong CP problem. The first is the axion, a dynamical mechanism that relaxes any initial value of the CP violating angle $\overline{θ}$ to zero. The second is the imposition of new symmetries that are believed to set $\overline{θ}$ to zero in the UV. The third is the acceptance of the fine tuning of parameters. We argue that the latter two solutions do not solve the strong CP problem. The θ term of QCD is not a parameter — it does not exist in the Hamiltonian. Rather, it is a property of the quantum state that our universe finds itself in, arising from the fact that there are CP violating states of a CP preserving Hamiltonian. It is not eliminated by imposing parity as a symmetry since the underlying theory is already parity symmetric and that does not preclude the existence of CP violating states. Moreover, since the value of θ realized in our universe is a consequence of measurement, it is inherently random and cannot be fine tuned by choice of parameters. Rather any fine tuning would require a tuning between parameters in the theory and the random outcome of measurement. Our results considerably strengthen the case for the existence of the axion and axion dark matter. The confusion around θ arises from the fact that unlike classical mechanics, the Hamiltonian and Lagrangian are not equivalent in quantum mechanics. The Hamiltonian defines the differential time evolution, whereas the Lagrangian is a solution to this evolution. Consequently, initial conditions could in principle appear in the Lagrangian but not in the Hamiltonian. This results in aspects of the initial condition such as θ misleadingly appearing in the Lagrangian as parameters. We comment on the similarity between the θ vacua and the violations of the constraint equations of classical gauge theories in quantum mechanics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Utilization Profiles of the CCSDS Unified Space Link Protocol (USLP)

The purpose of this paper is to identify the utilization profiles for interfacing the Data Protocol Sublayer using the Unified Space Link Protocols (USLP) (reference 1) with the space link coding procedures as specified in the CCSDS Coding & Synchronization Blue Books (references 2 through 5), used in both telecommand and telemetry applications. This paper describes how the USLP Protocol utilizes the coding and synchronization sublayer to support: a. Direct to Earth (DTE) telemetry links for engineering and science data b. Direct to Earth (DTE) telemetry links for very high rate science data c. Direct from Earth (DFE) command, sequencing and flight software loads d. Space to Space Links (Proximity) utilized by orbiters for data exchange to/from surface bound assets. The CCSDS has divided the functions of the Data Link Layer into two sublayers: the Data Link Protocol Sublayer (DLP-SL) and the Coding and Synchronization Sublayer (CS-SL). The Data Link Protocol Sublayer (DLP-SL) interfaces to the users, accepting the data that is to be transported, on the sending side of the link, and delivering that data on the receiving end. The Transfer Frame is the data unit that is transferred across the Data Link Protocol Sublayer and the Coding and Synchronization Sublayer boundary. The Coding and Synchronization Sublayer (CS-SL) provides the encoding, randomization, and frame synchronization functions that prepares the USLP Transfer Frame for transport across the space link. The CS-SL is divided into 2 processes: 1) The Frame Interface Processes (FIP) performs the interface functions required to prepare the data for delivery to the Coding/Decoding Process (CDP). This process includes prepending a Frame Start Marker to the provided frame, when management has designated that the frame is not to be aligned to the codeblock or when there is no block code used. 2) The Coding/Decoding Process (CDP) performs the forward error correction processes that are used to optimize the performance of the link and minimize the error rate. The CDP creates the symbol stream that is delivered to the Physical Layer. The transfer of the USLP transfer frames across different types of space links is the focus of this paper. The Protocol Data Unit (PDU) that is passed in both directions between the Data Link Protocol Sublayer (DLP-SL) and Coding and Synchronization Sublayer (CS-SL) is the transfer frame. The USLP frame structure provides flexibility that can be constrained by the functions utilized within the CS-SL that prepare the transfer frame for transit. For example, the USLP transfer frame contains a length field that enables the frame to be of variable length but CS-SL under certain conditions may constrain the frame to be fixed in length. This paper describes 5 operational modes available for use by the Data Link Layer to provide data exchange across the USLP space link. These modes are different because different operational requirements apply to vastly different types of space links and thus the communications implementation requirements differ. The environmental issues include the power or energy available, the distance between the end points of the link, the complexity of the equipment available at those end points, the atmospheric conditions and radiometric frequency selection. The CS-SL utilizes different forward error correcting codes supported by specific operational modes to configure the data for transit. This paper describes all of the operational modes in a series of data models which decompose the functionality between the Data Link Protocol Sublayer and the Coding and Synchronization sublayer. The operational modes described are: 1. Uncoded Mode: has been used for short links that contain significant available power to provide an acceptable frame error rate. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 2. Convolutional Only Mode: is currently the prime forward error correction coding used for the proximity links. The frames in this mode can be variable in length and typically use an error detection algorithm (i.e., CRC) to determine if there are errors in the received frame. 3. Variable Length Frame Aligned to Variable Length Codeblock (TC): is used for Direct from Earth links were power levels are high and the simple, least complex code i.e., the BCH code is used. This mode has been in use since the early 1970s. The BCH code is a short code and the decoder is easy to implement. 4. Fixed Length Frame Aligned to Fixed Length Codeblock (AOS/TM): was introduced when the concatenated Convolutional and Reed-Solomon Code was formulated to provide significant reduction in link data error rate and the ability to determine if there was an error in the decoded codeblock. The frame is aligned to the codeblock so that there is a one to one relationship of frame errors to codeblock errors without additional error detection coding being added. This mode requires the protocol frames to be the exact size of the message portion of the codeblock. 5. Frames Unaligned to Fixed Length Codeblocks (Currently used for very high rates and space to space links): This mode is currently used for missions that have a very high data rate that can be controlled adaptively as the environment changes and as the next generation operating mode for the proximity link. This mode from a coded data stream point of view is exactly like that described in 4. above, except that the frame need not be aligned to the codeblock. There is no requirement on frame length when using this mode. Thus when using USLP it can be used to support links that require short or long frames. There is also no mandatory requirement that frames cannot be separated by idle data reducing the tight data rate connection requirements between the data link protocol sublayer and the coding & synchronization sublayer. In conclusion, how these operational modes can be put to use in mission operational scenarios is described for Direct from Earth links (DFE), Direct to Earth links (DTE), and Proximity links.

Greenberg, E.↗

Supersonic/Hypersonic Laminar Heating Correlations for Rectangular and Impact-Induced Open and Closed Cavities

Impact and debris damage to the Space Shuttle Orbiter Thermal Protection System tiles is a random phenomenon, occurring at random locations on the vehicle surface, resulting in random geometrical shapes that are exposed to a definable range of surface flow conditions. In response to the 2003 Final Report of the Columbia Accident Investigation Board, wind tunnel aeroheating experiments approximating a wide range of possible damage scenarios covering both open and closed cavity flow conditions were systematically tested in hypersonic ground based facilities. These data were analyzed and engineering assessment tools for damage-induced fully-laminar heating were developed and exercised on orbit. These tools provide bounding approximations for the damaged-surface heating environment. This paper presents a further analysis of the baseline, zero-pressure-gradient, idealized, rectangular-geometry cavity heating data, yielding new laminar correlations for the floor-averaged heating, peak cavity endwall heating, and the downstream decay rate. Correlation parameters are derived in terms of cavity geometry and local flow conditions. Prediction Limit Uncertainty values are provided at the 95%, 99% and 99.9% levels of significance. Non-baseline conditions, including non-rectangular geometries and flows with known pressure gradients, are used to assess the range of applicability of the new correlations. All data variations fall within the 99% Prediction Limit Uncertainty bounds. Importantly, both open-flow and closed-flow cavity heating are combined into a single-curve parameterization of the heating predictions, and provide a concise mathematical model of the laminar cavity heating flow field with known uncertainty.

Everhart, Joel L.↗

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗