SEARCH · Search NASA

Results for “SLIP”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A model for the motion of the Philippine Sea plate consistent with NUVEL-1 and geological data

We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.

Seno, Tetsuzo↗

NASA NTRS1995

A source-sink model of the generation of plate tectonics from non-Newtonian mantle flow

A model of mantle convection which generates plate tectonics requires strain rate- or stress-dependent rheology in order to produce strong platelike flows with weak margins as well as strike-slip deformation and plate spin (i.e., toroidal motion). Here, we employ a simple model of source-sink driven surface flow to determine the form of such a rheology that is appropriate for Earth's present-day plate motions. In this model, lithospheric motion is treated as shallow layer flow driven by sources and sinks which correspond to spreading centers and subduction zones, respectively. Two plate motion models are used to derive the source sink field. As originally implied in the simpler Cartesian version of this model, the classical power law rheologies do not generate platelike flows as well as the hypothetical Whitehead-Gans stick-slip rheology (which incorporates a simple self-lubrication mechanism). None of the fluid rheologies examined, however, produce more than approximately 60% of the original maximum shear. For either plate model, the viscosity fields produced by the power law rheologies are diffuse, and the viscosity lows over strike-slip shear zones or pseudo-margins are not as small as over the prescribed convergent-divergent margins. In contrast, the stick-slip rheology generates very platelike viscosity fields, with sharp gradients at the plate boundaries, and margins with almost uniformly low viscosity. Power law rheologies with high viscosity contrasts, however, lead to almost equally favorable comparisons, though these also yield the least platelike viscosity fields. This implies that the magnitude of toroidal flow and platelike strength distributions are not necessarily related and thus may present independent constraints on the determination of a self-consistent plate-mantle rheology.

Bercovici, David↗

NASA NTRS1994

Multicomponent Gas Diffusion and an Appropriate Momentum Boundary Condition

Multicomponent gas diffusion is reviewed with particular emphasis on gas flows near solid boundaries-the so-called Kramers-Kistemaker effect. The aim is to derive an appropriate momentum boundary condition which governs many gaseous species diffusing together. The many species' generalization of the traditional single gas condition, either as slip or stick (no-slip), is not obvious, particularly for technologically important cases of lower gas pressures and very dissimilar molecular weight gases. No convincing theoretical case exists for why two gases should interact with solid boundaries equally but in opposite flow directions, such that the total gas flow exactly vanishes. ln this way, the multicomponent no-slip boundary requires careful treatment The approaches discussed here generally adopt a microscopic model for gas-solid contact. The method has the advantage that the mathematics remain tractable and hence experimentally testable. Two new proposals are put forward, the first building in some molecular collision physics, the second drawing on a detailed view of surface diffusion which does not unphysically extrapolate bulk gas properties to govern the adsorbed molecules. The outcome is a better accounting of previously anomalous experiments. Models predict novel slip conditions appearing even for the case of equal molecular weight components. These approaches become particularly significant in view of a conceptual contradiction found to arise in previous derivations of the appropriate boundary conditions. The analogous case of three gases, one of which is uniformly distributed and hence non-diffusing, presents a further refinement which gives unexpected flow reversals near solid boundaries. This case is investigated alone and for aggregating gas species near their condensation point. In addition to predicting new physics, this investigation carries practical implications for controlling vapor diffusion in the growth of crystals used in medical diagnosis (e.g. mercuric iodide) and semiconductors.

Noever, David A.↗

NASA NTRS2001

Evaluation of Ti-48Al-2Nb Under Fretting Conditions

An investigation was conducted to examine the fretting behavior of lambda-TiAl (Ti-48Al-2Cr-2Nb) in contact with a nickel-base superalloy (Inconel 718) in air at temperatures from 23 to 550 C. Fretting wear experiments were conducted with 9.4-mm-diameter hemispherical Inconel (IN) 718 pins in contact with Ti-48Al-2Cr-2Nb flats (and the reverse) at loads from 1 to 40 N and fretting frequencies from 50 to 160 Hz with slip amplitudes from 50 to 200 microns for 1 to 20 million fretting cycles. The results were similar for both combinations of pin and flat. Reference fretting wear experiments were also conducted with 9.4-mm-diameter hemispherical Ti-6Al-4V pins in contact with IN718 flats. The interfacial adhesive bonds between Ti-48Al-2Cr-2Nb and IN718 in contact were generally stronger than the cohesive bonds in the cohesively weaker Ti-48Al-2Cr-2Nb. The failed Ti-48Al-2Cr-2Nb subsequently transferred to the IN718 surface at any fretting condition. The wear scars produced on Ti-48Al-2Cr-2Nb contained metallic and oxide wear debris, scratches, plastically deformed asperities, cracks, and fracture pits. Oxide layers readily formed on the Ti-48Al-2Cr-2Nb surface at 550 C, but cracks easily occurred in the oxide layers. Factors including fretting frequency, temperature, slip amplitude, and load influenced the fretting behavior of Ti-48Al-2Cr-2Nb in contact with IN718. The wear volume loss of Ti-48Al-2Cr-2Nb generally decreased with increasing fretting frequency. The increasing rate of oxidation at elevated temperatures up to 200 C led to a drop in wear volume loss at 200 C. However, the fretting wear increased as the temperature was increased from 200 to 550 C. The highest temperatures of 450 and 550 C resulted in oxide film disruption with generation of cracks, loose wear debris, and pits on the Ti-48Al-2Cr-2Nb wear surface. The wear volume loss generally increased as the slip amplitude increased. The wear volume loss also generally increased as the load increased. Increasing slip amplitude and increasing load both tended to produce more metallic wear debris, causing severe abrasive wear in the contacting metals.

Miyoshi, Kazuhisa↗

NASA NTRS1999

Flow Trajectories in Friction Stir Welding

In the Friction Stir Welding (FSW) process, a rotating, shouldered tool with a threaded pin is inserted under pressure along the seam of two pieces of metal which are tightly clamped together, and secured against a rigid anvil underneath. The rotating pin travels along the seam and through a combination of pressure and friction heating produces a zone of plastic deformation around the pin within the workpiece on either side of the seam. As the pin is moved in the direction of welding, the plasticised material moves around the tool and bonds together behind it. The elements of the material flow behavior are a combination of three elements. There is a rotational transport of material being carried around the tool, extrusion of material being forced around the pin on both sides into the cavity created behind it, and a lifting and dropping of material as it is stirred and mixed by the rotating action of the pin. It was assumed that rotational motion of the plastic zone is the primary mechanism for transport of material around the welding tool. A kinematic mathematical model was used to compute trajectories of material movement for various distributions of rotational slip within the plastic zone. These trajectories were then compared with the results of an experiment that produced radiographs of markers embedded in a workpiece that was welded with the FSW process. It was assumed that the copper wire markers retained their original length as the aluminum material flowed around them. The kinematic model included a constraint so that the displacements were such that the total length of the wire markers did not increase after deformation. There was good agreement between the calculated trajectories for the case of localized slip at the outer surface of the plastic deformation zone and the radiographs of the copper wire markers. The trajectories differed markedly from the radiograph traces when a distributed slip zone was assumed. It was concluded that the flow field could be described generally by a rotating plug model with slip occurring on the surface at the outer boundary of the plastic zone.

Bernstein, Edward L.↗

NASA NTRS2002

On the Boundary Condition for Water at a Hydrophobic, Dense Surface

We study the no-slip boundary conditions for water at a hydrophobic (graphite) surface using non-equilibrium molecular-dynamics simulations. For the planar Couette flow, we find a slip length of 64 nm at 1 bar and 300 K, decreasing with increasing system pressure to a value of 31 nm at 1000 bar. Changing the properties of the interface to from hydrophobic to strongly hydrophilic reduces the slip to 14 nm. Finally, we study the flow of water past an array of carbon nanotubes mounted in an inline configuration with a spacing of 16.4 x 16.4 nm. For tube diameters of 1.25 and 2.50 nm we find drag coefficients in good agreement with the macroscopic, Navier-Stokes values. For carbon nanotubes, the no-slip condition is valid to within the definition of the position of the interface.

Walther, J. H.↗

NASA NTRS2004

The Multistage Compressor Facility

Research and developments of new aerospace technologies is one of Glenn Research Center's specialties. One facility that deals with the research of aerospace technologies is the High-speed Multistage Compressor Facility. This facility will be testing the performance and efficiency of an Ultra Efficient Engine Technology (UEET) two-stage compressor. There is a lot of preparation involved with testing something of this caliber. Before the test article can be installed into the test rig, the facility must be fully operational and ready to run. Meaning all the necessary instrumentation must be calibrated and installed in the facility. The test rig should also be in safe operating condition, and the proper safety permits obtained. In preparation for the test, the Multistage Compressor Facility went through a few changes. For instance the facility will now be utilizing slip rings, the gearbox went through some maintenance, new lubrications systems replaced the old ones, and special instrumentation needs to be fine tuned to achieve the maximum amount of accurate data. Slips rings help gather information off of a rotating device - in this case from a shaft - onto stationary contacts. The contacts (or brushes) need to be cooled to reduce the amount of frictional heat produced between the slip ring and brushes. The coolant being run through the slip ring is AK-225, a material hazardous to the ozone. To abide by the safety regulations the coolant must be run through a closed chiller system. A new chiller system was purchased but the reservoir that holds the coolant was ventilated which doesn t make the system truly closed and sealed. My task was to design and have a new reservoir built for the chiller system that complies with the safety guidelines. The gearbox had some safety issues also. Located in the back of the gearbox an inching drive was set up. When the inching drive is in use the gears and chain are bare and someone can easily get caught up in it. So to prevent anyone from getting hurt in the gears I designed a chain guard. Additional information is included in the original extended abstract.

Flegel, Ashlie↗

NASA NTRS2005

ICESat Observations of Topographic Change in the Northern Segment of the 2004 Sumatra-Andaman Islands Earthquake Rupture Zone

The Andaman Islands are located 120 km east of the Sunda trench in the northern quarter of the 1300 km long rupture zone of the 2004 Sumatra-Andaman Islands earthquake inferred from the distribution of aftershocks. Initial field reports indicate that several meters of uplift and up to a meter of submergence occurred on the western and eastern shorelines of the Andaman Islands, respectively, associated with the earthquake (Bilham, 2005). Satellite images also document uplift of western shoreline coral reef platforms above sea level. Body-wave (Ji, 2005; Yamamaka, 2005) and tide-gauge (Ortiz, 2005) slip inversions only resolve coseismic slip in the southern one-third to one-half of the rupture zone. The amount of coseismic slip in the Andaman Islands region is poorly constrained by these inversions. The Ice, Cloud, and land Elevation Satellite (ICESat), a part of the NASA Earth Observing System, is being used to document the spatial pattern of Andaman Islands vertical displacements in order to constrain models of slip distribution in the northern part of the rupture zone. ICESat carries the Geoscience Laser Altimeter System (GLAS) that obtains elevation measurements from 80 m diameter footprints spaced 175 m apart along profiles. For surfaces of low slope, single-footprint absolute elevation and horizontal accuracies of 10 cm and 6 m (1 sigma), respectively, referenced to the ITRF 2002 TOPEX/Poseidon ellipsoid are being obtained. Laser pulse backscatter waveforms enable separation of ground topography and overlying vegetation cover. During each 33-day observing period ICESat acquires three profiles crossing the Andaman Islands. A NNE-SSW oriented track consists of 1600 laser footprints along the western side of North, Middle, and South Andaman Islands and 240 laser footprints across the center of Great Andaman Island. Two NNW-SSE tracks consist of 440 footprints across Middle Andaman Island and 25 footprints across the west side of Sentinel Island. Cloud-free profiles were acquired in the fall of 2003 and 2004. During February-March, 2005 ICESat's precise pointing capability will be used to exactly repeat these three profiles, with a cross-track accuracy of better than 100 m, providing trench- parallel and -perpendicular observations of topographic change of the Andaman Islands that will compliment geodetic field surveys. The observed elevation changes will be compared to models of coseismic deformation associated with the mainshock and large aftershocks in the Andaman Islands region.

Harding, David↗

NASA NTRS2008

Visual Feedback for Rover-based Coring

Technology for coring from a low-mass rover has been developed to enable core sample acquisition where a planetary rover experiences moderate slip during the coring operation. A new stereo vision technique, Absolute Motion Visual Odometry, is used to measure rover slip during coring and the slip is accommodated through corresponding arm pose updating. Coring rate is controlled by feedback of themeasured force of the coring tool against the environment. Test results in the JPL Marsyard show for the first time that coring from a low-mass rover with slip is feasible.

sampling↗

NASA NTRS2005

Stroboscopic Goggles for Reduction of Motion Sickness

A device built around a pair of electronic shutters has been demonstrated to be effective as a prototype of stroboscopic goggles or eyeglasses for preventing or reducing motion sickness. The momentary opening of the shutters helps to suppress a phenomenon that is known in the art as retinal slip and is described more fully below. While a number of different environmental factors can induce motion sickness, a common factor associated with every known motion environment is sensory confusion or sensory mismatch. Motion sickness is a product of misinformation arriving at a central point in the nervous system from the senses from which one determines one s spatial orientation. When information from the eyes, ears, joints, and pressure receptors are all in agreement as to one s orientation, there is no motion sickness. When one or more sensory input(s) to the brain is not expected, or conflicts with what is anticipated, the end product is motion sickness. Normally, an observer s eye moves, compensating for the anticipated effect of motion, in such a manner that the image of an object moving relatively to an observer is held stationary on the retina. In almost every known environment that induces motion sickness, a change in the gain (in the signal-processing sense of gain ) of the vestibular system causes the motion of the eye to fail to hold images stationary on the retina, and the resulting motion of the images is termed retinal slip. The present concept of stroboscopic goggles or eyeglasses (see figure) is based on the proposition that prevention of retinal slip, and hence, the prevention of sensory mismatch, can be expected to reduce the tendency toward motion sickness. A device according to this concept helps to prevent retinal slip by providing snapshots of the visual environment through electronic shutters that are brief enough that each snapshot freezes the image on each retina. The exposure time for each snapshot is less than 5 ms. In the event that a higher rate of strobing is necessary for adequate viewing of the changing scene during rapid head movements, the rate of strobing (but not the exposure time) can be controlled in response to the readings of rate-of-rotation sensors attached to the device.

Reschke, M. F.↗

NASA NTRS2005

Tool for Viewing Faults Under Terrain

Multi Surface Light Table (MSLT) is an interactive software tool that was developed in support of the QuakeSim project, which has created an earthquake- fault database and a set of earthquake- simulation software tools. MSLT visualizes the three-dimensional geometries of faults embedded below the terrain and animates time-varying simulations of stress and slip. The fault segments, represented as rectangular surfaces at dip angles, are organized into collections, that is, faults. An interface built into MSLT queries and retrieves fault definitions from the QuakeSim fault database. MSLT also reads time-varying output from one of the QuakeSim simulation tools, called "Virtual California." Stress intensity is represented by variations in color. Slips are represented by directional indicators on the fault segments. The magnitudes of the slips are represented by the duration of the directional indicators in time. The interactive controls in MSLT provide a virtual track-ball, pan and zoom, translucency adjustment, simulation playback, and simulation movie capture. In addition, geographical information on the fault segments and faults is displayed on text windows. Because of the extensive viewing controls, faults can be seen in relation to one another, and to the terrain. These relations can be realized in simulations. Correlated slips in parallel faults are visible in the playback of Virtual California simulations.

Siegel, Herbert, L.↗

NASA NTRS2003

Microscope Cells Containing Multiple Micromachined Wells

Tech Briefs, May 2003 19 Manufacturing Microscope Cells Containing Multiple Micromachined Wells The cost per cell has been reduced substantially. John H. Glenn Research Center, Cleveland, Ohio An improved design for multiple-well microscope cells and an associated improved method of fabricating them have been devised. [As used here, "well" denotes a cavity that has a volume of about 1 or 2 L and that is used to hold a sample for examination under a microscope. As used here, "cell" denotes a laminate, based on a standard 1- by 3-in. (2.54- by 7.62-cm) microscope slide, that comprises (1) the slide as the lower layer, (2) an intermediate layer that contains holes that serve as the wells, and (3) a top layer that either consists of, or is similar to, a standard microscope-slide cover slip.] The improved design and method of fabrication make it possible to increase (relative to a prior design and method of fabrication) the number of wells per cell while reducing the fabrication loss and reducing the cost per cell to about one-tenth of the prior value. In the prior design and method, the slide, well, and cover-slip layers were made from silicate glass. The fabrication of each cell was a labor-intensive process that included precise cutting and grinding of the glass components, fusing of the glass components, and then more grinding and polishing to obtain desired dimensions. Cells of the prior design were expensive and fragile, the rate of loss in fabrication was high, and the nature of the glass made it difficult to increase the number of cells per well. Efforts to execute alternative prior designs in plastic have not yielded satisfactory results because, for typical applications, plastics are not sufficiently thermally or chemically stable, not sufficiently optically clear, and/or not hard enough to resist scratching. The figure depicts a cell of the present improved type. The slide and cover-slip layers are made of a low-thermal-expansion glass (Pyrex(TradeMark) or equivalent) and the intermediate (well layer) is made of SiO2 - a combination of materials that results in a laminate stronger than one made from layers of silicate glass. Before the layers are assembled into the laminate, the SiO2 layer is micromachined to form the wells plus shallow grooves that, when subsequently covered with the cover slip, become capillary channels that are used to fill the wells with samples. The micromachining is accomplished by use of the same patterning and etching techniques used to fabricate microelectromechanical systems (MEMS).

Turner, Walter↗

NASA NTRS2013

Soft Decision Analyzer

The Soft Decision Analyzer (SDA) is an instrument that combines hardware, firmware, and software to perform realtime closed-loop end-to-end statistical analysis of single- or dual- channel serial digital RF communications systems operating in very low signal-to-noise conditions. As an innovation, the unique SDA capabilities allow it to perform analysis of situations where the receiving communication system slips bits due to low signal-to-noise conditions or experiences constellation rotations resulting in channel polarity in versions or channel assignment swaps. SDA s closed-loop detection allows it to instrument a live system and correlate observations with frame, codeword, and packet losses, as well as Quality of Service (QoS) and Quality of Experience (QoE) events. The SDA s abilities are not confined to performing analysis in low signal-to-noise conditions. Its analysis provides in-depth insight of a communication system s receiver performance in a variety of operating conditions. The SDA incorporates two techniques for identifying slips. The first is an examination of content of the received data stream s relation to the transmitted data content and the second is a direct examination of the receiver s recovered clock signals relative to a reference. Both techniques provide benefits in different ways and allow the communication engineer evaluating test results increased confidence and understanding of receiver performance. Direct examination of data contents is performed by two different data techniques, power correlation or a modified Massey correlation, and can be applied to soft decision data widths 1 to 12 bits wide over a correlation depth ranging from 16 to 512 samples. The SDA detects receiver bit slips within a 4 bits window and can handle systems with up to four quadrants (QPSK, SQPSK, and BPSK systems). The SDA continuously monitors correlation results to characterize slips and quadrant change and is capable of performing analysis even when the receiver under test is subjected to conditions where its performance degrades to high error rates (30 percent or beyond). The design incorporates a number of features, such as watchdog triggers that permit the SDA system to recover from large receiver upsets automatically and continue accumulating performance analysis unaided by operator intervention. This accommodates tests that can last in the order of days in order to gain statistical confidence in results and is also useful for capturing snapshots of rare events.

Steele, Glen↗

NASA NTRS2016

Dynamic Visual Acuity and Landing Sickness in Crewmembers Returning from Long-Duration Spaceflight

Long-term exposure to microgravity causes sensorimotor adaptations that result in functional deficits upon returning to a gravitational environment. At landing the vestibular system and the central nervous system, responsible for coordinating head and eye movements, are adapted to microgravity and must re-adapt to the gravitational environment. This re-adaptation causes decrements in gaze control and dynamic visual acuity, with astronauts reporting oscillopsia and blurred vision. Dynamic visual acuity (DVA) is assessed using an oscillating chair developed in the Neuroscience Laboratory at JSC. This chair is lightweight and easily portable for quick deployment in the field. The base of the chair is spring-loaded and allows for manual oscillation of the subject. Using a metronome, the chair is vertically oscillated plus or minus 2 cm at 2 Hz by an operator, to simulate walking. While the subject is being oscillated, they are asked to discern the direction of Landolt-C optotypes of varying sizes and record their direction using a gamepad. The visual acuity thresholds are determined using an algorithm that alters the size of the optotype based on the previous response of the subject using a forced-choice best parameter estimation that is able to rapidly converge on the threshold value. Visual acuity thresholds were determined both for static (seated) and dynamic (oscillating) conditions. Dynamic visual acuity is defined as the difference between the dynamic and static conditions. Dynamic visual acuity measures will be taken prior to flight (typically L-180, L-90, and L-60) and up to eight times after landing, including up to 3 times on R plus 0. Follow up measurements will be taken at R plus 1 (approximately 36 hours after landing). Long-duration International Space Station crewmembers will be tested once at the refueling stop in Europe and once again upon return to Johnson Space Center. In addition to DVA, subjective ratings of motion sickness will be recorded throughout the testing. Using the chair as a portable and reliable way to test DVA, we aim to test returning astronauts to assess the amount of retinal slip that they experience. By comparing these measurements to their motion sickness scores (using a scale of 1 to 20 where 20 is vomiting), we will correlate the amount of retinal slip to the level of motion sickness experienced. In addition to testing this in returning astronauts, we will perform ground-based studies to determine the effectiveness of stroboscopic goggles in reducing retinal slip and improving DVA. Finally, we will employ stroboscopic goggles in the field to astronauts experiencing high levels of motion sickness to minimize retinal slip and reduce their symptoms.

Rosenberg, M.J.F↗

NASA NTRS2019

Continuously Improving Parametric Modeling with Historical Data on the ICESat-2 Mission

This paper delves into the details of the Joint Confidence Level (JCL) process performed for the Ice, Cloud, and Land Elevation Satellite (ICESat)-2 mission and how past performance was incorporated into subsequent JCL models to enable the project to continuously analyze potential slips to their launch readiness date (LRD). One year prior to the mission Preliminary Review (mPDR), the JCL model development process began. The first model was well received at the mPDR, held on October 10, 2012, and the input received by the Standing Review Board was incorporated into the model for the official data drop for key decision point (KDP)-C. The 70% JCL results of the October 2012 mPDR model forecast an LRD of February 2017 and associated cost of $830M. This result in 2012 immediately highlighted potential challenges with the project-planned LRD of July 2016. The year following the mPDR, the project had sustained a oneyear slip in the LRD due to problematic systems engineering requirement issues which impacted all project subsystems. This slip moved the project planned LRD from July 2016 to July 2017, an additional 5 months beyond the 2012 model’s 70% JCL result for the LRD of February 2017. As the project was quickly approaching the mission Critical Design Review (mCDR), the need for reliable JCL results increased significantly. The project held discussions on the JCL modeling process and focused on the input uncertainty distributions. Specifically, to identify the uncertainty distributions that the 2012 mPDR model would have needed to produce a 70% LRD result of July 2017. This led the project to compare multiple uncertainty distributions, and ultimately spurred the project to utilize uncertainty distributions that incorporated project past performance and historical data to forecast potential LRD slips. The revised results, created in 2014 and utilizing the new uncertainty distributions, showed that with 70% confidence, the ICESat-2 mission would launch in August 2018 with a cost of $1,044M. Today, ICESat-2 is scheduled to launch on September 15, 2018 with a project management (PM) agreement value of $1,056M. This illustrates how a JCL model can be continuously improved to produce valuable results for a project, even in cases of LRD delays. The primary reason for the ICESat-2 LRD delay is due to a laser failure on the primary instrument. Laser failure was one of the highest risk and uncertainty drivers within the JCL model. The project placed the most risk in this area of the model, and the model further identified the laser as the top risk driver and contributor to the LRD result. This further illustrates how a JCL can be used to predict and quantify possible issues on new technology missions.

Krygiel, Joseph↗

NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗

DOE OSTI2025

Comparative studies of cross-phase dynamics in turbulent momentum flux and particle flux at the tokamak edge

Turbulent transport events, including turbulent transport flux of momentum (i.e., turbulent momentum flux or Reynolds stress) and turbulent transport flux of particle (i.e., turbulent particle flux), have important effects on the confinement performance of magnetic confinement fusion devices. Poloidal Reynolds stress is the ensemble average of the product of radial velocity fluctuations and poloidal velocity fluctuations, i.e., $\langle {\widetilde{v}}_{r}{\widetilde{v}}_{\theta }\rangle$. Turbulent particle flux is the ensemble average of the product of radial velocity fluctuations and density fluctuations, i.e., $\langle \widetilde{n}{\widetilde{v}}_{r}\rangle$. Changes in either amplitude of fluctuations or cross phase between fluctuations can cause changes in turbulent transport. In this paper, cross-phase dynamics in the Reynolds stress and turbulent particle flux at the tokamak edge are studied in detail. Reynolds stress and turbulent particle flux are, respectively, written as the product of fluctuation amplitudes and an average cross-phase factor. The mathematical expressions of the average cross-phase factors are derived. The average cross-phase factors and the power spectra of cross phase are obtained using experimental measurement data. It is found that the cross-phase dynamics in Reynolds stress and particle flux are very different. Reynolds stress is found to be more sensitive to cross phase than particle flux is. In the strong $E\times B$ shear layer, spatial slips of cross phase lead to the obvious radial gradient of Reynolds stress. In the no/weak $E\times B$ shear region, the cross phase in Reynolds stress tends to lock. Here, phase locking refers to that the power spectra of phase tend to distribute around a fixed phase which does not change with radial position, while phase slip means that the power spectra of cross phase tend to distribute around a phase that varies with radial position. Phase slip or locking mainly describes the central phase weighted by the power spectra, while the phase scattering mainly describes the dispersion of the power spectrum distribution of the phase. The increased scattering of cross phase, which indicates the power spectra distribution of the phase is more dispersed, contributes to the decreased Reynolds stress for higher collisionality. The cross phase in particle flux tends to lock in both strong and no/weak shear regions. The degree of scattering of cross phase in the particle flux does not change obviously as collisionality increases. For higher collisionality, it is the increased density fluctuation amplitude rather than cross-phase dynamics that leads to the increased particle flux. The underlying physical mechanism that causes Reynolds stress and particle flux to exhibit different phase dynamics is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗