Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Unraveling Electronic Trends in O* and OH* Surface Adsorption in the MO 2 Transition-Metal Oxide Series

Understanding the bond strength of O* and OH* intermediates to metal-oxide surfaces is key to predicting the catalytic activity in oxygen-based electrochemistry. Here, we uncover highly non-linear trends in O* and OH* adsorption energies across the 3d, 4d, and 5d series of MO 2 transition-metal (TM) oxide surfaces computed within Hubbard- U corrected density functional theory (DFT + U ). Investigating the electronic structure with crystal orbital Hamiltonian populations (COHP) of the relevant metal–oxygen bonds reveals that the spin-dependent coupling strength between metal-d and oxygen-2p atomic orbitals together with the extent of filling of bonding and anti-bonding orbitals are the primary contributors to the adsorption energy. Importantly, we show that the integrated COHP obtained purely from bulk calculations is a highly accurate descriptor for surface adsorption energetics that captures trends across the group 5–12 TM oxide series within 0.19–0.36 eV. Our results suggest a pathway to prediction of adsorption energies for an arbitrary metal–ligand catalyst system.

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

On the Assessment of Acoustic Scattering and Shielding by Time Domain Boundary Integral Equation Solutions

Based on the time domain boundary integral equation formulation of the linear convective wave equation, a computational tool dubbed Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) has recently been under development. The time domain approach has a distinct advantage that the solutions at all frequencies are obtained in a single computation. In this paper, the formulation of the integral equation, as well as its stabilization by the Burton-Miller type reformulation, is extended to cases of a constant mean flow in an arbitrary direction. In addition, a "Source Surface" is also introduced in the formulation that can be employed to encapsulate regions of noise sources and to facilitate coupling with CFD simulations. This is particularly useful for applications where the noise sources are not easily described by analytical source terms. Numerical examples are presented to assess the accuracy of the formulation, including a computation of noise shielding by a thin barrier motivated by recent Historical Baseline F31A31 open rotor noise shielding experiments. Furthermore, spatial resolution requirements of the time domain boundary element method are also assessed using point per wavelength metrics. It is found that, using only constant basis functions and high-order quadrature for surface integration, relative errors of less than 2% may be obtained when the surface spatial resolution is 5 points-per-wavelength (PPW) or 25 points-per-wavelength squared (PPW2).

Hu, Fang Q.↗
DOE OSTI2025

Structure and Synthesizability of Iron–Sulfur Metal–Organic Frameworks

Sulfur-based metal–organic frameworks (MOFs) and coordination polymers (CPs) are an emerging class of hybrid materials that have received growing attention due to their magnetic, conductive, and catalytic properties with potential applications in electrocatalysis and energy storage. In this work, we report a high-throughput virtual screening protocol to predict the synthesizability of candidate metal–sulfur MOFs/CPs by computing the thermodynamically stable structures resulting from a particular combination of metal cluster, linker, cation, and synthetic conditions. Free energies are computed by using all-atom classical mechanical thermodynamic integration. Low-free-energy structures are refined using ab initio density functional theory, and pair distribution functions and powder X-ray diffraction patterns are calculated to complement and guide experimental structure determination. We validate the computational approach by retrospective predictions of the stable structure produced by experimental syntheses, and a subsequent screen predicts Fe 4 S 4 -BDT–TPP as a new thermodynamically stable one-dimensional (1D) CP comprising a redox-active Fe 4 S 4 cluster, a 1,4-benzenedithiolate (BDT) linker, and a tetraphenylphosphonium (TPP) countercation. Furthermore, this material is experimentally synthesized, and the 1D chain structure of the crystal is confirmed using microcrystal electron diffraction. The computational screening pipeline is generically transferable to neutral and ionic MOFs/CPs comprising arbitrary metal clusters, linkers, cations, and synthetic conditions, and we make it freely available as an open source tool to guide and accelerate the discovery and engineering of novel porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Social network structure and the spread of complex contagions from a population genetics perspective

Ideas, behaviors, and opinions spread through social networks. If the probability of spreading to a new individual is a non-linear function of the fraction of the individuals’ affected neighbors, such a spreading process becomes a “complex contagion”. This non-linearity does not typically appear with physically spreading infections, but instead can emerge when the concept that is spreading is subject to game theoretical considerations (e.g. for choices of strategy or behavior) or psychological effects such as social reinforcement and other forms of peer influence (e.g. for ideas, preferences, or opinions). Here we study how the stochastic dynamics of such complex contagions are affected by the underlying network structure. Motivated by simulations of complex contagions on real social networks, we present a framework for analyzing the statistics of contagions with arbitrary non-linear adoption probabilities based on the mathematical tools of population genetics. The central idea is to use an effective lower-dimensional diffusion process to approximate the statistics of the contagion. This leads to a tradeoff between the effects of ”selection” (microscopic tendencies for an idea to spread or die out), random drift, and network structure. Our framework illustrates intuitively several key properties of complex contagions: stronger community structure and network sparsity can significantly enhance the spread, while broad degree distributions dampen the effect of selection compared to random drift. Finally, we show that some structural features can exhibit critical values that demarcate regimes where global contagions become possible for networks of arbitrary size. Our results draw parallels between the competition of genes in a population and memes in a world of minds and ideas. Our tools provide insight into the spread of information, behaviors, and ideas via social influence, and highlight the role of macroscopic network structure in determining their fate.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1994

Near-surface ice on Mercury and the Moon: A topographic thermal model

A thermal model that can be easily adapted to craters of arbitrary shape is developed and applied to high-latitude impact craters on Mercury and the Moon, Chao Meng Fu crater at -87.5 deg L on Mercury, an unnamed bowl-shaped crater at 86.7 deg L on Mercury, and Peary crater at 88.6 deg L on the Moon. For an assumed input topography and grid of surface elements, the model computes for each element the irradiation from direct insolation and reflected and emitted radiation from other elements, taking into account shadowing by walls of the crater, partial obscuration of the solar disk near the poles and the diurnal, orbital, and seasonal cycles. Temperatures are computed over the surface grid as functions of depth and time from the surface to a specified depth and over the pertinent astronomical cycles, including the effects of direct and indirect surface irradiation, infrared radiation, heat conduction, and interior heating. Vapor fluxes and ice recession times are computed as functions of ice depth over the surface grid. Temperatures profiles, vapor fluxes, and ice recession times were computed for flat surfaces not associated with craters near the poles of Mercury and the Moon. It was found that water ice could have existed throughout geologic time within the maximum radar detection depth of recent observation of Mercury (J. K. Harmon and M. A. Slade, 1992, Science 258, 640-643) poleward of approximately 87 - 88 deg L on Mercury and poleward of approximately 73 deg L on the Moon. For Chao Meng Fu crater it was found that approximately 40% of the crater floor is permanently shadowed from direct solar insolation, while the remainder of the crater floor is periodically illuminated by a partially obscured Sun. Temperatures at the upper levels of the south wall can slightly exceed 550 K. Surface temperatures in the permanently shadowed region of the crater floor are under approximately 130 K, which could have allowed water ice to exist throughout geologic time within the radar detection depth of recent observation of Mercury. For small bowl-shaped crater on Mercury, it was found that most of the crater is permanently shadowed from direct solar radiation, except for a narrow semicircular band bordering the north rim. However, temperatures in the permanently shadowed region periodically reach a maximum near approximately 315 K due to efficient heating of the small crater by thermal emission and reflection from the small sunlit region, which periodically reaches temperatures exceeding 630 K. Water ice could not have existed throughout geologic time anywhere in this crater within the radar detection depth. For Peary crater on the Moon, the entire crater floor is permanently shadowed from direct solar insolation with maximum temperature under 120 K. The upper level of the north wall periodically reaches a maximum temperature near 310 K. The low temperatures on the crater floor would have allowed water ice to exist near the surface throughout geologic time, provided that the Moon's obliquity was always as low as it is at present.

Salvail, James R.↗
DOE OSTI2025

Investigating parameter trainability in the SNAP-displacement protocol of a qudit system

Here, in this study, we explore the universality of Selective Number-dependent Arbitrary Phase (SNAP) and Displacement gates for quantum control in qudit-based systems. However, optimizing the parameters of these gates poses a challenging task. Our main focus is to investigate the sensitivity of training any of the SNAP parameters in the SNAP-Displacement protocol. We analyze conditions that could potentially lead to the Barren Plateau problem in a qudit system and draw comparisons with multi-qubit systems. The parameterized ansatz we consider consists of blocks, where each block is composed of hardware operations, namely SNAP and Displacement gates [Fösel et al 2020 Efficient cavity control with snap gates arXiv:2004.14256]. Applying Variational Quantum algorithm (VQA) with observable and gate cost functions, we utilize techniques similar to those in [McClean et al 2018 Barren plateaus in quantum neural network training landscapes Nat. Commun.9 1–6] and [Cerezo et al 2021 Cost function dependent barren plateaus in shallow parametrized quantum circuits. Nat. Commun.12 1–12] along with the concept of t − design. Through this analysis, we make the following key observations: (a) The trainability of a SNAP-parameter does not exhibit a preference for any particular direction within our cost function landscape, (b) By leveraging the first and second moments properties of Haar measures, we establish new lemmas concerning the expectation of certain polynomial functions, and (c) utilizing these new lemmas, we identify a general condition that indicates an expected trainability advantage in a qudit system when compared to multi-qubit systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

A local dynamic model for large eddy simulation

The dynamic model is a method for computing the coefficient C in Smagorinsky's model for the subgrid-scale stress tensor as a function of position from the information already contained in the resolved velocity field rather than treating it as an adjustable parameter. A variational formulation of the dynamic model is described that removes the inconsistency associated with taking C out of the filtering operation. This model, however, is still unstable due to the negative eddy-viscosity. Next, three models are presented that are mathematically consistent as well as numerically stable. The first two are applicable to homogeneous flows and flows with at least one homogeneous direction, respectively, and are, in fact, a rigorous derivation of the ad hoc expressions used by previous authors. The third model in this set can be applied to arbitrary flows, and it is stable because the C it predicts is always positive. Finally, a model involving the subgrid-scale kinetic energy is presented which attempts to model backscatter. This last model has some desirable theoretical features. However, even though it gives results in LES that are qualitatively correct, it is outperformed by the simpler constrained variational models. It is suggested that one of the constrained variational models should be used for actual LES while theoretical investigation of the kinetic energy approach should be continued in an effort to improve its predictive power and to understand more about backscatter.

Ghosal, Sandip↗
NASA NTRS1991

An area model for on-chip memories and its application

An area model suitable for comparing data buffers of different organizations and arbitrary sizes is described. The area model considers the supplied bandwidth of a memory cell and includes such buffer overhead as control logic, driver logic, and tag storage. The model gave less than 10 percent error when verified against real caches and register files. It is shown that, comparing caches and register files in terms of area for the same storage capacity, caches generally occupy more area per bit than register files for small caches because the overhead dominates the cache area at these sizes. For larger caches, the smaller storage cells in the cache provide a smaller total cache area per bit than the register set. Studying cache performance (traffic ratio) as a function of area, it is shown that, for small caches, direct-mapped caches perform significantly better than four-way set-associative caches and, for caches of medium areas, both direct-mapped and set-associative caches perform better than fully associative caches.

Mulder, Johannes M.↗
NASA NTRS1994

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗
NASA NTRS1989

The Hole-Count Test Revisited: Effects of Test Specimen Thickness

For historical reasons the hole count, an important performance test for the Analytical Electron Microscope (AEM), is somewhat arbitrary yielding different numbers for different investigators. This was not a problem a decade ago when AEM specimens were often bathed with large fluxes of stray electrons and hard x rays. At that time the presence or absence of a thick Pt second condenser (C2) aperture could be detected by a simple comparison of the x-ray spectrum taken 'somewhere in the hole' with a spectrum collected on a 'typical thickness' of Mo or Ag foil. A high hole count of about 10-20% indicated that the electron column needed modifications; whereas a hole count of 1-2% was accepted for most AEM work. The absolute level of the hole count is a function of test specimen atomic number, overall specimen shape, and thin-foil thickness. In order that equivalent results may be obtained for any AEM in any laboratory in the world, this test must become standardized. The hole-count test we seek must be as simpl and as nonsubjective as the graphite 0.344nm lattice-line-resolution test. This lattice-resolution test spurred manufacturers to improve the image resolution of the TEM significantly in the 1970s and led to the even more stringent resolution tests of today. A similar phenomenon for AEM instruments would be welcome. The hole-count test can also indicate whether the spurious x-ray signal is generated by high-energy continuum x rays (bremsstrahlung) generated in the electron column (high K-line to L-line ratio) or uncollimated electrons passing through or around the C2 aperture (low K/L ratio).

Lyman, C. E.↗
DOE OSTI2026

Dimensional Reduction for Sampled Priors and Application to Photometric Redshift Distributions

A typical Bayesian inference on the values of some parameters of interest q from some data D involves running a Markov Chain (MC) to sample from the posterior $p$($q$,$n$|$D$) $\propto$ $\mathcal{L}$($D$|$q$,$n$)$p$(q)$p$($n$), where n are some nuisance parameters with a separable prior. In some cases, the nuisance parameters are high-dimensional, and their prior p(n) is itself defined only by a set of samples that have been drawn from some other MC. The MC for the posterior will typically require evaluation of p(n) at arbitrary values of n, i.e., one needs to provide a density estimator over the full n space from the provided samples. But the high dimensionality of n hinders both the density estimation and the efficiency of the MC for the posterior. We describe a solution to this problem: a linear compression of the n space into a much lower-dimensional space u, which projects away directions in n space that cannot appreciably alter $\mathcal{L}$. The algorithm for doing so is a slight modification to principal components analysis, and is less restrictive on p(n) than other proposed solutions to this issue. We demonstrate this “mode projection” technique using the analysis of 2-point correlation functions of weak lensing fields and galaxy density in the Dark Energy Survey, where n is a binned representation of the redshift distribution n(z) of the galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Apollo 9 multiband photography experiment S065

Fourier analysis was applied to microdensitometer scans of a selected region of one SO65 frame in each of the three black-and-white bands. The approach was unique because a somewhat arbitrary section of the image was used and not limited to available targets as in edge analysis. Comparison of duplicates and calculation of absolute SO65 MTF were done by applying linear systems theory to the spatial spectra of the image scans. It was found that the duplication process was nonlinear and resulted in general amplification of spatial frequency modulation. However, the increase in modulation was offset by a corresponding increase in the granularity of the copies. The amount of increase seemed to be related to the initial granularity, but a direct relationship was not verified. Band-to-band comparison of image quality was achieved in the form of signal-to-noise ratio curves as a function of spatial frequency for each band. From this standpoint the DD band was an order of magnitude better than the other two. These were several factors that restricted the analysis of the Apollo 9 imagery. Among these were the lack of precise sensitometric and optical system data on the high altitude photography. In addition it was determined that there was only one simultaneous pair of high altitude and SO65 frames. Finally, the original SO65 photography was not available for scanning (for obvious reasons), thus eliminating a reference base for granularity and SO65 MTF determination.

Schowengerdt, R. A.↗
NASA NTRS2011

A Self-Stabilizing Synchronization Protocol for Arbitrary Digraphs

This paper presents a self-stabilizing distributed clock synchronization protocol in the absence of faults in the system. It is focused on the distributed clock synchronization of an arbitrary, non-partitioned digraph ranging from fully connected to 1-connected networks of nodes while allowing for differences in the network elements. This protocol does not rely on assumptions about the initial state of the system, other than the presence of at least one node, and no central clock or a centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. There is no theoretical limit on the maximum number of participating nodes. The only constraint on the behavior of the node is that the interactions with other nodes are restricted to defined links and interfaces. This protocol deterministically converges within a time bound that is a linear function of the self-stabilization period. We present an outline of a deductive proof of the correctness of the protocol. A bounded model of the protocol was mechanically verified for a variety of topologies. Results of the mechanical proof of the correctness of the protocol are provided. The model checking results have verified the correctness of the protocol as they apply to the networks with unidirectional and bidirectional links. In addition, the results confirm the claims of determinism and linear convergence. As a result, we conjecture that the protocol solves the general case of this problem. We also present several variations of the protocol and discuss that this synchronization protocol is indeed an emergent system.

Malekpour, Mahyar R.↗
DOE OSTI2022

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

NASA Lewis Launch Collision Probability Model Developed and Analyzed

There are nearly 10,000 tracked objects orbiting the earth. These objects encompass manned objects, active and decommissioned satellites, spent rocket bodies, and debris. They range from a few centimeters across to the size of the MIR space station. Anytime a new satellite is launched, the launch vehicle with its payload attached passes through an area of space in which these objects orbit. Although the population density of these objects is low, there always is a small but finite probability of collision between the launch vehicle and one or more of these space objects. Even though the probability of collision is very low, for some payloads even this small risk is unacceptable. To mitigate the small risk of collision associated with launching at an arbitrary time within the daily launch window, NASA performs a prelaunch mission assurance Collision Avoidance Analysis (or COLA). For the COLA of the Cassini spacecraft, the NASA Lewis Research Center conducted an in-house development and analysis of a model for launch collision probability. The model allows a minimum clearance criteria to be used with the COLA analysis to ensure an acceptably low probability of collision. If, for any given liftoff time, the nominal launch vehicle trajectory would pass a space object with less than the minimum required clearance, launch would not be attempted at that time. The model assumes that the nominal positions of the orbiting objects and of the launch vehicle can be predicted as a function of time, and therefore, that any tracked object that comes within close proximity of the launch vehicle can be identified. For any such pair, these nominal positions can be used to calculate a nominal miss distance. The actual miss distances may differ substantially from the nominal miss distance, due, in part, to the statistical uncertainty of the knowledge of the objects positions. The model further assumes that these position uncertainties can be described with position covariance matrices. With these, and some additional simplifying assumptions, a closed-form solution for the probability of collision is obtained. This solution provides clear insights into how each of the independent parameters affects the probability of collision. It shows that for a given maximum probability of collision and prior knowledge of the objects position uncertainties and sizes, only knowledge of the nominal closest approach distance is required to make the launch/no launch decision. The model was completed and used for the mission assurance COLA analysis for the Cassini spacecraft, which was launched on a Titan IV/Centaur rocket on October 15, 1997. Although the model was specifically developed for the Cassini mission, it is clearly applicable for other launches as well. The effect of COLA closures on the launch window is shown. The bar represents the entire 140-min launch window on October 15, 1997; the blackened areas represent the loss of launch opportunities due to the potential for a collision with an orbiting object.

Bollenbacher, Gary↗
DOE OSTI2022

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗