SEARCH · Engineering Papers
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.
Shape preserving approximations are constructed by interpolating the data with polynomial splines of arbitrary degree. A regularity condition is formulated on the data which insures the existence of such a shape preserving spline, an algorithm is presented for its construction, and the uniform norm of the error is bound which results when the algorithm is used to produce an approximation to a given f epsilon Ca,b.
Molecular biomarkers are the most direct biosignatures of life on early Earth and a key target in the search for life on Mars. Lipid biomarkers are of particular interest given their ability to survive oxidative degradation and record microbial presence and activity of microorganisms that occurred billions of years ago (Eigenbrode, 2008). Environmental conditions that suspend biotic and abiotic degradative processes prior to lithification can lead to enhanced biomolecular preservation over geological time-scales. The hyperarid core of the Atacama Desert in northern Chile offers a unique environment to investigate lipid biomarker taphonomy under extreme and prolonged dryness. We investigated the accumulation and degree of preservation of lipid biomarkers in million-year-old hyperarid soils where primarily abiotic conditions influence their taphonomy. Soils were extracted and free and membrane bound lipids were analyzed across a vertical profile of 2.5 meters in the Yungay hyper-arid core of the Atacama Desert. Due to the extremely low inventory of biomass in Atacama soils, samples were collected by scientists wearing cleanroom suits to minimize anthropogenic contamination during sampling. Fatty acids were found to be well preserved in Yungay soils, and were most abundant in the clay-rich soils at approx.2 m depth (approx.750 ng of fatty acid methyl ester/g of soil). These buried clays layers were fluvially deposited approximately 2 million years ago, and have been excluded from exposure to rainwater and modern surficial processes since their emplacement (Ewing et al., 2008). Monocarboxylic fatty acid, monohydroxy fatty acid, glycerol tetraether, and n-alkane hydrocarbon content was found to change with depth. Lipid biomarker content in deeper soil layers is suggestive of soils having been formed at a time when environmental conditions were capable of supporting active microbial communities and plants. In short, total lipid extracts reveal a remarkable degree of lipid biomarker preservation even in the oldest soils analyzed (ca. 2 Myr) indicating that typical diagenetic processes of lipid destruction are arrested under extreme dryness. This result has implications for the search for molecular biomarkers on Mars, which could have experienced millions to billions of years of extreme hyperaridity.
Here a fully discrete surface finite element method is proposed for solving the viscous shallow water equations in a bounded Lipschitz domain on the sphere based on a general triangular mesh. The method consists of a modified Crank–Nicolson method in time and a Galerkin surface finite element method in space for the fluid thickness H and the fluid velocity u. A finite element space tangential to the sphere at all finite element nodes is proposed to approximate the fluid velocity u. The proposed method has second-order accuracy in time and first-order accuracy in space, and preserves mass conservation, well balancedness, tangency of velocity to the sphere, and energy decay. Numerical experiments are presented to illustrate the accuracy of the proposed method and the preservation of the physical properties, including mass conservation, well balancedness, and energy decay. A numerical simulation of ocean mesoscale activity on a circular basin with a continental shelf is provided.
Lossless compressors have very low compression ratios that do not meet the needs of today's large-scale scientific applications that produce vast volumes of data. Error-bounded lossy compression (EBLC) is considered a critical technique for the success of scientific research. Although EBLC allows users to set an error bound for the compression, users have been unable to specify the requirements on the compression quality, limiting practical use. Our contributions are: (1) We formulate the problem of configuring EBLC to preserve a user-defined metric as an optimization problem. This allows many classes of new metrics to be preserved, which improves over current practices. (2) We present a framework, OptZConfig, that can adapt to improvements in the search algorithm, compressor, and metrics with minimal changes, enabling future advancements in this area. (3) We demonstrate the advantages of our approach against the leading methods to configure compressors to preserve specific metrics. Here, our approach improves compression ratios against a specialized compressor by up to 3 x, has a 56x speedup over FRaZ, 1000x speedup over MGARD-QOI post tuning, and 110x speedup over systematic approaches which had not been bounded by compressors before.
Conjugate gradient type methods are considered for the solution of large linear systems Ax = b with complex coefficient matrices of the type A = T + i(sigma)I where T is Hermitian and sigma, a real scalar. Three different conjugate gradient type approaches with iterates defined by a minimal residual property, a Galerkin type condition, and an Euclidian error minimization, respectively, are investigated. In particular, numerically stable implementations based on the ideas behind Paige and Saunder's SYMMLQ and MINRES for real symmetric matrices are proposed. Error bounds for all three methods are derived. It is shown how the special shift structure of A can be preserved by using polynomial preconditioning. Results on the optimal choice of the polynomial preconditioner are given. Also, some numerical experiments for matrices arising from finite difference approximations to the complex Helmholtz equation are reported.
The early diagenetic interplay between reactive iron, sulfur, and organic matter in the bathymetrically isolated Santa Monica Basin (SMB) sediments are investigated in this study. Here, we explore solid-phase and porewater profiles from the basin, supplemented with a transect from 71 to 907 m water depth that includes oxygenated (>60 μM O 2 ) bottom waters near the coast and oxygen-deficient waters (~4 µM O 2 ) in the basin. The geochemical data of the basin sediments are further scrutinized by means of reactive transport modeling. The results show that the basin sediments do not follow the traditional geochemical signatures of oxygen-deficient settings. A lack of dissolved sulfide accumulation and sulfurized iron persists despite the sediments being deposited under reducing conditions (without bioturbation/bioirrigation), strong organic carbon input (TOC up to 5.0 wt%), and active dissimilatory sulfate reduction. Not only did we find an exceptional enrichment in highly reactive Fe in the surface sediments (~45 % of total Fe), but the enrichment of reactive Fe, including ferrihydrite, persists downcore and coexists with high levels of dissolved Fe. The enhanced preservation of Fe oxides and lack of iron-sulfide precipitation is in part explained by detection via Mössbauer spectra of iron oxides bounded to organic matter (Fe[III]-OM coprecipitates). The modeled Fe budget shows that most of the Fe oxides in the surface sediments are internally recycled by upward diffusion and subsequent oxidation of Fe 2+ . Sulfide oxidation coupled to Fe reduction effectively precludes sulfide accumulation while enhancing build-up of dissolved Fe, fueling the Fe cycle within the first 5 cm depth. Continuous reoxidation of Fe 2+ enhances the formation of Fe(III)-OM coprecipitates, limiting the amount of reactive organic matter. In the unavailability of labile organic matter, other than within the uppermost layers, the organic-rich sediment profiles are dominated by Fe cycling that limits the production and preservation of sulfides and enhances the preservation of Fe oxides and organic carbon. Finally, this study highlights key local controls on Fe availability in marginal basins and describes an intricate biogeochemical C-Fe-S cycling in modern and possibly ancient marine systems with important implications for Fe availability in the marine realm.
Time-domain modeling of the thermal deformation of crystal optics can help define acceptable operational ranges across the pulse-energy repetition-rate phase space. In this paper, we have studied the transient thermal deformation of a water-cooled diamond crystal for a cavity-based X-ray free-electron laser (CBXFEL), either an X-ray free-electron laser oscillator (XFELO) or a regenerative amplifier X-ray free-electron laser (RAFEL), by numerical simulations including finite-element analysis and advanced data processing. Pulse-by-pulse transient thermal deformation of a 50 µm-thick diamond crystal has been performed with X-ray pulse repetition rates between 50 kHz and 1 MHz. Results for temperature and thermal deformation have been compared with the results of transient analysis using a continuous wave (CW) power loading. Temperature and thermal deformation results from pulse-by-pulse transient analysis vary with time about the results for the CW case for the same average power. The variation amplitude increases with pulse energy and decreases with repetition rate. When the repetition rate increases to infinity, both temperature and thermal deformation converge to the results for the CW case. Two critical time scales for the operation of crystal optics in a CBXFEL are (1) first-turn time, i.e. the time for the XFEL pulse to complete the first turn around the cavity so that the crystal sees the recirculated XFEL pulse, and (2) period-end time, i.e. the time that the next electron bunch arrives for the amplification, so that the crystal outcouples the amplified FEL power. For the same average power, simulation results show that the crystal thermal deformation seen by the XFEL beam decreases with repetition rate at the first-turn time of a 300 m-long cavity and increases with repetition rate at the period-end time. For the wavefront preservation requirement of the crystal optics, a pulse-energy versus repetition-rate phase space has been established. The upper bounds of the pulse energy at both first-turn and period-end times decreases with repetition rate, especially at the period-end time. The upper bound of the thermal deformation of the crystal at the period-end time for any repetition frequency can be estimated from the CW case. For a water-cooled diamond crystal of dimension 5 mm × 5 mm × 0.05 mm, the time to reach a quasi steady-state is about 50 ms for temperature and 50 µs for thermal deformation.
An adaptive method in which cruise and nonlinear orbit determination problems can be solved using a single program is presented. It involves singular value decomposition augmented with an extended partial step algorithm. The extended partial step algorithm constrains the size of the correction to the spacecraft state and other solve-for parameters. The correction is controlled by an a priori covariance and a user-supplied bounds parameter. The extended partial step method is an extension of the update portion of the singular value decomposition algorithm. It thus preserves the numerical stability of the singular value decomposition method, while extending the region over which it converges. In linear cases, this method reduces to the singular value decomposition algorithm with the full rank solution. Two examples are presented to illustrate the method's utility.
Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.
The Distance Conjecture holds that any infinite-distance limit in the scalar field moduli space of a consistent theory of quantum gravity must be accompanied by a tower of light particles whose masses scale exponentially with proper field distance ||$\phi$|| as m ~ exp(–λ||$\phi$||), where λ is order-one in Planck units. While the evidence for this conjecture is formidable, there is at present no consensus on which values of λ are allowed. In this paper, we propose a sharp lower bound for the lightest tower in a given infinite-distance limit in d dimensions: λ ≥ $1/\sqrt{d-2}$. In support of this proposal, we show that (1) it is exactly preserved under dimensional reduction, (2) it is saturated in many examples of string/M-theory compactifications, including maximal supergravity in d = 4 – 10 dimensions, and (3) it is saturated in many examples of minimal supergravity in d = 4 – 10 dimensions, assuming appropriate versions of the Weak Gravity Conjecture. We argue that towers with λ < $1/\sqrt{d-2}$ discussed previously in the literature are always accompanied by even lighter towers with λ ≥ $1/\sqrt{d-2}$, thereby satisfying our proposed bound. We discuss connections with and implications for the Emergent String Conjecture, the Scalar Weak Gravity Conjecture, the Repulsive Force Conjecture, large-field inflation, and scalar field potentials in quantum gravity. In particular, we argue that if our proposed bound applies beyond massless moduli spaces to scalar fields with potentials, then accelerated cosmological expansion cannot occur in asymptotic regimes of scalar field space in quantum gravity.
Here, we study the problem of designing a distributed observer for an LTI system over a time-varying communication graph. The limited existing work on this topic imposes various restrictions either on the observation model or on the sequence of communication graphs. In contrast, we propose a single-time-scale distributed observer that works under mild assumptions. Specifically, our communication model only requires strong-connectivity to be preserved over non-overlapping, contiguous intervals that are even allowed to grow unbounded over time. We show that under suitable conditions that bound the growth of such intervals, joint observability is sufficient to track the state of any discrete-time LTI system exponentially fast, at any desired rate. We also develop a variant of our algorithm that is provably robust to worst-case adversarial attacks, provided the sequence of graphs is sufficiently connected over time. The key to our approach is the notion of a "freshness-index" that keeps track of the age-of-information being diffused across the network. Such indices enable nodes to reject stale estimates of the state, and, in turn, contribute to stability of the error dynamics.
Following on our previous work studying the orbits of quantum states under Clifford circuits via reachability graphs, we introduce contracted graphs whose vertices represent classes of quantum states with the same entropy vector. These contracted graphs represent the double cosets of the Clifford group, where the left cosets are built from the stabilizer subgroup of the starting state and the right cosets are built from the entropy-preserving operators. We study contracted graphs for stabilizer states, as well as 𝑊 states and Dicke states, discussing how the diameter of a state's contracted graph constrains the entropic diversity of its two-qubit Clifford orbit. We derive an upper bound on the number of entropy vectors that can be generated using any 𝑛-qubit Clifford circuit, for any quantum state. Here, we speculate on the holographic implications for the relative proximity of gravitational duals of states within the same Clifford orbit. Although we concentrate on how entropy evolves under the Clifford group, our double-coset formalism, and thus the contracted graph picture, is extendable to generic gate sets and generic state properties.
Abstract Kagome lattice hosts a plethora of quantum states arising from the interplay of topology, spin-orbit coupling, and electron correlations. Here, we report symmetry-breaking electronic orders tunable by an applied magnetic field in a model Kagome magnet FeSn consisting of alternating stacks of two-dimensional Fe 3 Sn Kagome and Sn 2 honeycomb layers. On the Fe 3 Sn layer terminated FeSn thin films epitaxially grown on SrTiO 3 (111) substrates, we observe trimerization of the Kagome lattice using scanning tunneling microscopy/spectroscopy, breaking its six-fold rotational symmetry while preserving the translational symmetry. Such a trimerized Kagome lattice shows an energy-dependent contrast reversal in dI/dV maps, which is significantly enhanced by bound states induced by Sn vacancy defects. This trimerized Kagome lattice also exhibits stripe modulations that are energy-dependent and tunable by an applied in-plane magnetic field, indicating symmetry-breaking nematicity from the entangled magnetic and charge degrees of freedom in antiferromagnet FeSn.
Clay minerals detected with orbital and in situ instruments in ancient Martian terrains constrain Mars' climate and aqueous alteration history. Early in its history, Mars experienced an atmospheric redox change and iron-bearing clay minerals may preserve the effects of that transition. Ferrous smectites, the thermodynamically predicted product of chemical weathering of basalts under anoxic conditions, may have undergone oxidation by exposure to chemical oxidants in the atmosphere or regolith, or by direct photooxidation at the surface. To assess these potential oxidation pathways, ferrous trioctahedral smectites of varying initial iron content were synthesized and subjected to oxidation by ultraviolet (UV) irradiation. Experimental UV irradiation under an anoxic atmosphere equivalent to approximately 7 years of flux on the Martian surface caused partial oxidation of smectite-bound Fe (Fe 3+ /ΣFe = 16–18%) and octahedral sheet contraction. Metal-OH vibrational bands in visible/near infrared (VNIR) reflectance spectra of oxidized smectites changed in band depth and asymmetry with higher iron content. X-ray diffraction patterns of UV irradiated samples indicate the formation of a mixed di- and trioctahedral smectite or a secondary nontronite phase, possibly on the surfaces of higher iron content smectites. Furthermore, these experiments suggest that UV irradiation is able to oxidize structurally bound iron in smectites without the presence of other chemical oxidants. Photooxidation may have influenced the mineralogy, both syndepositionally and postdepositionally, of Martian alteration assemblages formed near the surface and this process needs not be limited to one part of their formation history.
Clay minerals detected with orbital and in situ instruments in ancient Martian terrains constrain Mars' climate and aqueous alteration history. Early in its history, Mars experienced an atmospheric redox change and iron-bearing clay minerals may preserve the effects of that transition. Ferrous smectites, the thermodynamically predicted product of chemical weathering of basalts under anoxic conditions, may have undergone oxidation by exposure to chemical oxidants in the atmosphere or regolith, or by direct photooxidation at the surface. To assess these potential oxidation pathways, ferrous trioctahedral smectites of varying initial iron content were synthesized and subjected to oxidation by ultraviolet (UV) irradiation. Experimental UV irradiation under an anoxic atmosphere equivalent to approximately 7 years of flux on the Martian surface caused partial oxidation of smectite-bound Fe (Fe3+/ΣFe = 16–18%) and octahedral sheet contraction. Metal-OH vibrational bands in visible/near infrared (VNIR) reflectance spectra of oxidized smectites changed in band depth and asymmetry with higher iron content. X-ray diffraction patterns of UV irradiated samples indicate the formation of a mixed di- and trioctahedral smectite or a secondary nontronite phase, possibly on the surfaces of higher iron content smectites. These experiments suggest that UV irradiation is able to oxidize structurally bound iron in smectites without the presence of other chemical oxidants. Photooxidation may have influenced the mineralogy, both syndepositionally and postdepositionally, of Martian alteration assemblages formed near the surface and this process needs not be limited to one part of their formation history. Plain Language Summary: Martian mineral assemblages observed with orbiters and rover instruments allow us to understand the planet's past climate and aqueous activity. Hydrated minerals that contain Fe(II) are of particular interest as they would have been affected by the early atmosphere becoming more oxidizing overtime, changing the structural iron into Fe(III). In this study, we synthesized Fe(II) smectites with varying iron content and irradiated them with an ultraviolet (UV) light source to test whether they could be photooxidized. Iron in the smectite minerals was incompletely oxidized by UV radiation after being exposed for the equivalent of approximately 7 years of exposure on the early Mars surface. The smectites with higher iron content showed more changes with photooxidation in their visible/near infrared spectra as well as in their X-ray diffraction patterns, two datasets collected by instruments on currently active Martian missions. These results suggest that UV photooxidation is a plausible contributor to the mineralogy and redox state of clay mineral assemblages observed at the Martian surface.
Nongray gas radiation in a plane parallel slab bounded by gray, diffusely reflecting walls is studied using the discrete ordinates method. The spectral equation of transfer is averaged over a narrow wavenumber interval preserving the spectral correlation effect. The governing equations are derived by considering the history of multiple reflections between two reflecting wails. A closure approximation is applied so that only a finite number of reflections have to be explicitly included. The closure solutions express the physics of the problem to a very high degree and show relatively little error. Numerical solutions are obtained by applying a statistical narrow-band model for gas properties and a discrete ordinates code. The net radiative wail heat fluxes and the radiative source distributions are obtained for different temperature profiles. A zeroth-degree formulation, where no wall reflection is handled explicitly, is sufficient to predict the radiative transfer accurately for most cases considered, when compared with increasingly accurate solutions based on explicitly tracing a larger number of wail reflections without any closure approximation applied.
A series of self-semi-interpenetrating network (ssIPN) thin films based on PIM-1 structure were developed by end-cross-linking telechelic PIM-1 oligomers end-capped with curable carboxylic acid groups to form model networks, which are penetrated by linear high-molecular-weight PIM-1 chains. ssIPN films with systematically varied network content ranging from 10 to 30 wt % were comprehensively examined on their microstructure and gas permeation properties. Fresh PIM-1 ssIPN films exhibited gas separation performances close to those of as-cast linear PIM-1, where the films closely followed the upper bound trade-off line, gaining as much selectivity as they lose gas permeability as the network content increases. This indicates that gas permeability is largely preserved despite cross-linking. Wide-angle X-ray scattering supported this observation, with peaks shifting toward lower d-spacing as oligomer content increased, suggesting tighter chain packing with higher cross-linkable oligomer loading. In physical aging studies over two months, all ssIPN films outperform linear PIM-1 for H 2 /CH 4 and O 2 /N 2 separations. Notably, the PIM-1 ssIPN with 10 wt % network content showed the best physical aging resistance, with negligible permeability loss even after one month. Finally, these findings highlight a promising macromolecular strategy for enhancing the physical aging resistance of microporous polymer membranes for gas separation.
Data compression is becoming critical for data-intensive scientific applications. Scientists require compression techniques that accurately preserve derived quantities of interest (QoIs). Prior work has shown that a pipeline can be built to guarantee error on the primary data (PD) within user-defined bounds and achieve near-floating point QoI errors. In this paper, we present novel computational approaches for accelerating the pipeline and demonstrate results that enable concurrent execution of compression in parallel with the simulation nodes. This allows compression, including the writing of the required compression data, for the previous time step to be completed while the simulation proceeds with the current time step. Overall, the approach presented in this paper results in a 6–8 times improvement in computational overhead compared to previous work. These results were obtained using data generated by a large-scale fusion code called XGC, which produces hundreds of terabytes of data in a single day.