Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean current vector”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

MMS Extended Mission Design: Evaluation of a Lunar Gravity Assist Option

This paper will describe a study that was carried out on the design of a set of maneuvers that were considered for the later stages of an extended mission of the Magnetospheric Multiscale (MMS) mission. The goal of these maneuvers was to put MMS into a significantly different orbit from those flown heretofore, so allowing science collection in a different region of the magnetosphere. This study was made feasible by the fact that the rate at which fuel is being consumed to maintain small formations on the MMS high-apogee orbit is less than expected pre-flight: the current consumption rate is only about 2 kg/yr/spacecraft. In addition, the spacecraft finished the prime mission with a significant amount of fuel remaining: this was about 1-sigma above the mean when compared with pre-launch Monte Carlo simulations. The resulting situation is similar to that of a libration orbit mission, where station-keeping requires so little fuel that any margin at all will lead to an extensive mission lifetime. In the case of MMS, the spacecraft could, if desired, perform formation flying in the current orbit for several decades. Alternatively, the spacecraft could use a significant fraction of the remaining fuel to perform major orbit modifications, while still leaving enough to conduct formation flying for on the order of a decade. The extended mission maneuvers studied here are further apogee-raises, with the goal of setting up one or more lunar gravity assists. Geometry dictates that a lunar encounter is only achievable when the MMS apogee vector lies approximately in the lunar orbit plane: this limits the possible dates to mid-2021 or early 2027.

Design↗
DOE OSTI2022

Pseudogauge dependence of the spin polarization and of the axial vortical effect

The mean spin polarization vector of spin 1/2 particles in a relativistic fluid at local thermal equilibrium (LTE) in different pseudogauges (PGs), i.e., with different choices for the decomposition of orbital and spin angular momentum, is obtained. The spin polarization obtained in the canonical PG differs from the one obtained in the Belinfante PG. It is found that this difference can not be written by replacing the thermal vorticity with the spin potential in the usual polarization formula. Other PG choices affect the contribution of thermal shear to the spin polarization. So, the choice of a PG affects the predictions for the polarization measured in heavy-ion collisions. In general, it is shown that the Wigner function of a noninteracting Dirac field at LTE is affected by the PG transformations. Explicit expressions of the mean axial current are also calculated and it is found that even the axial vortical effect conductivity depends on the PG.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2001

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗
DOE OSTI2020

Incommensurate Spin Fluctuations in the Spin-Triplet Superconductor Candidate UTe 2

Spin-triplet superconductors are of extensive current interest because they can host topological state and Majorana fermions important for quantum computation. The uranium-based heavy-fermion superconductor UTe 2 has been argued as a spin-triplet superconductor similar to UGe2, URhGe, and UCoGe, where the superconducting phase is near (or coexists with) a ferromagnetic (FM) instability and spin-triplet electron pairing is driven by FM spin fluctuations. Here we use neutron scattering to show that, although UTe 2 exhibits no static magnetic order down to 0.3 K, its magnetism in the [0,K,L] plane is dominated by incommensurate spin fluctuations near an antiferromagnetic ordering wave vector and extends to at least 2.6 meV. We are able to understand the dominant incommensurate spin fluctuations of UTe 2 in terms of its electronic structure calculated using a combined density-functional and dynamic mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

CSM research: Methods and application studies

Computational mechanics is that discipline of applied science and engineering devoted to the study of physical phenomena by means of computational methods based on mathematical modeling and simulation, utilizing digital computers. The discipline combines theoretical and applied mechanics, approximation theory, numerical analysis, and computer science. Computational mechanics has had a major impact on engineering analysis and design. When applied to structural mechanics, the discipline is referred to herein as computational structural mechanics. Complex structures being considered by NASA for the 1990's include composite primary aircraft structures and the space station. These structures will be much more difficult to analyze than today's structures and necessitate a major upgrade in computerized structural analysis technology. NASA has initiated a research activity in structural analysis called Computational Structural Mechanics (CSM). The broad objective of the CSM activity is to develop advanced structural analysis technology that will exploit modern and emerging computers, such as those with vector and/or parallel processing capabilities. Here, the current research directions for the Methods and Application Studies Team of the Langley CSM activity are described.

Knight, Norman F., Jr.↗
DOE OSTI2024

Magnetic anisotropy in single-crystalline antiferromagnetic Mn 2 Au

Multiple recent studies have identified the metallic antiferromagnet Mn 2 ⁢Au to be a candidate for spintronic applications due to apparent in-plane anisotropy, preserved magnetic properties above room temperature, and current-induced Néel vector switching. Crystal growth is complicated by the fact that Mn 2 ⁢Au melts incongruently. We present a bismuth flux method to grow millimeter-scale bulk single crystals of Mn 2 ⁢Au in order to examine the intrinsic anisotropic electrical and magnetic properties. Flux quenching experiments reveal that the Mn 2 ⁢Au crystals precipitate below 550°⁢C, about 100⁢°⁢C below the decomposition temperature of Mn 2 ⁢Au. Bulk Mn 2 ⁢Au crystals have a room-temperature resistivity of 16–19 µ⁢Ωcm and a residual resistivity ratio of 41. Mn 2 ⁢Au crystals have a dimensionless susceptibility on the order of 10 –4 (SI units), comparable to calculated and experimental reports on powder samples. Single-crystal neutron diffraction confirms the in-plane magnetic structure. The tetragonal symmetry of Mn 2 ⁢Au constrains the ab-plane magnetic susceptibility to be constant, meaning that χ 100 =χ 110 in the low-field limit, below any spin-flop transition. We find that three measured magnetic susceptibilities χ 100 , χ 110 , and χ 001 are the same order of magnitude and agree with the calculated prediction, meaning the low-field susceptibility of Mn 2 ⁢Au is quite isotropic, despite clear differences in ab-plane and ac-plane magnetocrystalline anisotropy. Mn 2 ⁢Au is calculated to have an extremely high in-plane spin-flop field above 30 T, which is much larger than that of another in-plane antiferromagnet, Fe 2 ⁢As (less than 1 T). Finally, the subtle anisotropy of intrinsic susceptibilities may lead to dominating effects from shape, crystalline texture, strain, and defects in devices that attempt spin readout in Mn 2⁢ Au.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Experiments on near-wall structure of three-dimensional boundary layers

Investigations of three-dimensional turbulent boundary layers have shown basic differences between two- and three-dimensional flows. These differences can significantly impact the modeling of three-dimensional flows since many flow models are based on results from two-dimensional boundary layers. In many cases the shear stress vector direction has been shown to lag relative to the direction of the mean velocity gradient as the cross flow grows downstream. Coincidence of these vectors is necessary for a scalar eddy viscosity modeling assumption. A second effect is a reduction in magnitude of the shear stress and/or the shear stress to turbulence energy ratio, a(sub 1). This reduction has been observed in several experiments. Recent numerical simulations also indicate wall-layer structural differences between two- and three-dimensional boundary layers. The differences in structure between two- and three-dimensional boundary layers was also addressed in the experiment of Littell & Eaton. The experiment used two-point correlations to investigate the vortical structures in a three dimensional boundary layer on a spinning disk. It was found that each sign of longitudinal vortex is equally likely to exist, but one sign of vorticity is associated with a structure which is better at producing ejections. The goal of the current investigation is to study the structure of the inner layers. Among other questions, the differences between the effects deduced from the three-dimensional flow simulations and the effects seen in experiments can be examined. The research concentrates on the structure of the wall-layer through flow visualization and direct turbulence measurements down to y(+) = 5.

Flack, Karen A.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
NASA NTRSDate not supplied

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗
NASA NTRS1989

The use of vestibular models for design and evaluation of flight simulator motion

Quantitative models for the dynamics of the human vestibular system are applied to the design and evaluation of flight simulator platform motion. An optimal simulator motion control algorithm is generated to minimize the vector difference between perceived spatial orientation estimated in flight and in simulation. The motion controller has been implemented on the Vertical Motion Simulator at NASA Ames Research Center and evaluated experimentally through measurement of pilot performance and subjective rating during VTOL aircraft simulation. In general, pilot performance in a longitudinal tracking task (formation flight) did not appear to be sensitive to variations in platform motion condition as long as motion was present. However, pilot assessment of motion fidelity by means of a rating scale designed for this purpose, were sensitive to motion controller design. Platform motion generated with the optimal motion controller was found to be generally equivalent to that generated by conventional linear crossfeed washout. The vestibular models are used to evaluate the motion fidelity of transport category aircraft (Boeing 727) simulation in a pilot performance and simulator acceptability study at the Man-Vehicle Systems Research Facility at NASA Ames Research Center. Eighteen airline pilots, currently flying B-727, were given a series of flight scenarios in the simulator under various conditions of simulator motion. The scenarios were chosen to reflect the flight maneuvers that these pilots might expect to be given during a routine pilot proficiency check. Pilot performance and subjective rating of simulator fidelity was relatively insensitive to the motion condition, despite large differences in the amplitude of motion provided. This lack of sensitivity may be explained by means of the vestibular models, which predict little difference in the modeled motion sensations of the pilots when different motion conditions are imposed.

Bussolari, Steven R.↗
DOE OSTI2025

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗
NASA NTRS2002

Expression and affinity purification of recombinant proteins from plants

With recent advances in plant biotechnology, transgenic plants have been targeted as an inexpensive means for the mass production of proteins for biopharmaceutical and industrial uses. However, the current plant purification techniques lack a generally applicable, economic, large-scale strategy. In this study, we demonstrate the purification of a model protein, beta-glucuronidase (GUS), by employing the protein calmodulin (CaM) as an affinity tag. In the proposed system, CaM is fused to GUS. In the presence of calcium, the calmodulin fusion protein binds specifically to a phenothiazine-modified surface of an affinity column. When calcium is removed with a complexing agent, e.g., EDTA, calmodulin undergoes a conformational change allowing the dissociation of the calmodulin-phenothiazine complex and, therefore, permitting the elution of the GUS-CaM fusion protein. The advantages of this approach are the fast, efficient, and economical isolation of the target protein under mild elution conditions, thus preserving the activity of the target protein. Two types of transformation methods were used in this study, namely, the Agrobacterium-mediated system and the viral-vector-mediated transformation system. Copyright 2002 Elsevier Science (USA).

Non-NASA Center↗
NASA NTRS2008

Energetic Consistency and Coupling of the Mean and Covariance Dynamics

The dynamical state of the ocean and atmosphere is taken to be a large dimensional random vector in a range of large-scale computational applications, including data assimilation, ensemble prediction, sensitivity analysis, and predictability studies. In each of these applications, numerical evolution of the covariance matrix of the random state plays a central role, because this matrix is used to quantify uncertainty in the state of the dynamical system. Since atmospheric and ocean dynamics are nonlinear, there is no closed evolution equation for the covariance matrix, nor for the mean state. Therefore approximate evolution equations must be used. This article studies theoretical properties of the evolution equations for the mean state and covariance matrix that arise in the second-moment closure approximation (third- and higher-order moment discard). This approximation was introduced by EPSTEIN [1969] in an early effort to introduce a stochastic element into deterministic weather forecasting, and was studied further by FLEMING [1971a,b], EPSTEIN and PITCHER [1972], and PITCHER [1977], also in the context of atmospheric predictability. It has since fallen into disuse, with a simpler one being used in current large-scale applications. The theoretical results of this article make a case that this approximation should be reconsidered for use in large-scale applications, however, because the second moment closure equations possess a property of energetic consistency that the approximate equations now in common use do not possess. A number of properties of solutions of the second-moment closure equations that result from this energetic consistency will be established.

Cohn, Stephen E.↗
DOE OSTI2025

Superfluid weight of strongly inhomogeneous superconductors

In this study, we obtain the superfluid weight of strongly inhomogeneous superconductors. We find that, in general, the response of the superconductor’s pairing potential to the perturbing vector potential is essential to obtain the correct value of the superfluid weight. We consider two exemplary cases: the case when strong inhomogeneities in the pairing potential are induced by a periodic potential, and the case when superconducting vortices are induced by an external magnetic field. For both cases we show that the correction to the superfluid weight due to the response of the paring potential to the perturbing vector potential can be significant, it must be included to obtain quantitatively correct results, and that for the case when vortices are present the expression of the superfluid weight that does not include such correction returns qualitatively wrong results. The results are directly relevant for superconducting states in moiré multilayer where the almost inevitable spatial variations of the magic angle and the proximity of superconductivity to anisotropic phases induce inhomogeneities in the pairing potential, and to the efforts to optimally pin superconducting vortices to achieve large superconducting critical currents.

Bogoliubov-de Gennes equations↗
NASA NTRS2011

Influence of Vibrotactile Feedback on Controlling Tilt Motion After Spaceflight

We hypothesize that adaptive changes in how inertial cues from the vestibular system are integrated with other sensory information leads to perceptual disturbances and impaired manual control following transitions between gravity environments. The primary goals of this ongoing post-flight investigation are to quantify decrements in manual control of tilt motion following short-duration spaceflight and to evaluate vibrotactile feedback of tilt as a sensorimotor countermeasure. METHODS. Data is currently being collected on 9 astronaut subjects during 3 preflight sessions and during the first 8 days after Shuttle landings. Variable radius centrifugation (216 deg/s, <20 cm radius) in a darkened room is utilized to elicit otolith reflexes in the lateral plane without concordant canal or visual cues. A Tilt-Translation Sled (TTS) is capable of synchronizing pitch tilt with fore-aft translation to align the resultant gravitoinertial vector with the longitudinal body axis, thereby eliciting canal reflexes without concordant otolith or visual cues. A simple 4 tactor system was implemented to provide feedback when tilt position exceeded predetermined levels in either device. Closed-loop nulling tasks are performed during random tilt steps or sum-of-sines (TTS only) with and without vibrotactile feedback of chair position. RESULTS. On landing day the manual control performance without vibrotactile feedback was reduced by >30% based on the gain or the amount of tilt disturbance successfully nulled. Manual control performance tended to return to baseline levels within 1-2 days following landing. Root-mean-square position error and tilt velocity were significantly reduced with vibrotactile feedback. CONCLUSIONS. These preliminary results are consistent with our hypothesis that adaptive changes in vestibular processing corresponds to reduced manual control performance following G-transitions. A simple vibrotactile prosthesis improves the ability to null out tilt motion within a limited range of motion disturbances.

Wood, S. J.↗
NASA NTRS1989

An imaging vector magnetograph for the next solar maximum

Researchers describe the conceptual design of a new imaging vector magnetograph currently being constructed at the University of Hawaii. The instrument combines a modest solar telescope with a rotating quarter-wave plate, an acousto-optical tunable prefilter as a blocker for a servo-controlled Fabry-Perot etalon, CCD cameras, and on-line digital image processing. Its high spatial resolution (1/2 arcsec pixel size) over a large field of view (5 by 5 arcmin) will be sufficient to significantly measure, for the first time, the magnetic energy dissipated in major solar flares. Its millisecond tunability and wide spectral range (5000 to 7000 A) enable nearly simultaneous vector magnetic field measurements in the gas-pressure-dominated photosphere and magnetically-dominated chromosphere, as well as effective co-alignment with Solar-A's X ray images. Researchers expect to have the instrument in operation at Mees Solar Observatory (Haleakala) in early 1991. They have chosen to use tunable filters as wavelength-selection elements in order to emphasize the spatial relationships between magnetic field elements, and to permit construction of a compact, efficient instrument. This means that spectral information must be obtained from sequences of images, which can cause line profile distortions due to effects of atmospheric seeing.

Mickey, D. L.↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Classifying Agnostic Biosignatures using Raman, VNIR, and Elemental Data

How can we use our current wealth of terrestrial data, encompassing biogenic and abiogenic systems, to determine the distinguishing properties of life? SCOBI (Statistical Classification of Biosignature Information) uses machine learning techniques to algorithmically identify combinations of measurements that are “indicative of life”. A set of ~1000 observations, comprising elemental abundance, isotopic fractionation, VNIR reflectance, and (in progress) Raman spectra, have been assembled from existing literature and databases. The observations cover systems classified as “indicative alive” (e.g., cells, vegetation), “indicative non-alive” (e.g., fossils, teeth), “mixed indicative” (e.g., soil, pond water), or “non-indicative” (e.g., rocks, meteorites). VNIR data was preprocessed by linear interpolation from 400-2100 nm and smoothed with a Savitzky-Golay filter. To limit the amount of Earth-biochemistry-specific (non-agnostic) information included, the first five spectral features extracted were number of peaks, number of troughs, mean reflectance, mean peak width, and broadest peak width. To help further emphasize agnostic biosignatures, Earth-specific features such as chlorophylls have been manually flagged so that feature importance with and without them can be compared. Classifiers including k-nearest neighbors (KNN), Gaussian Naïve Bayes (GNB), logistic regression (LR), random forest (RF), and support vector machine (SVM) were implemented, as was a combination voting classifier. Performance metrics included false positive rates, false negative rates, and AUC with 50-50 test/train splits (Monte Carlo simulations). Key takeaways from this stage, prior to the inclusion of Raman spectra, are (1) the overall success rate of 0.933 AUC was most heavily influenced by the elemental abundance data; and (2) VNIR reflectance had the lowest classification performance with 0.52 AUC (58% of objects correctly classified). The next steps are to complete integration of Raman spectral data and to improve the approach to pre-processing and feature extraction for both types of spectral data, such as automated baseline removal, whole spectrum matching, and dimensionality reduction.

Biosignatures↗