Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random sampling”

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 847 records · Page 47

Defect Structure of Beta NiAl Using the BFS Method for Alloys

The semiempirical BFS method for alloys is generalized by replacing experimental input with first-principles results thus allowing for the study of complex systems. In order to examine trends and behavior of a system in the vicinity of a given point of the phase diagram a search procedure based on a sampling of selected configurations is employed. This new approach is applied to the study of the beta phase of the Ni-Al system, which exists over a range of composition from 45-60 at.% Ni. This methodology results in a straightforward and economical way of reproducing and understanding the basic features of this system. At the stoichiometric composition, NiAl should exist in a perfectly ordered B2 structure. Ni-rich alloys are characterized by antisite point defects (with Ni atoms in the Al sites) with a decrease in lattice parameters. On the Al-rich side of stoichiometry there is a steep decrease in lattice parameter and density with increasing Al content. The presence of vacancies in Ni sites would explain such behavior. Recent X-ray diffraction experiments suggest a richer structure: the evidence, while strongly favoring the presence of vacancies in Ni sites, also suggests the possibility of some vacancies in Al sites in a 3:1 ratio. Moreover, local ordering of vacant sites may be preferred over a random distribution of individual point defects.

Bozzolo, Guillermo↗

CMOS VLSI Active-Pixel Sensor for Tracking

An architecture for a proposed active-pixel sensor (APS) and a design to implement the architecture in a complementary metal oxide semiconductor (CMOS) very-large-scale integrated (VLSI) circuit provide for some advanced features that are expected to be especially desirable for tracking pointlike features of stars. The architecture would also make this APS suitable for robotic- vision and general pointing and tracking applications. CMOS imagers in general are well suited for pointing and tracking because they can be configured for random access to selected pixels and to provide readout from windows of interest within their fields of view. However, until now, the architectures of CMOS imagers have not supported multiwindow operation or low-noise data collection. Moreover, smearing and motion artifacts in collected images have made prior CMOS imagers unsuitable for tracking applications. The proposed CMOS imager (see figure) would include an array of 1,024 by 1,024 pixels containing high-performance photodiode-based APS circuitry. The pixel pitch would be 9 m. The operations of the pixel circuits would be sequenced and otherwise controlled by an on-chip timing and control block, which would enable the collection of image data, during a single frame period, from either the full frame (that is, all 1,024 ~ 1,024 pixels) or from within as many as 8 different arbitrarily placed windows as large as 8 by 8 pixels each. A typical prior CMOS APS operates in a row-at-a-time ( grolling-shutter h) readout mode, which gives rise to exposure skew. In contrast, the proposed APS would operate in a sample-first/readlater mode, suppressing rolling-shutter effects. In this mode, the analog readout signals from the pixels corresponding to the windows of the interest (which windows, in the star-tracking application, would presumably contain guide stars) would be sampled rapidly by routing them through a programmable diagonal switch array to an on-chip parallel analog memory array. The diagonal-switch and memory addresses would be generated by the on-chip controller. The memory array would be large enough to hold differential signals acquired from all 8 windows during a frame period. Following the rapid sampling from all the windows, the contents of the memory array would be read out sequentially by use of a capacitive transimpedance amplifier (CTIA) at a maximum data rate of 10 MHz. This data rate is compatible with an update rate of almost 10 Hz, even in full-frame operation

Pain, Bedabrata↗

State-space self-tuner for on-line adaptive control

Dynamic systems, such as flight vehicles, satellites and space stations, operating in real environments, constantly face parameter and/or structural variations owing to nonlinear behavior of actuators, failure of sensors, changes in operating conditions, disturbances acting on the system, etc. In the past three decades, adaptive control has been shown to be effective in dealing with dynamic systems in the presence of parameter uncertainties, structural perturbations, random disturbances and environmental variations. Among the existing adaptive control methodologies, the state-space self-tuning control methods, initially proposed by us, are shown to be effective in designing advanced adaptive controllers for multivariable systems. In our approaches, we have embedded the standard Kalman state-estimation algorithm into an online parameter estimation algorithm. Thus, the advanced state-feedback controllers can be easily established for digital adaptive control of continuous-time stochastic multivariable systems. A state-space self-tuner for a general multivariable stochastic system has been developed and successfully applied to the space station for on-line adaptive control. Also, a technique for multistage design of an optimal momentum management controller for the space station has been developed and reported in. Moreover, we have successfully developed various digital redesign techniques which can convert a continuous-time controller to an equivalent digital controller. As a result, the expensive and unreliable continuous-time controller can be implemented using low-cost and high performance microprocessors. Recently, we have developed a new hybrid state-space self tuner using a new dual-rate sampling scheme for on-line adaptive control of continuous-time uncertain systems.

Shieh, L. S.↗

Analysis of Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed ~ 1 ppm random errors for 8-10 km altitudes and ~30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2lidar on the NASA DC-8 and added an 02lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected ~linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. We demonstrated measurements over the California Central Valley, to stratus cloud tops over the Pacific Ocean, over mountain regions with snow, and over several areas with broken clouds. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗

Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption and Range During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar only on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear C02 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed in detail and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed -1 ppm random errors for 8-10 km altitudes and -30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2 lidar on the NASA DC-8 and added an O2 lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected -linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. The seven flights in the 2011 Ascends campaign were flown over a wide variety of surface and cloud conditions in the US, which produced a wide variety of lidar signal conditions. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗

Assessing the Effects of Various Surface Textures and Features on Turbulent Heat Transfer in Hypersonic Flight

Experiments were conducted in the NASA Ames Hypervelocity Free Flight Aerodynamic Facility (ballistic range) to quantify the effects on turbulent convective heat transfer of distributed surface roughness, and of isolated features, representative of thermal protection systems on atmospheric entry vehicles. The surface textures and features were applied on the conic frusta of 45o sphere-cone models having a nose-to-base radius ratio of 0.5, similar to the forebody geometry of the Galileo and Pioneer-Venus entry probes. Test conditions were selected to provide turbulent roughness Reynolds numbers, k+, in the ranges expected for outer planet entry missions. Turbulent flow on the conic frustum was achieved by tripping the flow on the sphere-segment nose cap with distributed surface roughness, created by sand-blasting the nose cap. Surface textures included distributed, acreage, roughness, as well as cavity and groove discrete features. The majority of the data to be presented are results for distributed roughness, which includes both random, sand-grain-like roughness, and regular pattern roughness. The pattern roughness was designed to represent the roughness on woven thermal protection system materials, such as NASA’s 3-D Medium Density Carbon Phenolic (3MDCP), also known as HEEET, developed by the Heatshield for Extreme Entry Environments Technology project. The pattern roughness tested is a 3-D wavy surface, and includes three different roughness element height-to-spacing ratios representative of two configurations of 3MDCP, and spanning ratios measured both before, and after, ablation in an arc jet test facility. The patterns tested in the ballistic range were laser-etched on metal models, and represented an idealized version of the real-world materials in which each roughness element was nearly identical. Additional tests were performed wherein the laser-etched patterns were degraded by sand-blasting with various sized grit media, to produce regular patterns with superimposed irregular roughness, more representative of flight materials. Results of each will be compared. The discrete features tested included cylindrical cavities and rectangular grooves of various width-to-depth ratios. Cavities represent either heatshield damage, such as from micro-meteoroid and orbital debris (MMOD) damage, or from designed penetrations, such as on the Genesis sample return capsule. The grooves were scaled representations of seams between segments of HEEET material in a notional tiled thermal protection system. The tests examined the effects on turbulent heating downstream of the isolated features. The tests were conducted at speeds between 3 km/s and 4 km/s in air between 0.15 atm and 0.25 atm (Mach numbers between 9 and 12). Roughness Reynolds numbers, k+, ranged from 12 to 70 for the sand roughness, and as high as 200 for the pattern roughness. Boundary-layer parameters required for calculating k+ were evaluated using computational fluid dynamics simulations using the DPLR (Data Parallel Line Relaxation) code. Each model included both rough- and smooth-wall segments, and heat transfer augmentation factors were determined as the ratio of the rough-wall to smooth-wall heat flux measured on each test.

hypersonic↗

Design and testing of a unique randomized gravity, continuous flow bioreactor

A rotating, null gravity simulator, or Couette bioreactor was successfully used for the culture of mammalian cells in a simulated microgravity environment. Two limited studies using Lipomyces starkeyi and Streptomyces clavuligerus were also conducted under conditions of simulated weightlessness. Although these studies with microorganisms showed promising preliminary results, oxygen limitations presented significant limitations in studying the biochemical and cultural characteristics of these cell types. Microbial cell systems such as bacteria and yeast promise significant potential as investigative models to study the effects of microgravity on membrane transport, as well as substrate induction of inactive enzyme systems. Additionally, the smaller size of the microorganisms should further reduce the gravity induced oscillatory particle motion and thereby improve the microgravity simulation on earth. Focus is on the unique conceptual design, and subsequent development of a rotating bioreactor that is compatible with the culture and investigation of microgravity effects on microbial systems. The new reactor design will allow testing of highly aerobic cell types under simulated microgravity conditions. The described reactor affords a mechanism for investigating the long term effects of reduced gravity on cellular respiration, membrane transfer, ion exchange, and substrate conversions. It offers the capability of dynamically altering nutrients, oxygenation, pH, carbon dioxide, and substrate concentration without disturbing the microgravity simulation, or Couette flow, of the reactor. All progeny of the original cell inoculum may be acclimated to the simulated microgravity in the absence of a substrate or nutrient. The reactor has the promise of allowing scientists to probe the long term effects of weightlessness on cell interactions in plants, bacteria, yeast, and fungi. The reactor is designed to have a flow field growth chamber with uniform shear stress, yet transfer high concentrations of oxygen into the culture medium. The system described allows for continuous, on line sampling for production of product without disturbing fluid and particle dynamics in the reaction chamber. It provides for the introduction of substrate, or control substances after cell adaptation to simulated microgravity has been accomplished. The reactor system provides for the nondisruptive, continuous flow replacement of nutrient and removal of product. On line monitoring and control of growth conditions such as pH and nutrient status are provided. A rotating distribution valve allows cessation of growth chamber rotation, thereby preserving the simulated microgravity conditions over longer periods of time.

Lassiter, Carroll B.↗

Bias properties of extragalactic distance indicators. 3: Analysis of Tully-Fisher distances for the Mathewson-Ford-Buchhorn sample of 1355 galaxies

The observational selection bias properties of the large Mathewson-Ford-Buchhorn (MFB) sample of axies are demonstrated by showing that the apparent Hubble constant incorrectly increases outward when determined using Tully-Fisher (TF) photometric distances that are uncorreted for bias. It is further shown that the value of H(sub 0) so determined is also multivlaued at a given redshift when it is calculated by the TF method using galaxies with differenct line widths. The method of removing this unphysical contradiction is developed following the model of the bias set out in Paper II. The model developed further here shows that the appropriate TF magnitude of a galaxy that is drawn from a flux-limited catalog not only is a function of line width but, even in the most idealistic cases, requires a triple-entry correction depending on line width, apparent magnitude, and catalog limit. Using the distance-limited subset of the data, it is shown that the mean intrinsic dispersion of a bias-free TF relation is high. The dispersion depends on line width, decreasing from sigma(M) = 0.7 mag for galaxies with rotational velocities less than 100 km s(exp-1) to sigma(M) = 0.4 mag for galaxies with rotational velocities greater than 250 km s(exp-1). These dispersions are so large that the random errors of the bias-free TF distances are too gross to detect any peculiar motions of individual galaxies, but taken together the data show again the offset of 500 km s(exp-1) fond both by Dressler & Faber and by MFB for galaxies in the direction of the putative Great Attractor but described now in a different way. The maximum amplitude of the bulk streaming motion at the Local Group is approximately 500 km s(exp-1) but the perturbation dies out, approaching the Machian frame defined by the CMB at a distance of approximately 80 Mpc (v is approximately 4000 km s(exp -1)). This decay to zero perturbation at v is approximately 4000 km s(exp -1) argues against existing models with a single attraction at approximately 4500 km s(exp -1) (the Great Attactor model) pulling the local region. Rather, the cause of the perturbation appears to be the well-known clumpy mass distribution within 4000 km s(exp -1) in the busy directions of Hydra, Centaurus, Antila and Dorado, as postulated earlier (Tammann & Sandage 1985).

Federspiel, Martin↗

Part 1: Aspects of lithospheric evolution on Venus. Part 2: Thermal and collisional histories of chondrite parent bodies

The geological evolution of distinctly different kinds of solar system objects is addressed. Venus has been observed over the past decade by orbital radars on both American and Soviet spacecraft. These surface measurements provide clues to the structure and evolution of the lithosphere. The parent bodies of chondritic meteorites, thought to resemble asteroids, represent the other end of the size spectrum of terrestrial objects. Their early thermal and collisional histories may be constrained by the chemical and textural record preserved in meteorite samples. Impact craters on Venus have been observed by the Soviet Venera 15/16 spacecraft. A formalism is presented by which the size-frequency distribution of impact craters may be used to estimate upper bounds on the mean global rates of volcanic resurfacing and lithospheric recycling on that planet over the past several hundred million years. The impact crater density reported from Venera observations, if valid for the entire Venus surface, indicates a mean volcanic flux no greater than 2 cu km/y, corresponding to a maximum average rate of resurfacing of about 4 km/b.y. For the lowest estimated mean crater retention age of the surface of Venus imaged by Venera 15/16, the rate of lithospheric recycling on Venus does not exceed 1.5 sq km/y. Ordinary chondrite meteorites show textural and chemical patterns indicative of varying intensities of thermal metamorphism. The conventional onion-shell model, which envisions highly metamorphosed material in the core and less intensely heated rocks near the surface, predicts an inverse relation between peak temperature and cooking rate, but none has been observed. A metamorphosed-planetesimal model is devised to explain this discrepancy, whereby heating occurs in planetesimals a few kilometers in radius which then accrete to form 100-km-radius parent bodies. Cooling rates are then randomly controlled by burial depth. Thermal and collisional constraints are examined, and the model is found to be applicable only to highly insulating Al-26-rich planetesimals that remain closely aggregated upon accretion. An alternative model is presented, in which onion-shell parent bodies are collisionally fragmented during metamorphism and then gravitationally reassembled. If reassembly times are short, then cooling rates would be determined by burial depth in the reaccreted parent body. This model, unlike previous ones, can explain both coherent and incoherent cooling of Breccia clasts by collisions during or after metamorphism, respectively.

Grimm, Robert E.↗

Uncovering Hazards Using a Multi-Objective Optimization to Explore the Faulty State-Space

Considering resilience when designing complex engineered systems is crucial to ensure the system is safe under unexpected hazardous scenarios. Traditional risk-based approaches, such as Failure Modes and Effects Analysis (FMEA) are useful for designing the system to mitigate hazardous scenarios that can be identified by the designer, but often require experience or prior knowledge of system failures to generate. More recently, researchers have developed simulation tools that enable the designer to model large sets of hazardous scenarios (driven by both internal faults and external factors) through simulation. While these tools enable a wider scope of fault modes to be evaluated (e.g., by injecting combined set of fault modes or injecting modes at different times), the resulting assessments (like FMEA) still require knowledge of the specific modes to be evaluated. However, failure to analyze a wide variety of fault scenarios can lead to an incomplete picture of the system resilience, especially to "surprise events'' which may be difficult for the designer to identify and predict beforehand. To overcome this challenge, previous work developed a fault sampling approach for resilience simulations which would procedurally-generate a wide variety of potential faults by systematically perturbing the health states of the system. While the resulting fault modes generated covered a much larger space hazards than would be otherwise considered (and identified many unique failure trajectories which would not have otherwise been identified), it also significantly increased the computational cost of the analysis and resulted in the simulation and analysis of a large set of essentially duplicate scenarios. Additionally, as the number of dimensions in the faulty state-space increases, the full elaboration of possible modes becomes computationally infeasible, justifying the use of a more targeted search. To resolve this limitation, this work proposes the use of a multiobjective optimization algorithm to search the health state space for potential fault modes that are both (1) hazardous and (2) unique. To solve this type of problem, this work proposes the use of a cooperative co-evolutionary algorithm. To demonstrate this approach, it will be applied to a model of an autonomous rover which uses line markings to navigate, focusing on potential hazards in the drive system which could cause the rover to crash. To determine the merit of the approach, it will further be compared with the previously-presented range elaboration approach and a random mode generation approach on the basis of computational efficiency and found modes.

Resilience↗

Hypercortisolemia alters muscle protein anabolism following ingestion of essential amino acids

Debilitating injury is accompanied by hypercortisolemia, muscle wasting, and disruption of the normal anabolic response to food. We sought to determine whether acute hypercortisolemia alters muscle protein metabolism following ingestion of a potent anabolic stimulus: essential amino acids (EAA). A 27-h infusion (80 microg. kg(-1). h(-1)) of hydrocortisone sodium succinate mimicked cortisol (C) levels accompanying severe injury (>30 microg/dl), (C + AA; n = 6). The control group (AA) received intravenous saline (n = 6). Femoral arteriovenous blood samples and muscle biopsies were obtained during a primed (2.0 micromol/kg) constant infusion (0.05 micromol. kg(-1). min(-1)) of l-[ring-(2)H(5)]phenylalanine before and after ingestion of 15 g of EAA. Hypercortisolemia [36.5 +/- 2.1 (C + AA) vs. 9.0 +/- 1.0 microg/dl (AA)] increased postabsorptive arterial, venous, and muscle intracellular phenylalanine concentrations. Hypercortisolemia also increased postabsorptive and post-EAA insulin concentrations. Net protein balance was blunted (40% lower) following EAA ingestion but remained positive for a greater period of time (60 vs. 180 min) in the C + AA group. Thus, although differences in protein metabolism were evident, EAA ingestion improved muscle protein anabolism during acute hypercortisolemia and may help minimize muscle loss following debilitating injury.

Clinical Trial↗

Epstein-Barr virus shedding by astronauts during space flight

Patterns of Epstein-Barr virus (EBV) reactivation in 32 astronauts and 18 healthy age-matched control subjects were characterized by quantifying EBV shedding. Saliva samples were collected from astronauts before, during, and after 10 space shuttle missions of 5-14 days duration. At one time point or another, EBV was detected in saliva from each of the astronauts. Of 1398 saliva specimens from 32 astronauts, polymerase chain reaction analysis showed that 314 (23%) were positive for EBV DNA. Examination by flight phase showed that 29% of the saliva specimens collected from 28 astronauts before flight were positive for EBV DNA, as were 16% of those collected from 25 astronauts during flight and 16% of those collected after flight from 23 astronauts. The mean number of EBV copies from samples taken during the flights was 417 per mL, significantly greater (p<.05) than the number of viral copies from the preflight (40) and postflight (44) phases. In contrast, the control subjects shed EBV DNA with a frequency of 3.7% and mean number of EBV copies of 40 per mL of saliva. Ten days before flight and on landing day, titers of antibody to EBV viral capsid antigen were significantly (p<.05) greater than baseline levels. On landing day, urinary levels of cortisol and catecholamines were greater than their preflight values. In a limited study (n=5), plasma levels of substance P and other neuropeptides were also greater on landing day. Increases in the number of viral copies and in the amount of EBV-specific antibody were consistent with EBV reactivation before, during, and after space flight.

Randomized Controlled Trial↗

A model for interphase chromosomes and evaluation of radiation-induced aberrations

We have developed a theoretical model for evaluating radiation-induced chromosomal exchanges by explicitly taking into account interphase (G(0)/G(1)) chromosome structure, nuclear organization of chromosomes, the production of double-strand breaks (DSBs), and the subsequent rejoinings in a faithful or unfaithful manner. Each of the 46 chromosomes for human lymphocytes (40 chromosomes for mouse lymphocytes) is modeled as a random polymer inside a spherical volume. The chromosome spheres are packed randomly inside a spherical nucleus with an allowed overlap controlled by a parameter Omega. The rejoining of DSBs is determined by a Monte Carlo procedure using a Gaussian proximity function with an interaction range parameter sigma. Values of Omega and sigma have been found which yield calculated results of interchromosomal aberration frequencies that agree with a wide range of experimental data. Our preferred solution is one with an interaction range of 0.5 microm coupled with a relatively small overlap parameter of 0.675 microm, which more or less confirms previous estimates. We have used our model with these parameter values and with resolution or detectability limits to calculate yields of translocations and dicentrics for human lymphocytes exposed to low-LET radiation that agree with experiments in the dose range 0.09 to 4 Gy. Five different experimental data sets have been compared with the theoretical results. Essentially all of the experimental data fall between theoretical curves corresponding to resolution limits of 1 Mbp and 20 Mbp, which may reflect the fact that different investigators use different limits for sensitivity or detectability. Translocation yields for mouse lymphocytes have also been calculated and are in good agreement with experimental data from 1 cGy to 10 cGy. There is also good agreement with recent data on complex aberrations. Our model is expected to be applicable to both low- and high-LET radiation, and we include a sample prediction of the yield of interchromosomal rejoining in the dose range 0.22 Gy to 2 Gy of 1000 MeV/nucleon iron particles. This dose range corresponds to average particle traversals per nucleus ranging from 1.0 to 9.12.

Non-NASA Center↗

Active Suppression of Instabilities in Engine Combustors

A method of feedback control has been proposed as a means of suppressing thermo-acoustic instabilities in a liquid- fueled combustor of a type used in an aircraft engine. The basic principle of the method is one of (1) sensing combustor pressure oscillations associated with instabilities and (2) modulating the rate of flow of fuel to the combustor with a control phase that is chosen adaptively so that the pressure oscillations caused by the modulation oppose the sensed pressure oscillations. The need for this method arises because of the planned introduction of advanced, lean-burning aircraft gas turbine engines, which promise to operate with higher efficiencies and to emit smaller quantities of nitrogen oxides, relative to those of present aircraft engines. Unfortunately, the advanced engines are more susceptible to thermoacoustic instabilities. These instabilities are hard to control because they include large dead-time phase shifts, wide-band noise characterized by amplitudes that are large relative to those of the instabilities, exponential growth of the instabilities, random net phase walks, and amplitude fluctuations. In this method (see figure), the output of a combustion-pressure sensor would be wide-band-pass filtered and then further processed to generate a control signal that would be applied to a fast-actuation valve to modulate the flow of fuel. Initially, the controller would rapidly take large phase steps in order to home in, within a fraction of a second, to a favorable phase region within which the instability would be reduced. Then the controller would restrict itself to operate within this phase region and would further restrict itself to operate within a region of stability, as long as the power in the instability signal was decreasing. In the phase-shifting scheme of this method, the phase of the control vector would be made to continuously bounce back and forth from one boundary of an effective stability region to the other. Computationally, this scheme would be implemented by the adaptive sliding phaser averaged control (ASPAC) algorithm, which requires very little detailed knowledge of the combustor dynamics. In the ASPAC algorithm, the power of the instability signal would be calculated from the wide-bandpass- filtered combustion-pressure signal and averaged over a period of time (typically of the order of a few hundredths of a second) corresponding to the controller updating cycle [not to be confused with the controller sampling cycle, which would be much shorter (typically of the order of 10(exp -4) second)].

Kopasakis, George↗

Large-Eddy Simulation Code Developed for Propulsion Applications

A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of Smagorinsky s model for the subgrid-scale turbulence, and it uses generalized curvilinear coordinates to allow analysis of a wide range of geometries. The code runs in parallel on shared memory multiprocessor computers and is written in Fortran 90 with dynamic memory allocation. A sample result for a Mach-1.4 round jet is presented in the figure. Instantaneous Mach number contours in several cross-planes downstream of the nozzle exit are shown, illustrating how an LES captures the large unsteady three-dimensional turbulent structures present in the jet.

DeBonis, James R.↗

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗

Portable Electromyograph

A portable electronic apparatus records electromyographic (EMG) signals in as many as 16 channels at a sampling rate of 1,024 Hz in each channel. The apparatus (see figure) includes 16 differential EMG electrodes (each electrode corresponding to one channel) with cables and attachment hardware, reference electrodes, an input/output-and-power-adapter unit, a 16-bit analog-to-digital converter, and a hand-held computer that contains a removable 256-MB flash memory card. When all 16 EMG electrodes are in use, full-bandwidth data can be recorded in each channel for as long as 8 hours. The apparatus is powered by a battery and is small enough that it can be carried in a waist pouch. The computer is equipped with a small screen that can be used to display the incoming signals on each channel. Amplitude and time adjustments of this display can be made easily by use of touch buttons on the screen. The user can also set up a data-acquisition schedule to conform to experimental protocols or to manage battery energy and memory efficiently. Once the EMG data have been recorded, the flash memory card is removed from the EMG apparatus and placed in a flash-memory- card-reading external drive unit connected to a personal computer (PC). The PC can then read the data recorded in the 16 channels. Preferably, before further analysis, the data should be stored in the hard drive of the PC. The data files are opened and viewed on the PC by use of special- purpose software. The software for operation of the apparatus resides in a random-access memory (RAM), with backup power supplied by a small internal lithium cell. A backup copy of this software resides on the flash memory card. In the event of loss of both main and backup battery power and consequent loss of this software, the backup copy can be used to restore the RAM copy after power has been restored. Accessories for this device are also available. These include goniometers, accelerometers, foot switches, and force gauges.

De Luca, Gianluca↗

Evaluating Meteorological Dust Events and Machine-Learning Based Dust Identification in Geostationary Satellite Imagery

NASA scientists in the Short-term Prediction Research and Transition Center (SPoRT) developed a physically-based machine learning approach to identify dust in satellite imagery with a focus on night-time dust detection (Berndt et al. 201; DustTracker-AI). NASA/NOAA Geostationary Environmental Operational Satellite-16 (GOES-16) imagery was used for training and model inputs. The training, testing and validation data set consists of 28 events in the Southwest United States, capturing dust and null events in the region from 2018-2020.With 83 distinct images and millions of pixels a random forest model was trained and validated, correctly labeling 85% of dust pixels.For the first time, the model was run in near-real time production during the spring of 2022 and dust probability visualizations were made available to NOAA National Weather Service (NWS) forecasters to assess its utility for dust forecasting. Results indicated the model helped increase the confidence in the presence of dust and enabled dust tracking for a longer period of time into the night-time hours. Forecaster assessment and running the model in near real-time allowed for the team to determine the types of events missed, captured, and false alarms. To gain additional context on model performance,the SPoRT team sought to gather more detailed information on the training database(e.g., meteorological characteristics and drivers). The goal of this project was to identify the meteorological drivers for the dust events and create a database which synthesized information from observations, forecaster discussions, and analyses pertaining to the dust events to understand the types of events currently used to train the model. A more detailed meteorological synopsis was created for each dust event in the training, testing, and validation datasets. Following the completion of the database and documentation, the classification details revealed that 88% of the dust events were synoptically driven while mesoscale events were less prevalent in model datasets. Meteorological conditions found such as mixing layer depth and wind velocity had mean values of 645mb and 21kt respectively.With conditions of deep mixed layers and moderate to strong surface winds a mesoscale thunderstorm outflow event was considered and subsequently added to the model training data set to test the impact of additional mesoscale training data. The model was retrained and then qualitatively tested on a sample thunderstorm outflow case that the original model was unable to identify. Preliminary results showed potential that the addition of more mesoscale events included in the training data could help to better identify indistinct and localized dust events.

Connor Welch↗