Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Insights into Copper Sulfide Formation from Cu and S K edge XAS and DFT studies

An understanding of the fundamentals of the reaction between CuO with trace amounts of H 2 S to form CuS products is critical for the optimal utilization of this process in sulfur removal applications. Unfortunately, CuS is a complex material, featuring various Cu 2-x S compounds (with 0 ≤ x ≤ 1), distorted crystal phases, and varying electronic structures and coordination environments of Cu and S ions. In this work, we combine ex situ and in situ X-ray absorption spectroscopy (XAS) at S and Cu K edges, fixed bed sorption experiments, DFT simulations, and other characterization techniques to speciate the CuS products formed at different temperatures (298–383 K) and from CuO sorbents with different crystallite sizes (2.8–40 nm). The results of our analysis identify the formation of a distorted CuS layer at the surface of CuO crystals with disulfide groups with shorter Cu–S bonds and higher delocalization of the positive charge of the Cu center into (S 1– )2. This distorted CuS layer dominates the XAS signal at lower temperatures (298–323 K) and at the initial stages of sulfidation at higher temperatures (353 and 383 K) where conversion is low (<40%). First-principles atomistic simulations confirm the thermodynamic favorability of the formation of surface (S 1– )2 on both CuO (111) and ($\bar{1}11$) surfaces, providing further support for our experimental observations. Furthermore, these simulations reveal that the presence of disulfide bonds stabilized surface hydroxyl groups, leading to lower Gibbs Free Energies of their surface migration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Applicability of Loads Estimation Techniques Using Sparse Acceleration Sensor Data to Spacecraft Structural Health Monitoring

The use of structural health monitoring systems on spacecraft structures can play a crucial role in ensuring the safety, reliability, and longevity of the structure by gathering and analyzing onboard sensor data. Of specific importance is monitoring for excessive loading at critical interfaces as any off-nominal structural excitations experienced by spacecraft structures can cause early unpredicted high structural life consumption or damage. The availability and cost of flight-certified sensors along with the size of spacecraft structures and allowable payload mass drives the need for a method to estimate loads using sparsely-located sensors. Numerous approaches such as physics-based, statistical learning, and physics-enhanced statistical learning algorithms have gained popularity among structural prognostics applications. However, developing noise-robust prediction models to assess loads and structural life predictions from a sparse multi-sensor data acquisition system can be a challenging task. This paper discusses the evaluation of physics-based versus machine-learning algorithms for predicting loads and structural life at mission critical locations on the spacecraft structure using a finite element loads analysis with the application of simulated noise and noise reduction techniques. To estimate the loads from accelerations, the physics-based algorithm leverages a loads transformation matrix from a Craig-Bampton reduced finite element model. A System Equivalent Reduction Expansion Process (SEREP) and a pseudo-inverse approach are considered to expand from the onboard sensor degrees of freedom to the Craig-Bampton model degrees of freedom. The machine learning algorithm provides a data driven solution/mapping of the sensor accelerations to the loads at the mission critical locations using a high dimensionality analysis. Although these strategies produce comparable loads prediction without noise, the limitations of these strategies with incorporating simulated noise and noise reduction techniques with low signal to noise ratio signals are evaluated. The study demonstrates the immense potential of statistical learning algorithms for sparse structural prognostic models and enhancing signal denoising techniques. These findings also highlight the need for noise-resilient prognostic models and low-noise data acquisition systems onboard spacecraft structures.

Spacecraft Structural Health Monitoring↗
DOE OSTI2021

Artificial intelligence techniques for integrative structural biology of intrinsically disordered proteins

Here, we outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self-organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure-function relationships - however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Harnessing Uncertainty through Functional Data Analysis in Gas Breakthrough Data

Detecting subsurface explosions from radionuclide gas migration through rock fractures is an effective way to identify nuclear activity. Los Alamos National Laboratory (LANL) has developed simulation methods, based on data from the 1962 Hardhat underground nuclear test, to predict gas breakthrough times at the surface. However, these methods rely on an imperfect understanding of the relationship between rock damage and fracture permeability. Our clinic project studies methods for predicting breakthrough curves that characterize total mass produced as a function of time, as well as quantifying the uncertainty associated with these predictions. The model that is currently employed to relate damage to permeability uses an empirically motivated power-law expression, with a range of parameter values that are compatible with the experimental Hardhat data. We develop emulators, built from functional data analysis techniques and trained on simulation data, that rapidly predict the gas breakthrough curve given a damage field and given the parameter values of the power-law equation. Using Bayesian regression, we address the problem of uncertainty quantification in our emulators. Finally, in order to test the robustness of the model, we further validate it on a damage field representing different physical conditions.

58 GEOSCIENCES↗
NASA NTRS2004

Statistical Comparison and Improvement of Methods for Combining Random and Harmonic Loads

Structures in many environments experience both random and harmonic excitation. A variety of closed-form techniques has been used in the aerospace industry to combine the loads resulting from the two sources. The resulting combined loads are then used to design for both yield ultimate strength and high cycle fatigue capability. This paper examines the cumulative distribution function (CDF) percentiles obtained using each method by integrating the joint probability density function of the sine and random components. A new Microsoft Excel spreadsheet macro that links with the software program Mathematics is then used to calculate the combined value corresponding to any desired percentile along with a curve fit to this value. Another Excel macro is used to calculate the combination using a Monte Carlo simulation. Unlike the traditional techniques, these methods quantify the calculated load value with a Consistent percentile. Using either of the presented methods can be extremely valuable in probabilistic design, which requires a statistical characterization of the loading. Also, since the CDF at high probability levels is very flat, the design value is extremely sensitive to the predetermined percentile; therefore, applying the new techniques can lower the design loading substantially without losing any of the identified structural reliability.

Brown, Andrew M.↗
NASA NTRS1977

Computer considerations for real time simulation of a generalized rotor model

Scaled equations were developed to meet requirements for real time computer simulation of the rotor system research aircraft. These equations form the basis for consideration of both digital and hybrid mechanization for real time simulation. For all digital simulation estimates of the required speed in terms of equivalent operations per second are developed based on the complexity of the equations and the required intergration frame rates. For both conventional hybrid simulation and hybrid simulation using time-shared analog elements the amount of required equipment is estimated along with a consideration of the dynamic errors. Conventional hybrid mechanization using analog simulation of those rotor equations which involve rotor-spin frequencies (this consititutes the bulk of the equations) requires too much analog equipment. Hybrid simulation using time-sharing techniques for the analog elements appears possible with a reasonable amount of analog equipment. All-digital simulation with affordable general-purpose computers is not possible because of speed limitations, but specially configured digital computers do have the required speed and consitute the recommended approach.

Howe, R. M.↗
NASA NTRS1973

Measurement techniques for antennas in dissipative media.

Methods used to simulate dissipative media environments for antennas are reported, and special techniques are developed to determine the electrical properties of antennas in such media. A new device for measuring antenna input admittance is presented which permits a simple measurement of admittance over a broad frequency range where conventional methods are cumbersome in field use. An apparatus for measuring antenna current and charge distributions is described, and the associated probing errors are discussed. Measured admittances and current and charge distributions are compared with current theories for thin wire linear and loop antennas in the dissipative media and are found to be in good agreement.

Scott, L. D.↗
DOE OSTI2021

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗
DOE OSTI2017

System Modeling and Simulation of In-Stream Hydrokinetic Turbines for Power Management and Control

Electricity generation from moving currents without using dams (i.e., in-stream hydrokinetic electricity) has the potential to introduce multiple GW of renewable power to U.S. grids. Here this study evaluates a control system designed to regulate the generator rotor rate (rpm) to improve power production from in-stream hydrokinetic turbines. The control algorithm is evaluated using numerical models of both a rigidly mounted tidal turbine (TT) and a moored ocean current turbine (OCT) coupled to an induction electric machine model. The moored simulation utilizes an innovative approach for coupling a multiple degrees-of-freedom (DOF) nonlinear hydrodynamic/mechanical turbine model with a nonlinear electromechanical generator model. Based on the turbine torque-speed characteristic, as well as the asynchronous machine features, a proportional–integral (PI) controller is used to generate a correction term for the frequency of the three-phase sinusoidal voltages that are supplied to the asynchronous generator. The speed control of the induction generator through the supply frequency is accomplished by a simplified voltage source inverter (VSI). The simplified VSI consists of control voltage sources (CVSs), while the comparison with a real VSI using diodes and transistors, which are controlled by pulse width modulation (PWM) technique, is also presented. Simulations are used to evaluate the developed algorithms showing that rpm fluctuations are around 0.02 for a tidal turbine operating in a wave field with a 6 m significant wave height and around 0.005 for a moored ocean current turbine operating in a wave field with a 2 m significant wave height.

16 TIDAL AND WAVE POWER↗
NASA NTRS1980

A new method of measuring radio source parameters of a partially polarized distributed source from spacecraft observations

An experimental system is analyzed which is designed to determine the angular and polarization properties of low frequency radio sources from measurements made on a spinning spacecraft. The system has been optimized to provide high accuracies for sources at high as well as low elevation angles. Theoretical expressions are derived for the response of this system to a partially polarized point source. Integrations are then carried out to get the system response to a uniform circular distributed source. Data processing techniques are derived so that computer simulations can be carried out to investigate the accuracy of this technique. It is shown that using 24 measurements of a partially polarized source (with Q = U = V = 0.5), taken in one spacecraft rotation, the resulting rms errors in angular position are less than one degree and the errors in determining the Stokes parameters are generally 1-10% for a wide range of source elevations.

Manning, R.↗
NASA NTRS1986

Robust detection, isolation and accommodation for sensor failures

The objective is to extend the recent advances in robust control system design of multivariable systems to sensor failure detection, isolation, and accommodation (DIA), and estimator design. This effort provides analysis tools to quantify the trade-off between performance robustness and DIA sensitivity, which are to be used to achieve higher levels of performance robustness for given levels of DIA sensitivity. An innovations-based DIA scheme is used. Estimators, which depend upon a model of the process and process inputs and outputs, are used to generate these innovations. Thresholds used to determine failure detection are computed based on bounds on modeling errors, noise properties, and the class of failures. The applicability of the newly developed tools are demonstrated on a multivariable aircraft turbojet engine example. A new concept call the threshold selector was developed. It represents a significant and innovative tool for the analysis and synthesis of DiA algorithms. The estimators were made robust by introduction of an internal model and by frequency shaping. The internal mode provides asymptotically unbiased filter estimates.The incorporation of frequency shaping of the Linear Quadratic Gaussian cost functional modifies the estimator design to make it suitable for sensor failure DIA. The results are compared with previous studies which used thresholds that were selcted empirically. Comparison of these two techniques on a nonlinear dynamic engine simulation shows improved performance of the new method compared to previous techniques

Emami-Naeini, A.↗
NASA NTRS1991

A trajectory planning scheme for spacecraft in the space station environment

Simulated annealing is used to solve a minimum fuel trajectory problem in the space station environment. The environment is special because the space station will define a multivehicle environment in space. The optimization surface is a complex nonlinear function of the initial conditions of the chase and target crafts. Small permutations in the input conditions can result in abrupt changes to the optimization surface. Since no prior knowledge about the number or location of local minima on the surface is available, the optimization must be capable of functioning on a multimodal surface. It was reported in the literature that the simulated annealing algorithm is more effective on such surfaces than descent techniques using random starting points. The simulated annealing optimization was found to be capable of identifying a minimum fuel, two-burn trajectory subject to four constraints which are integrated into the optimization using a barrier method. The computations required to solve the optimization are fast enough that missions could be planned on board the space station. Potential applications for on board planning of missions are numerous. Future research topics may include optimal planning of multi-waypoint maneuvers using a knowledge base to guide the optimization, and a study aimed at developing robust annealing schedules for potential on board missions.

Soller, Jeffrey Alan↗
NASA NTRS1999

Virtual Petaflop Simulation: Parallel Potential Solvers and New Integrators for Gravitational Systems

"The orbit of any one planet depends on the combined motion of all the planets, not to mention the actions of all these on each other. To consider simultaneously all these causes of motion and to define these motions by exact laws allowing of convenient calculation exceeds, unless I am mistaken, the forces of the entire human intellect" -Isaac Newton 1687. Epochal surveys are throwing down the gauntlet for cosmological simulation. We describe three keys to meeting the challenge of N-body simulation: adaptive potential solvers, adaptive integrators and volume renormalization. With these techniques and a dedicated Teraflop facility, simulation can stay even with observation of the Universe. We also describe some problems in the formation and stability of planetary systems. Here, the challenge is to perform accurate integrations that retain Hamiltonian properties for 10(exp 13) timesteps.

Lake, George↗
NASA NTRS2007

Potential of VIIRS Time Series Data for Aiding the USDA Forest Service Early Warning System for Forest Health Threats: A Gypsy Moth Defoliation Case Study

This report details one of three experiments performed during FY 2007 for the NASA RPC (Rapid Prototyping Capability) at Stennis Space Center. This RPC experiment assesses the potential of VIIRS (Visible/Infrared Imager/Radiometer Suite) and MODIS (Moderate Resolution Imaging Spectroradiometer) data for detecting and monitoring forest defoliation from the non-native Eurasian gypsy moth (Lymantria dispar). The intent of the RPC experiment was to assess the degree to which VIIRS data can provide forest disturbance monitoring information as an input to a forest threat EWS (Early Warning System) as compared to the level of information that can be obtained from MODIS data. The USDA Forest Service (USFS) plans to use MODIS products for generating broad-scaled, regional monitoring products as input to an EWS for forest health threat assessment. NASA SSC is helping the USFS to evaluate and integrate currently available satellite remote sensing technologies and data products for the EWS, including the use of MODIS products for regional monitoring of forest disturbance. Gypsy moth defoliation of the mid-Appalachian highland region was selected as a case study. Gypsy moth is one of eight major forest insect threats listed in the Healthy Forest Restoration Act (HFRA) of 2003; the gypsy moth threatens eastern U.S. hardwood forests, which are also a concern highlighted in the HFRA of 2003. This region was selected for the project because extensive gypsy moth defoliation occurred there over multiple years during the MODIS operational period. This RPC experiment is relevant to several nationally important mapping applications, including agricultural efficiency, coastal management, ecological forecasting, disaster management, and carbon management. In this experiment, MODIS data and VIIRS data simulated from MODIS were assessed for their ability to contribute broad, regional geospatial information on gypsy moth defoliation. Landsat and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data were used to assess the quality of gypsy moth defoliation mapping products derived from MODIS data and from simulated VIIRS data. The project focused on use of data from MODIS Terra as opposed to MODIS Aqua mainly because only MODIS Terra data was collected during 2000 and 2001-years with comparatively high amounts of gypsy moth defoliation within the study area. The project assessed the quality of VIIRS data simulation products. Hyperion data was employed to assess the quality of MODIS-based VIIRS simulation datasets using image correlation analysis techniques. The ART (Application Research Toolbox) software was used for data simulation. Correlation analysis between MODIS-simulated VIIRS data and Hyperion-simulated VIIRS data for red, NIR (near-infrared), and NDVI (Normalized Difference Vegetation Index) image data products collectively indicate that useful, effective VIIRS simulations can be produced using Hyperion and MODIS data sources. The r(exp 2) for red, NIR, and NDVI products were 0.56, 0.63, and 0.62, respectively, indicating a moderately high correlation between the 2 data sources. Temporal decorrelation from different data acquisition times and image misregistration may have lowered correlation results. The RPC experiment also generated MODIS-based time series data products using the TSPT (Time Series Product Tool) software. Time series of simulated VIIRS NDVI products were produced at approximately 400-meter resolution GSD (Ground Sampling Distance) at nadir for comparison to MODIS NDVI products at either 250- or 500-meter GSD. The project also computed MODIS (MOD02) NDMI (Normalized Difference Moisture Index) products at 500-meter GSD for comparison to NDVI-based products. For each year during 2000-2006, MODIS and VIIRS (simulated from MOD02) time series were computed during the peak gypsy moth defoliation time frame in the study area (approximately June 10 through July 27). Gypsy moth defoliation mapping products from simated VIIRS and MOD02 time series were produced using multiple methods, including image classification and change detection via image differencing. The latter enabled an automated defoliation detection product computed using percent change in maximum NDVI for a peak defoliation period during 2001 compared to maximum NDVI across the entire 2000-2006 time frame. Final gypsy moth defoliation mapping products were assessed for accuracy using randomly sampled locations found on available geospatial reference data (Landsat and ASTER data in conjunction with defoliation map data from the USFS). Extensive gypsy moth defoliation patches were evident on screen displays of multitemporal color composites derived from MODIS data and from simulated VIIRS vegetation index data. Such defoliation was particularly evident for 2001, although widespread denuded forests were also seen for 2000 and 2003. These visualizations were validated using aforementioned reference data. Defoliation patches were visible on displays of MODIS-based NDVI and NDMI data. The viewing of apparent defoliation patches on all of these products necessitated adoption of a specialized temporal data processing method (e.g., maximum NDVI during the peak defoliation time frame). The frequency of cloud cover necessitated this approach. Multitemporal simulated VIIRS and MODIS Terra data both produced effective general classifications of defoliated forest versus other land cover. For 2001, the MOD02-simulated VIIRS 400-meter NDVI classification produced a similar yet slightly lower overall accuracy (87.28 percent with 0.72 Kappa) than the MOD02 250-meter NDVI classification (88.44 percent with 0.75 Kappa). The MOD13 250-meter NDVI classification had a lower overall accuracy (79.13 percent) and a much lower Kappa (0.46). The report discusses accuracy assessment results in much more detail, comparing overall classification and individual class accuracy statistics for simulated VIIRS 400-meter NDVI, MOD02 250-meter NDVI, MOD02-500 meter NDVI, MOD13 250-meter NDVI, and MOD02 500-meter NDMI classifications. Automated defoliation detection products from simulated VIIRS and MOD02 data for 2001 also yielded similar, relatively high overall classification accuracy (85.55 percent for the VIIRS 400-meter NDVI versus 87.28 percent for the MOD02 250-meter NDVI). In contrast, the USFS aerial sketch map of gypsy moth defoliation showed a lower overall classification accuracy at 73.64 percent. The overall classification Kappa values were also similar for the VIIRS (approximately 0.67 Kappa) versus the MOD02 (approximately 0.72 Kappa) automated defoliation detection product, which were much higher than the values exhibited by the USFS sketch map product (overall Kappa of approximately 0.47). The report provides additional details on the accuracy of automated gypsy moth defoliation detection products compared with USFS sketch maps. The results suggest that VIIRS data can be effectively simulated from MODIS data and that VIIRS data will produce gypsy moth defoliation mapping products that are similar to MODIS-based products. The results of the RPC experiment indicate that VIIRS and MODIS data products have good potential for integration into the forest threat EWS. The accuracy assessment was performed only for 2001 because of time constraints and a relative scarcity of cloud-free Landsat and ASTER data for the peak defoliation period of the other years in the 2000-2006 time series. Additional work should be performed to assess the accuracy of gypsy moth defoliation detection products for additional years.The study area (mid-Appalachian highlands) and application (gypsy moth forest defoliation) are not necessarily representative of all forested regions and of all forest threat disturbance agents. Additional work should be performed on other inland and coastal regions as well as for other major forest threats.

Spruce, Joseph P.↗
NASA NTRSDate not supplied

Electroluminescence Imaging: A Quantitative Characterization Technique to Measure Dust Occlusion of Solar Cells

Electroluminescence (EL) imaging is a qualitative characterization technique that is typically used to identify cracks, corrosion, and other defects in solar cells. It consists of imaging a cell under forward bias, where the solar cell emits photons due to radiative electron-hole pair recombination. Over the past few years, our team has expanded this into a quantitative technique with the help of image processing. We have primarily used this method to investigate lunar dust occlusion of solar cells and arrays. Lunar dust accumulation on solar cells is a major concern because it directly limits light accessible to the cell, decreasing power output. Many teams are working to develop dust mitigation technology to protect these arrays on the surface of the Moon, but thoroughly characterizing their efficacy is important to ensure their success prior to launch. Here, we present EL as a quantitative characterization technique to observe dust coverage on solar cells that, when coupled with IV performance measurements, can offer unique insights into how dust coverage impacts power output. Quantitative electroluminescence imaging works by running EL images through an image processing script that first grayscales the image then plots a histogram of the brightness of each pixel. On its own, it does not offer much insight into a solar cell’s performance. However, when comparing images to a baseline pristine, undamaged, or uncoated solar cell, it can quickly provide information about the impacts of surface contaminants or damage to the cell. Here, we present a case study where quantitative EL is used to measure the efficacy of dust mitigation technology for flexible solar arrays and discuss the lessons learned about dust mitigation, testing with lunar simulant, and this characterization technique.

photovoltaics↗
DOE OSTI2021

Mesoscale Convective Systems Dominate the Energetics of the South Asian Summer Monsoon Onset

The onset of the South Asian summer monsoon (SASM) is featured as an abrupt increase of precipitation and atmospheric column-integrated moist static energy (MSE). Using convection permitting simulations, a novel cloud tracking technique, and an MSE budget, the relative roles of mesoscale convective systems (MCSs) and non-MCS isolated deep convection in the SASM onset are studied. The analysis shows MCSs play a crucial role in the sudden increases of rainfall, diabatic MSE sources, and the subsequent enhanced export of column-integrated MSE during the onset. A “diabatic sources-energy export” plane analysis delineates the different roles of MCSs and non-MCS deep convection in the SASM onset. Although deep convection always leads to a net energy export, mesoscale organization of deep convection dominates the convection associated positive feedbacks that accelerate the energy cycle of the SASM onset. Sensitivity experiment shows that the results are insensitive to the model resolutions.

Chen, Xingchao↗
NASA NTRS1996

Spectral Irradiance Measurements of Simulated Lightning in Planetary Atmospheres

Optical emissions from lightning provide information to estimate the altitude, latitude and longitude of lightning storms, the energy of the discharges, and the production of important trace species such as HCN and C2H2. Knowledge of the spectra of planetary lightning is needed to deduce the total energy dissipated by lightning (and thereby, the production of trace gases) and to help design experiments to detect and track lightning storms. Measurements of the spectral irradiance from approximately 380 to 820 nm are reported for laboratory simulations of lightning in the atmospheres of Venus, Jupiter, and Titan. In our laboratory, laser-induced plasmas (LIP) are used to simulate lightning discharges. This technique avoids contamination of the spectra by elect-ode material and maintains a safe environment while allowing the use of flammable gases, such as hydrogen and methane, found in outer planet atmospheres. The observations were made at I and 5 bars pressure for Venus and Jupiter and at 1 bar for the Titan mixture. At a pressure of one bar, our results show prominent lines from H(sub alpha), H(sub beta), H(sub gamma), and H(sub delta) lines of the Balmer Series of atomic hydrogen, a single line from Helium at 588 nm and strong continuum radiation. At pressures of 5 bars, the H(sub alpha) and H(sub beta) lines are wider, the H(sub gamma), and H(sub delta) lines merge into the continuum because of pressure broadening, and the helium line at 588 nm is no longer visible. The observed spectra of simulated lightning in the venusian atmosphere at 1 and 5 bars shows that the OI multiplet at 777.7 nm dominates the spectra, but weak features due to atomic carbon and singly excited and singly ionized oxygen atoms are also visible. Although lightning has not yet been observed on Titan, it conceivable that some form of lightning discharge could be occurring. Therefore experiments on a Titan atmosphere mixture were conducted. The most prominent features seen in the simulated spectra are the H(sub alpha) at 656 nm. from the dissociation of the methane and two strong features at 740 and 820 nm from singly excited nitrogen atoms. Most of the features between 400 and 656 nm are due to ionized nitrogen atoms. A broad, highly structured feature centered near 389 nm is molecular band radiation from CN molecules formed in the cooling plasma of atomic nitrogen, carbon, and hydrogen. It is observed in the lightning simulations of both the venusian and Titan atmospheres, For Venus, Jupiter, and Titan, narrow spectral features that are ideal for detecting lightning in the presence of sunlight are present in the spectra. These features could allow lightning to be detected and tracked by Earth-based instruments.

Borucki, W. J.↗
DOE OSTI2025

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗