Search NASA⌕ Search

SEARCH · Search NASA

Results for “residual function”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Implicit application of polynomial filters in a k-step Arnoldi method

The Arnoldi process is a well known technique for approximating a few eigenvalues and corresponding eigenvectors of a general square matrix. Numerical difficulties such as loss of orthogonality and assessment of the numerical quality of the approximations as well as a potential for unbounded growth in storage have limited the applicability of the method. These issues are addressed by fixing the number of steps in the Arnoldi process at a prescribed value k and then treating the residual vector as a function of the initial Arnoldi vector. This starting vector is then updated through an iterative scheme that is designed to force convergence of the residual to zero. The iterative scheme is shown to be a truncation of the standard implicitly shifted QR-iteration for dense problems and it avoids the need to explicitly restart the Arnoldi sequence. The main emphasis of this paper is on the derivation and analysis of this scheme. However, there are obvious ways to exploit parallelism through the matrix-vector operations that comprise the majority of the work in the algorithm. Preliminary computational results are given for a few problems on some parallel and vector computers.

Sorensen, D. C.↗
DOE OSTI2022

Comparison of PsbQ and Psb27 in photosystem II provides insight into their roles

Photosystem II (PSII) catalyzes the oxidation of water at its active site that harbors a high-valent inorganic Mn 4 CaO x cluster called the oxygen-evolving complex (OEC). Extrinsic subunits generally serve to protect the OEC from reductants and stabilize the structure, but diversity in the extrinsic subunits exists between phototrophs. Recent cryo-electron microscopy experiments have provided new molecular structures of PSII with varied extrinsic subunits. We focus on the extrinsic subunit PsbQ, that binds to the mature PSII complex, and on Psb27, an extrinsic subunit involved in PSII biogenesis. PsbQ and Psb27 share a similar binding site and have a four-helix bundle tertiary structure, suggesting they are related. Here, we use sequence alignments, structural analyses, and binding simulations to compare PsbQ and Psb27 from different organisms. We find no evidence that PsbQ and Psb27 are related despite their similar structures and binding sites. Evolutionary divergence within PsbQ homologs from different lineages is high, probably due to their interactions with other extrinsic subunits that themselves exhibit vast diversity between lineages. This may result in functional variation as exemplified by large differences in their calculated binding energies. Psb27 homologs generally exhibit less divergence, which may be due to stronger evolutionary selection for certain residues that maintain its function during PSII biogenesis which is consistent with their more similar calculated binding energies between organisms. Previous experimental inconsistencies, low confidence binding simulations, and recent structural data suggest that Psb27 is likely to exhibit flexibility that may be an important characteristic of its activity. Furthermore, the analysis provides insight into the functions and evolution of PsbQ and Psb27, and an unusual example of proteins with similar tertiary structures and binding sites that probably serve different roles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Mapping the functional landscape of the receptor binding domain of T7 bacteriophage by deep mutational scanning

The interaction between a bacteriophage and its host is mediated by the phage's receptor binding protein (RBP). Despite its fundamental role in governing phage activity and host range, molecular rules of RBP function remain a mystery. Here, we systematically dissect the functional role of every residue in the tip domain of T7 phage RBP (1660 variants) by developing a high-throughput, locus-specific, phage engineering method. This rich dataset allowed us to cross compare functional profiles across hosts to precisely identify regions of functional importance, many of which were previously unknown. Substitution patterns showed host-specific differences in position and physicochemical properties of mutations, revealing molecular adaptation to individual hosts. We discovered gain-of-function variants against resistant hosts and host-constricting variants that eliminated certain hosts. To demonstrate therapeutic utility, we engineered highly active T7 variants against a urinary tract pathogen. Our approach presents a generalized framework for characterizing sequence–function relationships in many phage–bacterial systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Studying the effects of overlapping objects in dark energy

Observing the clustering of galaxies allows us to calculate cosmological parameters necessary for understanding dark energy. However, as the density of observed objects increases, the probability of these objects blending likewise increases, causing multiple galaxies to be observed as one. This affects the inferred values of parameters such as the galaxy bias (b) and the matter energy density (ΩM). To see whether the bias from incorrectly inferring the galaxy count is significant, we compare the correlation functions in simulated data for “true” and “observed” data sets with one-to-one and multiple-to-one correspondences, respectively. For each data set, we create two correlation functions: one “measured” function directly relying on the galaxies’ positions using the TreeCorr python library, and one “model” derived mathematically from the galaxies’ power spectrum using the Cosmological Core Library (CCL). By minimizing the residual between these two functions, we compute the ideal values for b and ΩM across the various possible redshifts that position the galaxies in three dimensional space. This minimization is done with an Markov chain Monte Carlo (MCMC) estimate that finds one value of ΩM and ten values for b corresponding to the ten redshift bins ranging from z = 0.2 to z = 1.2. We find that neither b nor ΩM is particularly affected by inclusion of blended galaxies. Though there is room for improvement, the data suggests that the fluctuations we found are a result of noise or limitations on the modeling rather than blending explicitly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Effect of Process History on Grain Boundaries and Dislocation Substructures on Functional Properties of Nb for SRF Cavities: Plastic Formability and Microstructural Evolution. Final report

While much progress has been made in improving the performance of cavities in the past decade, reproducible performance remains elusive. This renewal of “The cost of grain boundaries on the performance of superconducting cavities” examines in further depth the interactions between forming and annealing processes and functional performance at the scale of dislocation substructure and grain boundaries. There are five tasks identified in the renewal, which are: Characterization / analysis of current cavity technology: Effect of deformation, heat treatment, surface treatments, and residual damage on cavity function (Lee and Bieler) Characterization of dislocation mechanisms in single and bicrystals (Bieler) Quantifying effects of impurities and deformation mechanisms in Nb using density functional theory calculations (Solanki) Crystal plasticity constitutive model development (Pourboghrat and Eisenlohr) Modeling and measuring thermal conductivity in thin layers (Wright)

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of Polymer Rheology on Gel Treatment Performance of Horizontal Wells with Severe Channeling

Summary Gel treatment has been a cost-effective method to control the conformance of a reservoir with severe heterogeneity problems. The water channels in such reservoirs can be classified as open fractures or high permeability porous media with pore-throat network. Many simulation studies have been conducted to discuss gel treatment performance for conformance control. However, nobody considered the polymer rheology difference in open fractures and porous media in simulation. Previous simulation studies also ignored the residual resistance factor as a function of rock permeability rather than a constant parameter. In this study, a conceptual simulation model was established to simulate the linear flow system for the reservoir with horizontal wells considering the two factors mentioned above. The results demonstrate that the gel treatment always provides the better placement results in the open fracture type channel than pore-throat network type channel. Moreover, it is very necessary to consider residual resistance factor as a function of permeability, which is based on the experimental results and can provide much greater plugging efficiency in the higher permeable channels than constant residual resistance factor. Sensitivity analysis studies and multifactor analysis indicate that increasing oil viscosity and permeability ratio has a strong positive influence on conformance control results, which indicate in-situ gel treatment can be better applied in heavy or viscous oil reservoirs with fracture-like channels. Besides, the results also indicate that in reservoirs with severe channeling problem where channel velocity was high, the differences of gelant placement and profile improvement in models with two different types of channels could be enlarged greatly.

Engineering↗
DOE OSTI2023

Targeting Aspergillus allergen oryzin with a chemical probe at atomic precision

We report the molecular basis of Aspergillus fumigatus oryzin, allergen Asp f 13, or alkaline proteinase ALP1, containing the sequence motif His–Asp–Ser of the subtilisin family, structure, and function at atomic detail. Given the resolution of the data (1.06 Å), we use fragment molecular replacement with ideal polyalanine α-helices to determine the first crystal structure of oryzin. We probe the catalytic serine through formation of an irreversible bond to a small molecule compound, specifically labeling it, describing the amino acid residues performing the catalytic function. Defined by a self-processed pro-peptide, the active site architecture shapes up pocket-like subsites that bind to and unveil the S1′–S4′ substrate binding preferences. We use molecular modeling to dock a model of the pro-peptide in the S1–S4 region and to dock collagen along the active site cleft. Opposite to the face harboring the catalytic serine, the enzyme binds to a calcium ion in a binding site created by backbone flipping. We use thermal unfolding to show that this metal ion provides structural stability. With no known host inhibitor identified thus far, this structure may hasten the progress of developing new therapeutic agents for diseases caused by pathogenic fungi.

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

Axial residual stresses in boron fibers

The axial residual stress distribution as a function of radius was determined from the fiber surface to the core including the average residual stress in the core. Such measurements on boron on tungsten (B/W) fibers show that the residual stresses for 102, 142, 203, and 366 micron diameter fibers were similar, being compressive at the surface and changing monotonically to a region of tensile within the boron. At approximately 25 percent of the original radius, the stress reaches a maximum tensile stress of about 860 mn/sq.m and then decreases to a compressive stress near the tungsten boride core. Data were presented for 203 micron diameter B/W fibers that show annealing above 900 C reduces the residual stresses. A comparison between 102 micron diameter B/W and boron on carbon (b/C) shows that the residual stresses were similar in the outer regions of the fibers, but that large differences near and in the core were observed. The effects of these residual stresses on the fracture of boron fibers were discussed.

Behrendt, D. R.↗
NASA NTRS1978

Axial residual stresses in boron fibers

A method of measuring axial residual stresses in boron fibers is presented. With this method, the axial residual stress distribution as a function of radius is determined from the fiber surface to the core including the average residual stress in the core. Such measurements on boron on tungsten (B/W) fibers show that the residual stresses for 102, 142, 203, and 366 micron diam fibers are similar, being compressive at the surface and changing monotonically to a region of tensile stress within the boron. At approximately 25% of the original radius, the stress reaches a maximum tensile stress of about 860 MN sq m and then decreases to a compressive stress near the tungsten boride core. Data are presented for 203-micron diam B/W fibers that show annealing above 900 C reduces the residual stresses. A comparison between 102-micron diam B/W and boron on carbon (B/C) show that the residual stresses are similar in the outer regions of the fibers, but that large differences near and in the core are observed. Fracture of boron fibers is discussed.

Behrendt, D. R.↗
NASA NTRS1993

Fabrication and properties of functionally graded NiAl/Al2O3 composites

A modified sedimentation process was used in the production of a functionally gradient material (FGM), NiAl/Al2O3. A simple finite element model was used to guide our design and fabrication efforts by estimating residual stress states as a function of composite structure. This approach could lead to tailored designs that enhance or avoid specific residual stress states. Thermal cycling tests were factored into the model to predict time dependent or steady-state internal temperature and stress profiles. Four-point bend tests were conducted to establish the mechanical load-displacement behavior of a single interlayer FGM at room temperature, 800 and 1000 K. Room temperature bend strength of the FGM was 3-4 times that of the base NiAl. At elevated temperatures, composite fracture occurred in a gradual, noncatastrophic mode involving NiAl retardation of a succession of cracks originating in the alumina face.

Miller, D. P.↗
NASA NTRS2017

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗
DOE OSTI2025

Functional and Structural Studies on the Esperamicin Thioesterase and Progress toward Understanding Enediyne Core Biosynthesis

Enediynes are among the most potent antitumor and antibacterial natural products. Studies on their biosynthetic pathways have identified a shared, linear polyene precursor generated from an iterative type I polyketide synthase (PKSE) as the source of the enediyne warhead. A key step is the release of this polyene from the PKSE by a discrete thioesterase (TE). Here, in this study, we used X-ray crystallography, site-directed mutagenesis, and heterologous coexpression of PKSEs and TEs to elucidate how enediyne TEs mediate the production of the polyene. We solved the structure of wild-type EspE7 from esperamicin producer Actinomodura verrucosospora. The substrate binding pocket was also defined upon serendipitous cocrystallization of an EspE7 mutant with a fatty acyl-CoA ligand. Structural data and in vitro activity assays with EspE7 mutants provide strong evidence that Glu68 in EspE7 and the analogous Glu residue in other enediyne TEs functions as a key catalytic residue, thus supporting a hydrolysis mechanism for enediyne TEs that aligns with that of Pseudomonas sp. 4-HB-CoA TE. Furthermore, combinations of 9- and 10-membered enediyne PKSEs and TEs produced 1,3,5,7,9,11,13-pentadecaheptaene (1) as the major product. Thus, the data further support previous conclusions that 1 serves as the sole precursor for the biosynthesis of all enediyne cores.

Enediynes↗
NASA NTRSDate not supplied

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↗
DOE OSTI2025

Computational screening of fly ash zeolite sorbents for boric acid removal

In the United States, many impoundments at coal-fired power plants contain elevated contaminants like arsenic, boron, barium, and selenium. Zeolites synthesized from fly ash show promise as sorbents for these contaminants. However, optimizing sorption capacity is challenging due to numerous possible topologies, silicon to aluminum (Si/Al) ratios, and cation types. In this study, molecular simulations are used to design cationic zeolites for boric acid adsorption. Force field models based on quantum mechanical calculations (PBE + D2) for Na-, Ca-, Mn-, and Fe-exchanged chabazite and LTA are presented. The new D2FF force fields reproduce DFT energies with about half the error of UFF. Zeolite performance depends on Si/Al ratio and cation type, with low Si/Al ratio chabazite (CHA) and phillipsite (PHI) zeolite frameworks exchanged with Ca 2+ or Na + /Ca 2+ mixtures showing the highest adsorption. In conclusion, these findings suggest tailored fly ash-derived zeolites could provide effective boron removal from leachate ponds.

CCR impoundment↗
DOE OSTI2024

Double winding condensate CFT

We conjecture a two-dimensional conformal field theory built out of a linear dilaton and a compact β γ system deformed by winding condensates on each of the compact cycles. In particular, explicit expressions of the residues of the correlation functions are provided. As a worldsheet theory, it describes a stringy black hole in AdS 3 (supported by the Kalb-Ramond flux) at the Hawking-Page temperature. It is expected to be connected to the Horowitz-Polchinski-like solution near the Hagedorn temperature of thermal AdS 3 and Euclidean BTZ blackhole as we vary the temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Deformation-enhanced hierarchical multiscale structure heterogeneity in a Pd-Si bulk metallic glass

Here, the multiscale structures in a Pd 82 Si 18 binary bulk metallic glass before and after deformation were studied using electron microscopies, high-energy synchrotron X-ray diffraction, and small-angle scattering techniques. The experimental results revealed an enhancement of hierarchical structure heterogeneities on multiple length scales after deformation. Hierarchical multiple shear bands of high number density were observed after bending, introducing complex but periodically distributed residual strain. Pair distribution function analysis revealed that the connectivity of the short-range clusters on the medium-range scale determines the packing density difference between the tension side and the compression side in the sample after bending. In-situ synchrotron X-ray diffraction study also revealed a transformation of connection modes among short-range clusters under uniaxial tension and compression, which is consistent with those of triaxial tension/compression parts upon bending in Pd 82 Si 18 glassy alloys. The nanoscale heterogeneities for metallic glasses after deformation observed by small-angle scattering and transmission electron microscopy may be attributed to the nanoscale amorphous phase separation and interacting multiple shear bands enhanced by plastic deformation. Our findings suggested that the enhancement of hierarchical heterogeneous structure on multiple length scales may explain the excellent plasticity of Pd-Si glassy alloys, deepening the understanding of structure-property relation during plastic deformation in metallic glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stellarator optimization for good magnetic surfaces at the same time as quasisymmetry

A method is demonstrated to optimize a stellarator's geometry to eliminate magnetic islands and achieve other desired physics properties at the same time. For many physics quantities that have been used in stellarator optimization, including quasisymmetry, neoclassical transport, and magnetohydrodynamic stability, it is convenient to use a magnetic equilibrium representation that assures the existence of magnetic surfaces. However, this representation hides the possible presence of magnetic islands, which are typically undesirable. To include both surface-based objectives and island widths in a single optimization, two fixed-boundary equilibrium calculations are run at each iteration of the optimization: one that enforces the existence of magnetic surfaces (the Variational Moments Equilibrium Code) [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)] and one that does not (the Stepped Pressure Equilibrium Code) [Hudson et al., Phys. Plasmas 19, 112502 (2012)]. By penalizing the island residues in the objective function, the two magnetic field representations are brought into agreement during the optimization. Here, an example is presented in which, particularly on the surface where quasisymmetry was targeted, quasisymmetry is achieved more accurately than in previously published examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗