Sensitivity analysis for Monte Carlo simulations of integro-partial differential equations.
Abstract not provided.
Bond, Stephen D↗
SEARCH · Engineering Papers
Results for “differential analysis”
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.
Precise positioning capabilities for TOPEX using differential GPS
NASA's Ocean Topographic Experiment (TOPEX), to be launched in 1991, is the first mission designed to reach the decimeter accuracy needed for the solution of the general mean circulation problem. An experimental tracking capability for TOPEX is studied using differential measurements with satellites of the U.S. DOD's Global Positioning System (GPS). Two data types are studied: (1) integrated Doppler from GPS carrier phase, and (2) GPS P-code pseudo-range. Results of covariance analysis predict that with differential GPS techniques, 5-10 cm average TOPEX altitude accuracies can be achieved over data arcs of two hours.
Lichten, S. M.↗
The stability of coupled renewal-differential equations with econometric applications
Concepts and results are presented in the fields of mathematical modeling, economics, and stability analysis. A coupled renewal-differential equation structure is presented as a modeling form for systems possessing hereditary characteristics, and this structure is applied to a model of the Austrian theory of business cycles. For realistic conditions, the system is shown to have an infinite number of poles, and conditions are presented which are both necessary and sufficient for all poles to lie strictly in the left half plane.
Rhoten, R. P.↗
Omics-Lethal Human Viruses, Interferon-Stimulated Response Experiment IFNaIHH001
The purpose of this experiment was to evaluate the host interferon-stimulated cellular response to interferon alpha (IFNα) treatment. Sample data was obtained from human immortalized human hepatocyte cells (IHH) treated with or without IFNα for mRNA and miRNA transcriptome expression analysis.
59 BASIC BIOLOGICAL SCIENCES↗
Lithium loss, resistance growth, electrode expansion, gas evolution, and Li plating: Analyzing performance and failure of commercial large-format NMC-Gr lithium-ion pouch cells
Not Available
Gasper, Paul (ORCID:0000000188349458)↗
Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints
The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.
97 MATHEMATICS AND COMPUTING↗
Differentiable Simulation of a Liquid Argon Time Projection Chamber
Liquid argon time projection chambers (LArTPCs) are widely used in particle detection for their tracking and calorimetric capabilities. The particle physics community actively builds and improves high-quality simulators for such detectors in order to develop physics analyses in a realistic setting. The ability of these simulators to mimic real, measured data is limited by the modeling of the physical detectors used for data collection. This modeling can be improved by performing dedicated calibration measurements. Conventional approaches calibrate individual detector parameters or processes one at a time. However, the impact of detector processes is entangled, making this a poor description of the underlying physics. We introduce a differentiable simulator that enables a gradient-based optimization, allowing for the first time a simultaneous calibration of all detector parameters. We describe the procedure of making a differentiable simulator, highlighting the challenges of retaining the physics quality of the standard, non-differentiable version while providing meaningful gradient information. We further discuss the advantages and drawbacks of using our differentiable simulator for calibration. Finally, we provide a starting point for extensions to our approach, including applications of the differentiable simulator to physics analysis pipelines.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Direct preemptive stabilization of $m,n$ = 2,1 neoclassical tearing modes by electron cyclotron current drive in the DIII-D low-torque ITER baseline scenario
We report the first preemptive direct stabilization of $m,n$ = 2,1 neoclassical tearing modes (NTMs) in DIII-D low-torque ITER baseline scenario plasmas by electron cyclotron waves. These experiments strongly support that the observed stabilization is achieved by replacing the missing bootstrap current in the island by electron cyclotron current drive (ECCD). These plasmas (with and without ECCD) are stable to 2,1 NTMs when the differential rotation between the q = 1 and q = 2 surfaces (Δf 1,2 ) is sustained above a critical value (Δf CRIT ≈ 1 kHz), but those evolving to low differential rotation (Δf 1,2 < Δf CRIT ) routinely develop a disruptive 2,1 island at the time and frequency of a sawtooth precursor. Preemptive, local and direct stabilization by ECCD is tested by scanning the current drive amplitude and location inside and outside of the q = 2 rational surface shot by shot. Analysis of only low differential rotation time windows, i.e. when stabilization from Δf 1,2 is absent, shows focused ECCD at q = 2 prevents the onset of 2,1 NTMs. These results are important, as they give the first demonstration of disruptive NTM control by ECCD in the plasma scenario planned for ITER.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods
Background: Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. Methods: In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, the median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.
Finsam Samson↗
The role of stiffness in training and generalization of ResNets
Neural ordinary differential equations (NODEs) have recently regained popularity as large-depth limits of a large class of neural networks. In particular, residual neural networks (ResNets) are equivalent to an explicit Euler discretization of an underlying NODE, where the transition from one layer to the next is one time step of the discretization. The relationship between continuous and discrete neural networks has been of particular interest. Notably, analysis from the ordinary differential equation viewpoint can potentially lead to new insights for understanding the behavior of neural networks in general. In this work, we take inspiration from differential equations to define the concept of stiffness for a ResNet via the interpretation of a ResNet as the discretization of a NODE. Here, we then examine the effects of stiffness on the ability of a ResNet to generalize, via computational studies on example problems coming from climate and chemistry models. We find that penalizing stiffness does have a unique regularizing effect, but we see no benefit to penalizing stiffness over L 2 regularization (penalization of network parameter norms) in terms of predictive performance.
97 MATHEMATICS AND COMPUTING↗
An Empirical Model of Coronal Heating Deduced from SOHO Observations of Transition Region Dynamics
The research carried out under the above award was concerned with the analysis of "explosive events" in the solar atmosphere. Explosive events are small energy releases that reveal themselves as skewed, broadened line profiles in the 10(exp 5) K degree material of the solar transition region. Using a data set of such events from the Solar Ultraviolet Measurements of Emitted Radiation (SUMER) instrument on the NASA SoHo satellite, with unprecedented temporal, spatial and spectral resolution, we have evaluated (Winebarger, Emslie, Mariska, & Warren 1999, Ap. J., 526, 471) the energetics of these events through application of the Velocity Differential Emission Measure (VDEM) analysis technique introduced by Newton, Emslie & Mariska, (1995, Ap. J, 447, 915). This work has allowed us to evaluate the contribution of various terms (such as enthalpy flux, free-streaming of high-energy electrons, etc.) to the energy budget of these events.
Emslie, A. Gordon↗
First Extraction of Single Differential Cross-Sections on Argon for CC1$\mu$2p0$\pi$ Event Topologies in the MicroBooNE Detector
The precise measurement of cross-sections for a variety of interactions is critical to the success of forthcoming flagship neutrino experiments. Of special interest are neutrino interactions that leave the nucleus in a 2-particle 2-hole state (2p2h). This note will present the first cross-section measurement for the production of 2p2h states on argon. Using MicroBooNE data collected from the first three years of operation, we select events corresponding to a charged-current $ν_µ$ interaction that left the argon nucleus in a 2p2h state. These interactions produce a topology with one muon and two protons in the final state (CC1$\mu$2p0$\pi$). This analysis targets a single differential cross-section measurement for CC1$\mu$2p0$\pi$ event topologies in a variety of kinematic variables. Comparisons are made to a set of theoretical models that explore different cross-section modeling configurations.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Identification of LiH and nanocrystalline LiF in the solid–electrolyte interphase of lithium metal anodes
A comprehensive understanding of the solid–electrolyte interphase (SEI) composition is crucial to developing high-energy batteries based on lithium metal anodes. A particularly contentious issue concerns the presence of LiH in the SEI. In this work, we report on the use of synchrotron-based X-ray diffraction and pair distribution function analysis to identify and differentiate two elusive components, LiH and LiF, in the SEI of lithium metal anodes. LiH is identified as a component of the SEI in high abundance, and the possibility of its misidentification as LiF in the literature is discussed. LiF in the SEI is found to have different structural features from LiF in the bulk phase, including a larger lattice parameter and a smaller grain size (<3 nm). These characteristics favour Li + transport and explain why an ionic insulator, like LiF, has been found to be a favoured component for the SEI. Finally, pair distribution function analysis reveals key amorphous components in the SEI.
36 MATERIALS SCIENCE↗
Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories
Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.
entry↗
Multi-Body Analysis of a Tiltrotor Configuration
The paper describes the aeroelastic analysis of a tiltrotor configuration. The 1/5 scale wind tunnel semispan model of the V-22 tiltrotor aircraft is considered. The analysis is performed by means of a multi-body code, based on an original formulation. The differential equilibrium problem is stated in terms of first order differential equations. The equilibrium equations of every rigid body are written, together with the definitions of the momenta. The bodies are connected by kinematic constraints, applied in form of Lagrangian multipliers. Deformable components are mainly modelled by means of beam elements, based on an original finite volume formulation. Multi-disciplinar problems can be solved by adding user-defined differential equations. In the presented analysis the equations related to the control of the swash-plate of the model are considered. Advantages of a multi-body aeroelastic code over existing comprehensive rotorcraft codes include the exact modelling of the kinematics of the hub, the detailed modelling of the flexibility of critical hub components, and the possibility to simulate steady flight conditions as well as wind-up and maneuvers. The simulations described in the paper include: 1) the analysis of the aeroelastic stability, with particular regard to the proprotor/pylon instability that is peculiar to tiltrotors, 2) the determination of the dynamic behavior of the system and of the loads due to typical maneuvers, with particular regard to the conversion from helicopter to airplane mode, and 3) the stress evaluation in critical components, such as the pitch links and the conversion downstop spring.
Ghiringhelli, G. L.↗
Observations of Mantle Seismic Anisotropy Using Array Techniques: Shear-Wave Splitting of Beamformed SmKS Phases
Shear-wave splitting measurements are commonly used to resolve seismic anisotropy in both the upper and lowermost mantle. Typically, such techniques are applied to SmKS phases that have reflected (m-1) times off the underside of the core-mantle boundary before being recorded. Practical constraints for shear-wave splitting studies include the limited number of suitable phases as well as the large fraction of available data discarded because of poor signal-to-noise ratios (SNRs) or large measurement uncertainties. Array techniques such as beamforming are commonly used in observational seismology to enhance SNRs, but have not been applied before to improve SmKS signal strength and coherency for shear wave splitting studies. Here, we investigate how a beamforming methodology, based on slowness and backazimuth vespagrams to determine the most coherent incoming wave direction, can improve shear-wave splitting measurement confidence intervals. Through the analysis of real and synthetic seismograms, we show that (a) the splitting measurements obtained from the beamformed seismograms (beams) reflect an average of the single-station splitting parameters that contribute to the beam; (b) the beams have (on average) more than twice as large SNRs than the single-station seismograms that contribute to the beam; (c) the increased SNRs allow the reliable measurement of shear wave splitting parameters from beams down to average single-station SNRs of 1.3. Beamforming may thus be helpful to more reliably measure splitting due to upper mantle anisotropy. Moreover, here we show that beamforming holds potential to greatly improve detection of lowermost mantle anisotropy by demonstrating differential SKS–SKKS splitting analysis using beamformed USArray data.
54 ENVIRONMENTAL SCIENCES↗
Omics-Lethal Human Viruses, Ebola Experiment EHUH001
The purpose of this experiment was to evaluate the human host response to Zaire Ebola wild-type virus and mutant virus infection. Samples were obtained from human hepatoma carcinoma cells (HUH-7) infected with Zaire Ebola ΔVP30-WT background, in addition to mutant viruses Ebola-ΔsGP (eliminating expression of soluble glycoprotein ssGP) and Ebola-Δmucin (encoding a glycoprotein lacking the mucin domain) for mRNA and miRNA transcriptome expression analysis. Secondary host-associated viral dataset downloads contain one or more statistically processed (normalization data transformation) quantitative dataset collections resulting in qualitative expression analyses of primary host-pathogen experimental study designs. Transcriptomics dataset downloads have a direct relationship to a primary sample submission corresponding to a specific Ebola virus infection.
59 BASIC BIOLOGICAL SCIENCES↗
Omics-Lethal Human Viruses, Ebola Experiment EHUH003
The purpose of this experiment was to evaluate the human host response to Zaire Ebola wild-type virus infection. Samples were obtained from human hepatoma carcinoma cells (HUH-7) infected with Zaire Ebola ΔVP30-WT background for mRNA and miRNA transcriptome expression analysis. Secondary host-associated viral dataset downloads contain one or more statistically processed (normalization data transformation) quantitative dataset collections resulting in qualitative expression analyses of primary host-pathogen experimental study designs. Transcriptomics dataset downloads have a direct relationship to a primary sample submission corresponding to a specific Ebola virus infection.
59 BASIC BIOLOGICAL SCIENCES↗