Search NASA⌕ Search

SEARCH · Search NASA

Results for “Scaling and squaring method”

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.

393 records · Page 22

U.S. Coast Guard Boat Recovery Simulation at the NASA Ames Vertical Motion Simulator

The U.S. Coast Guard routinely uses the Over-The-Horizon (OTH-IV), a cutter deployed, rigid-hulled inflatable boat for rescue and law enforcement operations. The term "cutter" refers to a Coast Guard vessel 65 feet in length or greater with accommodations for crew to live aboard and the ability to deploy smaller boats including the OTH-IV. A deployment method employed by some cutter class vessels is a davit crane system and specialized hook mechanism to connect the smaller boat. The launch and recovery operation requires several crewmembers, where one is tasked with releasing and connecting the hook to the small boat. Manipulating the hook can pose a risk to the crewmember depending on sea conditions combined with fatigue level or task expertise. Equipment and crewmembers are tested extensively due to the inherent risk of the task, but the sea conditions cannot be controlled. To address the sea condition variable, the U.S. Coast Guard and NASA leveraged the Vertical Motion Simulator (VMS) as a platform to test new equipment or train crewmembers in varying sea conditions. With 6 degrees of freedom and a vertical displacement of +/- 22ft, the VMS is capable of simulating conditions up to sea-state 5. The proof of concept experiment took place in September 2018 and demonstrated that the VMS can accurately simulate varying sea-states, collect performance data, and design a reliable and safe system for participants. The U.S. Coast Guard provided boat displacement data from the Large Amplitude Motion Program (LAMP). The program generates multi-directional waves with a cosine squared spreading function which produced a time history response of displacement data (surge, sway, heave, roll, pitch, yaw) of the OTH-IV's center of gravity (C.G). The displacement data was differentiated into acceleration evaluated at the crewmember's position at the front (bow) of the boat. Those accelerations were transformed to the VMS coordinate system and used as command inputs to the motion system's washout filters. The VMS uses Interchangeable Cabs (ICab) to perform flight simulations on a variety of aircraft or spacecraft configurations. An ICab typically contains seats, restraints, pilot controls, instrument panel, and out-the-window projectors all enclosed in a familiar cockpit configuration. To replicate the launch and recovery of a small boat, all the conventional features of a flight simulator and canopy were removed to be replaced with a to-scale OTH-IV bow mockup. The bow section was constructed from aluminum, and a connection point known as the davit ring. Additionally, a hook replica was suspended by a cable from the ceiling as if it were hanging from a davit crane aboard a cutter vessel. The last step was providing an auxiliary control to close the distance between suspended hook and the davit ring on the bow mockup. Without control the boat would simulate the wave motion, but with no guarantee that the crewmember will be in reach of the hook. Since the wave motion was repeatable, an additional acceleration command sent to the washout filters to move the simulator in the surge, sway, and heave directions, independently. The commands were rate limited such that the additional movement was not detectable with respect to the wave motion. The results showed that for each wave profile, a unique set of acceleration commands provided opportunities for crewmembers to gain control of the hook and connect it to the davit ring.

Riccobono, Nicholas S.↗

Redshift evolution and covariances for joint lensing and clustering studies with DESI Y1

ABSTRACT Galaxy–galaxy lensing (GGL) and clustering measurements from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) data set promise to yield unprecedented combined-probe tests of cosmology and the galaxy–halo connection. In such analyses, it is essential to identify and characterize all relevant statistical and systematic errors. We forecast the covariances of DESI Y1 GGL + clustering measurements and the systematic bias due to redshift evolution in the lens samples. Focusing on the projected clustering and GGL correlations, we compute a Gaussian analytical covariance, using a suite of N-body and lognormal simulations to characterize the effect of the survey footprint. Using the DESI one percent survey data, we measure the evolution of galaxy bias parameters for the DESI luminous red galaxy (LRG) and bright galaxy survey (BGS) samples. We find mild evolution in the LRGs in $0.4 < z < 0.8$, subdominant to the expected statistical errors. For BGS, we find less evolution for brighter absolute magnitude cuts, at the cost of reduced sample size. We find that for a redshift bin width $\Delta z = 0.1$, evolution effects on DESI Y1 GGL is negligible across all scales, all fiducial selection cuts, all fiducial redshift bins. Galaxy clustering is more sensitive to evolution due to the bias squared scaling. Nevertheless the redshift evolution effect is insignificant for clustering above the 1-halo scale of $0.1h^{-1}$ Mpc. For studies that wish to reliably access smaller scales, additional treatment of redshift evolution is likely needed. This study serves as a reference for GGL and clustering studies using the DESI Y1 sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗

A fast and accurate domain decomposition nonlinear manifold reduced order model

Here, this paper integrates nonlinear-manifold reduced order models (NM-ROMs) with domain decomposition (DD). NM ROMs approximate the full order model (FOM) state in a nonlinear-manifold by training a shallow, sparse autoencoder using FOM snapshot data. These NM-ROMs can be advantageous over linear-subspace ROMs (LS-ROMs) for problems with slowly decaying Kolmogorov n-width. However, the number of NM-ROM parameters that need to be trained scales with the size of the FOM. Moreover, for “extreme-scale” problems, the storage of high-dimensional FOM snapshots alone can make ROM training expensive. To alleviate the training cost, this paper applies DD to the FOM, computes NM-ROMs on each subdomain, and couples them to obtain a global NM-ROM. This approach has several advantages: Subdomain NM-ROMs can be trained in parallel, involve fewer parameters to be trained than global NM-ROMs, require smaller subdomain FOM dimensional training data, and can be tailored to subdomain specific features of the FOM. The shallow, sparse architecture of the autoencoder used in each subdomain NM-ROM allows application of hyper-reduction (HR), reducing the complexity caused by nonlinearity and yielding computational speedup of the NM-ROM. This paper provides the first application of NM-ROM (with HR) to a DD problem. In particular, this paper details an algebraic DD reformulation of the FOM, training a NM-ROM with HR for each sub domain, and a sequential quadratic programming (SQP) solver to evaluate the coupled global NM-ROM. Theoretical convergence results for the SQP method and a priori and a posteriori error estimates for the DD NM-ROM with HR are provided. The proposed DD NM-ROM with HR approach is numerically compared to a DD LS-ROM with HR on the 2D steady-state Burgers’ equation, showing an order of magnitude improvement in accuracy of the proposed DD NM-ROM over the DD LS-ROM.

97 MATHEMATICS AND COMPUTING↗

Fabrication Method for LOBSTER-Eye Optics in <110> Silicon

Soft x-ray optics can use narrow slots to direct x-rays into a desirable pattern on a focal plane. While square-pack, square-pore, slumped optics exist for this purpose, they are costly. Silicon (Si) is being examined as a possible low-cost replacement. A fabrication method was developed for narrow slots in <110> Si demonstrating the feasibility of stacked slot optics to replace micropores. Current micropore optics exist that have 20-micron-square pores on 26-micron pitch in glass with a depth of 1 mm and an extent of several square centimeters. Among several proposals to emulate the square pore optics are stacked slot chips with etched vertical slots. When the slots in the stack are positioned orthogonally to each other, the component will approach the soft x-ray focusing observed in the micropore optics. A specific improvement Si provides is that it can have narrower sidewalls between slots to permit greater throughput of x-rays through the optics. In general, Si can have more variation in slot geometry (width, length). Further, the sidewalls can be coated with high-Z materials to enhance reflection and potentially reduce the surface roughness of the reflecting surface. Narrow, close-packed deep slots in <110> Si have been produced using potassium hydroxide (KOH) etching and a patterned silicon nitride (SiN) mask. The achieved slot geometries have sufficient wall smoothness, as observed through scanning electron microscope (SEM) imaging, to enable evaluation of these slot plates as an optical element for soft x-rays. Etches of different angles to the crystal plane of Si were evaluated to identify a specific range of etch angles that will enable low undercut slots in the Si <110> material. These slots with the narrow sidewalls are demonstrated to several hundred microns in depth, and a technical path to 500-micron deep slots in a precision geometry of narrow, closepacked slots is feasible. Although intrinsic stress in ultrathin wall Si is observed, slots with walls approaching 1.5 microns can be achieved (a significant improvement over the 6-micron walls in micro - pore optics). The major advantages of this technique are the potential for higher x-ray throughout (due to narrow slot walls) and lower cost over the existing slumped micropore glass plates. KOH etching of smooth sidewalls has been demonstrated for many applications, suggesting its feasibility for implementation in x-ray optics. Si cannot be slumped like the micropore optics, so the focusing will be achieved with millimeter-scale slot plates that populate a spherical dome. The possibility for large-scale production exists for Si parts that is more difficult to achieve in micropore parts.

Chervenak, James↗

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗

Contaminant Transport Modeling for Technology Evaluation and Long-Term Monitoring in the Tims Branch Testbed, SC - 20343

Studies conducted at U.S. DOE sites have shown the presence of heavy metals, radionuclides, and volatile organic compounds in surface water, groundwater, and soil as a result of nuclear activity in the Cold War era. Since the 1990's, innovative cleanup methods have been implemented in the Tims Branch watershed at Savannah River Site (SRS) to limit the contaminant flux to the stream that have reduced the contaminant concentrations to acceptable regulatory levels in the dissolved phase. A tin-based treatment which effectively eliminated all local anthropogenic mercury inputs to this ecosystem resulted in a known step function addition of inert tin oxide particles which now serve as a potential tracer for sedimentation and particle transport processes in the stream. The long-term effectiveness of this and other remediation techniques and the potential for remobilization of adsorbed contaminant in sediment during extreme hydrologic conditions however remains unclear. It is therefore important to understand not only the fate and transport of dissolved contaminants, but also the movement of sediment and the relevant interactions with dissolved contaminant. To narrow this knowledge gap, a study is being conducted using the Tims Branch watershed as a stream-scale ecosystem test-bed to identify the primary transport processes of major contaminants of concern (such as mercury, nickel and uranium) with an emphasis on interactions with sediment transport. This involves the development of a fully distributed hydrologic watershed model of the Tims Branch watershed to predict streamflow under extreme weather conditions, as well as the development of a comprehensive contaminant transport model that can properly account for coupled contaminant and sediment transport. Review of relevant research reports and peer-reviewed journals revealed that aside from advection-dispersion transport of dissolved contaminants, adsorption and desorption with suspended solids and bed sediment also play an important role in the transport of those contaminants of concern. To develop the fully distributed hydrologic watershed model, the MIKE SHE 2-dimensional (2D) land surface/3D groundwater model that simulates surface/subsurface hydrologic processes (such as overland flow, evapotranspiration, and infiltration) was coupled with a 1D streamflow model that accounts for stream water hydraulics (such as hydraulic structures, cross-sections, and network). To model contaminant transport, the MIKE 11 streamflow component was coupled with the MIKE 11 AD module that simulates solute transport through advection and dispersion, and MIKE ECO Lab module that accounts for both sediment transport and interactions with dissolved contaminant. At this stage, the development and optimization of the fully distributed hydrologic model has been completed, achieving satisfactory statistical results between observed and predicted discharge as indicated by a root mean square error (RMSE) of 0.039 cms and a Nash-Sutcliffe efficiency coefficient (NSE) of 0.764. The ongoing development of the contaminant transport model has also yielded realistic results from preliminary tests. Results from this study are a key to evaluating the effectiveness of tin (II)-based mercury treatment of wetlands at the SRS site, and are also relevant to evaluating the potential of using this type of novel remediation technology in other mercury-contaminated stream systems. Knowledge acquired from this research will also support interpretation of historical data on the trends of contaminant concentration distribution in Tims Branch, particularly considering the effect of extreme hydrological events on the stream flow and pollutant transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Monte Carlo predictions of DNA fragment-size distributions for large sizes after HZE particle irradiation

DSBs (double-strand breaks) produced by densely ionizing space radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. DSB clustering at large scales, from >100 Mbp down to approximately 2 kbp, is modeled using a Monte-Carlo algorithm. A random-walk model of chromatin is combined with a track model, that predicts the radial distribution of energy from an ion, and the RLC (randomly-located-clusters) formalism, in software called DNAbreak. This model generalizes the random-breakage model, whose broken-stick fragment-size distribution is applicable to low-LET radiation. DSB induction due to track interaction with the DNA volume depends on the radiation quality parameter Q. This dose-independent parameter depends only weakly on LET. Multi-track, high-dose effects depend on the cluster intensity parameter lambda, proportional to fluence as defined by the RLC formalism. After lambda is determined by a numerical experiment, the model reduces to one adjustable parameter Q. The best numerical fits to the experimental data, determining Q, are obtained. The knowledge of lambda and Q allows us to give biophysically based extrapolations of high-dose DNA fragment-size data to low doses or to high LETs.

NASA Center JSC↗

Robust In-Situ Strain Measurements to Monitor CO 2 Storage

The goal of this project was to develop and demonstrate robust instrumentation to monitor the in-situ strain tensor in order to improve the reliability and security of CO 2 storage in geologic formations. We met the original goals of the project and the major overarching accomplishment is the advancement of strain tensor monitoring from an intriguing concept to a commercially available technology with a solid foundation of novel instruments supported by theoretical analyses and validation experiments. The main accomplishments of the project are summarized below. We designed, built and evaluated nine new optical fiber strainmeters and tiltmeters using Michelson interferometers to measure deformation with ultra-high resolution at both shallow and deep point locations in the subsurface. These are the robust strainmeters that motivated the title of the project. We designed, built and evaluated a novel method of measuring distributed strain in optical fibers with nanostrain resolution, and cm-scale location, and sampling into the seismic band. The new method is called Coherence-length-gated Microwave Photonics Interfereometry (CMPI). CMPI technology has advantages over existing commercial DAS and DSS methods. We developed and demonstrated capabilities to deploy instruments in the field and used them to measure strain caused by ambient signals like barometric pressure and tides, as well as induced signals like surface loading and pore pressure changes from pumping tests. We deployed a working strainmeter at 1,700 ft depth, slightly above an active reservoir. This is to our knowledge the greatest depth a strainmeter has been deployed and the techniques we used can readily be extended to greater depths. Optical fiber borehole tensor strainmeter techology was advanced from a TRL 4 at the start, to a TRL of 7 at the conclusion of the project. The project included advances in simulations and theoretical analyses. We developed and demonstrated a computational workflow that uses machine learning to reduce the computational requirements and make it practical to use Bayesian inversion to solve large numerical poroelastic analyses needed to interpret strain tensor field data. We evaluated the strain tensor fields and time series that would be caused by leaks of CO 2 or other fluids from reservoirs. These simulations demonstrated that signals from leaks could be measured with instruments developed for the project, opening a potentially new method for ensuring storage security. We showed that strains in caprock can be used to estimate pressure in a reservoir. This avoids the need to drill monitoring wells into the reservoir, and it expands the capabilities of monitoring in the caprock. The project includes a derivation and application of a novel analytical solution to the strains in the vicinity of a pressurized poroelastic inclusion. This solution explains field data measured during injeciton tests at the North Avant Field, and it will simplify future interpretation of strain tensor data. The project included a broad range of experiments, and of the most significant is the characterization of the strain tensor at an array three strainmeters during six injection tests at the North Avant Field, Oklahoma. This demonstrated repeatability of the strain signal measured by the new instruments developed for the project, and it showed similarities between the strain signal at shallow depths and pressure in the underlying reservoir. We also demonstrated that useful strain data can be measured at reservoir depths. This confirms that strain tensor data can be measured throughout the caprock over a reservoir. The project demonstrated the feasibility of using the strain tensor and distributed strain measured in caprock during a variety of different well tests where the pumping rate was constant, sinusoidal and positive, or a periodic square wave with zero net rate. This further strengthens the validity of using strain data to characterize reservoirs and aquifers. We also demonstrated that strain caused be fluctuations of air pressure and water pressure in the vadose zone can be measured and interpreted, suggesting that high resolution distributed strain measurements hold promise for monitoring the vadose zone. The project partially supported nine graduate students in the Environmental Engineering, Hydrogeology, Electrical Engineering programs at Clemson University. The research was described in nine journal papers, 23 talks and conference abstracts. Additional journal papers are in preparation. A new company called Tensora was started to provide strainmeter technology for commercial applications.

01 COAL, LIGNITE, AND PEAT↗

Investigation of Spiral Bevel Gear Condition Indicator Validation Via AC-29-2C Using Damage Progression Tests

This report documents the results of spiral bevel gear rig tests performed under a NASA Space Act Agreement with the Federal Aviation Administration (FAA) to support validation and demonstration of rotorcraft Health and Usage Monitoring Systems (HUMS) for maintenance credits via FAA Advisory Circular (AC) 29-2C, Section MG-15, Airworthiness Approval of Rotorcraft (HUMS) (Ref. 1). The overarching goal of this work was to determine a method to validate condition indicators in the lab that better represent their response to faults in the field. Using existing in-service helicopter HUMS flight data from faulted spiral bevel gears as a "Case Study," to better understand the differences between both systems, and the availability of the NASA Glenn Spiral Bevel Gear Fatigue Rig, a plan was put in place to design, fabricate and test comparable gear sets with comparable failure modes within the constraints of the test rig. The research objectives of the rig tests were to evaluate the capability of detecting gear surface pitting fatigue and other generated failure modes on spiral bevel gear teeth using gear condition indicators currently used in fielded HUMS. Nineteen final design gear sets were tested. Tables were generated for each test, summarizing the failure modes observed on the gear teeth for each test during each inspection interval and color coded based on damage mode per inspection photos. Gear condition indicators (CI) Figure of Merit 4 (FM4), Root Mean Square (RMS), +/- 1 Sideband Index (SI1) and +/- 3 Sideband Index (SI3) were plotted along with rig operational parameters. Statistical tables of the means and standard deviations were calculated within inspection intervals for each CI. As testing progressed, it became clear that certain condition indicators were more sensitive to a specific component and failure mode. These tests were clustered together for further analysis. Maintenance actions during testing were also documented. Correlation coefficients were calculated between each CI, component, damage state and torque. Results found test rig and gear design, type of fault and data acquisition can affect CI performance. Results found FM4, SI1 and SI3 can be used to detect macro pitting on two more gear or pinion teeth as long as it is detected prior to progressing to other components or transitioning to another failure mode. The sensitivity of RMS to system and operational conditions limit its reliability for systems that are not maintained at steady state. Failure modes that occurred due to scuffing or fretting were challenging to detect with current gear diagnostic tools, since the damage is distributed across all the gear and pinion teeth, smearing the impacting signatures typically used to differentiate between a healthy and damaged tooth contact. This is one of three final reports published on the results of this project. In the second report, damage modes experienced in the field will be mapped to the failure modes created in the test rig. The helicopter CI data will then be re-processed with the same analysis techniques applied to spiral bevel rig test data. In the third report, results from the rig and helicopter data analysis will be correlated. Observations, findings and lessons learned using sub-scale rig failure progression tests to validate helicopter gear condition indicators will be presented.

gears↗

CODEX clusters

Context: Large area catalogs of galaxy clusters constructed from ROSAT All-Sky Survey provide the basis for our knowledge of the population of clusters thanks to long-term multiwavelength efforts to follow up observations of these clusters. Aims: The advent of large area photometric surveys superseding previous, in-depth all-sky data allows us to revisit the construction of X-ray cluster catalogs, extending the study to lower cluster masses and higher redshifts and providing modeling of the selection function. Methods: We performed a wavelet detection of X-ray sources and made extensive simulations of the detection of clusters in the RASS data. We assigned an optical richness to each of the 24 788 detected X-ray sources in the 10 382 square degrees of the Baryon Oscillation Spectroscopic Survey area using red sequence cluster finder redMaPPer version 5.2 run on Sloan Digital Sky Survey photometry. We named this survey COnstrain Dark Energy with X-ray (CODEX) clusters. Results: We show that there is no obvious separation of sources on galaxy clusters and active galactic nuclei (AGN) based on the distribution of systems on their richness. This is a combination of an increasing number of galaxy groups and their selection via the identification of X-ray sources either by chance or by groups hosting an AGN. To clean the sample, we use a cut on the optical richness at the level corresponding to the 10% completeness of the survey and include it in the modeling of the cluster selection function. We present the X-ray catalog extending to a redshift of 0.6. Conclusions: The CODEX suvey is the first large area X-ray selected catalog of northern clusters reaching fluxes of 10 –13 ergs s –1 cm –2 . We provide modeling of the sample selection and discuss the redshift evolution of the high end of the X-ray luminosity function (XLF). Our results on z < 0.3 XLF agree with previous studies, while we provide new constraints on the 0.3 < z < 0.6 XLF. We find a lack of strong redshift evolution of the XLF, provide exact modeling of the effect of low number statistics and AGN contamination, and present the resulting constraints on the flat ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗

Roll-to-Roll Advanced Materials Manufacturing DOE Laboratory Collaboration (FY2020 Final Report)

R2R processing is used to manufacture a wide range of products for various applications which span many industrial business sectors. The overall R2R methodology has been in use for decades and this continuous technique traditionally involves deposition of material(s) onto moving webs, carriers or other continuous belt-fed or conveyor-based processes that enable successive steps to build a final version which serves to support the deposited materials. Established methods that typify R2R processing include tape casting, silk-screen printing, reel-to-reel vacuum deposition/coating, and R2R lithography. Products supported by R2R manufacturing include micro-electronics, electro-chromic window films, PVs, fuel cells for energy conversion, battery electrodes for energy storage, and barrier and membrane materials. Due to innovation in materials and process equipment, high-quality yet very low-cost multilayer technologies have the potential to be manufactured on a very cost-competitive basis. To move energy-related products from high-cost niche applications to the commercial sector, the means must be available to enable manufacture of these products in a cost-competitive manner that is affordable. Fortunately, products such as fuel cells, thin- and mid-film PVs, batteries, electrochromic and piezoelectric films, water separation membranes, and other energy saving technologies readily lend themselves to manufacture using R2R approaches. However, more early-stage research is needed to solve the challenge of linking the materials (particles, polymers, solvents, additives) used in ink and slurry formulations and the coating and drying processes to the ultimate performance of the final R2R product, especially for a process that uses multiple layers of deposition to achieve the end product. To solve the problems associated with these challenges, the R2R Collaboration is executing a research program with outcomes that will ultimately link modeling, processing, metrology and defect detection tools, thereby directly relating the properties of constituent particles and processing conditions to the performance of final devices. This collaborative approach was designed to foster identification and development of materials and processes related to R2R for clean-energy materials development. Using computational and experimental capabilities by acknowledged subject matter experts within the supported National Laboratory system, this project leverages the capabilities and expertise at each of five National Laboratories to further the development of multilayer technologies that will enable high-volume, cost-competitive platforms. A typical R2R process has three steps: (1) mixing of particles and various constituents in a slurry, (2) coating of the ink/slurry mixture on a substrate, and (3) drying/curing and processing of the coating. Final performance of devices made via R2R processes is dependent on the active materials (e.g., electrochemical particles in battery or fuel cell electrodes) and the device structure that stems from the governing component interactions within the various steps. However, a fundamental understanding of the underlying mechanisms and phenomena is still lacking, which is why industrial-scale R2R process development and manufacturing is still largely empirical in nature. The FY 2019 through FY 2021 program addresses aspects of the following two targets from the AMO Multi-Year Program Plan: (1) Target 8.1 Develop technologies to reduce the cost per manufactured throughput of continuous R2R manufacturing processes. (A) Increasing throughput of R2R processes by 5 times for batteries (to 50 square feet per minute (50 ft 2 /min)) and capacitors and 10 times for printed electronics and the manufacture of other substrates and MEs used in support of these products. (B) Developing resolution capabilities to enable registration and alignment that will detect, align, and co-deposit multiple layers of coatings and print < 1-micron (1 µm) features using continuous process scalable for commercial production. (C) Developing scalable and reliable R2R processes for solution deposition of ultra-thin (<10 nm) films for active and passive materials. (D) Develop in-line multilayer coating technology on thin films with yields greater than 95%. (2) Target 8.2 Develop in-line instrumentation tools that will evaluate the quality of single and multilayer materials in-process. (A) Developing in-line QC technologies and methodologies for real-time identification of defects and expected product properties “in-use/application” during continuous processing at all size-scales with a focus on the “micro” and “nano” scale traces, lines, and devices, i.e., <1 μm at 300 ft./min for R2R processing in air and <10 nm at 20 ft./min for vacuum (B) Developing technologies to increase the measurement frequency of surface rheology without significant cost increases with a goal of a 10-nanometer in-line profilometry at a production rate of 100,000 square millimeters per minute (100,000 mm 2 /min).

42 ENGINEERING↗

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗