SEARCH · Engineering Papers

Results for “inversion theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗

DOE OSTI2025

Evaluating design safety margins in the American Society of Mechanical Engineers graphite core components design-by-analysis assessments

Graphite is an important material being used for core components in next-generation high-temperature gas-cooled nuclear reactors. The selection of graphite grade for a specific Designer is a complex task, dependent on reactor conditions, component functionality, and required reliability. The American Society of Mechanical Engineers (ASME) provides two semi-probabilistic design-by-analysis assessments to evaluate graphite core components against design reliability targets. The simplified assessment uses a 2-parameter Weibull distribution to describe the graphite grade’s tensile-strength distribution to establish component stress limits. The full assessment uses the 3-parameter Weibull distribution and a modified Weakest-Link Theory approach to calculate a component design probability of failure. The paper defines recommended assessment rules, which are the as-written simplified assessment and the full assessment with parameter lower bounds, the modulus update with threshold reduction, and the 2027 grouping rules. Code rules are applied to three grades: 2114, IG-110, and NBG-18. The baseline margin calculation is developed using the experimental tensile dogbone specimen. Percent margin is defined as the percent reduction in the median experimental load to obtain the allowable load per ASME assessments. Under the recommended rules, the SRC–1 margin in the simplified assessment ranged from 40.2 % to 52.7 % among the grades in this study and from 36.1 % to 49.8 % in the full assessment. The full assessment only decreases the margins by 2.5–4.5 % for the SRC-1 components and 0–1.5 % for the SRC-2 components for this baseline case. Margin is inversely related to material median strength (i.e., the strongest grade, 2114, has the lowest margin).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Securely Aggregated Coded Matrix Inversion

Coded computing is a method for mitigating straggling workers in a centralized computing network, by using erasure-coding techniques. Federated learning is a decentralized model for training data distributed across client devices. In this work we propose approximating the inverse of an aggregated data matrix, where the data is generated by clients; similar to the federated learning paradigm, while also being resilient to stragglers. To do so, we propose a coded computing method based on gradient coding. We modify this method so that the coordinator does not access the local data at any point; while the clients access the aggregated matrix in order to complete their tasks. Here, the network we consider is not centrally administrated, and the communications which take place are secure against potential eavesdroppers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Shell-structure and asymmetry effects in level densities

Level density [Formula: see text] is derived for a nuclear system with a given energy [Formula: see text], neutron [Formula: see text], and proton [Formula: see text] particle numbers, within the semiclassical extended Thomas–Fermi and periodic-orbit theory beyond the Fermi-gas saddle-point method. We obtain [Formula: see text], where [Formula: see text] is the modified Bessel function of the entropy [Formula: see text], and [Formula: see text] is related to the number of integrals of motion, except for the energy [Formula: see text]. For small shell structure contribution one obtains within the micro–macroscopic approximation (MMA) the value of [Formula: see text] for [Formula: see text]. In the opposite case of much larger shell structure contributions one finds a larger value of [Formula: see text]. The MMA level density [Formula: see text] reaches the well-known Fermi gas asymptote for large excitation energies, and the finite micro-canonical limit for low excitation energies. Fitting the MMA [Formula: see text] to experimental data on a long isotope chain for low excitation energies, due mainly to the shell effects, one obtains results for the inverse level density parameter [Formula: see text], which differs significantly from that of neutron resonances.

Physics↗

DOE OSTI2020

On negative induced polarization in frequency domain measurements

SUMMARY Induced polarization (IP) has been widely used to non-invasively characterize electrical conduction and polarization in the subsurface resulting from an applied electric field. Earth materials exhibit a lossy capacitance defined by a negative intrinsic phase in frequency-domain IP (FDIP) or a positive intrinsic chargeability in time-domain IP (TDIP). However, error-free positive apparent phase or negative apparent chargeability (i.e. negative IP effects) can occur in IP measurements over heterogeneous media. While negative IP effects in TDIP data sets have been discussed, no studies have addressed this topic in detail for FDIP measurements. We describe theory and numerical modelling to explain the origin of negative IP effects in FDIP measurements. A positive apparent phase may occur when a relatively high polarizability feature falls into negative sensitivity zones of complex resistivity measurements. The polarity of the apparent phase is determined by the distribution of subsurface intrinsic phase and resistivity, with the resistivity impacting the apparent phase polarity via its control on the sensitivity distribution. A physical explanation for the occurrence of positive apparent phase data is provided by an electric circuit model representing a four-electrode measurement. We also show that the apparent phase polarity will be frequency dependent when resistivity changes significantly with frequency (i.e. in the presence of significant IP effects). Consequently, negative IP effects manifest themselves in the shape of apparent phase spectra recorded with multifrequency (spectral IP) data sets. Our results imply that positive apparent phase measurements should be anticipated and should be retained during inversion and interpretation of single frequency and spectral IP data sets.

Geochemistry & Geophysics↗

DOE OSTI2024

Time-reversal invariant topological moiré flat band: A platform for the fractional quantum spin Hall effect

Motivated by recent observation of the quantum spin Hall effect in monolayer germanene and twisted bilayer transition-metal-dichalcogenides (TMDs), we study the topological phases of moir twisted bilayers with time-reversal symmetry and spin sz conservation. By using a continuum model description which can be applied to both germanene and TMD bilayers, we show that at small twist angles the emergent moir flat bands can be topologically nontrivial due to inversion symmetry breaking. Each of these flat bands admits a lowest-Landau-level description for each spin projection in the chiral limit and at magic twist angle. Furthermore, this allows for the construction of a many-body Laughlin state with time-reversal symmetry which can be stabilized by a short-range pseudopotential, and therefore serves as an ideal platform for realizing the so-far elusive fractional quantum spin Hall effect with emergent spin-1/2 U(1) symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2021

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗

NASA NTRS1991

Transition operators in acoustic-wave diffraction theory. I - General theory. II - Short-wavelength behavior, dominant singularities of Zk0 and Zk0 exp -1

A formal theory of the scattering of time-harmonic acoustic scalar waves from impenetrable, immobile obstacles is established. The time-independent formal scattering theory of nonrelativistic quantum mechanics, in particular the theory of the complete Green's function and the transition (T) operator, provides the model. The quantum-mechanical approach is modified to allow the treatment of acoustic-wave scattering with imposed boundary conditions of impedance type on the surface (delta-Omega) of an impenetrable obstacle. With k0 as the free-space wavenumber of the signal, a simplified expression is obtained for the k0-dependent T operator for a general case of homogeneous impedance boundary conditions for the acoustic wave on delta-Omega. All the nonelementary operators entering the expression for the T operator are formally simple rational algebraic functions of a certain invertible linear radiation impedance operator which maps any sufficiently well-behaved complex-valued function on delta-Omega into another such function on delta-Omega. In the subsequent study, the short-wavelength and the long-wavelength behavior of the radiation impedance operator and its inverse (the 'radiation admittance' operator) as two-point kernels on a smooth delta-Omega are studied for pairs of points that are close together.

Hahne, G. E.↗

DOE OSTI2021

The diversity of three-dimensional photonic crystals

Many butterflies, birds, beetles, and chameleons owe their spectacular colors to the microscopic patterns within their wings, feathers, or skin. When these patterns, or photonic crystals, result in the omnidirectional reflection of commensurate wavelengths of light, it is due to a complete photonic band gap (PBG). The number of natural crystal structures known to have a PBG is relatively small, and those within the even smaller subset of notoriety, including diamond and inverse opal, have proven difficult to synthesize. Here, we report more than 150,000 photonic band calculations for thousands of natural crystal templates from which we predict 351 photonic crystal templates – including nearly 300 previously-unreported structures – that can potentially be realized for a multitude of applications and length scales, including several in the visible range via colloidal self-assembly. With this large variety of 3D photonic crystals, we also revisit and discuss oft-used primary design heuristics for PBG materials.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Turbulence link to L-mode, I-mode, and H-mode confinement in the DIII-D tokamak

Understanding the physics of low-confinement (L-), improved-confinement (I-), and high-confinement (H-) modes is critical for fusion reactors. The finding herein reports observations of two types of turbulence coexisting near the L-mode edge, one magnetohydrodynamic (MHD)-like and another micro-tearing mode (MTM)-like, linked to the H-mode and I-mode confinement in the DIII-D tokamak. Ion-scale magnetic and density turbulence is measured using a Faraday-effect radial-interferometer-polarimeter and beam-emission-spectroscopy (BES). Broadband turbulence spectra of up to ~600 kHz are observed in two discharges where transitions between L-mode, I-mode, and H-mode occurs. Turbulence is found to be inversely correlated with confinement, meaning lower turbulence power at higher confinement. Distinctively, the high-frequency (HF, >~100 kHz) magnetic turbulence power changes by the most (55%) during transitions primarily involving energy confinement change, whereas the low-frequency (LF, <~100 kHz) magnetic and density turbulence power changes by the most (80%) during transitions primarily involving particle confinement change. The LF turbulence amplitude oscillates with and leads to deuterium-alpha emission oscillations before an H-mode. These results imply that HF turbulence mainly affects energy confinement whereas LF turbulence can affect particle confinement. The magnetic and density turbulence exhibits coherence up to 0.6 and cross-phase magnitude close to π/2 in most cases, suggesting they have a common origin in both the LF and HF ranges. BES suggests that LF turbulence resides at the edge (ρ = 0.95) and HF turbulence can be at the outer core (ρ = 0.8) or edge (ρ = 0.95). Comparisons of measurements, theory, and gyrokinetic simulations suggest that HF turbulence is MTM-like in all cases, whereas LF turbulence is more consistent with MHD-like modes and the exact instability might change during transitions—except that a drift-wave origin is possible in a low collisionality H-mode. These results suggest that the H-mode involves suppressed MHD-like turbulence, whereas the I-mode mitigates MTM-like turbulence along with largely unchanged MHD-like turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Site Mixing for Engineering Magnetic Topological Insulators

The van der Waals compound, MnBi 2 Te 4 , is the first intrinsic magnetic topological insulator, providing a materials platform for exploring exotic quantum phenomena such as the axion insulator state and the quantum anomalous Hall effect. However, intrinsic structural imperfections lead to bulk conductivity, and the roles of magnetic defects are still unknown. With higher concentrations of the same types of magnetic defects, the isostructural compound MnSb 2 Te 4 is a better model system for a systematic investigation of the connections among magnetism, topology, and lattice defects. In this work, the impact of antisite defects on the magnetism and electronic structure is studied in MnSb 2 Te 4 . Mn-Sb site mixing leads to complex magnetic structures and tunes the interlayer magnetic coupling between antiferromagnetic and ferromagnetic. The detailed nonstoichiometry and site mixing of MnSb2Te4 crystals depend on the growth parameters, which can lead to ≈40% of Mn sites occupied by Sb and ≈15% of Sb sites by Mn in as-grown crystals. Single-crystal neutron diffraction and electron microscopy studies show nearly random distribution of the antisite defects. Band structure calculations suggest that the Mn-Sb site mixing favors a ferromagnetic interlayer coupling, consistent with experimental observation, but is detrimental to the band inversion required for a nontrivial topology. Overall, our results suggest a long-range magnetic order of Mn ions sitting on Bi sites in MnBi 2 Te 4 . The effects of site mixing should be considered in all layered heterostructures that consist of alternating magnetic and topological layers, including the entire family of MnTe(Bi 2 Te 3 )n, its Sb analogs, and their solid solution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS2004

High-Energy Emission From Millisecond Pulsars

The X-ray and gamma-ray spectrum of rotation-powered millisecond pulsars is investigated in a model for acceleration and pair cascades on open field lines above the polar caps. Although these pulsars have low surface magnetic fields, their short periods allow them to have large magnetospheric potential drops, but the majority do not produce sufficient pairs to completely screen the accelerating electric field. In these sources, the primary and secondary electrons continue to accelerate to high altitude and their Lorentz factors are limited by curvature and synchrotron radiation reaction. The accelerating particles maintain high Lorentz factors and undergo cyclotron resonant absorption of radio emission, that produces and maintains a large pitch angle, resulting in a strong synchrotron component. The resulting spectra consist of several distinct components: curvature radiation from primary electrons dominating from 1 - 100 GeV, synchrotron radiation from primary and secondary electrons dominating up to about 100 MeV, and much weaker inverse-Compton radiation from primary electrons a t 0.1 - 1 TeV. We find that the relative size of these components depends on pulsar period, period derivative, and neutron star mass and radius with the level of the synchrotron component also depending sensitively on the radio emission properties. This model is successful in describing the observed X-ray and gamma-ray spectrum of PSR J0218+4232 as synchrotron radiation, peaking around 100 MeV and extending up to a turnover around several GeV. The predicted curvature radiation components from a number of millisecond pulsars, as well as the collective emission from the millisecond pulsars in globular clusters, should be detectable with AGILE and GLAST. We also discuss a hidden population of X-ray-quiet and radio-quiet millisecond pulsars which have evolved below the pair death line, some of which may be detectable by telescopes sensitive above 1 GeV. Subject headings: pulsars: general - radiation mechanisms: nonthermal - stars: neutron - gamma rays: theory

Harding, Alice K.↗

DOE OSTI2022

Spectral density reconstruction with Chebyshev polynomials

Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semiexclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. Here, in the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to obtain the local density of states for graphene in a magnetic field as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Artificial intelligence-driven approaches for materials design and discovery

Materials design is an important component of modern science and technology, yet traditional approaches rely heavily on trial and error and can be inefficient. Computational techniques, enhanced by modern artificial intelligence, have reshaped the landscape of designing new materials. Among these approaches, inverse design has shown great promise in designing materials that meet specific property requirements. Here, in this Review, we present key computational advances in materials design over the past few decades. We follow the evolution of relevant materials design techniques, from high-throughput forward machine learning methods and evolutionary algorithms, to advanced artificial intelligence strategies such as reinforcement learning and deep generative models. We highlight the paradigm shift from conventional screening approaches to inverse generation driven by deep generative models. Finally, we discuss current challenges and future perspectives of materials inverse design. This Review may serve as a brief guide to the approaches, progress and outlook of designing future functional materials with technological relevance.

computational methods↗

DOE OSTI2022

Mapping the N = 40 island of inversion: Precision mass measurements of neutron-rich Fe isotopes

Nuclear properties across the chart of nuclides are key to improving and validating our understanding of the strong interaction in nuclear physics. Here we present high-precision mass measurements of neutron-rich Fe isotopes performed at the TITAN facility. The multiple-reflection time-of-flight mass spectrometer (MR-ToF-MS), achieving a resolving power greater than 600000 for the first time, enabled the measurement of 63–70 Fe, including first-time high-precision direct measurements (δm/m ≈ 10 -7 ) of 68–70 Fe, as well as the discovery of a long-lived isomeric state in 69 Fe. These measurements are accompanied by both mean-field and ab initio calculations using the most recent realizations which enable theoretical assignment of the spin-parities of the 69 Fe ground and isomeric states. Together with mean-field calculations of quadrupole deformation parameters for the Fe isotope chain, these results benchmark a maximum of deformation in the N = 40 island of inversion in Fe and shed light on trends in level densities indicated in the newly refined mass surface.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Grey Galaxies’ as an endpoint of the Kerr-AdS superradiant instability

Kerr-AdS d+1 black holes for d ≥ 3 suffer from classical superradiant instabilities over a range of masses above extremality. We conjecture that these instabilities settle down into Grey Galaxies (GGs) — a new class of coarse-grained solutions to Einstein’s equations which we construct in d = 3. Grey Galaxies are made up of a black hole with critical angular velocity ω = 1 in the ‘centre’ of AdS, surrounded by a large flat disk of thermal bulk gas that revolves around the centre of AdS at the speed of light. The gas carries a finite fraction of the total energy, as its parametrically low energy density and large radius are inversely related. GGs exist at masses that extend all the way down to the unitarity bound. Their thermodynamics is that of a weakly interacting mix of Kerr-AdS black holes and the bulk gas. Their boundary stress tensor is the sum of a smooth ‘black hole’ contribution and a peaked gas contribution that is delta function localized around the equator of the boundary sphere in the large N limit. We also construct another class of solutions with the same charges; ‘Revolving Black Holes (RBHs)’. RBHs are macroscopically charged SO(d, 2) descendants of AdS-Kerr solutions, and consist of ω = 1 black holes revolving around the centre of AdS at a fixed radial location but in a quantum wave function in the angular directions. RBH solutions are marginally entropically subdominant to GG solutions and do not constitute the endpoint of the superradiant instability. Nonetheless, we argue that supersymmetric versions of these solutions have interesting implications for the spectrum of supersymmetric states in, e.g. N = 4 Yang-Mills theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗