Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical 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 937 records · Page 52

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Testing convolutional neural network based deep learning systems: a statistical metamorphic approach

Machine learning technology spans many areas and today plays a significant role in addressing a wide range of problems in critical domains,i.e., healthcare, autonomous driving, finance, manufacturing, cybersecurity,etc. Metamorphic testing (MT) is considered a simple but very powerful approach in testing such computationally complex systems for which either an oracle is not available or is available but difficult to apply. Conventional metamorphic testing techniques have certain limitations in verifying deep learning-based models (i.e., convolutional neural networks (CNNs)) that have a stochastic nature (because of randomly initializing the network weights) in their training. In this article, we attempt to address this problem by using a statistical metamorphic testing (SMT) technique that does not require software testers to worry about fixing the random seeds (to get deterministic results) to verify the metamorphic relations (MRs). We propose seven MRs combined with different statistical methods to statistically verify whether the program under test adheres to the relation(s) specified in the MR(s). We further use mutation testing techniques to show the usefulness of the proposed approach in the healthcare space and test two CNN-based deep learning models (used for pneumonia detection among patients). The empirical results show that our proposed approach uncovers 85.71% of the implementation faults in the classifiers under test (CUT). Furthermore, we also propose an MRs minimization algorithm for the CUT, thus saving computational costs and organizational testing resources.

Computer Science↗

Spectral eclipse mapping of the accretion disk in the nova-like variable UX Ursae Majoris

We analyze narrow-band eclipse light curves of the nova-like cataclysmic variable UX UMa, obtained from low-resolution spectra spanning lambda lambda 3600-9800 A . The light curves for narrow bands in the continuum as well as those for individual spectral lines are treated independently, and are used to construct images of the accretion disk's brightness distribution using the maximum-entropy eclipse-mapping technique. Particular attention is paid to the propagation of statistical uncertainties in the data and to how the analysis may introduce systematic errors in the final result. From the many narrrow band images we have reconstructed the spectra from isolated parts of the accretion disk. These spectra reveal that the inner disk radiates a continuum spectrum which peaks in the near UV and has the hydrogen Balmer lines in absorption (with the exception of H-alpha), whereas the outer disk is much fainter, has a much redder spectrum, and has Balmer emission lines. Our analysis reveals the presence of an uneclipsed component of the total light, whose spectrum is very red and has Balmer lines in emission. This unexpected feature of the eclipse mapping technique offers a new tool for an independent assessment of the secondary star's spectrum in eclipsing cataclysmic variables.

Rutten, R. G. M.↗

Problems in Recovering Crater Production Functions

Two techniques proposed for the recovery of the crater production functions on densely cratered surfaces were appraised. The first technique classifies craters according to their degree of rim preservation, then interprets the statistics of the freshest craters as closely mimicing the full population's production function. This method is misleading for densely cratered surfaces. The second technique measures and weighs each crater diameter by the fraction of the crater that remains unoverlapped by subsequent impacts. The method is theoretically sound even at saturation crater densities, though it too fails at high densities.

Woronow, A.↗

Extensions of the space trajectories error analysis programs

A generalized covariance analysis technique which permits the study of the sensitivity of linear estimation algorithms to errors in a priori statistics has been developed and programed. Several sample cases are presented to illustrate the use of this technique. Modifications to the Simulated Trajectories Error Analysis Program (STEAP) to enable targeting a multiprobe mission of the Planetary Explorer type are discussed. The logic for the mini-probe targeting is presented. Finally, the initial phases of the conversion of the Viking mission Lander Trajectory Reconstruction (LTR) program for use on Venus missions is discussed. An integrator instability problem is discussed and a solution proposed.

Adams, G. L.↗

Microbial Risk Assessment

Historically, microbiological spaceflight requirements have been established in a subjective manner based upon expert opinion of both environmental and clinical monitoring results and the incidence of disease. The limited amount of data, especially from long-duration missions, has created very conservative requirements based primarily on the concentration of microorganisms. Periodic reevaluations of new data from later missions have allowed some relaxation of these stringent requirements. However, the requirements remain very conservative and subjective in nature, and the risk of crew illness due to infectious microorganisms is not well defined. The use of modeling techniques for microbial risk has been applied in the food and potable water industries and has exceptional potential for spaceflight applications. From a productivity standpoint, this type of modeling can (1) decrease unnecessary costs and resource usage and (2) prevent inadequate or inappropriate data for health assessment. In addition, a quantitative model has several advantages for risk management and communication. By identifying the variable components of the model and the knowledge associated with each component, this type of modeling can: (1) Systematically identify and close knowledge gaps, (2) Systematically identify acceptable and unacceptable risks, (3) Improve communication with stakeholders as to the reasons for resource use, and (4) Facilitate external scientific approval of the NASA requirements. The modeling of microbial risk involves the evaluation of several key factors including hazard identification, crew exposure assessment, dose-response assessment, and risk characterization. Many of these factors are similar to conditions found on Earth; however, the spaceflight environment is very specialized as the inhabitants live in a small, semi-closed environment that is often dependent on regenerative life support systems. To further complicate modeling efforts, microbial dose-response characteristics may be affected by a potentially dysfunctional crew immune system during a mission. In addition, microbial virulence has been shown to change under certain conditions during spaceflight, further complicating dose-response characterization. An initial study of the applicability of microbial risk assessment techniques was performed using Crew Health Care System (CHeCS) operational data from the International Space Station potable water systems. The risk of infection from potable water was selected as the flight systems and microbial ecology are well defined. This initial study confirmed the feasibility of using microbial risk assessment modeling for spaceflight systems. While no immediate threat was detected, the study identified several medically significant microorganisms that could pose a health risk if uncontrolled. The study also identified several specific knowledge gaps in making a risk assessment and noted that filling these knowledge gaps is essential as the risk estimates may change by orders of magnitude depending on the answers. The current phase of the microbial risk assessment studies focuses on the dose-response relationship of specific infectious agents, focusing on Salmonella enterica Typhimurium, Pseudomonas spp., and Escherichia coli, as their evaluation will provide a better baseline for determining the overall hazard characterization. The organisms were chosen as they either have been isolated on spacecraft or have an identified route of infection during a mission. The characterization will utilize dose-response models selected either from the peer-reviewed literature and/or by using statistical approaches. Development of these modeling and risk assessment techniques will help to optimize flight requirements and to protect the safety, health, and performance of the crew.

Ott, C. M.↗

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Introduction to Data Assimilation

Atmospheric data assimilation is a class of techniques used for producing descriptions of fields of air temperature, pressure, humidity, wind, etc. on a spatial grid or in terms of a finite functional representation. These are then used to initialize numerical weather forecasts or to analyze the atmosphere for other purposes. The techniques combine past, present, and even future observations in an approximate statistically optimal way. Various types of statistical or physically-based models and their corresponding adjoints are employed to relate diverse fields in both time and space and to relate what is observed to what is being analyzed. Computationally, the problem is very demanding and onstraining on the techniques that can be employed on a routine basis.

Atmospheric data assimilation↗

High mass and halo resolution from fast low resolution simulations

Generating mocks for future sky surveys requires large volumes and high resolutions, which is computationally expensive even for fast simulations. Here we try to develop numerical schemes to calibrate various halo and matter statistics in fast low resolution simulations compared to high resolution N-body and hydrodynamic simulations. For the halos, we improve the initial condition resolution and develop a halo finder "relaxed-FoF", where we allow different linking lengths for different halo mass and velocity dispersions. We show that our relaxed-FoF halo finder improves the common statistics, such as halo bias, halo mass function, halo auto power spectrum, cross correlation coefficient with the reference halo catalog, and halo-matter cross power spectrum. We also calibrate small-scale velocities of small halos to improve the power spectrum in redshift space. For the matter statistics, we incorporate the potential gradient descent (PGD) method into fast simulations to improve the matter distribution at nonlinear scales. By building a lightcone output, we show that the PGD method significantly improves the weak lensing convergence tomographic power spectrum. With these improvements FastPM is comparable to the high resolution full N-body simulation of the same mass resolution, with two orders of magnitude fewer time steps. These techniques can be used to improve the halo and matter statistics of FastPM simulations for mock catalogs of future surveys such as DESI and LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Modern estimation techniques applied to microwave sensing of the marine boundary layer

Previous efforts in the area of satellite microwave sensing of the marine boundary layer have relied upon linear regression techniques to extract geophysical parameters from the microwave measurement data. The approach used in the present paper shifts emphasis away from the generation of regression weighting matrices which implicitly assume that the data are linear in the parameters to be determined and that the problem is statistically stationary. The idea is simply to employ modern computational estimation techniques to obtain parameter estimates from nonlinear noisy measurements. The approach is limited only to the region of validity of Grody's (1976) model. Attention is focused on documenting how estimation techniques, in particular the square root information filter (SRIF), are used to solve a nonlinear function optimization problem.

Bierman, G. J.↗

Shielding Effectiveness in a Two-Dimensional Reverberation Chamber Using Finite-Element Techniques

Reverberation chambers are attaining an increased importance in determination of electromagnetic susceptibility of avionics equipment. Given the nature of the variable boundary condition, the ability of a given source to couple energy into certain modes and the passband characteristic due the chamber Q, the fields are typically characterized by statistical means. The emphasis of this work is to apply finite-element techniques at cutoff to the analysis of a two-dimensional structure to examine the notion of shielding-effectiveness issues in a reverberating environment. Simulated mechanical stirring will be used to obtain the appropriate statistical field distribution. The shielding effectiveness (SE) in a simulated reverberating environment is compared to measurements in a reverberation chamber. A log-normal distribution for the SE is observed with implications for system designers. The work is intended to provide further refinement in the consideration of SE in a complex electromagnetic environment.

Bunting, Charles F.↗

MAG4 Versus Alternative Techniques for Forecasting Active-Region Flare Productivity

MAG4 is a technique of forecasting an active region's rate of production of major flares in the coming few days from a free-magnetic-energy proxy. We present a statistical method of measuring the difference in performance between MAG4 and comparable alternative techniques that forecast an active region's major-flare productivity from alternative observed aspects of the active region. We demonstrate the method by measuring the difference in performance between the "Present MAG4" technique and each of three alternative techniques, called "McIntosh Active-Region Class," "Total Magnetic Flux," and "Next MAG4." We do this by using (1) the MAG4 database of magnetograms and major-flare histories of sunspot active regions, (2) the NOAA table of the major-flare productivity of each of 60 McIntosh active-region classes of sunspot active regions, and (3) five technique-performance metrics (Heidke Skill Score, True Skill Score, Percent Correct, Probability of Detection, and False Alarm Rate) evaluated from 2000 random two-by-two contingency tables obtained from the databases. We find that (1) Present MAG4 far outperforms both McIntosh Active-Region Class and Total Magnetic Flux, (2) Next MAG4 significantly outperforms Present MAG4, (3) the performance of Next MAG4 is insensitive to the forward and backward temporal windows used, in the range of one to a few days, and (4) forecasting from the free-energy proxy in combination with either any broad category of McIntosh active-region classes or any Mount Wilson active-region class gives no significant performance improvement over forecasting from the free-energy proxy alone (Present MAG4).

Falconer, David A.↗

MAG4 Versus Alternative Techniques for Forecasting Active-Region Flare Productivity

MAG4 is a technique of forecasting an active region's rate of production of major flares in the coming few days from a free-magnetic-energy proxy. We present a statistical method of measuring the difference in performance between MAG4 and comparable alternative techniques that forecast an active region's major-flare productivity from alternative observed aspects of the active region. We demonstrate the method by measuring the difference in performance between the "Present MAG4" technique and each of three alternative techniques, called "McIntosh Active-Region Class," "Total Magnetic Flux," and "Next MAG4." We do this by using (1) the MAG4 database of magnetograms and major-flare histories of sunspot active regions, (2) the NOAA table of the major-flare productivity of each of 60 McIntosh active-region classes of sunspot active regions, and (3) five technique-performance metrics (Heidke Skill Score, True Skill Score, Percent Correct, Probability of Detection, and False Alarm Rate) evaluated from 2000 random two-by-two contingency tables obtained from the databases. We find that (1) Present MAG4 far outperforms both McIntosh Active-Region Class and Total Magnetic Flux, (2) Next MAG4 significantly outperforms Present MAG4, (3) the performance of Next MAG4 is insensitive to the forward and backward temporal windows used, in the range of one to a few days, and (4) forecasting from the free-energy proxy in combination with either any broad category of McIntosh active-region classes or any Mount Wilson active-region class gives no significant performance improvement over forecasting from the free-energy proxy alone (Present MAG4).

Falconer, David A.↗

Gas-Phase Study of the Elementary Reaction of the D1-Ethynyl Radical (C 2 D; X 2 Σ + ) with Propylene (C 3 H 6 ; X 1 A') under Single-Collision Conditions

The bimolecular gas-phase reactions of the D1-ethynyl radical (C 2 D; X 2 Σ + ) with propylene (C 3 H 6 ; X 1 A’) and partially substituted D3-3,3,3-propylene (C 2 H 3 CD 3 ; X 1 A’) were studied under single collision conditions utilizing the crossed molecular beams technique. Combining our laboratory data with electronic structure and statistical calculations, the D1-ethynyl radical is found to add without barrier to the C1 and C2 carbons of the propylene reactant, resulting in doublet C 5 H 6 D intermediate(s) with lifetime(s) longer than their rotational period(s). These intermediates undergo isomerization and unimolecular decomposition via atomic hydrogen loss through tight exit transition states forming predominantly cis/trans-3-penten-1-yne ((HCC)CH=CH(CH 3 )) and to a minor amount 3-methyl-3-buten-1-yne ((HCC)C(CH 3 )=CH 2 ) via overall exoergic reactions. Although the title reaction does not lead to the cyclopentadiene molecule (c-C 5 H 6 , X 1 A 1 ), high temperature environments can convert the identified acyclic C 5 H 6 isomers through hydrogen atom assisted isomerization to cyclopentadiene (c-C 5 H 6 , X 1 A 1 ). Since both the ethynyl radical and propylene reactants have been observed in cold interstellar environments such as TMC-1 and the reaction is exoergic and all barriers lie below the energy of the separated reactants, the these C 5 H 6 product isomers are predicted to form in those low temperature regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

2023 Roadmap on molecular modelling of electrochemical energy materials

New materials for electrochemical energy storage and conversion are the key to the electrification and sustainable development of our modern societies. Molecular modelling based on the principles of quantum mechanics and statistical mechanics as well as empowered by machine learning techniques can help us to understand, control and design electrochemical energy materials at atomistic precision. Therefore, this roadmap, which is a collection of authoritative opinions, serves as a gateway for both the experts and the beginners to have a quick overview of the current status and corresponding challenges in molecular modelling of electrochemical energy materials for batteries, supercapacitors, CO 2 reduction reaction, and fuel cell applications.

25 ENERGY STORAGE↗

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗

Framework for an adaptive integrated observation system using a hierarchy of machine learning approaches

Focal Area(s): 1. Data acquisition enabled by machine learning, AI, and advanced methods including experimental/network design/optimization, and hardware-related efforts involving AI. 2. Insight gleaned from complex measurements using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI Science Challenge and Rationale: Atmospheric processes are stochastic, occur at scales from the micrometer to many kilometers, and are constantly changing over time. Characterizing these interactions and associated environmental conditions using traditional measurement techniques is difficult and can take years to build statistics on atmospheric phenomena that occurs episodically. Developing new innovative approaches to modify sampling strategies in real-time would enable the routine collection of targeted measurements focused on a specific set of science questions.

54 ENVIRONMENTAL SCIENCES↗