Search NASA⌕ Search

SEARCH · Search NASA

Results for “Differentiability”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Physics-based compact modeling of electro-thermal memristors: Negative differential resistance, local activity, and non-local dynamical bifurcations

Leon Chua's Local Activity theory quantitatively relates the compact model of an isolated nonlinear circuit element, such as a memristor, to its potential for desired dynamical behaviors when externally coupled to passive elements in a circuit. However, the theory's use has often been limited to potentially unphysical toy models and analyses of small-signal linear circuits containing pseudo-elements (resistors, capacitors, and inductors), which provide little insight into required physical, material, and device properties. Furthermore, the Local Activity concept relies on a local analysis and must be complemented by examining dynamical behavior far away from the steady-states of a circuit. In this work, we review and study a class of generic and extended one-dimensional electro-thermal memristors (i.e., temperature is the sole state variable), re-framing the analysis in terms of physically motivated definitions and visualizations to derive intuitive compact models and simulate their dynamical behavior in terms of experimentally measurable properties, such as electrical and thermal conductance and capacitance and their derivatives with respect to voltage and temperature. Within this unified framework, we connect steady-state phenomena, such as negative differential resistance, and dynamical behaviors, such as instability, oscillations, and bifurcations, through a set of dimensionless nonlinearity parameters. In particular, we reveal that the reactance associated with electro-thermal memristors is the result of a phase shift between oscillating current and voltage induced by the dynamical delay and coupling between the electrical and thermal variables. We thus, demonstrate both the utility and limitations of local analyses to understand non-local dynamical behavior. Critically for future experimentation, the analyses show that external coupling of a memristor to impedances within modern sourcing and measurement instruments can dominate the response of the total circuit, making it impossible to characterize the response of an uncoupled circuit element for which a compact model is desired. However, these effects can be minimized by proper understanding of the Local Activity theory to design and utilize purpose-built instruments.

Brown, Timothy D. (ORCID:0000000347068355)↗
DOE OSTI2023

Structural inference of networked dynamical systems with universal differential equations

Networked dynamical systems are common throughout science in engineering; e.g., biological networks, reaction networks, power systems, and the like. For many such systems, nonlinearity drives populations of identical (or near-identical) units to exhibit a wide range of nontrivial behaviors, such as the emergence of coherent structures (e.g., waves and patterns) or otherwise notable dynamics (e.g., synchrony and chaos). Here, we seek to infer (i) the intrinsic physics of a base unit of a population, (ii) the underlying graphical structure shared between units, and (iii) the coupling physics of a given networked dynamical system given observations of nodal states. These tasks are formulated around the notion of the Universal Differential Equation, whereby unknown dynamical systems can be approximated with neural networks, mathematical terms known a priori (albeit with unknown parameterizations), or combinations of the two. We demonstrate the value of these inference tasks by investigating not only future state predictions but also the inference of system behavior on varied network topologies. The effectiveness and utility of these methods are shown with their application to canonical networked nonlinear coupled oscillators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Non-destructive depth-resolved characterization of residual strain fields in high electron mobility transistors using differential aperture x-ray microscopy

Localized residual stress and elastic strain concentrations in microelectronic devices often affect the electronic performance, resistance to thermomechanical damage, and, likely, radiation tolerance. A primary challenge for the characterization of these concentrations is that they exist over sub-μm length-scales, precluding their characterization by more traditional residual stress measurement techniques. Here, in this study, we demonstrate the use of synchrotron x-ray-based differential aperture x-ray microscopy (DAXM) as a viable, non-destructive means to characterize these stress and strain concentrations in a depth-resolved manner. DAXM is used to map two-dimensional strain fields between the source and the drain in a gallium nitride (GaN) layer within high electron mobility transistors (HEMTs) with sub-μm spatial resolution. Strain fields at various positions in both pristine and irradiated HEMT specimens are presented in addition to a preliminary stress analysis to estimate the distribution of various stress components within the GaN layer. γ-irradiation is found to significantly reduce the lattice plane spacing in the GaN along the sample normal direction, which is attributed to radiation damage in transistor components bonded to the GaN during irradiation.

42 ENGINEERING↗
DOE OSTI2023

Memory-based parameterization with differentiable solver: Application to Lorenz ’96

Physical parameterizations (or closures) are used as representations of unresolved subgrid processes within weather and global climate models or coarse-scale turbulent models, whose resolutions are too coarse to resolve small-scale processes. These parameterizations are typically grounded on physically based, yet empirical, representations of the underlying small-scale processes. Machine learning-based parameterizations have recently been proposed as an alternative solution and have shown great promise to reduce uncertainties associated with the parameterization of small-scale processes. Yet, those approaches still show some important mismatches that are often attributed to the stochasticity of the considered process. This stochasticity can be due to coarse temporal resolution, unresolved variables, or simply to the inherent chaotic nature of the process. To address these issues, we propose a new type of parameterization (closure), which is built using memory-based neural networks, to account for the non-instantaneous response of the closure and to enhance its stability and prediction accuracy. We apply the proposed memory-based parameterization, with differentiable solver, to the Lorenz ’96 model in the presence of a coarse temporal resolution and show its capacity to predict skillful forecasts over a long time horizon of the resolved variables compared to instantaneous parameterizations. This approach paves the way for the use of memory-based parameterizations for closure problems.

Mathematics↗
DOE OSTI2025

X-ray spectroscopy of multi-temperature plasmas using the differential emission measure formalism

We present a theoretical construct that nominally underlies spectroscopic data analysis of multi-temperature plasmas, known as the differential emission measure (DEM). From a data analytic perspective, the DEM formalism is used to derive temperature distributions from line spectra that are formed in the presence of temperature gradients and by time integrations of evolving plasmas. From a modeling perspective, DEMs are convenient intermediaries between radiation hydrodynamics simulations and spectroscopic measurements acquired in the laboratory. The DEM concept and its associated methodologies were originally developed by spectroscopists working with astrophysical data. We borrow from these earlier investigations. In this manuscript, intended primarily as a tutorial, we discuss the basic concepts, but also augment various aspects of the theory by the way of extension and example, including a detailed treatment of various weighting and averaging schemes, intended to mitigate ambiguities that often arise when reporting temperature information. We focus on high-temperature plasmas that are not in local thermodynamic equilibrium and the x-ray spectra that they produce, although the core ideas presented here are applicable to spectroscopy in other energy bands. A few examples involving the derivation and manipulation of model DEMs in simple geometries are provided.

Liedahl, Duane A. [Lawrence Livermore National Lab↗
DOE OSTI2026

A flexible and differentiable coil proxy for stellarator equilibrium optimization

Balancing plasma performance and coil cost is a significant challenge when designing a stellarator power plant. Most current stellarator designs are produced through two-stage optimization: stage-1 for the equilibrium and stage-2 for a coil design that reproduces its magnetic configuration. Because few proxies connect both stages, two-stage optimization can produce plasmas that have high-quality physical properties but overly complex coils. In recent years, single-stage optimization has increasingly been used to optimize the plasma and coils simultaneously in order to improve the plasma–coil balance. However, all existing single-stage tools are specialized for filament coils, cannot model coil systems containing permanent magnets or dipole arrays, and continue to be challenged by numerical problems. The quasi-single-stage optimization finds a middle ground by integrating a coil optimization subproblem into stage-1 optimization. We present a flexible, differentiable coil complexity proxy based on the newly developed QUADCOIL coil optimization code. QUADCOIL is fast and can target realistic coil metrics and constraints that are unavailable to codes with comparable speed. We demonstrate the effectiveness and flexibility of the QUADCOIL proxy by presenting two QSS optimization studies. The first study produces a permanent magnet solution for the MUSE stellarator with 34% fewer magnets than previous solutions. The second study produces a coil solution for the ARIES-CS stellarator with 30% reductions in both peak and root-mean-square force.

Fu, Lanke [Princeton Plasma Physics Laboratory (PP↗
DOE OSTI2016

Quantization of the Sobolev space of half-differentiable functions

A quantization of the Sobolev space V=H{sub 0}{sup 1/2}(S{sup 1},R) of half- differentiable functions on the circle, which is closely connected with string theory, is constructed. The group QS(S{sup 1}) of quasisymmetric circle homeomorphisms acts on V by reparametrizations, but this action is not smooth. Nevertheless, a quantum infinitesimal action of QS(S{sup 1}) on V can be defined, which enables one to construct a quantum algebra of observables which is associated with the system (V,QS(S{sup 1})). Bibliography: 7 titles. (paper)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Adaptive differentiation and rapid evolution of a soil bacterium along a climate gradient

Significance Increasing evidence suggests that evolutionary processes frequently shape ecological patterns; however, most microbiome studies thus far have focused on only the ecological responses of these communities. By using parallel field experiments and focusing in on a model soil bacterium, we showed that bacterial “species” are differentially adapted to local climates, leading to changes in their composition. Furthermore, we detected strain-level evolution, providing direct evidence that both ecological and evolutionary processes operate on annual timescales. The consideration of eco-evolutionary dynamics may therefore be important to understand the response of soil microbiomes to future environmental change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Crystal structure of the WFIKKN2 follistatin domain reveals insight into how it inhibits growth differentiation factor 8 (GDF8) and GDF11

Growth differentiation factor 8 (GDF8; also known as myostatin) and GDF11 are closely related members of the transforming growth factor β (TGF-β) family. GDF8 strongly and negatively regulates skeletal muscle growth, and GDF11 has been implicated in various age-related pathologies such as cardiac hypertrophy. GDF8 and GDF11 signaling activities are controlled by the extracellular protein antagonists follistatin; follistatin-like 3 (FSTL3); and WAP, follistatin/kazal, immunoglobulin, Kunitz, and netrin domain–containing (WFIKKN). All of these proteins contain a follistatin domain (FSD) important for ligand binding and antagonism. Here, we investigated the structure and function of the FSD from murine WFIKKN2 and compared it with the FSDs of follistatin and FSTL3. Using native gel shift and surface plasmon resonance analyses, we determined that the WFIKKN2 FSD can interact with both GDF8 and GDF11 and block their interactions with the type II receptor activin A receptor type 2B (ActRIIB). Further, we solved the crystal structure of the WFIKKN2 FSD to 1.39 Å resolution and identified surface-exposed residues that, when substituted with alanine, reduce antagonism of GDF8 in full-length WFIKKN2. Furthermore, comparison of the WFIKKN2 FSD with those of follistatin and FSTL3 revealed differences in both the FSD structure and position of residues within the domain that are important for ligand antagonism. Taken together, our results indicate that both WFIKKN and follistatin utilize their FSDs to block the type II receptor but do so via different binding interactions.

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

Stable calcium-free myocilin olfactomedin domain variants reveal challenges in differentiating between benign and glaucoma-causing mutations

Nonsynonymous gene mutations can be beneficial, neutral, or detrimental to the stability, structure, and biological function of the encoded protein, but the effects of these mutations are often not readily predictable. For example, the β-propeller olfactomedin domain of myocilin (mOLF) exhibits a complex interrelationship among structure(s), stability, and aggregation. Numerous mutations within mOLF are linked to glaucoma; the resulting variants are less stable, aggregation-prone, and sequestered intracellularly, causing cytotoxicity. In this study, we report the first stable mOLF variants carrying substitutions in the calcium-binding site that exhibit solution characteristics indistinguishable from those of glaucoma variants. Crystal structures of these stable variants at 1.8–2.0-Å resolution revealed features that we could not predict by molecular dynamics simulations, including loss of loop structure, helix unwinding, and a blade shift. Double mutants that combined a stabilizing substitution and a selected glaucoma-causing single-point mutant rescued in vitro folding and stability defects. In the context of full-length myocilin, secretion of stable single variants was indistinguishable from that of the WT protein, and the double mutants were secreted to varying extents. In summary, our finding that mOLF can tolerate particular substitutions that render the protein stable despite a conformational switch emphasizes the complexities in differentiating between benign and glaucoma-causing variants and provides new insight into the possible biological function of myocilin.

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

On-the-Fly Generation of Differential Resonance Scattering Probability Distribution Functions for Monte Carlo Codes

Current Monte Carlo codes use one of three models: (1) the asymptotic scattering model, (2) the free gas scattering model, or (3) the S(α,β) model, depending on the neutron energy and the specific Monte Carlo code. This thesis addresses the consequences of using the free gas scattering model, which assumes that the neutron interacts with atoms in thermal motion in a monatomic gas in thermal equilibrium at material temperature, T. Most importantly, the free gas model assumes the scattering cross section is constant over the neutron energy range, which is usually a good approximation for light nuclei, but not for heavy nuclei where the scattering cross section may have several resonances in the epithermal region. Several researchers in the field have shown that the exact resonance scattering model is temperaturedependent, and neglecting the resonances in the lower epithermal range can under-predict resonance absorption due to the upscattering phenomenon mentioned above, leading to an over-prediction of keff by several hundred pcm. Existing methods to address this issue involve changing the neutron weights or implementing an extra rejection scheme in the free gas sampling scheme, and these all involve performing the collision analysis in the center-of-mass frame, followed by a conversion back to the laboratory frame to continue the random walk of the neutron. The goal of this paper was to develop a sampling methodology that (1) accounted for the energydependent scattering cross sections in the collision analysis and (2) was performed in the laboratory frame,avoiding the conversion to the center-of-mass frame. The energy dependence of the scattering cross section was modeled with even-ordered polynomials (2nd and 4th order) to approximate the scattering cross section in Blackshaw’s equations for the moments of the differential scattering PDFs. These moments were used to sample the outgoing neutron speed and angle in the laboratory frame on-the-fly during the random walk of the neutron. Results for criticality studies on fuel pin and fuel assembly calculations using methods developed in this dissertation showed very close comparison to results using the reference Dopplerbroadened rejection correction (DBRC) scheme.

42 ENGINEERING↗
DOE OSTI2021

Measurements of Neutron Scattering from a Copper Sample Using a Quasi-Differential Method in the Region from 2 keV to 20 MeV

Neutron scattering from a copper sample was measured at Rensselaer Polytechnic Institute utilizing the quasi-differential method. The measurement spanned the energy range from 0.5 to 20 MeV using the high-energy scattering system and from 2 keV to 0.5 MeV using the new mid-energy scattering system. Copper was selected as a material of interest to measure due to large discrepancies between experiments and simulations of the Zeus benchmark. The Zeus benchmark consists of a copper reflected highly enriched uranium system, and the angular distribution of copper scattering was thought to potentially be the cause of the discrepancy. The copper measurements found differences in the scattering response particularly in the incident energy region from 1 to 2 MeV for the high-energy measurement and from 2 to 4 keV in the mid-energy system. These differences are particularly noticeable at angles near 90 deg in the high-energy system and back angles in the mid-energy system. Additionally, for ENDF/B-VIII.0 there is a large discrepancy at the forward angle in the energy range around 0.5 MeV. For these reasons, a new evaluation of copper scattering utilizing these results is recommended and perhaps could help to improve the agreement with the Zeus benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Element differentiation with a Hartmann- based X-ray phase imaging system

Significant efforts are currently ongoing in X-ray imaging to provide multimodal imaging systems, targeting better sensitivity and specificity for both biomedical or non-destructive testing applications. Knowing the elemental composition of specific structures, such as breast microcalcifications in mammography, would help to differentiate malign and benign tumours. Standard X-ray Phase Contrast Imaging techniques provide only qualitative information on elements with similar absorption properties. However, their chemical composition can be determined from the measurement of the phase as it is directly related to the optical index of elemental materials. We will show new experimental results obtained with an X-ray phase imaging system based on a Hartmann mask. Early data treatment succeeded in retrieving both the real and imaginary parts of the refractive index. The system demonstrates the capability to discriminate materials based on elemental composition.

36 MATERIALS SCIENCE↗
DOE OSTI2015

TESTING MICROVARIABILITY IN QUASAR DIFFERENTIAL LIGHT CURVES USING SEVERAL FIELD STARS

Microvariability consists of small timescale variations of low amplitude in the photometric light curves of quasars and represents an important tool to investigate their inner core. Detection of quasar microvariations is challenging because of their non-periodicity, as well as the need for high monitoring frequency and a high signal-to-noise ratio. Statistical tests developed for the analysis of quasar differential light curves usually show either low power or low reliability, or both. In this paper we compare two statistical procedures to perform tests on several stars with enhanced power and high reliability. We perform light curve simulations of variable quasars and non-variable stars and analyze them with statistical procedures developed from the F-test and the analysis of variance. The results show a large improvement in the power of both statistical probes and a larger reliability when several stars are included in the analysis. The results from the simulations agree with those obtained from observations of real quasars. The high power and high reliability of the tests discussed in this paper improve the results that can be obtained from short and long timescale variability studies. These techniques are not limited to quasar variability; on the contrary, they can be easily implemented for other sources, such as variable stars. Their applications to future research and to the analysis of large-field photometric monitoring archives could reveal new variable sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Negative differential mobility and trapping in active matter systems

Using simulations, we explore the average velocity as a function of applied drift force for active matter particles moving through a random obstacle array. We find that for low drift force, there is an initial flow regime where the mobility increases linearly with drive, while for higher drift forces a regime of negative differential mobility appears in which the velocity decreases with increasing drive due to the trapping of active particles behind obstacles. A fully clogged regime exists at very high drift forces when all the particles are permanently trapped behind obstacles. We discover for increasing activity that the overall mobility is nonmonotonic, with an enhancement of the mobility for small levels of activity and a decrease in mobility for large activity levels. We show how these effects evolve as a function of disk and obstacle density, active run length, drift force, and motor force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

Harnessing mass differential confinement effects in magnetized rotating plasmas to address new separation needs

The ability to separate large volumes of mixed species based on atomic mass appears desirable for a variety of emerging applications with high societal impact. One possibility to meet this objective consists in leveraging mass differential effects in rotating plasmas. Beyond conventional centrifugation, rotating plasmas offer in principle additional ways to separate elements based on mass. Single ion orbits show that ion radial mass separation in a uniform magnetized plasma column can be achieved by applying a tailored electric potential profile across the column, or by driving a rotating magnetic field within the column. Furthermore, magnetic pressure and centrifugal effects can be combined in a non-uniform geometry to separate ions based on mass along the field lines. In conclusion, practical application of these separation schemes hinges on the ability to produce the desirable electric and magnetic field configuration within the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗