Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic gradients”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sampling in Long-Screened Wells: Issues, Misconceptions, and Solutions

The issues associated with long-screened wells (LSWs) (and open boreholes) at contaminated sites are well documented in the groundwater literature but are still not fully appreciated in practice. As established in seminal and review papers going back over three decades, the interpretation of sampling results from LSWs is challenging in the presence of vertical hydraulic gradients and borehole flow; furthermore, LSWs allow for vertical redistribution of contamination between aquifer layers. Acknowledgment of these issues has led to the development of new technologies and well designs to enable discrete-zone monitoring (DZM), yet LSWs remain common for many reasons, for example, as multipurpose wells, for geophysical logging, and (or) as legacy installations. Despite the literature on LSWs and despite the adoption of DZM at many sites, the use of LSWs persists and the challenges of interpreting sampling results from LSWs remain. In this issue paper, we provide a conceptual overview of the problems posed by LSWs and review existing literature and past work to improve the interpretation of sampling in LSWs. We draw on experience from previous studies at the Hanford Site in eastern WA, USA, and use synthetic examples to illustrate key concepts and challenges for interpretation. A recently published analytical modeling framework is used to develop illustrative synthetic examples and demonstrate a workflow for building scientific intuition to understand issues around interpreting samples from LSWs, which is critical to effective characterization and groundwater remediation at sites with LSWs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

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

Carbon and Sulfur Isotopic Signatures of Ancient Life and Environment at the Microbial Scale: Neoarchean Shales and Carbonates

An approach to coordinated, spatially resolved, in situ carbon isotope analysis of organic matter and carbonate minerals, and sulfur three- and four-isotope analysis of pyrite with an unprecedented combination of spatial resolution, precision, and accuracy is described. Organic matter and pyrite from eleven rock samples of Neoarchean drill core express nearly the entire range of delta(sup 13)C, delta(sup 34)S, Delta(sup 33)S, and Delta(sup 36)S known from the geologic record, commonly in correlation with morphology, mineralogy, and elemental composition. A new analytical approach (including a set of organic calibration standards) to account for a strong correlation between H/C and instrumental bias in SIMS delta(sup 13)C measurement of organic matter is identified. Small (2-3 microns) organic domains in carbonate matrices are analyzed with sub-permil accuracy and precision. Separate 20- to 50-micron domains of kerogen in a single approx. 0.5 cu cm sample of the approx. 2.7 Ga Tumbiana Formation have delta(sup 13)C = −52.3 +/- 0.1per mille and −34.4 +/- 0.1per mille, likely preserving distinct signatures of methanotrophy and photoautotrophy. Pyrobitumen in the approx. 2.6 Ga Jeerinah Formation and the approx. 2.5 Ga Mount McRae Shale is systematically 13C-enriched relative to co-occurring kerogen, and associations with uraniferous mineral grains suggest radiolytic alteration. A large range in sulfur isotopic compositions (including higher Delta(sup 33)S and more extreme spatial gradients in Delta(sup 33)S and Delta(sup 36)S than any previously reported) are observed in correlation with morphology and associated mineralogy. Changing systematics of delta(sup 34)S, Delta(sup 33)S, and Delta(sup 36)S, previously investigated at the millimeter to centimeter scale using bulk analysis, are shown to occur at the micrometer scale of individual pyrite grains. These results support the emerging view that the dampened signature of mass-independent sulfur isotope fractionation (S-MIF) associated with the Mesoarchean continued into the early Neoarchean, and that the connections between methane and sulfur metabolism affected the production and preservation of S-MIF during the first half of the planet's history.

isotopic signatures↗
DOE OSTI2022

RASER MRI: Magnetic resonance images formed spontaneously exploiting cooperative nonlinear interaction

The spatial resolution of magnetic resonance imaging (MRI) is limited by the width of Lorentzian point spread functions associated with the transverse relaxation rate 1/T 2 * . Here, we show a different contrast mechanism in MRI by establishing RASER (radio-frequency amplification by stimulated emission of radiation) in imaged media. RASER imaging bursts emerge out of noise and without applying radio-frequency pulses when placing spins with sufficient population inversion in a weak magnetic field gradient. Small local differences in initial population inversion density can create stronger image contrast than conventional MRI. This different contrast mechanism is based on the cooperative nonlinear interaction between all slices. On the other hand, the cooperative nonlinear interaction gives rise to imaging artifacts, such as amplitude distortions and side lobes outside of the imaging domain. Contrast mechanism and artifacts are explored experimentally and predicted by simulations on the basis of a proposed RASER MRI theory.

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

Analysis and testing of two-dimensional vented Coanda ejectors with asymmetric variable area mixing sections

The analysis of asymmetric, curved (Coanda) ejector flow has been completed using a finite difference technique and a quasi-orthogonal streamline coordinate system. The boundary layer type jet mixing analysis accounts for the effect of streamline curvature in pressure gradients normal to the streamlines and on eddy viscosities. The analysis assured perfect gases, free of pressure discontinuities and flow separation and treated three compound flows of supersonic and subsonic streams. Flow parameters and ejector performance were measured in a vented Coanda flow geometry for the verification of the computer analysis. A primary converging nozzle with a discharge geometry of 0.003175 m x 0.2032 m was supplied with 0.283 cu m/sec of air at about 241.3 KPa absolute stagnation pressure and 82 C stagnation temperature. One mixing section geometry was used with a 0.127 m constant radius Coanda surface. Eight tests were run at spacing between the Coanda surface and primary nozzle 0.01915 m and 0.318 m and at three angles of Coanda turning: 22.5 deg, 45.0 deg, and 75.0 deg. The wall static pressures, the loci of maximum stagnation pressures, and the stagnation pressure profiles agree well between analytical and experimental results.

Maroti, L. A.↗
DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

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

Transient development of perturbations in stratified shear flow

Transient development of perturbations in inviscid stratified shear flow is investigated. Use is made of closed form analytic solutions that allow concise identification of optimally growing plane-wave solutions for the case of an unbounded flow with constant shear and stratification. For the case of channel flow, variational techniques are employed to determine the optimally growing disturbances. The maximum energy growth attained over a specific time interval decreases continuously with increasing stratification, and no special significance attaches to Ri = 0.25. Indeed, transient growth can be substantial even for Ri = O(1). A general lower bound on the energy growth attained by an optimal perturbation in a stratified flow over a given time interval is the square root of the growth attained by the corresponding perturbation in unstratified flow. Enhanced perturbation persistence is found for mean-flow stratification lying in the range Ri between 0.1 and 0.3. Small but finite perturbations in mean flow with Ri less than 0.4 produce regions with locally negative total density gradient, which are expected to overturn. Although the perturbations are of wave form, buoyancy fluxes mediate transfer between perturbation kinetic and potential energy during transient development, thus implying that buoyancy flux is not a determinative diagnostic for distinguishing between waves and turbulence in stratified flows.

Farrell, Brian F.↗
NASA NTRS1997

The Effect of Temperature Dependent Material Nonlinearities on the Response of Piezoelectric Composite Plates

Previously developed analytical formulations for piezoelectric composite plates are extended to account for the nonlinear effects of temperature on material properties. The temperature dependence of the composite and piezoelectric properties are represented at the material level through the thermopiezoelectric constitutive equations. In addition to capturing thermal effects from temperature dependent material properties, this formulation also accounts for thermal effects arising from: (1) coefficient of thermal expansion mismatch between the various composite and piezoelectric plies and (2) pyroelectric effects on the piezoelectric material. The constitutive equations are incorporated into a layerwise laminate theory to provide a unified representation of the coupled mechanical, electrical, and thermal behavior of smart structures. Corresponding finite element equations are derived and implemented for a bilinear plate element with the inherent capability to model both the active and sensory response of piezoelectric composite laminates. Numerical studies are conducted on a simply supported composite plate with attached piezoceramic patches under thermal gradients to investigate the nonlinear effects of material property temperature dependence on the displacements, sensory voltages, active voltages required to minimize thermal deflections, and the resultant stress states.

Lee, Ho-Jun↗
DOE OSTI2022

Predicting solid state material platforms for quantum technologies

Semiconductor materials provide a compelling platform for quantum technologies (QT). However, identifying promising material hosts among the plethora of candidates is a major challenge. Therefore, we have developed a framework for the automated discovery of semiconductor platforms for QT using material informatics and machine learning methods. Different approaches were implemented to label data for training the supervised machine learning (ML) algorithms logistic regression, decision trees, random forests and gradient boosting. We find that an empirical approach relying exclusively on findings from the literature yields a clear separation between predicted suitable and unsuitable candidates. In contrast to expectations from the literature focusing on band gap and ionic character as important properties for QT compatibility, the ML methods highlight features related to symmetry and crystal structure, including bond length, orientation and radial distribution, as influential when predicting a material as suitable for QT.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonlinear Shear Rheology of Entangled Polymer Rings

Steady-state shear viscosity η($\bar{γ}$) of unconcatenated ring polymer melts as a function of the shear rate $\bar{γ}$ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity η ~ $\bar{γ}$ –0.56 ± 0.02 , independent of chain length, for Weissenberg numbers up to about 10 2 . Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with η ~ $\bar{γ}$ -0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high $\bar{γ}$. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ~ 10 4 , the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime η ~ $\bar{γ}$ –0.57 , we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large $\bar{γ}$. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/$\bar{γ}$. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox–Merz rule for ring polymers.

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

A New Correction Technique for Strain-Gage Measurements Acquired in Transient-Temperature Environments

Significant strain-gage errors may exist in measurements acquired in transient-temperature environments if conventional correction methods are applied. As heating or cooling rates increase, temperature gradients between the strain-gage sensor and substrate surface increase proportionally. These temperature gradients introduce strain-measurement errors that are currently neglected in both conventional strain-correction theory and practice. Therefore, the conventional correction theory has been modified to account for these errors. A new experimental method has been developed to correct strain-gage measurements acquired in environments experiencing significant temperature transients. The new correction technique has been demonstrated through a series of tests in which strain measurements were acquired for temperature-rise rates ranging from 1 to greater than 100 degrees F/sec. Strain-gage data from these tests have been corrected with both the new and conventional methods and then compared with an analysis. Results show that, for temperature-rise rates greater than 10 degrees F/sec, the strain measurements corrected with the conventional technique produced strain errors that deviated from analysis by as much as 45 percent, whereas results corrected with the new technique were in good agreement with analytical results.

Richards, W. Lance↗
DOE OSTI2022

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗
DOE OSTI2025

Introducing GPU Acceleration into the Python-Based Simulations of Chemistry Framework

We introduce the first version of GPU4P Y SCF, a module that provides GPU acceleration of methods in P Y SCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using the Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree–Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of 2 orders of magnitude with respect to the multithreaded CPU Hartree–Fock code of P Y SCF and the performance comparable to other open-source GPU-accelerated quantum chemical packages, including GAMESS and QUICK, on a single NVIDIA A100 GPU.

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

An Integrated Tool for the Coupled Thermal and Mechanical Analysis of Pyrolyzing Heatshield Materials

Materials that pyrolyze at elevated temperature have been commonly used as thermal protection materials in hypersonic flight, and advanced pyrolyzing materials for this purpose continue to be developed. Because of the large temperature gradients that can arise in thermal protection materials, significant thermal stresses can develop. Advanced applications of pyrolytic materials are calling for more complex heatshield configurations, making accurate thermal stress analysis more important, and more challenging. For non-pyrolyzing materials, many finite element codes are available and capable of performing coupled thermal-mechanical analyses. These codes do not, however, have a built-in capability to perform analyses that include pyrolysis effects. When a pyrolyzing material is heated, one or more components of the original virgin material pyrolyze and create a gas. This gas flows away from the pyrolysis zone to the surface, resulting in a reduction in surface heating. A porous residue, referred to as char, remains in place of the virgin material. While the processes involved can be complex, it has been found that a simple physical model in which virgin material reacts to form char and pyrolysis gas, will yield satisfactory analytical results. Specifically, the effects that must be modeled include: (1) Variation of thermal properties (density, specific heat, thermal conductivity) as the material composition changes; (2) Energy released or absorbed by the pyrolysis reactions; (3) Energy convected by the flow of pyrolysis gas from the interior to the surface; (4) The reduction in surface heating due to surface blowing; and (5) Chemical and mass diffusion effects at the surface between the pyrolysis gas and edge gas Computational tools for the one-dimensional thermal analysis these materials exist and have proven to be reliable design tools. The objective of the present work is to extend the analysis capabilities of pyrolyzing materials to axisymmetric configurations, and to couple thermal and mechanical analyses so that thermal stresses may be efficiently and accurately calculated.

Pronchick, Stephen W.↗
DOE OSTI2021

Solvent-induced electrochemistry at an electrically asymmetric carbon Janus particle

Chemical doping through heteroatom substitution is often used to control the Fermi level of semiconductor materials. Doping also occurs when surface adsorbed molecules modify the Fermi level of low dimensional materials such as carbon nanotubes. A gradient in dopant concentration, and hence the chemical potential, across such a material generates usable electrical current. This opens up the possibility of creating asymmetric catalytic particles capable of generating voltage from a surrounding solvent that imposes such a gradient, enabling electrochemical transformations. In this work, we report that symmetry-broken carbon particles comprised of high surface area single-walled carbon nanotube networks can effectively convert exothermic solvent adsorption into usable electrical potential, turning over electrochemical redox processes in situ with no external power supply. The results from ferrocene oxidation and the selective electro-oxidation of alcohols underscore the potential of solvent powered electrocatalytic particles to extend electrochemical transformation to various environments.

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

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗
DOE OSTI2021

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

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

Microfluidic Investigation of Salinity-Induced Oil Recovery in Porous Media during Chemical Flooding

High and low salinity water flooding are common oil recovery processes performed in the oil fields for extracting crude oil from the reservoir. These processes are often performed sequentially, naturally establishing non-uniform salinity in the porous subsurface. In this article, we investigate oil transport in porous media induced by salinity change upon flooding with high and low salinity water. As we observe a large number of impervious dead-ends from three-dimensional imaging of the actual reservoir, we identify that these areas play an important role in oil recovery where the oil transport is governed by the salinity change rather than hydrodynamics. The salinity gradients induced upon high salinity water flooding provide pathways to enhance the transport of oil drops trapped in the dead-end regions via non-equilibrium effects. However, above a critical salinity, we observe a rapid aggregation of drops that lead to the complete blockage of the pore space, thereby inhibiting oil recovery. We also find that, at an intermediate salinity where the drop aggregation is modest, the aggregation rather promotes the oil recovery. Furthermore, our observations suggest that there exist optimal salinity conditions for maximizing oil recovery during chemical flooding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗