Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ordering transformation”

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 487 records · Page 27

Gas chromatographic concepts for the analysis of planetary atmospheres

Over the last few years, new gas chromatographic (GC) concepts were developed for use on board spacecraft or any other restricted environments for determining the chemical composition of the atmosphere and surface material of various planetary bodies. Future NASA Missions include an entry probe that will be sent to Titan and various spacecraft that will land on Mars. In order to be able to properly respond to the mission science requirements and physical restrictions imposed on the instruments by these missions, GC analytical techniques are being developed. Some of these techniques include hardware and mathematical techniques that will improve GC sensitivity and increase the sampling rate of a GC descending through a planetary atmosphere. The technique of Multiplex Gas Chromatography (MGC) is an example of a technique that was studied in a simulated Titan atmosphere. In such an environment, the atmospheric pressure at instrument deployment is estimated to be a few torr. Thus, at such pressures, the small amount of sample that is acquired might not be enough to satisfy the detection requirements of the gas chromatograph. In MGC, many samples are pseudo-randomly introduced to the chromatograph without regard to elution of preceding components. The resulting data is then reduced using mathematical techniques such as cross-correlation of Fourier Transforms. Advantages realized from this technique include: improvement in detection limits of several orders of magnitude and increase in the number of analyses that can be conducted in a given period of time. Results proving the application of MGC at very low pressures emulating the same atmospheric pressures that a Titan Probe will encounter when the instruments are deployed are presented. The sample used contained hydrocarbons that are expected to be found in Titan's atmosphere. In addition, a new selective modulator was developed to monitor water under Martian atmospheric conditions. Since this modulator is selective only to water, the need for a GC column is eliminated. This results in further simplification of the instrument.

Valentin, J. R.↗

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trisodium Phosphate Phases and Solubility in Alkaline Solutions Relevant to Radioactive Waste Processing

The U.S. Department of Energy’s Hanford Site faces significant challenges in managing millions of gallons of legacy radioactive waste, where phosphate precipitation can obstruct pipelines during retrieval and processing. To resolve long-standing inconsistencies in reported solubility and clarify factors governing trisodium phosphate hydrates, we examined the Na3PO4:NaOH:H2O system. Powder and single-crystal X-ray diffraction revealed that commercial precursors undergo transformations that produce multiple hydrates, including three previously unreported phases comprising an ordered polymorph of Na3(PO4)·12H2O·1/6NaOH, Na3PO4·5H2O, and Na3PO4·9H2O. Computational modeling indicated that the ordered dodecahydrate is more stable than its disordered counterpart, suggesting kinetic persistence of structural disorder. Solubility measurements were conducted with solutions prepared either from anhydrous Na3PO4 or from Na3(PO4)·12H2O·1/6NaOH and revealed that release of interstitial NaOH from the hydrate precursor elevated solution alkalinity, thereby significantly reducing phosphate solubility relative to solutions prepared with anhydrous Na3PO4 across 20–44 °C. These results begin to reconcile inconsistencies in prior solubility data and clarify how phase composition dictates phosphate precipitation under alkaline conditions.

Graham, Trenton R. (ORCID:0000000189078004)↗

Thermodynamics of high-pressure ice phases explored with atomistic simulations

Most experimentally known high-pressure ice phases have a body-centred cubic (bcc) oxygen lattice. Our large-scale molecular-dynamics simulations with a machine-learning potential indicate that, amongst these bcc ice phases, ices VII, VII' and X are the same thermodynamic phase under different conditions, whereas superionic ice VII" has a first-order phase boundary with ice VII'. Moreover, at about 300 GPa, the transformation between ice X and the Pbcm phase has a sharp structural change but no apparent activation barrier, whilst at higher pressures the barrier gradually increases. Our study thus clarifies the phase behaviour of the high-pressure ices and reveals peculiar solid–solid transition mechanisms not known in other systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Tuning of electronic properties in highly lattice-mismatched epitaxial SmN

Here, we establish the relationship between native N vacancies, introduced through varying growth parameters, and the structural and transport properties in SmN thin films grown via molecular beam epitaxy grown on MgO(001). The varying levels of N vacancies introduced by varying the ratio of Sm to N atoms during deposition creates excess carriers that radically transform the electrical behavior of the film, over a range of five orders of magnitude, from highly resistive to highly conductive, and unlocking a phase transition evidenced by the presence of a ferromagnetic feature in resistivity. X-ray photoelectron spectroscopy results show that this effect is much less pronounced when varying the available nitrogen species. These samples retain a highly crystal quality despite being grown on a substrate with a lattice mismatch of 20%, alleviating the strain by forming a highly strained SmO oxide layer. The integration between SmN and several transition metal nitride compounds has the potential to unlock new architectures for Josephson junction devices.

36 - MATERIALS SCIENCE↗

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. Here, these methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

16 O Electroweak Response Functions from First Principles

We present calculations of various electroweak response functions for the 16 O nucleus obtained using coupled-cluster theory in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-leading order and next-to-next-to-leading order in chiral effective field theory and perform a Bayesian analysis to assess uncertainties. Our results are in good agreement with available electron-scattering data at |𝐪|≈326 MeV/c. Additionally, we provide several predictions for the weak response functions in the quasielastic peak region at |𝐪| =300 and 400 MeV/c, which are critical for long-baseline neutrino experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Metallization of colloidal crystals

Colloidal crystals formed by size-asymmetric binary particles coassemble into a wide variety of colloidal compounds with lattices akin to ionic crystals. Recently, a transition from a compound phase with a sublattice of small particles to a metal-like phase in which the small particles are delocalized has been predicted computationally and observed experimentally. In this colloidal metallic phase, the small particles roam the crystal maintaining the integrity of the lattice of large particles, as electrons do in metals. A similar transition also occurs in superionic crystals, termed sublattice melting. Here, we use energetic principles and a generalized molecular dynamics model of a binary system of functionalized nanoparticles to analyze the transition to sublattice delocalization in different coassembled crystal phases as a function of temperature (T), number of grafted chains on the small particles, and number ratio between the small and large particles n s : n l . Here we find that n s : n l is the primary determinant of crystal type due to energetic interactions and interstitial site filling, while the number of grafted chains per small particle determines the stability of these crystals. We observe first-order sublattice delocalization transitions as T increases, in which the host lattice transforms from low- to high-symmetry crystal structures, including A20 → bct → bcc, A d → bct → bcc, and bcc → bcc/fcc → fcc transitions and lattices. Analogous sublattice transitions driven primarily by lattice vibrations have been seen in some atomic materials exhibiting an insulator-metal transition also referred to as metallization. We also find minima in the lattice vibrations and diffusion coefficient of small particles as a function of n s : n l , indicating enhanced stability of certain crystal structures for n s : n l values that form compounds.

36 MATERIALS SCIENCE↗

On the Impact of High-Order Harmonic Generation in Electrical Distribution Systems

The modern power grid has seen a rise in the integration of non-linear loads, presenting a significant concern for operators. These loads introduce unwanted harmonics, leading to potential issues such as overheating and improper functioning of circuit breakers. In pursuing a more sustainable grid, the adoption of electric vehicles (EVs) and photovoltaic (PV) systems in residential networks has increased. Understanding and examining the effects of high-order harmonic frequencies beyond $1.5$ kHz is crucial to understanding their impact on the operation and planning of electrical distribution systems under varying nonlinear loading conditions. This study investigates a diverse set of critical power electronic loads within a household modeled using PSCAD/EMTdc, analyzing their unique harmonic spectra. This information is utilized to run the time-series harmonic analysis program in OpenDSS on a modified IEEE 34 bus test system model. The impact of high-order harmonics is quantified using metrics that evaluate total harmonic distortion (THD), transformer harmonic-driven eddy current loss component, and propagation of harmonics from the source to the substation transformer.

Peerzada, Aaqib A. [BATTELLE (PACIFIC NW LAB)]↗

Distribution Market-Clearing and Pricing Considering Coordination of DSOs and ISO: An EPEC Approach

Distribution-level electricity market provides a platform for trading energy and grid services from large-amount of small-scale distributed energy resources (DERs) located on distribution grids. The behavior of the DERs in the distribution electricity market may ultimately impact the market-clearing and locational marginal prices (LMPs) in the wholesale market. This paper proposes a bi-level optimization model for distribution market clearing and distribution locational marginal pricing (DLMP) considering the interactions between distribution and transmission wholesale markets. In the proposed model, the upper-level model represents the distribution system operator (DSO) market-clearing and the lower-level model represents the wholesale market-clearing by the independent system operator (ISO). The LMP at the substation will impact the DER dispatch and power demands of the DSO as well as the DLMP. In turn, the power demands of the DSO will further impact the ISO market and its LMPs. The equilibrium problem with equilibrium constraints (EPEC) approach is applied to find the equilibria of multiple DSOs and the ISO. The EPEC problem is transformed into a single-level mixed-integer convex problem in order to allow for efficient solving. Finally, the effectiveness of the proposed model and solution method are demonstrated through case studies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Powder Diffraction Crystallography : 98 Years as Plan B?

A historical essay written in 2014 to accompany a talk in the Transactions Symposium "100 Years of Crystallography" at the 2014 American Crystallographic Association Annual Meeting (talk at https://history.amercrystalassn.org/h-brian-toby-videos). Abstract: Powder diffraction patterns were used for structure solution almost contemporaneously with the Bragg’s initial single-crystal determinations. Structures have continued to be determined from powders, in the subsequent century. There has been a nearly explosive growth in this technique since the 1970’s. This was initially driven by neutron diffraction, which prompted the invention of Rietveld analysis, but synchrotrons have gone on to have at least as large an impact. There is no apparent upper limit to the scale of systems that can be structurally studied with powder diffraction; novel protein structures have been determined. When short-range ordering or disordered regions of a material are of interest, the Fourier transform of a powder diffraction pattern, the pair distribution function can be modeled.

Toby, Brian H.↗

Uncertainties in ab initio nuclear structure calculations with chiral interactions

We present theoretical ground state energies and their uncertainties for p -shell nuclei obtained from chiral effective field theory internucleon interactions as a function of chiral order, fitted to two- and three-body data only. We apply a Similary Renormalization Group transformation to improve the numerical convergence of the many-body calculations, and discuss both the numerical uncertainties arising from basis truncations and those from omitted induced many-body forces, as well as chiral truncation uncertainties. With complete Next-to-Next-to-Leading (N 2 LO) order two- and three-body interactions, we find significant overbinding for the ground states in the upper p -shell, but using higher-order two-body potentials, in combination with N 2 LO three-body forces, our predictions agree with experiment throughout the p -shell to within our combined estimated uncertainties. The uncertainties due to chiral order truncation are noticeably larger than the numerical uncertainties, but they are expected to become comparable to the numerical uncertainties at complete N 3 LO.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of the use of a singularity element in finite element analysis of center-cracked plates

Two different methods are applied to the analyses of finite width linear elastic plates with central cracks. Both methods give displacements as a primary part of the solution. One method makes use of Fourier transforms. The second method employs a coarse mesh of triangular second-order finite elements in conjunction with a single singularity element subjected to appropriate additional constraints. The displacements obtained by these two methods are in very good agreement. The results suggest considerable potential for the use of a cracked element for related crack problems, particularly in connection with the extension to nonlinear material behavior.

Mendelson, A.↗

Frequency domain stability criteria for distributed nonlinear sampled-data systems.

Popov-type sufficient conditions are presented for the absolute stability of a class of closed loop sampled-data systems. The closed loop system consists of a linear distributed element, a feedback controller (linear or nonlinear), a sampler, and possibly a zero-order hold circuit. The distributed element is assumed to be finite Hankel transformable, and such that its dynamics can be represented by a transfer function which is the ratio of the multiple (Laplace and finite Hankel) transforms of output and input. The input is assumed to be distributed. The stability criterion presented here parallels the criterion for distributed systems that are finite Fourier transformable. An example is given to illustrate the applications of the stability criterion.

Nagarajan, N.↗

Signature extension through the application of cluster matching algorithms to determine appropriate signature transformations

Signature extension is intended to increase the space-time range over which a set of training statistics can be used to classify data without significant loss of recognition accuracy. A first cluster matching algorithm MASC (Multiplicative and Additive Signature Correction) was developed at the Environmental Research Institute of Michigan to test the concept of using associations between training and recognition area cluster statistics to define an average signature transformation. A more recent signature extension module CROP-A (Cluster Regression Ordered on Principal Axis) has shown evidence of making significant associations between training and recognition area cluster statistics, with the clusters to be matched being selected automatically by the algorithm.

Lambeck, P. F.↗

Numerical solution of three-dimensional free turbulent shear flows

Numerical results are presented for three-dimensional compressible turbulent jet and wake flows. An alternating direction implicit numerical procedure is used to solve the finite-difference form of the parabolic elliptic Navier-Stokes equations. A coordinate transformation maps the boundaries at infinity into a finite computational domain in order to properly specify infinity boundary conditions as well as contain the downstream growth of the viscous flow field in a fixed computational grid. Turbulence closure is achieved through an algebraic mixing length eddy viscosity model. Numerical results for supersonic flow are presented for an axisymmetric jet, an elliptical jet, an elliptical wake, and two interacting rectangular jets. Experimental data were not available for comparison with the numerical results. However, the results compare well with empirical results for free shear flows.

Oh, Y. H.↗

Development of a three-dimensional Navier-Stokes code on CDC star-100 computer

A three-dimensional code in body-fitted coordinates was developed using MacCormack's algorithm. The code is structured to be compatible with any general configuration, provided that the metric coefficients for the transformation are available. The governing equations are developed in primitive variables in order to facilitate the incorporation of physical boundary conditions and turbulence-closure models. MacCormack's two-step, unsplit, time-marching algorithm is used to solve the unsteady Navier-Stokes equations until steady-state solution is achieved. Cases discussed include (1) flat plate in supersonic free stream; (2) supersonic flow along an axial corner; (3) subsonic flow in an axial corner at M infinity = 0.95; and (4) supersonic flow in an axial corner at M infinity 1.5.

Vatsa, V. N.↗