Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “metric geometry”

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

Sources of Uncertainty in Ultrasonic Phase-Based Assessment of Interfacial Stiffness in Adhesive Joints

Quantitative nondestructive evaluation (NDE) of adhesive joint strength is critical for certifying fastener-free bonded aerospace structures. Unfortunately, available adhesive NDE approaches are only suitable for detecting gross defects and delamination and cannot reliably detect variations in strength in intact bonds. Intact (i.e., nominally defect free) bonds are challenging to inspect because they, by definition, contain no macroscopic defects that would produce a strong NDE signature. Even in the absence of obvious defects, however, bond strength can vary substantially due to factors such as contamination, variations in the adhesive cure processes, and environmental exposure. For ultrasonic bond inspections, the effect of interactions from an imperfect adhesive interface are subtle and may include changes in the frequency content or phase of the ultrasonic signal. The bond condition can be deduced from these subtle changes in the acoustic response using an accompanying physics-based model of adhesive joint interaction. For example, researchers have used the spring network adhesion model of Baik and Thompson [1] to determine the interfacial stiffness of the joint from its ultrasonic response (e.g., [2]). This interfacial stiffness metric describes the bonds resistance to separation between the adhesive and adherent and has shown a strong correlation with joint strength. In this study, interfacial stiffness is assessed using the ultrasonic phase change associated with bond interaction in varying bond conditions. The phase change is measured precisely using a phase-locked-loop instrument over a swept range of frequencies. The phase transition at resonance, where the bond thickness is half the acoustic wavelength, produces a unique signature with which to fit a physics-based model. Interfacial stiffness is determined by minimizing the difference between the theoretical and experimental phase response. Experimental results are presented for glass specimens bonded with aerospace adhesives and a sensitivity analysis is performed to examine the effects of parameters that may be uncertain in the industrial inspection environment – namely bond thickness, adhesive and adherent properties, and porosity. The results are a first step toward an uncertainty quantified version of the measurement technique, show correlation between interfacial stiffness and bond strength, and are another key step toward maturing this approach for use in complex materials and joint geometries. [1] Baik, J.M. and Thompson, R.B., 1984. “Ultrasonic scattering from imperfect interfaces: a quasi-static model”. Journal of Nondestructive Evaluation, 4. [2] Haldren, H.A., Perey, D.F., Yost, W.T., Cramer, K.E., and Gupta, M.C., 2019. “Swept-frequency ultrasonic phase evaluation of adhesive bonding in tri-layer structures”. Journal of the Acoustical Society of America, 145.

Matt Webster↗

Fusing GeoNEX and VIIRS Surface BRDF Retrievals: Exploring a GEO-LEO Synergy

The Bidirectional Reflectance Distribution Function or BRDF, which describes the dependency of surface reflectance on the illumination-view geometries, are the foundation of many high-level satellite products for terrestrial and aquatic system monitoring. The latest geostationary sensors like GOES ABI provide high frequent (~10 minutes) observations of the Earth surface that feature continuously changing sun angles, allowing us to retrieve surface BRDF with dedicated atmospheric correction algorithms like MAIAC (Multi-Angle Implementation of Atmospheric Correction). For mid-latitude locations, because geostationary satellites have fixed view angles in the back-scattering directions, the angular sampling of surface BRDF by GEO sensors is not comprehensive. This study explores a GEO-LEO synergy to address this issue. We first extract concurrent GeoNEX and VIIRS BRDF data with the best quality (cloud-free and low aerosol loading) at chosen AERONET sites. We then compare the magnitude and the shape factors of the two set of BRDF parameters as well as their variations through the season. We calculate the “distances” between the GeoNEX and VIIRS BRDF by using them to cross-predict the top-of-atmosphere reflectance measured by their counterpart and evaluating the corresponding prediction errors. This metric allows us to derive a set of optimized BRDF parameters that minimize such distances or prediction errors, which are considered as the fused BRDF result. We validate the algorithm with reserved AERONET data and then apply it to generate the GEO-LEO BRDF synergy over CONUS. We expect the fused BRDF to have reduced uncertainties as compared to the source GeoNEX or VIIRS data and may find broadly application in deriving other high-level satellite products.

Geostationary satellite↗

ISRU Pilot Excavator: Bucket Drum Scaling Experimental Results

NASA’s Space Technology Mission Directorate (STMD) is funding the development of a robotic excavator called the “ISRU Pilot Excavator” which will be a technology demonstration of excavating and transporting 10 metric tons of lunar regolith on the surface of the moon with a 30k-class robotic excavator. ISRU Pilot Excavator will be the next generation of robotic excavator to use bucket drums as excavation tools. This is an evolution of the Regolith Advanced Surface Systems Operations Robot (RASSOR) developed at NASA’s Kennedy Space Center. Bucket drums are hollow cylinders with regularly spaced scoops around the perimeter. The drums rotate in one direction to collect regolith with the scoops. The regolith slides down an internal baffling system inside the drum which prevents the regolith from falling back out of the scoops. The captured regolith can then be transported while held in the drum and then deposited by rotating the drum in the opposite direction allowing the regolith to slide back down the baffling and out of the excavation scoops. Bucket drums were developed by Lockheed Martin in 2008 and used on multiple robotic excavator prototypes ever since. However the forces on a bucket drum and considerations for scaling have not been measured in detail. Bucket drums are challenging to model using classical blade\bucket equations because of their unique geometry. Therefore this experiment was performed to measure the forces on three bucket drums of the same geometry at different scales. Small: 9.4” (239mm) dia. x 8.1” (206mm) width, Medium: 11.6” (294mm) dia. x 10” (254mm) width, and Large: 17” (432mm) dia. x 14.1”(358mm) width. The test stand consisted of an actuated gantry with controlled motion in the vertical (Z) and horizontal (X) axes and a single rotation axis (R). The bucket drums were individually mounted to the rotary axis of the test stand and translated across a prepared bed of BP-1 lunar regolith simulant at a specified linear speed and cutting depth. The test stand was outfitted with a torque sensor in line with the rotation of the drum (R) and a 3 axis (X, Y, and Z) load cell. In addition to the three sizes of bucket drums the linear excavation speed and cutting depth were test variables. The results of these experiments show the relationship between the three scales of bucket drums for factors such as: excavation force, torque due to regolith rotation inside the drum, excavation energy, time to fill, etc. and will be discussed in detail in this paper. This fundamental data will be used in the design of the ISRU Pilot Excavator and can inform the design of future bucket drum excavators.

Jason Michael Schuler↗

ISRU Pilot Excavator: Bucket Drum Scaling Experimental Results

NASA’s Space Technology Mission Directorate (STMD) is funding the development of a robotic excavator called the “ISRU Pilot Excavator” (IPEx) which will be a technology demonstration of excavating and transporting 10 metric tons of lunar regolith on the surface of the moon with a 30kg-class robotic excavator. IPEx will be the next generation of robotic excavators to use bucket drums as excavation tools. This is an evolution of the Regolith Advanced Surface Systems Operations Robot (RASSOR) developed at NASA’s Kennedy Space Center (KSC). Bucket drums are hollow cylinders with regularly spaced scoops around the perimeter. The drums rotate in one direction to collect regolith with the scoops. The regolith slides down an internal baffling system inside the drum which prevents the regolith from falling back out of the scoops (see Figure 1). The captured regolith can then be transported while held in the drum and then deposited by rotating the drum in the opposite direction allowing the regolith to slide back down the baffling and out of the excavation scoops. Bucket drums were developed by Lockheed Martin in 2008 and used on multiple robotic excavator prototypes ever since. However, the forces on a bucket drum and considerations for scaling have not been measured in detail. Bucket drums are challenging to model using classical blade\bucket equations because of their unique geometry. Therefore, this experiment was performed to measure the forces on three bucket drums of the same geometry at different scales. Small: 9.4” (239mm) dia. x 8.1” (206mm) width, Medium: 11.6” (294mm) dia. x 10” (254mm) width, and Large: 17” (432mm) dia. x 14.1” (358mm) width. The test stand consisted of an actuated gantry with controlled motion in the vertical (Z) and horizontal (X) axes and a single rotation axis (R). The bucket drums were individually mounted to the rotary axis of the test stand and translated across a prepared bed of BP-1 lunar regolith simulant at a specified linear speed and cutting depth. The test stand was outfitted with a torque sensor in line with the rotation of the drum (R) and a 3 axis (X, Y, and Z) load cell. In addition to the three sizes of bucket drums the linear excavation speed and cutting depth were test variables. The results of these experiments show the relationship between the three scales of bucket drums for factors such as: excavation force, torque due to regolith rotation inside the drum, excavation energy, time to fill, etc. and will be discussed in detail in this paper. This fundamental data will be used in the design of IPEx and can inform the design of future bucket drum excavators.

RASSOR↗

Geometric origin of the energy-momentum tensor improvement terms

In a flat background, the canonical energy momentum tensor of Lorentz and conformally invariant matter field theories can be improved to a symmetric and traceless tensor that gives the same conserved charges. We argue that the geometric origin of this improvement process is unveiled when the matter theory is coupled to metric-affine gravity. In particular, we show that the Belinfante-Rosenfeld improvement terms correspond to the matter theory’s hypermomentum. The improvement terms in conformally invariant matter theories are also related to the hypermomentum; however, a general proof would require an extended investigation. We demonstrate our results through various examples, such as the free massless scalar, the Maxwell field, Abelian p-forms, the Dirac field, and a nonunitary massless scalar field. Possible applications of our method for theories that break Lorentz or special conformal invariance are briefly discussed.

classical solutions in field theory↗

Experimental Study of the Invariance of Pressure Gain with Respect to the Dynamics of Multiple Competing Waves in a Rotating Detonation Combustor

Changes in the overall performance of a rotating detonation combustor with respect to changes in operation mode and wave dynamics arising by operation with fixed inlet/exit geometry but at different combustor, lengths are investigated experimentally. The air inlet, fuel injection, and exit constriction geometry are held constant while only the length of the detonation channel is varied from 71 to 137 mm (which corresponds to about 10 to 20 channel widths). Operation of H2/air over a range of air mass flow rates and equivalence ratios are considered for every chamber length. The number and speed of (primary) detonation and secondary waves are characterized through high-speed pressure measurements in the detonation channel and aft chemiluminescence videos. The number of waves is found to increase with length while detonation wave speed decreases significantly. Particular emphasis is given to characterize a phenomenon that is observed at operation with longer combustor lengths and higher mass flow rates. The phenomenon manifests as a super-cycle behavior with a period equal to many detonation wave rotational periods and is characterized by a periodic and structured ascending/descending sequence of the number, speed, and direction of both (primary) detonation and secondary waves. This phenomenon is likely a manifestation of the system failing in achieving operation with a higher number of detonation waves as length and/or mass flow rate are increased. The performance of the device is quantified in terms of measured thrust and pressure gain (through the use of the equivalent available pressure). Both metrics are essentially found to be invariant with respect to combustor length and most importantly, mode of operation. Surprisingly, even operation with complex wave dynamics arising from transitions between multiple competing wave systems does not appear to alter the overall global performance of the device but rather, it remains defined by the total (capture) air mass flow rate, equivalence ratio, and inlet/outlet areas.

33 ADVANCED PROPULSION SYSTEMS↗

The turbomachine blading design using S2-S1 approach

The boundary conditions corresponding to the design problem when the blades being simulated by the bound vorticity distribution are presented. The 3D flow is analyzed by the two steps S2 - S1 approach. In the first step, the number of blades is supposed to be infinite, the vortex distribution is transformed into an axisymmetric one, so that the flow field can be analyzed in a meridional plane. The thickness distribution of the blade producing the flow channel striction is taken into account by the modification of metric tensor in the continuity equation. Using the meridional stream function to define the flow field, the mass conservation is satisfied automatically. The governing equation is deduced from the relation between the azimuthal component of the vorticity and the meridional velocity. The value of the azimuthal component of the vorticity is provided by the hub to shroud equilibrium condition. This step leads to the determination of the axisymmetric stream sheets as well as the approximate camber surface of the blade. In the second step, the finite number of blades is taken into account, the inverse problem corresponding to the blade to blade flow confined in each stream sheet is analyzed. The momentum equation implies that the free vortex of the absolute velocity must be tangential to the stream sheet. The governing equation for the blade to blade flow stream function is deduced from this condition. At the beginning, the upper and the lower surfaces of the blades are created from the camber surface obtained from the first step with the assigned thickness distribution. The bound vorticity distribution and the penetrating flux conservation applied on the presumed blade surface constitute the boundary conditions of the inverse problem. The detection of this flux leads to the rectification of the geometry of the blades.

Luu, T. S.↗

Covariance Shaping Over Riemannian Manifolds for Massive MIMO Communication

Acquiring accurate instantaneous channel state information (CSI) is a challenging aspect of massive multi-input multi-output (MIMO) communication. Utilizing statistical information, such as channel covariance matrix, to design statistical beamforming vectors is robust when compared to instantaneous CSI. In this paper, we propose a novel MIMO covariance shaping scheme over Riemannian manifolds. It serves as an effective statistical beamforming solution to a number of close proximity user equipment (UE) that are undergoing substantial channel correlation. Proposed algorithm exploits the Hermitian positive definite nature of covariance matrices lying over Riemannian manifold. We introduce Wasserstein distance function as a Riemannian metric to measure distances between channel covariance matrices. Furthermore, K-means clustering technique is utilized to effectively identify the optimal shape of effective optimal covariance matrices. Our findings suggest that maximizing the geodesic distance between covariance matrices ultimately leads to a corresponding increase in the network throughput, as determined by the beamforming vector used to shape the covariance matrices. Simulation results validate that the proposed solution converges faster than Euclidean-based state-of-the-art, while maintaining the same computational complexity. Finally, the sum rate performance asymptotically achieves full capacity for two-UE case and more than 96% of the upper bound exhaustive search benchmark for multi-UE scenario.

42 ENGINEERING↗

Image Segmentation, Registration, Compression, and Matching

A novel computational framework was developed of a 2D affine invariant matching exploiting a parameter space. Named as affine invariant parameter space (AIPS), the technique can be applied to many image-processing and computer-vision problems, including image registration, template matching, and object tracking from image sequence. The AIPS is formed by the parameters in an affine combination of a set of feature points in the image plane. In cases where the entire image can be assumed to have undergone a single affine transformation, the new AIPS match metric and matching framework becomes very effective (compared with the state-of-the-art methods at the time of this reporting). No knowledge about scaling or any other transformation parameters need to be known a priori to apply the AIPS framework. An automated suite of software tools has been created to provide accurate image segmentation (for data cleaning) and high-quality 2D image and 3D surface registration (for fusing multi-resolution terrain, image, and map data). These tools are capable of supporting existing GIS toolkits already in the marketplace, and will also be usable in a stand-alone fashion. The toolkit applies novel algorithmic approaches for image segmentation, feature extraction, and registration of 2D imagery and 3D surface data, which supports first-pass, batched, fully automatic feature extraction (for segmentation), and registration. A hierarchical and adaptive approach is taken for achieving automatic feature extraction, segmentation, and registration. Surface registration is the process of aligning two (or more) data sets to a common coordinate system, during which the transformation between their different coordinate systems is determined. Also developed here are a novel, volumetric surface modeling and compression technique that provide both quality-guaranteed mesh surface approximations and compaction of the model sizes by efficiently coding the geometry and connectivity/topology components of the generated models. The highly efficient triangular mesh compression compacts the connectivity information at the rate of 1.5-4 bits per vertex (on average for triangle meshes), while reducing the 3D geometry by 40-50 percent. Finally, taking into consideration the characteristics of 3D terrain data, and using the innovative, regularized binary decomposition mesh modeling, a multistage, pattern-drive modeling, and compression technique has been developed to provide an effective framework for compressing digital elevation model (DEM) surfaces, high-resolution aerial imagery, and other types of NASA data.

Yadegar, Jacob↗

LatticeAnalytics: Strut-Level Visualization and Inspection of Additively Manufactured Lattice Structures

Additive manufacturing (AM) is revolutionizing the production of custom components with complex internal geometries, essential for high-performance applications in diverse fields such as medicine and defense. These AM parts optimize strength while minimizing weight by utilizing internal lattice structures consisting of large quantities of small interconnected struts. However, the complexity of these structures, combined with the challenges of using X-ray Computed Tomography (XCT) data, makes validation of part reliability difficult. This ultimately inhibits the development of novel parts for our collaborating material scientists. Here, we introduce LatticeAnalytics, a novel framework specifically designed for visual inspection of defects in these lattice structures. Our framework offers an end-to-end solution that includes the data management of XCT scans, enables remote access for geographically dispersed teams through a web-based dashboard, and incorporates novel visualizations. Our analysis is facilitated by a coarse alignment between the lattice’s nominal model, a spatial graph, and the XCT data. We employ a simple VR-based approach for fast and rough alignment, followed by an offline registration and identification of the struts. With the nodes and struts aligned and identified in the volume, our framework allows querying of subvolumes containing a single strut at multiple resolutions. This avoids computation over the entire lattice and also allow for easy parallelization of down-stream computations, such as strut-specific metrics. To depict a fast overview of the strut quality, we introduce two innovative visual encodings, crucial for our collaborators’ research in creating novel AM parts: the Contour View and the Roughness Map, which depict critical geometrical and surface features of individual struts in standardized two 2D views. We evaluated the integrated system through expert interviews. The feedback confirms the framework’s practicality and its effectiveness in enhancing current inspection workflows. It solves major bottlenecks for our collaborators, ultimately helping them create novel parts with advanced properties.

Miao, Haichao [Lawrence Livermore National Laborat↗

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

Annihilation-to-nothing: a quantum gravitational boundary condition for the Schwarzschild black hole

The interior of a static Schwarzschild metric can be written in terms of two functions, similar to some models of anisotropic cosmology. With a suitable choice of canonical variables, we solve the Wheeler-DeWitt equation (WDW) inside the horizon of a Schwarzschild black hole. By imposing classicality near the horizon, and requiring boundedness of the wave function, we get a rather generic solution of the WDW equation, whose steepest-descent solution, i.e., the ridge of the wave function, coincides nicely with the classical trajectory. However, there is an ambiguity in defining the arrow of time which leads to two possible interpretations—(i) if there is only one arrow of time, one can infer that the steepest-descent of the wave function follows the classical trajectory throughout: coming from the event horizon and going all the way down to the singularity, while (ii) if there are two different arrows of time in two separate regimes, it can be inferred that the steepest-descent of the wave function comes inwards from the event horizon in one region while it moves outwards from the singularity in the other region, and there exists an annihilation process of these two parts of the wave function inside the horizon. Additionally, adopting the second interpretation could shed light on the information loss paradox: as time goes on, probabilities for histories that include black holes and singularities decay to zero and eventually only trivial geometries dominate.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

ACES M and S: Unmitigated Factorial Encounter Study on DAA/TCAS Interoperability

Realization of the expected proliferation of Unmanned Aircraft System (UAS) operations in the National Airspace System (NAS) depends on the development and validation of standards for UAS Detect and Avoid (DAA) Systems. The RTCA Special Committee 228 is charged with leading the development of draft Minimum Operational Performance Standards (MOPS) for UAS DAA Systems. NASA, as a participating member of RTCA SC-228 is committed to supporting the development and validation of draft requirements for DAA alerting and guidance systems. This presentation contains the results of two combinatorial encounter analysis studies using NASA's SAA Control fast-time simulation capability for this purpose. In these studies, encounters between two aircraft were simulated one at a time for the full factorial combination of encounter geometries (e.g., encounter angle, CPA offset) and aircraft performance (e.g., ownership and intruder ground speeds and vertical rates). The first study analyzes the relationships (e.g., timeline) between the different alerting-safety regions in the SC-228 MOPS (in order of increasing severity): 1) DAA warning alert, 2) well clear recovery (WCR) guidance, 3) DAA-Collision Avoidance (CA), and 4) TCAS RA. This study will focus primarily on encounter situations in which TCAS RA occurs prior to any of the other alerting-safety boundaries. In particular, this study will investigate whether using vertical distance or vertical distance at closest point of approach (i.e., vertical miss distance or VMD) is more appropriate for the definition of the DAA-CA region. In addition, cases where transitions between different regions skip an intermediate region will be analyzed. The second study in this presentation explores a proposal to use an altitude rate error threshold to determine if vertical maneuvers are acceptable for DAA WCR guidance against non-cooperative intruders. This study incorporates the radar from the Honeywell sensor model and examines a series of pairwise encounters between a non-cooperative intruder and a UAS ownship, with different combinations of intruder states and ownship performance levels. The study uses SAA Control as a simulation platform and pilot model, and Omnibands to provide DWC recovery guidance. Two simulation sets, one that allows vertical DWC recovery guidance and one that does not, are compared to determine if encounters with altitude rate errors above 250 feet-per-minute are more likely to have more severe losses of well clear, as determined by the Loss of Well-Clear Severity metric.

SaaControl↗

Entanglement types for two-qubit states with real amplitudes

We study the set of two-qubit pure states with real amplitudes and their geometrical representation in the three-dimensional sphere. In this representation, we show that the maximally entangled states—those locally equivalent to the Bell states—form two disjoint circles perpendicular to each other. We also show that taking the natural Riemannian metric on the sphere, the set of states connected by local gates are equidistant to this pair of circles. Moreover, the unentangled or so-called product states are π/4 units away to the maximally entangled states. This is, the unentangled states are the farthest away to the maximally entangled states. In this way, if we define two states to be equivalent if they are connected by local gates, we have that there are as many equivalent classes as points in the interval [0,π/4] with the point 0 corresponding to the maximally entangled states. The point π/4 corresponds to the unentangled states which geometrically are described by a torus. Finally, for every 0<d<π/4 the point d corresponds to a disjoint pair of torus. Finally, we also show how this geometrical interpretation allows to clearly see that any pair of two-qubit states with real amplitudes can be connected with a circuit that only has single-qubit gates and one controlled-Z gate.

97 MATHEMATICS AND COMPUTING↗

Inflatable Space Structures Technology Development for Large Radar Antennas

There has been recent interest in inflatable space-structures technology for possible applications on U.S. Department of Defense (DOD) missions because of the technology's potential for high mechanical-packaging efficiency, variable stowed geometry, and deployment reliability. In recent years, the DOD sponsored Large Radar Antenna (LRA) Program applied this new technology to a baseline concept: a rigidizable/inflatable (RI) perimeter-truss structure supporting a mesh/net parabolic reflector antenna. The program addressed: (a) truss concept development, (b) regidizable materials concepts assessment, (c) mesh/net concept selection and integration, and (d) developed potential mechanical-system performance estimates. Critical and enabling technologies were validated, most notably the orbital radiation durable regidized materials and the high modulus, inflatable-deployable truss members. These results in conjunction with conclusions from previous mechanical-packaging studies by the U.S. Defense Advanced Research Projects Agency (DARPA) Special Program Office (SPO) were the impetus for the initiation of the DARPA/SPO Innovative Space-based Antenna Technology (ISAT) Program. The sponsor's baseline concept consisted of an inflatable-deployable truss structure for support of a large number of rigid, active radar panels. The program's goal was to determine the risk associated with the application of these new RI structures to the latest in radar technologies. The approach used to define the technology maturity level of critical structural elements was to: (a) develop truss concept baseline configurations (s), (b) assess specific inflatable-rigidizable materials technologies, and (c) estimate potential mechanical performance. The results of the structures portion of the program indicated there was high risk without the essential materials technology flight experiments, but only moderate risk if the appropriate on-orbit demonstrations were performed. This paper covers both programs (LRA and ISAT) in two sections, Parts 1 and 2 respectively. Please note that the terms strut, tube, and column are all used interchangeably and refer to the basic strut element of a truss. Also, the paper contains a mix of English and metric dimensional descriptions that reflect prevailing technical discipline conventions and common usage.

Freeland, R. E.↗

Geophysical monitoring using active seismic techniques at the Citronelle Alabama CO 2 storage demonstration site

Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.

58 GEOSCIENCES↗