Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient structure”

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 883 records · Page 49

The physical basis for solvent flow in organic solvent nanofiltration

Organic solvent nanofiltration (OSN) is an emerging membrane technology that could revolutionize chemical separations in numerous vital industries. Despite its significance, there remains a lack of fundamental understanding of solvent transport mechanisms in OSN membranes. Here, we use an extended Flory-Rehner theory, nonequilibrium molecular dynamic simulations, and organic solvent transport experiments to demonstrate that solvent flow in OSN membranes is driven by a pressure gradient. We show that solvent molecules migrate as clusters through interconnected pathways within the membrane pore structure, challenging the widely accepted diffusion-based view of solvent transport in OSN. We further reveal that solvent permeance is dependent on solvent affinity to the OSN membrane, which, in turn, controls the membrane pore structure. Our fundamental insights lay the scientific groundwork for the development of next-generation OSN membranes.

42 ENGINEERING↗

Controlling Molecular Structure and Spin with Multiconfigurational Quantum Chemistry (Final Technical Report)

For many first-row transition metal complexes, structure-property relationships can be obtained from high-level molecular geometry optimizations and subsequent electronic structure studies. The project funded under this award utilized newly implemented fully internally contracted (FIC) nuclear gradients for extended multi-state (XMS) complete active space second-order multireference perturbation theory (CASPT2) to explore geometry changes in first-row transition metal coordination complexes. At the start of the project, only one full geometry optimization using FIC-CASPT2 analytical gradients had been reported (J. Chem. Theory Comput., 2016, 12 (8), 3781), and many open-questions regarding the performance and achievable accuracy in applying such computations to larger complexes persisted. Key deliverables in the report include the implementation of the numerical Hessian and subsequent vibrational analysis in the BAGEL program package. This includes both full Hessian vibrational analysis (FHVA) and a partial Hessian vibrational analysis (PHVA). The latter of which allows us to reduce the significant cost of computing the vibrational frequencies in a large molecule by focusing on the modes of highest interest. The gradient code also proved useful in the development of an approach to improve the Hubbard-U correction by removing bias in predicting spin-splitting in Fe(II) complexes via plane-wave density functional theory (DFT). Finally, we showed for three families of complexes (spin-crossover (SCO) complexes, metallocorroles complexes, and those with metal-metal bonds) that CASPT2 can result in good molecular geometries and established best practices in undertaking this work. We are now using this approach in collaborative projects with experimental groups, which would not have been possible at the start of this project. A common theme also arose through this work that static and dynamic correlation must be recovered in a balanced way to yield quantitative results. By systematically varying how both effects were included, we were able to make key chemical insights. This work supported the training of two postdoctoral scholars, one graduate student, and one undergraduate student in applying high-level wavefunction based methods to challenging systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revisiting Dropout: Escaping Pressure for Training Neural Networks with Multiple Costs

A common approach to jointly learn multiple tasks with a shared structure is to optimize the model with a combined landscape of multiple sub-costs. However, gradients derived from each sub-cost often conflicts in cost plateaus, resulting in a subpar optimum. In this work, we shed light on such gradient conflict challenges and suggest a solution named Cost-Out, which randomly drops the sub-costs for each iteration. We provide the theoretical and empirical evidence of the existence of escaping pressure induced by the Cost-Out mechanism. While simple, the empirical results indicate that the proposed method can enhance the performance of multi-task learning problems, including two-digit image classification sampled from MNIST dataset and machine translation tasks for English from and to French, Spanish, and German WMT14 datasets.

97 MATHEMATICS AND COMPUTING↗

Thermal static bending of deployable interlocked booms

Metal ribbons processed with a heat-forming treatment are enabled to form tubelike structures when deployed from a roll. Deployable booms of this have been utilized for gravity-gradient stabilization on the RAE, ATS, and Nimbus D satellites. An experimental thermal-mechanics test apparatus was developed to measure the thermal static bending and twist of booms up to 3 meters long. The apparatus was calibrated by using the correlation between calculated and observed thermal bending of a seamless tube. Thermal static bending values of 16 interlocked deployable booms were observed to be within a factor of 2.5 of the values calculated from seamless-tube theory. Out-of-Sun-plane thermal bending was caused by complex heat transfer across the interlocked seam. Significant thermal static twisting was not observed.

Staugaitis, C. L.↗

A kinetic-theory approach to turbulent chemically reacting flows

The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

Chung, P. M.↗

Galactic cosmic ray intensity from 1 to 9 A.U

The paper is concerned with comparative observations of galactic cosmic ray intensity from 1 to 9 AU during the period April 1973-January 1976. The dominant intensity variation is found to be a coherent temporal one with a peak-to-peak amplitude of about 35% of the mean, while a further cyclic variation of about 8% peak-to-peak was found to have a 26-day period and is attributed to solar-corotating interplanetary magnetic field structures. Detailed time lag analysis of the intensity ratios for the two spacecraft yields a radial gradient of +2.0 (plus or minus 0.5) percent per AU, integral above a particle energy of 80 MeV/nucleon.

Van Allen, J. A.↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

Redistribution of pulmonary blood flow during unilateral hypoxia in prone and supine dogs

We used fluorescent-labeled microspheres in pentobarbital-anesthetized dogs to study the effects of unilateral alveolar hypoxia on the pulmonary blood flow distribution. The left lung was ventilated with inspired O2 fraction of 1.0, 0.09, or 0.03 in random order; the right lung was ventilated with inspired O2 fraction of 1.0. The lungs were removed, cleared of blood, dried at total lung capacity, then cubed to obtain approximately 1,500 small pieces of lung ( approximately 1.7 cm3). The coefficient of variation of flow increased (P < 0.001) in the hypoxic lung but was unchanged in the hyperoxic lung. Most (70-80%) variance in flow in the hyperoxic lung was attributable to structure, in contrast to only 30-40% of the variance in flow in the hypoxic lung (P < 0.001). When adjusted for the change in total flow to each lung, 90-95% of the variance in the hyperoxic lung was attributable to structure compared with 70-80% in the hypoxic lung (P < 0.001). The hilar-to-peripheral gradient, adjusted for change in total flow, decreased in the hypoxic lung (P = 0.005) but did not change in the hyperoxic lung. We conclude that hypoxic vasoconstriction alters the regional distribution of flow in the hypoxic, but not in the hyperoxic, lung.

NASA Discipline Cardiopulmonary↗

Gradient-Based Optimization of the Common Research Model Wing Subject to CFD-Based Gust and Flutter Constraints

The linearized frequency-domain method was recently implemented in the stabilized finite element solver in NASA’s FUN3D code. Previous work by the authors used this method for enforcing flutter constraints during gradient-based optimizations. More recently, the solver was expanded to account for continuous (also known as stochastic) gust responses. This paper expands on recent Common Research Model wing optimization work, which demonstrated gradient-based optimization with flutter and stochastic gust constraints, among others. While that work utilized FUN3D for static aeroelastic solutions but relied on doublet lattice aerodynamics for gust and flutter responses, the present work replaces these unsteady aerodynamic analyses with those of FUN3D’s linearized frequency-domain solver. With analytic derivatives available, gradient-based optimization is performed through the use of the OpenMDAO/MPhys libraries with over 700 shape, structural, and aerodynamic design variables and over 10 nonlinear constraints. Comparisons of analysis results and optimized designs are made between doublet lattice and linearized frequency-domain solutions.

aeroelasticity↗

Modeling-Based Design and Optimization of a Gradient Composite Transition Joint

An innovative additively manufactured gradient composite transition joint (AM-GCTJ) has been designed to join dissimilar metals, to address the pressing issue of premature failure observed in conventional dissimilar metal welds (DMWs) when subjected to increased cyclic operating conditions of fossil fuel power plants. The transition design, guided by computational modeling, developed a gradient composite material distribution, facilitating a smooth transition in material volume fraction and physical properties between different alloys. This innovative design seeks to alleviate structural challenges arising from distinct material properties, including high thermal stress and potential cracking issues resulting from the thermal expansion mismatch typically observed in conventional DMWs. In this study, we investigated the creep properties of transition joints comprising Grade 91 steel and 304 stainless steel through a combination of simulations and creep testing experiments. The implementation of a gradient composite design in the plate transition joint resulted in a significant enhancement of creep resistance when compared to the baseline conventional DMW. For instance, the creep rupture life of the transition joint was improved by > 400% in a wide range of temperature and stress testing conditions. Meanwhile, the failure location shifted to the base material of Grade 91 steel. Such enhancement can be primarily attributed to the strong mechanical constraint facilitated by the gradient composite design, which effectively reduced the stresses on the less creep-resistant alloy in the transition zone. Beyond examining plate joints, it is crucial to assess the deformation response of tubular transition joints under pressure loading and transient temperature conditions to substantiate and demonstrate the effectiveness of the design. The simulation results affirm that the tubular transition joint demonstrates superior resistance compared to its counterpart DMW when subjected to multiaxial stresses in tubular structures. In addition, optimization of the transition joint’s geometry dimensions has been conducted to diminish the accumulated deformation and enhance the service life. Lastly, the scalability and potential of the innovative transition joints for large-diameter pipe applications are addressed.

Zhang, Wei↗

The effect of secondary phases on microstructure and irradiation damage in an as-built additively manufactured 316 L stainless steel with a hafnium compositional gradient

Additive manufacturing (AM) or rapid prototyping has become a crucial tool for reducing both cost and time while increasing efficiency in qualifying structural materials for reactor use. In this study, directed energy deposition (DED) was used to develop an as-built 316 L stainless steel sample with three regions of increasing Hf-dopant to study the effects of Hf on the irradiation response of the material. Morphological and microstructural changes were analyzed before and after 2 MeV proton irradiation at 360 °C to a damage of 2.5 dpa at ~ 10 µm below the surface. The addition of Hf effectively suppressed radiation-induced damage (dislocation loops, radiation-induced segregation) due to enhanced point defect recombination. The radiation damage seen in the as-built sample was further compared to a thermo-mechanically treated counterpart of the same fabrication and irradiation parameters which was found to behave superiorly. In conclusion, the increased radiation resistance of this material may be attributed to the as-built microstructure, which includes undissolved Hf particles, delta ferrite grains and cellular sub-grain boundaries that can hinder defect motion.

316L Stainless Steel↗

Mo-100 Disk Structural Integrity Tests I

The NorthStar pressed and sintered 29 mm diameter, 0.5 mm thick Mo-100 disks will be subjected to large temperature gradients when in beam during production. The beam preferentially heats the disk center, while the aggressive cooling by the helium keeps areas outside the beam spot relatively cool. The result is a temperature gradient of up to 900°C over the radial distance of 14.5 mm. To test the disk integrity with this high temperature gradient, a special holder for a single disk was made such that proper cooling can be applied as in the design target and the disk heating can be replicated with induction heating. NorthStar press and sintered disks were used for these tests.

36 MATERIALS SCIENCE↗

Structural and electronic properties of rare-earth chromites: A computational and experimental study

Here, in this work, the structural, optical, and electronic properties of rare-earth perovskites of the general formula R⁢ CrO 3 , where R⁢ represents the rare-earth Gd, Tb, Dy, Ho, Er, and Tm, have been studied in detail. These compounds were synthesized through a facile citrate route. X-ray diffraction, Raman spectroscopy, and UV-Visible spectroscopy were utilized to reveal the structural evolutions in R⁢ ⁢CrO 3 . The lattice parameters, Cr 3+ –O 2 –Cr 3+ bond angle, and CrO 6 octahedral distortions were found to strongly depend on the ionic radii of R⁢ . First-principles calculations based on density-functional theory within the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) and strongly constrained and appropriately normed (SCAN) meta-GGA were also employed to calculate the structural and electronic properties of R⁢ ⁢CrO 3 . The ground-state energy, lattice constants, electronic structures, and density of states of R⁢ ⁢CrO 3 were calculated. These provide some insights into the electronic characteristics of the R⁢ ⁢CrO 3 compounds. The calculated values of lattice parameters and band gaps with Hubbard U correction (SCAN+ U ) agree well with values measured experimentally and show more accuracy in predicting the ground-state crystal structure and band structure compared to PBE+ U approximation. The band gap of R⁢ ⁢CrO 3 is found to be independent of the ionic radii of R⁢ from both experiments and calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗

Micro-structural features and material properties impact on adhesive metal joints via computational modeling and machine learning

The quality of structural bonding in practical applications depends on various factors arising from materials, pre-processing conditions, and manufacturing. Understanding how these factors influence bonding performance and determining their relative importance are of significant interest. Thus, this study evaluates the effects of microstructural features and material properties on the structural strength of adhesively-bonded metal joints at the submillimeter scale, utilizing a combination of Finite Element Modeling (FEM) and Machine Learning (ML) with Gradient Boosting Regression (GBR). The microstructural features include adhesive thickness, internal voids within the adhesive, adherend-adhesive interfacial voids, void size and volume fraction, and surface roughness. The material properties include the constitutive behavior of the adhesive, as well as the adherend-adhesive interfacial strength and fracture energy. The changes in structural strength and morphologies of the bonded metal structures with respect to different microstructural features and material properties were clarified by FEM. By further leveraging ML-GBR, the sequence of importance of these factors affecting bonding performance across various scenarios was summarized. This work provides valuable insights into the development of improved structural bonding for adhesive joints in industries such as automotive , aerospace, and beyond.

36 MATERIALS SCIENCE↗

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Imaging X-ray observations and the temperature structure of the Cygnus Loop

Imaging X-ray observations of the northeast portion of the Cygnus Loop (NGC 6992) reveal that X-ray emission extends 5 arc min beyond the optical filaments. Comparison with Fe XIV 5303 forbidden line measurements across the emitting shell is consistent with a temperature gradient across the shell from 10 to the 6.1 K (outside) to 10 to the 6.3 K (inside). Both X-ray and Fe XIV forbidden line intensity and spatial structure can be explained by a Sedov adiabatic blast wave with a shock velocity of 290 plus or minus 15 km/s.

Tuohy, I. R.↗