Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Thoughts Without Content are Empty, Intuitions Without Concepts are Blind - Determinism and Contingency Revisited

Was the emergence of life a predictable outcome of chemical evolution on earth? Could evolution produce life very different from ours? These are one of the oldest questions in the field of the origin of life that not only have broad philosophical implications but also impact how we approach the problem from the methodological standpoint. Framing the issue in terms of the dichotomy between contingency and determinism is not a fortunate because these two terms in their conventional meaning are neither mutually exclusive nor jointly exhaustive. Determinism, represented in natural sciences by Newtonian physics, relies on the assumption that every event is causally determined by a chain of previous events. In the context of the origin of life it means that once the initial conditions on the early earth have been specified further evolution follows inevitably. Considering uncertainties about conditions on the prebiotic earth, many plausible sets of initial conditions can be defined, each followed by a separate deterministic trajectory. This conventional understanding of determinism does not admit contingency. Further, it has no implications for evaluating how many sets of initial conditions lead to the emergence of life. It appears that a better framing of the problem is as follows: given plausible sets of initial conditions on the early earth how probable and broadly spread are evolutionary trajectories that lead to life? Instead of undertaking an impossible task of specifying microscopic initial conditions for all components of the system one uses a reduced representation of this system and specify only a small set of essential macroscopic parameters, values (or ranges of values) of which can be identified, inferred or estimated from experiment, theory or historical record. The following evolutionary trajectories are still governed by laws of physics and chemistry but become probabilistic and "contingency" is admitted as variations in other variables in the system. A similar reasoning is common in other fields of science, for example in statistical mechanics. Some trajectories lead to life, perhaps in different forms, whereas others do not. Of our true interest is the ratio of these two outcomes. The issue of determinism does not directly enter the picture. The debate about the likelihood of the emergence of life is quite old. One view holds that the origin of life is an event governed by chance, and the result of so many random events (contingencies) is unpredictable. This view was eloquently expressed by Monod. In his book "Chance or Necessity" he argued that life was a product of "nature's roulette." In an alternative view, expressed in particular by deDuve and Morowitz, the origin of life is considered a highly probable or even inevitable event (although its details need not be determined in every respect). Only in this sense the origin of life can be considered a "deterministic event".

Pohorille, Andrew↗
DOE OSTI2016

Clinical Experience and Evaluation of Patient Treatment Verification With a Transit Dosimeter

Purpose: To prospectively evaluate a protocol for transit dosimetry on a patient population undergoing intensity modulated radiation therapy (IMRT) and to assess the issues in clinical implementation of electronic portal imaging devices (EPIDs) for treatment verification. Methods and Materials: Fifty-eight patients were enrolled in the study. Amorphous silicon EPIDs were calibrated for dose and used to acquire images of delivered fields. Measured EPID dose maps were back-projected using the planning computed tomographic (CT) images to calculate dose at prespecified points within the patient and compared with treatment planning system dose offline using point dose difference and point γ analysis. The deviation of the results was used to inform future action levels. Results: Two hundred twenty-five transit images were analyzed, composed of breast, prostate, and head and neck IMRT fields. Patient measurements demonstrated the potential of the dose verification protocol to model dose well under complex conditions: 83.8% of all delivered beams achieved the initial set tolerance level of Δ{sub D} of 0 ± 5 cGy or %Δ{sub D} of 0% ± 5%. Importantly, the protocol was also sensitive to anatomic changes and spotted that 3 patients from 20 measured prostate patients had undergone anatomic change in comparison with the planning CT. Patient data suggested an EPID-reconstructed versus treatment planning system dose difference action level of 0% ± 7% for breast fields. Asymmetric action levels were more appropriate for inversed IMRT fields, using absolute dose difference (−2 ± 5 cGy) or summed field percentage dose difference (−6% ± 7%). Conclusions: The in vivo dose verification method was easy to use and simple to implement, and it could detect patient anatomic changes that impacted dose delivery. The system required no extra dose to the patient or treatment time delay and so could be used throughout the course of treatment to identify and limit systematic and random errors in dose delivery for patient groups.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2010

A Model for Backscattering from Quasi Periodic Corn Canopies at L-Band

In this study, a model for backscattering at L-band from a corn canopy is proposed. The canopy consists of a quasi-periodic distribution of stalks and a random distribution of leaves. The Distorted Born Approximation (DBA) is employed to calculate the single scattered return from the corn field. The new feature of the method is that the coherence of the stalks in the row direction is incorporated in the model in a systematic fashion. Since the wavelength is on the order of the distance between corn stalks in a row, grating lobe behavior is observed at certain azimuth angles of incidence. The results are compared with experimental values measured in Huntsville, Alabama in 1998. The mean field and the effective dielectric constant of the canopy are obtained by using the Foldy approximation. The stalks are placed in the effective medium in a two dimensional lattice to simulate the row structure of a corn field. In order to mimic a real corn field, a quasi-periodic stalk distribution is assumed where the stalks are given small random perturbations about their lattice locations. Corn leaves are also embedded in the effective medium and the backscattered field from the stalks and the leaves is computed. The backscattering coefficient is calculated and averaged over successive stalk position perturbations. It is assumed that soil erosion has smoothed the soil sufficiently so that it can be assumed flat. Corn field backscatter data was collected from cornfields during the Huntsville 98 experimental campaign held at Alabama A&M University Research Station, Huntsville, Alabama in 1998 using the NASA/GW truck mounted radar. Extensive ground truth data was collected. This included soil moisture measurements and corn plant architectural data to be used in the model. In particular, the distances between the stalks in a single row have been measured. The L-band radar backscatter data was collected for both H and V polarizations and for look angles of 15o and 45o over a two week period under varying soil moisture conditions. These measured backscattering values will be compared with the model backscattering values and a discussion of the results will be presented.

Lang, R.↗
NASA NTRS2018

Characterization Testing of Lockheed Martin Micro1-2 Cryocoolers Optimized for 220 K Environment

The Mapping Imaging Spectrometer for Europa (MISE) Instrument on the Europa Clipper mission has baselined a Lockheed Martin high-power Micro1-2 pulse tube cryocooler operating at 135 Hz with a 220 K heat rejection temperature. This paper describes the testing and results of two Lockheed Martin Micro1-2 coolers optimized for these conditions. The thermal performance of the microcoolers was measured in vacuum for heat reject temperatures between 220 and 230 K. The coolers were operated with input powers ranging from 5 W to 40 W and drive frequency between 125 Hz and 145 Hz. The optimal drive frequency was strongly dependent on heat reject temperature. In addition, the exported forces and torques of the coolers were measured at 300 K heat rejection for input powers ranging from 5 W to 60 W and drive frequency between 120 Hz and 160 Hz. The exported forces were dependent on both piston amplitude and drive frequency. Moreover, the following were measured on one of the coolers: DC and AC magnetic fields at various locations, the off-state conductance of the pulse tube, and the effect of inclination angle of the pulse tube relative to gravity on the performance of the cooler. Finally, one of the coolers was subjected to random vibration testing.

Rodriguez, J.I.↗
NASA NTRS2014

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using MBH have become well developed during the past decade. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing RVs from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by Englander significantly improves MBH performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness, where efficiency is finding better solutions in less time, and robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive RWs originally developed in the field of statistical physics.

mission design↗
NASA NTRS2014

An Investigation of Low Earth Orbit Internal Charging

Low Earth orbit is usually considered a relatively benign environment for internal charging threats due to the low flux of penetrating electrons with energies of a few MeV that are encountered over an orbit. There are configurations, however, where insulators and ungrounded conductors used on the outside of a spacecraft hull may charge when exposed to much lower energy electrons of some 100's keV in a process that is better characterized as internal charging than surface charging. For example, the minimal radiation shielding afforded by thin thermal control materials such as metalized polymer sheets (e.g., aluminized Kapton or Mylar) and multilayer insulation may allow electrons of 100's of keV to charge underlying materials. Yet these same thermal control materials protect the underlying insulators and ungrounded conductors from surface charging currents due to electrons and ions at energies less than a few keV as well as suppress the photoemission, secondary electron, and backscattered electron processes associated with surface charging. We investigate the conditions required for this low Earth orbit "internal charging" to occur and evaluate the environments for which the process may be a threat to spacecraft. First, we describe a simple one-dimensional internal charging model that is used to compute the charge accumulation on materials under thin shielding. Only the electron flux that penetrates exposed surface shielding material is considered and we treat the charge balance in underlying insulation as a parallel plate capacitor accumulating charge from the penetrating electron flux and losing charge due to conduction to a ground plane. Charge dissipation due to conduction can be neglected to consider the effects of charging an ungrounded conductor. In both cases, the potential and electric field is computed as a function of time. An additional charge loss process is introduced due to an electrostatic discharge current when the electric field reaches a prescribed breakdown strength. For simplicity, the amount of charge lost in the discharge is treated as a random percentage of the total charge between a set maximum and minimum amount so a user can consider partial discharges of insulating materials (small loss of charge) or arcing from a conductor (large loss of charge). We apply the model to electron flux measurements from the NOAA-19 spacecraft to demonstrate that charging can reach levels where electrostatic discharges occur and estimate the magnitude of the discharge.

NeergaardParker, Linda↗
NASA NTRS2003

Sheet Flows, Avalanches, and Dune Evolution on Earth and Mars

This investigation is a collaboration between researchers at Cornell University, the University of Florida, and the University of Rennes 1, France. Flow modeling at Cornell University focused on mechanisms for the suspension and transport of wind-blown sand that are important in both terrestrial and Martian environments. These mechanisms include the saltation (or jumping) of grains, collisions between grains, and the interaction of grains with the velocity fluctuations of the turbulent wind. Of particular interest are sheet flows; these are relatively thin, highly concentrated regions of grains flowing near the ground under the influence of a strong turbulent wind. In them, the grains are suspended by interparticle collisions. Sheet flows may be relatively rare events, but they have the capacity to move great amounts of sand. In order to describe sheet flows, a turbulent mixture theory was formulated for particles in a fluid in which fluctuations in the volume fiaction of the particles take place on the scale of the turbulent eddies. Ensemble averaged equations for particle and fluid mass, momentum, and energy and fluid rate of dissipation were expressed in terms of Farve (concentration) averaged velocities and the associated velocity fluctuations. Correlations that describe the turbulent suspension of particles and dissipation of turbulent energy of both phases due to fluid particle interactions were modeled and boundary conditions at the bed and at the upper surface of the collisional flow were formulated. The boundary conditions at the upper surface were tested in a numerical simulation developed at the University of Florida. Steady and unsteady solutions for steady and unsteady fully-developed flows were obtained over a range of wind speeds fiom the lowest for which collisional between particles occurred to at which turbulent suspension is found to dominate collisional suspension. Below the value of the wind speed at which collisions between particles were unimportant, numerical solutions were obtained for the velocity distribution function and the resulting fields of concentration, particle and gas mean velocity, and particle shear stress for the steady two-dimensional saltation of spherical sand particles driven by a turbulent wind over a bed characterized by a simple relationship (the splash function) between the properties of incoming particles and those of the rebounding particles and other particles ejected fiom the bed. At the University of Rennes 1, experiments devoted to the characterization of the splash function for beds consisting of either random or ordered arrays of spheres in two- dimensions were completed. These indicated the role played by the packing geometry in the rebound and ejection of grains. Preliminary experiments on response of a three- dimensional collision bed to a collision with a single particle were performed. Data was taken with a single camera focused on the plane of collision. Here, for example, the decrease of the effective coefficient of restitution of the bed with an increase of the angle of incidence of the incoming particle has been measured. Other experiments on avalanches at Rennes studied the properties of the flows of particles that are responsible for the motion of the leeward side of a dune. In these, the dependence of the initiation of avalanches on the packing and depth of the particles was measured. Particle migration was studied in inclined flows of a binary mixture of disks and the mechanisms of diffision and segregation were isolated and characterized. The influence of side wall on dense, rapid inclined flows was measured and shown to be the reason why the angle of the free surface in such flows can exceed the static angle of repose. Future research will be devoted to a better understanding the transition between saltating (collisionless) and collisional flows as the wind speed the increases. This will involve the understanding of the evolution of the splash function as clisions with the bed become more numerous, more frequent, and more violent.

Source record↗
NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗
NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗
DOE OSTI2026

Measurement of Local Impedance Characteristics of In-Plane PEM Electrolyzer Component Features via Segmented Cell

Segmented cells are a class of advanced diagnostic devices that enable the measurement of current distribution within the active area of an electrochemical cell, allowing for the evaluation of localized impacts from operating conditions, flow fields and component inhomogeneities. Segmented cell devices have been applied to support R&D efforts for the advancement of polymer electrolyte membrane fuel cell technologies and more recently proton exchange membrane water electrolysis (PEMWE). Beyond current density distributions, segmented cells have been successfully deployed to collect the distribution of other relevant parameters such as high frequency resistance and electrode potentials. Very recently, a promising impedance method has been reported for water electrolysis where a one-dimensional segmented cell is coupled to a multichannel potentiostat, to enable the measurement of local cell impedances, further expanding the diagnostic capabilities of the segmented cell. In this study, we will present recent efforts to accurately measure local impedance characteristics for a PEMWE cell using a two-dimensional segmented cell device coupled with a multichannel potentiostat. The distribution of local resistances, i.e. high frequency, charge transfer, mass transport and catalyst layer, and iR-free overpotentials were successfully measured after eliminating sources of random and systematic error. We will discuss the application of this diagnostic on a cell containing a 1 cm feature simulating an ionomer skin. The results from these tests give insights if such an irregularity needs to be classified as a defect, and highlight the usefulness of coupling impedance spectroscopy with segmented cells for holistic spatial diagnostics of electrochemical devices.

08 HYDROGEN↗
DOE OSTI2019

1D PIC-DSMC analysis of a high-pressure nanosecond pulse discharge breakdown in helium

Kinetic simulations of plasma phenomena during and after formation of the conductive plasma channel of a nanosecond pulse discharge are analyzed and compared to existing experimental measurements. Particle-in-cell with direct simulation Monte Carlo collisions (PIC-DSMC) modeling is used to analyze a discharge in helium at 200 Torr and 300 K over a 1 cm gap. The analysis focuses on physics that would not be reproduced by fluid models commonly used at this high number density and collisionality, specifically non-local and stochastic phenomena. Similar analysis could be used to improve the predictive capability of lower fidelity or reduced order models. First, the modeling results compare favorably with experimental measurements of electron number density, temperature, and 1D electron energy distribution function at the same conditions. Second, it is shown that the ionization wave propagates in a stochastic, stepwise manner, dependent on rare, random ionization events ahead of the ionization wave when the ionization fraction in front of the ionization wave is very low, analagous to the stochastic branching of streamers in 3D. Third, analysis shows high-energy runaway electrons accelerated in the cathode layer produce electron densities in the negative glow region over an order of magnitude above those in the positive column. Finally, future work to develop reduced order models of these two phenomena would improve the accuracy of fluid plasma models without the cost of PIC-DSMC simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

An improved method to access initial states in relativistic heavy-ion collisions

Abstract Observables in heavy-ion collisions are generally categorized into centralities, which reflect an average over events within a range of impact parameters including a wide variety of initial-state configurations. A multiple binning method using spectator neutrons within each centrality has been shown to provide access to events with rare initial-state conditions. This work suggests an improvement in quantifying the difference between standard centrality and spectator neutron binning towards accessing the initial-state properties. A selection of events with higher initial-state density at a fixed participating nucleon number was observed to result in larger final-state particle production and smaller elliptic flow. The relative difference between observables in centrality and spectator binning shows reduced sensitivity for the observables dominated by impact parameter fluctuations in the initial state, such as triangular flow. This property renders the spectator binning method a good candidate for separating geometric contributions from random fluctuations in the initial state towards final-state observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Binary Black Hole Mergers from Field Triples: Properties, Rates, and the Impact of Stellar Evolution

We consider the formation of binary black hole (BH) mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a BH binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov–Kozai cycles) with stellar evolution. After a BH triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most BH mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a BH merger rate in the range (0.3–1.3) Gpc{sup −3} yr{sup −1}, or up to ≈2.5 Gpc{sup −3} yr{sup −1} if the BH orbital planes have initially random orientation. Finally, we show that BH mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ≈10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Adiabatic phase-conserving processes for executing quantum operations with ultracold atoms

We have studied the regimes of deterministic single-atom Rydberg excitation in the conditions of Rydberg blockade and the methods of compensation for the dynamic phase of the wave function during the adiabatic passage. Using these methods, we have proposed schemes of single-qubit and two-qubit quantum states with mesoscopic atomic ensembles containing a random number of atoms, considred as quibits. The double adiabatic passage of the Förster resonance for two interacting atoms with a deterministic phase shift can be used for the implementation of two-qubit gates with reduced sensitivity of the gate fidelity to the fluctuations of the interatomic distance. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion

We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

High Pressure Hydrogen Pressure Relief Devices: Accelerated Life Testing and Application Best Practices

Pressure relief devices (PRDs ) are used to protect high pressure systems from burst failure caused by overpressurization. Codes and standards require the use of PRDs for the safe design of many pressurized systems. These systems require high reliability due to the risks associated with a burst failure. Hydrogen service can increase the risk of PRD failure due to material property degradation caused by hydrogen attack. The National Renewable Energy Laboratory (NREL) has conducted an accelerated life test on a conventional spring loaded PRD. Based on previous failures in the field, the nozzles specific to these PRDs are of particular interest. A nozzle in a PRD is a small part that directs the flow of fluid toward the sealing surface to maintain the open state of the valve once the spring force is overcome. The nozzle in this specific PRD is subjected to the full tensile force of the fluid pressure. These nozzles are made from 440C material, which is a type of hardened steel that is commonly chosen for high pressure applications because of its high strength properties. In a hydrogen environment, however, 440C is considered a worst case material since hydrogen attack results in a loss of almost all ductility and thus 440C is prone to fatigue and material failure. Accordingly, 440C is not recommended for hydrogen service. Conducting an accelerated life test on a PRD with 440C material provides information on necessary and sufficient conditions required to produce crack initiation and failure. The accelerated life test also provides information on other PRD failure modes that are somewhat statistically random in nature.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1999

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗
NASA NTRS2010

Scaled CMOS Technology Reliability Users Guide

The desire to assess the reliability of emerging scaled microelectronics technologies through faster reliability trials and more accurate acceleration models is the precursor for further research and experimentation in this relevant field. The effect of semiconductor scaling on microelectronics product reliability is an important aspect to the high reliability application user. From the perspective of a customer or user, who in many cases must deal with very limited, if any, manufacturer's reliability data to assess the product for a highly-reliable application, product-level testing is critical in the characterization and reliability assessment of advanced nanometer semiconductor scaling effects on microelectronics reliability. A methodology on how to accomplish this and techniques for deriving the expected product-level reliability on commercial memory products are provided.Competing mechanism theory and the multiple failure mechanism model are applied to the experimental results of scaled SDRAM products. Accelerated stress testing at multiple conditions is applied at the product level of several scaled memory products to assess the performance degradation and product reliability. Acceleration models are derived for each case. For several scaled SDRAM products, retention time degradation is studied and two distinct soft error populations are observed with each technology generation: early breakdown, characterized by randomly distributed weak bits with Weibull slope (beta)=1, and a main population breakdown with an increasing failure rate. Retention time soft error rates are calculated and a multiple failure mechanism acceleration model with parameters is derived for each technology. Defect densities are calculated and reflect a decreasing trend in the percentage of random defective bits for each successive product generation. A normalized soft error failure rate of the memory data retention time in FIT/Gb and FIT/cm2 for several scaled SDRAM generations is presented revealing a power relationship. General models describing the soft error rates across scaled product generations are presented. The analysis methodology may be applied to other scaled microelectronic products and their key parameters.

Microelectronics Reliability↗