Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered statistics”

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,351 records · Page 75

U-net architected deep material network training with microstructure local field information

The Deep Material Network (DMN) has recently emerged as a powerful reduced-order modeling framework for simulating the mechanical response of heterogeneous materials such as composites. Unlike most data-driven approaches that directly learn a material’s response under prescribed loading, the DMN acts as a homogenization operator, learning the kinematic constraints and mechanical interactions of the underlying microstructure. However, traditional DMN training relies exclusively on homogenized effective properties derived from Direct Numerical Simulations (DNS), discarding the rich local field data that govern microstructural interactions. In this work, we extend the DMN framework to incorporate such local field information into the offline training process. Utilizing a U-Net architecture, we augment the DMN training objective to include the first and second statistical moments of the local stress fields obtained from linear DNS. This ensures that the learned network topology not only fits the effective stiffness but also accurately reflects the internal local stress and strain partitioning of the microstructure. The results confirm that supervising the localization process during training yields a superior surrogate model, reducing local prediction errors by an order of magnitude and significantly improving generalization to unseen nonlinear constitutive behaviors compared to traditional DMNs.

36 MATERIALS SCIENCE↗

Thermodynamic consistency and fluctuations in mesoscopic stochastic simulations of reactive gas mixtures

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

Polimeno, Matteo [University of California, Merced↗

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗

Bayesian recursive image estimation.

The enhancement of images that are characterized only by statistical data where the picture contains additive noise is considered. The random process representing the output of the scanner is characterized by the output of a dynamic system with white noise input. The dynamic system describes a first-order vector-Markov process. The procedure of Kalman filtering is then utilized to recursively determine the minimum mean-square error estimate of the image. The result is then extended to obtain the smoothing of the data.

Nahi, N. E.↗

Linear and circular polarization of the cosmic background radiation

New data which consist of continued measurements of the linear polarization of the cosmic background radiation as well as the first measurement of the circular polarization are examined. Eleven declinations have been surveyed for linear polarization and one declination for circular polarization, all at 9 mm wavelength. No evidence is found for either a significant linear or circular component with statistical errors on the linear component of 20-60 microK for various models. For linear polarization, a 95 percent confidence level limit of 0.1 mK (0.00003) for an axisymmetric anisotropic model is achieved, while for spherical harmonics through third order, a corresponding limit of 0.2 mK is achieved. For a declination of 37 deg, a limit of 12 mK is placed on the time-varying component and 20 mK on the dc component of the circular polarization at the 95 percent confidence level. At 37 percent declination, the sensitivity per beam patch (7 deg) is 0.2 mK.

Lubin, P.↗

Coronal propagation of flare associated electrons and protons

A statistical study of characteristic times and intensities of 36 solar particle events observed between 1977 and 1979 by the Kiel Cosmic Ray Experiment on board HELIOS-1 and -2 has been carried out. For approx. 0.5 MeV electrons we order the times of maximum and the absolute intensities with respect to angular distance from the parent flare. Discussion of coronal parameters in terms of Reid's model leads to typical time constants for coronal diffusion and escape.

Schellert, G.↗

Error modeling for differential GPS

Differential Global Positioning System (DGPS) positioning is used to accurately locate a GPS receiver based upon the well-known position of a reference site. In utilizing this technique, several error sources contribute to position inaccuracy. This thesis investigates the error in DGPS operation and attempts to develop a statistical model for the behavior of this error. The model for DGPS error is developed using GPS data collected by Draper Laboratory. The Marquardt method for nonlinear curve-fitting is used to find the parameters of a first order Markov process that models the average errors from the collected data. The results show that a first order Markov process can be used to model the DGPS error as a function of baseline distance and time delay. The model's time correlation constant is 3847.1 seconds (1.07 hours) for the mean square error. The distance correlation constant is 122.8 kilometers. The total process variance for the DGPS model is 3.73 sq meters.

Blerman, Gregory S.↗

Attitude Estimation Signal Processing: A First Report on Possible Algorithms and Their Utility

In this brief effort, time has been of the essence. The data had to be acquired from APL/Lincoln Labs, stored, and sorted out to obtain the pertinent streams. This has been a significant part of this effort and hardware and software problems have been addressed with the appropriate solutions to accomplish this part of the task. Passed this, some basic and important algorithms are utilized to improve the performance of the attitude estimation systems. These algorithms are an essential part of the signal processing for the attitude estimation problem as they are utilized to reduce the amount of the additive/multiplicative noise that in general may or may not change its structure and probability density function, pdf, in time. These algorithms are not currently utilized in the processing of the data, at least, we are not aware of their use in this attitude estimation problem. Some of these algorithms, like the variable thresholding, are new conjectures, but one would expect that someone somewhere must have utilized this kind of scheme before. The variable thresholding idea is a straightforward scheme to use in case of a slowly varying pdf, or statistical moments of the unwanted random process. The algorithms here are kept simple but yet effective for processing the data and removing the unwanted noise. For the most part, these algorithms can be arranged so that their consecutive and orderly execution would complement the preceding algorithm and improve the overall performance of the signal processing chain.

Riasati, Vahid R.↗

Improved Modeling in a Matlab-Based Navigation System

An innovative approach to autonomous navigation is available for low earth orbit satellites. The system is developed in Matlab and utilizes an Extended Kalman Filter (EKF) to estimate the attitude and trajectory based on spacecraft magnetometer and gyro data. Preliminary tests of the system with real spacecraft data from the Rossi X-Ray Timing Explorer Satellite (RXTE) indicate the existence of unmodeled errors in the magnetometer data. Incorporating into the EKF a statistical model that describes the colored component of the effective measurement of the magnetic field vector could improve the accuracy of the trajectory and attitude estimates and also improve the convergence time. This model is identified as a first order Markov process. With the addition of the model, the EKF attempts to identify the non-white components of the noise allowing for more accurate estimation of the original state vector, i.e. the orbital elements and the attitude. Working in Matlab allows for easy incorporation of new models into the EKF and the resulting navigation system is generic and can easily be applied to future missions resulting in an alternative in onboard or ground-based navigation.

Deutschmann, Julie↗

Alternative Test Methods for Electronic Parts

It is common practice within NASA to test electronic parts at the manufacturing lot level to demonstrate, statistically, that parts from the lot tested will not fail in service using generic application conditions. The test methods and the generic application conditions used have been developed over the years through cooperation between NASA, DoD, and industry in order to establish a common set of standard practices. These common practices, found in MIL-STD-883, MIL-STD-750, military part specifications, EEE-INST-002, and other guidelines are preferred because they are considered to be effective and repeatable and their results are usually straightforward to interpret. These practices can sometimes be unavailable to some NASA projects due to special application conditions that must be addressed, such as schedule constraints, cost constraints, logistical constraints, or advances in the technology that make the historical standards an inappropriate choice for establishing part performance and reliability. Alternate methods have begun to emerge and to be used by NASA programs to test parts individually or as part of a system, especially when standard lot tests cannot be applied. Four alternate screening methods will be discussed in this paper: Highly accelerated life test (HALT), forward voltage drop tests for evaluating wire-bond integrity, burn-in options during or after highly accelerated stress test (HAST), and board-level qualification.

Plante, Jeannette↗

Entry, Descent, and Landing technological barriers and crewed MARSvehicle performance analysis

Mars has been explored historically only by robotic crafts, but a crewed mission encompasses several new engineering challenges – high ballistic coefficient entry, hypersonic decelerators, guided entry for reaching intended destinations within acceptable margins for error in the landing ellipse, and payload mass are all critical factors for evaluation. A comprehensive EDL parametric analysis has been conducted in support of a high mass landing architecture by evaluating three types of vehicles −70° Sphere Cone, Ellipsled and SpaceX hybrid architecture called Red Dragon as potential candidate options for crewed entry vehicles. Aerocapture at the Martian orbit of about 400 km and subsequent Entry-from-orbit scenarios were investigated at velocities of 6.75 km/s and 4 km/s respectively. A study on aerocapture corridor over a range of entry velocities (6–9 km/s) suggests that a hypersonic L/D of 0.3 is sufficient for a Martian aerocapture. Parametric studies conducted by varying aeroshell diameters from 10 m to 15 m for several entry masses up to 150 mt are summarized and results reveal that vehicles with entry masses in the range of about 40–80 mt are capable of delivering cargo with a mass on the order of 5–20 mt. For vehicles with an entry mass of 20 mt to 80 mt, probabilistic Monte Carlo analysis of 5000 cases for each vehicle were run to determine the final landing ellipse and to quantify the statistical uncertainties associated with the trajectory and attitude conditions during atmospheric entry. Strategies and current technological challenges for a human rated Entry, Descent, and Landing to the Martian surface are presented in this study.

Prabhakar Subrahmanyam↗

System Engineers and Decisions: It?s All about Knowledge

In order to guarantee that a system meets adequate levels of reliability and availability, system performances are continuously monitored and analyzed thanks to the technological advancements driving the Industry 4.0 revolution. An Industry 4.0 approach is typically based on advanced statistical, big data mining, machine learning, and internet-of-things methods designed to detect anomalies in the behavior of system, detect the most likely failure modes, and provide indications to system engineers on when maintenance activities should be performed before system performance are deemed unacceptable (which can be generated by diagnostic and prognostic methods). However, these analyses, which are designed to automatize and increase the efficacy of the system maintenance program, require large amount of data which can come in various forms: numeric, textual, images, sounds etc. Such data constitutes the historic knowledge benchmark to track system performances and support system engineer decisions. Here we claim that data is not sufficient to support this kind of analyses when applied to systems characterized by complex architectures and behaviors. Robust system engineer decisions require the ability to understand the system operational context that lies behind the observed data elements. In this respect, system models are in fact necessary to “put data in context” and capture relationships between data elements. Industry 4.0 methods require in fact contextual knowledge as a basis upon which hypotheses can be generated and assumptions tested. In our view, for complex systems, model-based system engineering (MBSE) models can afford this contextual knowledge, as they are typically used to describe systems architecture and dynamic behaviors. System knowledge is here intended as the blending of collected data and system architecture which takes the form of a “knowledge graph”. A knowledge graph is a database which consists of a large set of nodes (in our case an entity can be either a data or an MBSE element) which are linked to each other. The types of nodes and links follow a pre-defined topology, sometimes also refers as an ontology, that is designed to fit the actual decisions that needs to be performed. We show here how a knowledge graph can be defined to support system engineer maintenance decisions and how the same graph can be built based on system MBSE models and pre-processed data from numeric (through anomaly detections and diagnostic methods) and textual elements (through technical language processing TLP).

97 - MATHEMATICS AND COMPUTING↗

Statistical Properties of SGR 1900+14 Bursts

We study the statistics of soft gamma repeater (SGR) bursts using a database of 187 events detected with BATSE and 837 events detected with the Rossi X-Ray Timing Explorer Proportional Counter Array: all events are from SGR 1900+14 during its 1998-1999 active phase. We find that the fluence or energy distribution of bursts is consistent with a power law of index 1.66, over 4 orders of magnitude. This scale-free distribution resembles the Gutenberg-Richter law for earthquakes and gives evidence for self-organized criticality in SGRS. The distribution of time intervals between successive bursts from SGR 1900+14 is consistent with a lognormal distribution. There is no correlation between burst intensity and the waiting times till the next burst, but there is some evidence for a correlation between burst intensity and the time elapsed since the previous burst. We also find a correlation between the duration and the energy of the bursts, but with significant scatter. In all these statistical properties, SGR bursts resemble earthquakes and solar flares more closely than they resemble any known accretion-powered or nuclear-powered phenomena. Thus, our analysis lends support to the hypothesis that the energy source for SGR bursts is internal to the neutron star and plausibly magnetic.

Gogus, Ersin↗

Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Accurate detection of successive breakup and coalescence events requires that the snapshot interval be an order of magnitude smaller than the characteristic breakup and coalescence times to capture these successive events while minimizing the identification of repeated confounding events. Together, these algorithms serve as a toolbox for detailed analysis of two-phase simulations, and enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Statistics of pressure fluctuations in turbulent kinetic plasmas

ABSTRACT In this study, we explore the statistics of pressure fluctuations in kinetic collisionless turbulence. A 2.5D kinetic particle-in-cell simulation of decaying turbulence is used to investigate pressure balance via the evolution of thermal and magnetic pressure in a plasma with β of order unity. We also discuss the behaviour of thermal, magnetic, and total pressure structure functions and their corresponding wavenumber spectra. The total pressure spectrum exhibits a slope of −7/3 extending for about a decade in the ion-inertial range. In contrast, shallower −5/3 spectra are characteristic of the magnetic pressure and thermal pressure. The steeper total pressure spectrum is a consequence of cancellation caused by density-magnetic field magnitude anti-correlation. Further, we evaluate higher order total pressure structure functions in an effort to discuss intermittency and compare the power exponents with higher order structure functions of velocity and magnetic fluctuations. Finally, applications to astrophysical systems are also discussed.

Astronomy & Astrophysics↗

Galaxy Zoo DESI: large-scale bars as a secular mechanism for triggering AGNs

ABSTRACT Despite the evidence that supermassive black holes (SMBHs) co-evolve with their host galaxy, and that most of the growth of these SMBHs occurs via merger-free processes, the underlying mechanisms which drive this secular co-evolution are poorly understood. We investigate the role that both strong and weak large-scale galactic bars play in mediating this relationship. Using 48 871 disc galaxies in a volume-limited sample from Galaxy Zoo DESI, we analyse the active galactic nucleus (AGN) fraction in strongly barred, weakly barred, and unbarred galaxies up to $z = 0.1$ over a range of stellar masses and colours. After controlling for stellar mass and colour, we find that the optically selected AGN fraction is $31.6 \pm 0.9$ per cent in strongly barred galaxies, $23.3 \pm 0.8$ per cent in weakly barred galaxies, and $14.2 \pm 0.6$ per cent in unbarred disc galaxies. These are highly statistically robust results, strengthening the tantalizing results in earlier works. Strongly barred galaxies have a higher fraction of AGNs than weakly barred galaxies, which in turn have a higher fraction than unbarred galaxies. Thus, while bars are not required in order to grow an SMBH in a disc galaxy, large-scale galactic bars appear to facilitate AGN fuelling, and the presence of a strong bar makes a disc galaxy more than twice as likely to host an AGN than an unbarred galaxy at all galaxy stellar masses and colours.

Astronomy & Astrophysics↗

Measurement of the production and elliptic flow of (anti)nuclei in Xe-Xe collisions at √sNN = 5.44 TeV

Measurements of (anti)deuteron and (anti) 3He production in the rapidity range |𝑦|<0.5 as a function of the transverse momentum and event multiplicity in Xe-Xe collisions at a center-of-mass energy per nucleon-nucleon pair of √𝑠𝑁⁢𝑁=5.44 TeV are presented. The coalescence parameters 𝐵2 and 𝐵3 are measured as a function of the transverse momentum per nucleon. The ratios between (anti)deuteron and (anti) 3He yields and those of (anti)protons and pions are reported as a function of the mean charged-particle multiplicity density and compared with two implementations of the statistical hadronization model and with coalescence predictions. The elliptic flow of (anti)deuterons is measured for the first time in Xe-Xe collisions and shows features similar to those already observed in Pb-Pb collisions, i.e., the mass ordering at low transverse momentum and the meson-baryon grouping at intermediate transverse momentum. The production of nuclei is particularly sensitive to the chemical freeze-out temperature of the system created in the collision, which is extracted from a grand-canonical-ensemble-based thermal fit, performed for the first time including light nuclei along with light-flavor hadrons in Xe-Xe collisions. The extracted chemical freeze-out temperature 𝑇chem=(154.2±1.1) MeV in Xe-Xe collisions is similar to that observed in Pb-Pb collisions and close to the crossover temperature predicted by lattice quantum chromodynamics calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Diffusion flame in homologous turbulent shear flows.

The simplifed statistical theory developed previously is employed to analyze the equilibrium and near-equilibrium combustion of initially unmixed reactants. It is found that the flame zone in the limit of large Damkohler numbers is very thick and is of the order of the local integral scale of turbulence. This is in contrast to the existing phenomenological theories which predict the infinitesimally thin flame sheet, in the same limit, as it is with the laminar diffusion flame. Qualitative agreements with the available experimental results are shown. It is found that singularities exist at the edges of the flame which are removed as Damkohler number is reduced. Also, it is found that the heat transfer may take place against the local mean temperature gradient in certain regions within the flame.

Chung, P. M.↗