Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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,639 records · Page 91

Direct Acceleration of Pickup Ions at The Solar Wind Termination Shock: The Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

Ellison, Donald C.↗

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗

Modification at the C2'-O-Position with 2-Methylbenzothiophene Induces Unique Structural Changes and Thermal Transitions on Duplexes of RNA and DNA

Oligonucleotides can be chemically modified for a variety of applications that include their use as biomaterials, in therapeutics, or as tools to understand biochemical processes, among others. This work focuses on the functionalization of oligonucleotides of RNA and DNA (12- or 14-nucleotides long) with methylbenzothiophene (BT), at the C2'-O-position, which led to unique structural features. Circular dichroism (CD) analyses showed that positioning the BT units on one strand led to significant thermal destabilization, while duplexes where each strand contained 4-BT rings formed a distinct arrangement with cooperativity/interactions among the modifications (evidenced from the appearance of a band with positive ellipticity at 235 nm). Interestingly, the structural arrays displayed increased duplex stabilization (>10 °C higher than the canonical analogue) as a function of [Na + ] with an unexpected structural rearrangement at temperatures above 50 °C. Density functional theory–polarizable continuum model (DFT-PCM) calculations were carried out, and the analyses were in agreement with induced structural changes as a function of salt content. A model was proposed where the hydrophobic surface allows for an internal nucleobase rearrangement into a more thermodynamically stable structure, before undergoing full denaturation, with increased heat. While this behavior is not common, B- to Z-form duplex transitions can occur and are dependent on parameters that were probed in this work, i.e., temperature, nature of modification, or ionic content. To take advantage of this phenomenon, we probed the ability of the modified duplexes to be recognized by Zα (an RNA binding protein that targets Z-form RNA) via electrophoretic analysis and CD. Interestingly, the protein did not bind to canonical duplexes of DNA or RNA; however, it recognized the modified duplexes, in a [monovalent/divalent salt] dependent manner. Overall, the findings describe methodology to attain unique structural motifs of modified duplexes of DNA or RNA, and control their behavior as a function of salt concentration. While their affinity to RNA binding proteins, and the corresponding mechanism of action, requires further exploration, the tunable properties can be of potential use to study this, and other, types of modifications. The novel arrays that formed, under the conditions described herein, provide a useful way to explore the structure and behavior of modified oligonucleotides, in general.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laboratory modeling and analysis of aircraft-lightning interactions

Modeling studies of the interaction of a delta wing aircraft with direct lightning strikes were carried out using an approximate scale model of an F-106B. The model, which is three feet in length, is subjected to direct injection of fast current pulses supplied by wires, which simulate the lightning channel and are attached at various locations on the model. Measurements are made of the resulting transient electromagnetic fields using time derivative sensors. The sensor outputs are sampled and digitized by computer. The noise level is reduced by averaging the sensor output from ten input pulses at each sample time. Computer analysis of the measured fields includes Fourier transformation and the computation of transfer functions for the model. Prony analysis is also used to determine the natural frequencies of the model. Comparisons of model natural frequencies extracted by Prony analysis with those for in flight direct strike data usually show lower damping in the in flight case. This is indicative of either a lightning channel with a higher impedance than the wires on the model, only one attachment point, or short streamers instead of a long channel.

Turner, C. D.↗

Climate Drives Modeled Forest Carbon Cycling Resistance and Resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond (resistance) and recover from (resilience) to rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the ED-2.2 model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900, and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes, and interacted with one another. We then calculated metrics of functional ecosystem stability. Modeled gross primary production exhibited a mean resistance of -0.13, -0.23, and -0.42 (unitless) in response to increasing disturbance severity of 45%, 65%, and 85% mortality, respectively; and resilience of 0.002, 0.004, 0.008 (yr-1) respectively. This pattern held true for other carbon fluxes, meaning that in general simulated GPP declined more initially but recovered more quickly from higher mortality. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

54 ENVIRONMENTAL SCIENCES↗

2‐Aminobenzenethiol‐Functionalized Silver‐Decorated Nanoporous Silicon Photoelectrodes for Selective CO 2 Reduction

Abstract A molecularly thin layer of 2‐aminobenzenethiol (2‐ABT) was adsorbed onto nanoporous p‐type silicon (b‐Si) photocathodes decorated with Ag nanoparticles (Ag NPs). The addition of 2‐ABT alters the balance of the CO 2 reduction and hydrogen evolution reactions, resulting in more selective and efficient reduction of CO 2 to CO. The 2‐ABT adsorbate layer was characterized by Fourier transform infrared (FTIR) spectroscopy and modeled by density functional theory calculations. Ex situ X‐ray photoelectron spectroscopy (XPS) of the 2‐ABT modified electrodes suggests that surface Ag atoms are in the +1 oxidation state and coordinated to 2‐ABT via Ag−S bonds. Under visible light illumination, the onset potential for CO 2 reduction was −50 mV vs. RHE, an anodic shift of about 150 mV relative to a sample without 2‐ABT. The adsorption of 2‐ABT lowers the overpotentials for both CO 2 reduction and hydrogen evolution. A comparison of electrodes functionalized with different aromatic thiols and amines suggests that the primary role of the thiol group in 2‐ABT is to anchor the NH 2 group near the Ag surface, where it serves to bind CO 2 and also to assist in proton transfer.

Kan, Miao↗

Addressing bias in bagging and boosting regression models

As artificial intelligence (AI) becomes widespread, there is increasing attention on investigating bias in machine learning (ML) models. Previous research concentrated on classification problems, with little emphasis on regression models. This paper presents an easy-to-apply and effective methodology for mitigating bias in bagging and boosting regression models, that is also applicable to any model trained through minimizing a differentiable loss function. Our methodology measures bias rigorously and extends the ML model's loss function with a regularization term to penalize high correlations between model errors and protected attributes. We applied our approach to three popular tree-based ensemble models: a random forest model (RF), a gradient-boosted model (GBT), and an extreme gradient boosting model (XGBoost). We implemented our methodology on a case study for predicting road-level traffic volume, where RF, GBT, and XGBoost models were shown to have high accuracy. Despite high accuracy, the ML models were shown to perform poorly on roads in minority-populated areas. Our bias mitigation approach reduced minority-related bias by over 50%.

97 MATHEMATICS AND COMPUTING↗

Capacitor Voltage Ripple Estimation and Optimal Sizing of Modular Multi-Level Converters for Variable-Speed Drives

This article deals with the sizing of submodule capacitors for modular multilevel converters (MMCs) used for variable speed drives for medium voltage and high voltage electric machines. Variable speed drives are required to run over a wide frequency range from zero to a few hundred Hz or even over 1 kHz. Under low frequency operation (<30 Hz), a large current ripple flowing through the MMC submodule capacitors will cause significant capacitor voltage fluctuation. This article proposes an optimal submodule capacitor selection method based on the common-mode voltage injection. The proposed sizing model includes two functions: 1) accurately estimating the capacitor voltage ripple based on the drive system parameters; and 2) identifying the minimum capacitance of the submodule capacitors based on the allowed maximum capacitor voltage fluctuations. Simulation and experimental results show that the error between the actual and estimated values of the voltage fluctuation is within 3% using the proposed model.

42 ENGINEERING↗

Reuseable Objects Software Environment (ROSE): Introduction to Air Force Software Reuse Workshop

The Reusable Objects Software Environment (ROSE) is a common, consistent, consolidated implementation of software functionality using modern object oriented software engineering including designed-in reuse and adaptable requirements. ROSE is designed to minimize abstraction and reduce complexity. A planning model for the reverse engineering of selected objects through object oriented analysis is depicted. Dynamic and functional modeling are used to develop a system design, the object design, the language, and a database management system. The return on investment for a ROSE pilot program and timelines are charted.

Cottrell, William L.↗

Opacity probability distribution functions for electronic systems of CN and C2 molecules including their stellar isotopic forms.

The basis and techniques are presented for generating opacity probability distribution functions for the CN molecule (red and violet systems) and the C2 molecule (Swan, Phillips, Ballik-Ramsay systems), two of the more important diatomic molecules in the spectra of carbon stars, with a view to including these distribution functions in equilibrium model atmosphere calculations. Comparisons to the CO molecule are also shown. T he computation of the monochromatic absorption coefficient uses the most recent molecular data with revision of the oscillator strengths for some of the band systems. The total molecular stellar mass absorption coefficient is established through fifteen equations of molecular dissociation equilibrium to relate the distribution functions to each other on a per gram of stellar material basis.

Querci, F.↗

X-Ray Constraints on the Warm-Hot Intergalactic Medium

Three observational constraints can be placed on a warm-hot intergalactic medium (WHIM) using ROSAT Position Sensitive Proportional Counter (PSPC) pointed and survey data, the emission strength, the energy spectrum, and the fluctuation spectrum. The upper limit to the emission strength of the WHIM is 7.5 +/- 1.0 keV/(s*sq cm*sr*keV) in the 3/4 keV band, an unknown portion of which value may be due to our own Galactic halo. The spectral stape of the WHIM emission can be described as thermal emission with logT = 6.42, although the true spectrum is more likely to come from a range of temperatures. The values of emission strength and spectral shape are in reasonable agreement with hydrodynamical cosmological models. The autocorrelation function in the 0.44 keV < E < 1.21 keV band range, w(theta), for the extragalactic soft X-ray background (SXRB) which includes both the WHIM and contributions due to point sources, is approx. < 0.002 for 10 min < 0 < 20 min in the 3/4 keV band. This value is lower than the Croft et al. (2000) cosmological model by a factor of approx. 5, but is still not inconsistent with cosmological models. It is also found that the normalization of the extragalactic power law component of the soft X-ray background spectrum must be 9.5 +/- 0.9 keV/(s*sq cm*sr*keV) to be consistent with the ROSAT All-Sky Survey.

Kuntz, K. D.↗

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗

Identification of aerodynamic indicial functions using flight data

It is pointed out that the use of indicial function representation provides a model superior to the aerodynamic derivative model. Specific derivatives can be approximated from the indicial models. The model can also be used to compute equivalent stability and control parameters not usually available from flight data. It is shown that derivatives regarding the angle-of-attack and the side slip angle can be derived directly from the indicial functions without any identifiability problem. Attention is given to the pitch moment coefficient, linear indicial function representation, the identification problem for the pitch moment equation, the identifiability of linear systems, parametric representations of the indicial functions, an identification technique, angle-of-attack and pitch rate dynamics in the pitch plane, multivariate linear models, nonlinear aerodynamic indicial functions, measurement system accuracy, and poststall and spin-entry data from a scaled research vehicle.

Gupta, N. K.↗

Quantifying the effect of 3D models on moment tensor results using synthetic data

Moment tensors provide vital information on seismic source properties for seismic events. Moment tensors require seismic wavespeed models to compute Green’s functions, which measure the impulse response between a given source and receiver. Traditionally, researchers have used one-dimensional velocity models to calculate Green’s functions since 1D Green’s functions are computationally cheap to compute. Local 1D velocity models can also accurately model waveforms at short distances (< 500 km). However, 1D velocity models do not account for lateral heterogeneity, which can cause significant misfit in tectonically complex regions such as the Middle East (Covellone and Savage, 2012). Green’s functions calculated using 3D seismic wavespeed models have been shown to perform better in tectonically complex regions (Covellone and Savage, 2012; Kintner and Modrak, 2022), so we are interested in quantifying the effect of considering 3D structure on moment tensor inversion results. The Middle East is an ideal study area for a synthetic moment tensor test for two reasons. Firstly, the Middle East is a tectonically complex region that has been heavily studied. Secondly, the Middle East has significant seismic activity throughout the region, but imperfect station coverage due to limited open data through large swaths of the domain. The tectonic complexity and uneven station coverage will test real-world performance even in a synthetic experiment.

58 GEOSCIENCES↗

The computation of generalized cross-validation functions through householder tridiagonalization with applications to the fitting of interaction spline models

An efficient algorithm for computing the generalized cross-validation function for the general cross-validated regularization/smoothing problem is provided. This algorithm is appropriate for problems where no natural structure is available, and the regularization/smoothing problem is solved (exactly) in a reproducing kernel Hilbert space. It is particularly appropriate for certain multivariate smoothing problems with irregularly spaced data, and certain remote sensing problems, such as those that occur in meteorology, where the sensors are arranged irregularly. The algorithm is applied to the fitting of interaction spline models with irregularly spaced data and two smoothing parameters; favorable timing results are presented. The algorithm may be extended to the computation of certain generalized maximum likelihood (GML) functions. Application of the GML algorithm to a problem in numerical weather forecasting, and to a broad class of hypothesis testing problems, is noted.

Gu, Chong↗

Turbulence and the Stabilization Principle

Further results of research, reported in several previous NASA Tech Briefs articles, were obtained on a mathematical formalism for postinstability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence). To recapitulate: Fictitious control forces are introduced to couple the dynamical equations with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in ordinary perceived three-dimensional space is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. Con sequently, the postinstability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable. The previously reported findings are analyzed from the perspective of the authors Stabilization Principle, according to which (1) stability is recognized as an attribute of mathematical formalism rather than of underlying physics and (2) a dynamical system that appears unstable when modeled by differentiable functions only can be rendered stable by modifying the dynamical equations to incorporate intrinsic stochasticity.

Zak, Michail↗