Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Predictable third harmonic generation in GaAs metasurfaces through group theory inverse design of meta-atoms

Here we report the group theory-based inverse design of meta-atoms for a dielectric metasurface in GaAs with predictable optical linear response and third harmonic generation (THG). Six sharp Fano resonances have been observed with a corresponding polarization dependence as predicted by group theory and the meta-atom’s symmetry in the D2h point group. THG has been observed for two modes under x-polarization excitation and one mode for y-polarization, in agreement with theoretical symmetry predictions. The polarization-dependent THG aspect ratio was observed to reach as high as 108. Through strategic structural or symmetry-preserving perturbations, it was shown that the THG can be enhanced or reduced by a factor of 7. The highest THG conversion efficiency was estimated to be 3.1×10-7 at the pump intensity of 1.51 MW/cm2. This high THG conversion efficiency indicates that our group theory approach to modal engineering opens a new path towards optical nonlinearity tuning in dielectric metasurfaces.

Aoueille, Andrew↗
NASA NTRS1979

FORTRAN subroutines for out-of-core solutions of large complex linear systems

The design and usage of two main subprograms using direct methods to solve large linear complex systems, of the form Ax = b, whose coeffficient matrices are too large to be stored in core are described. The first main subprogram is for systems whose coefficient matrices are of a particular sparse structure, namely, the matrix A can be written in the form B + D, where B is a block-banded system, and D has only a few columns of nonzeros. Key elements of the algorithms used in the subprograms include: the data structure, the strategy for preserving numerical stability, the adaptability of the algorithms for dense systems as well as for block-profile systems.

Yip, E. L.↗
DOE OSTI2020

Overset meshes for incompressible flows: On preserving accuracy of underlying discretizations

This study on overset meshes for incompressible-flow simulations is motivated by accurate prediction of wind farm aerodynamics involving large motions and deformations of components with complex geometry. Specifically, using first-order hyperbolic and elliptic equation proxies for the incompressible Navier-Stokes (NS) equations, we investigate the influence of information exchange between overset meshes on numerical performance where the underlying discretization is second-order accurate. The first aspect of information exchange surrounds interpolation of solution where we examine Lagrange and point-cloud-based interpolation for creating constraint equations between overset meshes. To maintain overall second-order accuracy, higher-order interpolation is required for elliptic problems, but linear interpolation is sufficient for hyperbolic problems in first-order form. Higher-order point-cloud-based interpolation provides a pathway to maintaining accuracy in unstructured meshes, but at higher complexity. The second aspect of information exchange focuses on comparing the approaches of overset single system (OSS) and overset Additive Schwarz (OAS) for coupling the linear systems of the overlapping meshes. While the former involves a single linear system, in the latter the discrete linear systems are solved separately, and solving the global system is accomplished through outer iterations and sequential information exchange in a Jacobi fashion. For the test cases studied, accuracy for hyperbolic problems is maintained by performing two outer iterations, whereas many outer iterations are required for elliptic systems. The order-of-accuracy studies explored here are critical for verifying the overset-mesh coupling algorithms used in engineering simulations. Accuracy of these simulations themselves is, however, quantified using engineering quantities of interest such as drag, power, etc. Consequently, we conclude with numerical experiments using NS equations for incompressible flows where we show that linear interpolation and few outer iterations are sufficient for achieving asymptotic convergence of engineering quantities of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

FedOSAA: Improving Federated Learning with One-Step Anderson Acceleration

Federated learning (FL) is a distributed machine learning approach that enables multiple local clients and a central server to collaboratively train a model while keeping the data on their own devices. First-order methods, particularly those incorporating variance reduction techniques, are the most widely used FL algorithms due to their simple implementation and stable performance. However, these methods tend to be slow and require a large number of communication rounds to reach the global minimizer. We propose FedOSAA, a novel approach that preserves the simplicity of first-order methods while achieving the rapid convergence typically associated with second-order methods. Our approach applies one Anderson acceleration (AA) step following classical local updates based on first-order methods with variance reduction, such as FedSVRG and SCAFFOLD, during local training. This AA step is able to leverage curvature information from the history points and gives a new update that approximates the Newton-GMRES direction, thereby significantly improving the convergence. We establish a local linear convergence rate to the global minimizer of FedOSAA for smooth and strongly convex loss functions. Numerical comparisons show that FedOSAA substantially improves the communication and computation efficiency of the original first-order methods, achieving performance comparable to second-order methods like GIANT.

Feng, Xue [University of California, Davis]↗
NASA NTRS2015

In Situ Detection of Organic Molecules on the Martian Surface With the Mars Organic Molecule Analyzer (MOMA) on Exomars 2018

The Mars Organic Molecule Analyzer (MOMA) investigation on the 2018 ExoMars rover will examine the chemical composition of samples acquired from depths of up to two meters below the martian surface, where organics may be protected from radiative and oxidative degradation. The MOMA instrument is centered around a miniaturized linear ion trap (LIT) that facilitates two modes of operation: i) pyrolysisgas chromatography mass spectrometry (pyrGC-MS); and, ii) laser desorptionionization mass spectrometry (LDI-MS) at ambient Mars pressures. The LIT also enables the structural characterization of complex molecules via complementary analytical capabilities, such as multi-frequency waveforms (i.e., SWIFT) and tandem mass spectrometry (MSMS). When combined with the complement of instruments in the rovers Pasteur Payload, MOMA has the potential to reveal the presence of a wide range of organics preserved in a variety of mineralogical environments, and to begin to understand the structural character and potential origin of those compounds.

ExoMars2018↗
NASA NTRS1992

Comparison of Predictable Smooth Ocular and Combined Eye-Head Tracking Behaviour in Patients with Lesions Affecting the Brainstem and Cerebellum

We compared the ability of eight normal subjects and 15 patients with brainstem or cerebellar disease to follow a moving visual stimulus smoothly with either the eyes alone or with combined eye-head tracking. The visual stimulus was either a laser spot (horizontal and vertical planes) or a large rotating disc (torsional plane), which moved at one sinusoidal frequency for each subject. The visually enhanced Vestibulo-Ocular Reflex (VOR) was also measured in each plane. In the horizontal and vertical planes, we found that if tracking gain (gaze velocity/target velocity) for smooth pursuit was close to 1, the gain of combined eye-hand tracking was similar. If the tracking gain during smooth pursuit was less than about 0.7, combined eye-head tracking was usually superior. Most patients, irrespective of diagnosis, showed combined eye-head tracking that was superior to smooth pursuit; only two patients showed the converse. In the torsional plane, in which optokinetic responses were weak, combined eye-head tracking was much superior, and this was the case in both subjects and patients. We found that a linear model, in which an internal ocular tracking signal cancelled the VOR, could account for our findings in most normal subjects in the horizontal and vertical planes, but not in the torsional plane. The model failed to account for tracking behaviour in most patients in any plane, and suggested that the brain may use additional mechanisms to reduce the internal gain of the VOR during combined eye-head tracking. Our results confirm that certain patients who show impairment of smooth-pursuit eye movements preserve their ability to smoothly track a moving target with combined eye-head tracking.

Grant, Michael P.↗
DOE OSTI2024

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
NASA NTRS1993

Variational data assimilation with a semi-Lagrangian semi-implicit global shallow-water equation model and its adjoint

An adjoint model is developed for variational data assimilation using the 2D semi-Lagrangian semi-implicit (SLSI) shallow-water equation global model of Bates et al. with special attention being paid to the linearization of the interpolation routines. It is demonstrated that with larger time steps the limit of the validity of the tangent linear model will be curtailed due to the interpolations, especially in regions where sharp gradients in the interpolated variables coupled with strong advective wind occur, a synoptic situation common in the high latitudes. This effect is particularly evident near the pole in the Northern Hemisphere during the winter season. Variational data assimilation experiments of 'identical twin' type with observations available only at the end of the assimilation period perform well with this adjoint model. It is confirmed that the computational efficiency of the semi-Lagrangian scheme is preserved during the minimization process, related to the variational data assimilation procedure.

Li, Y.↗
DOE OSTI2024

AMReX v2024

The software framework, AMReX, supports the development of block-structured adaptive mesh refinement (AMR) algorithms for solving systems of partial differential equations. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving the essential local descriptions of different physical processes in complex multiphysics algorithms. AMR uses a hierarchical representation of the solution at multiple levels of resolution where the solution on each level is defined on the union of data containers at that resolution. These data containers, which represent the solution over a logically rectangular subregion of the domain, can contain field data defined on a mesh, Lagrangian particles or combinations of both. In addition to these basic data types, AMReX supports a multilevel embedded boundary representation of complex geometry; linear solvers for cell-centered and nodal data; asynchronous I/O in a native format readable by ParaView, VisIt and yt; and interfaces to hypre and PETSc solvers. AMReX enables applications to run on distributed memory architectures with multicore CPUs and with GPU accelerators. AMReX uses a lightweight abstraction layer that effectively hides the details of the architecture from the application. The framework currently supports CUDA, HIP and SYCL for GPU acceleration and OpenMP for multi-core CPU architectures.

Almgren, Ann↗
NASA NTRS1995

Improved reduced-resolution satellite imagery

The resolution of satellite imagery is often traded-off to satisfy transmission time and bandwidth, memory, and display limitations. Although there are many ways to achieve the same reduction in resolution, algorithms vary in their ability to preserve the visual quality of the original imagery. These issues are investigated in the context of the Landsat browse system, which permits the user to preview a reduced resolution version of a Landsat image. Wavelets-based techniques for resolution reduction are proposed as alternatives to subsampling used in the current system. Experts judged imagery generated by the wavelets-based methods visually superior, confirming initial quantitative results. In particular, compared to subsampling, the wavelets-based techniques were much less likely to obscure roads, transmission lines, and other linear features present in the original image, introduce artifacts and noise, and otherwise reduce the usefulness of the image. The wavelets-based techniques afford multiple levels of resolution reduction and computational speed. This study is applicable to a wide range of reduced resolution applications in satellite imaging systems, including low resolution display, spaceborne browse, emergency image transmission, and real-time video downlinking.

Ellison, James↗
NASA NTRS2014

Elastic Model Transitions: a Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes both quadratically constrained least squares (LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical displacements. This approach is applicable to the simulation of the elastic behavior of launch vehicles and other structures that utilize multiple LTI finite element model (FEM) derived mode sets that are propagated throughout time. The time invariant nature of the elastic data for discrete segments of the launch vehicle trajectory presents a problem of how to properly transition between models while preserving motion across the transition. In addition, energy may vary between flex models when using a truncated mode set. The LSQI-DSM algorithm can accommodate significant changes in energy between FEM models and carries elastic motion across FEM model transitions. Compared with previous approaches, the LSQI-DSM algorithm shows improvements ranging from a significant reduction to a complete removal of transients across FEM model transitions as well as maintaining elastic motion from the prior state.

Jurenko, Robert J.↗
NASA NTRS2006

Chemical, Physical, and Mechanical Characterization of Isocyanate Cross-linked Amine-Modified Silica Aerogels

We describe a new mechanically strong lightweight porous composite material obtained by encapsulating the skeletal framework of amine-modified silica aerogels with polyurea. The conformal polymer coating preserves the mesoporous structure of the underlying silica framework and the thermal conductivity remains low at 0.041 plus or minus 0.001 W m(sup -1 K(sup -1). The potential of the new cross-linked silica aerogels for load-carrying applications was determined through characterization of their mechanical behavior under compression, three-point bending, and dynamic mechanical analysis (DMA). A primary glass transition temperature of 130 C was identified through DMA. At room temperature, results indicate a hyperfoam behavior where in compression cross-linked aerogels are linearly elastic under small strains (less than 4%) and then exhibit yield behavior (until 40% strain), followed by densification and inelastic hardening. At room temperature the compressive Young's modulus and the Poisson's ratio were determined to be 129 plus or minus 8 MPa and 0.18, respectively, while the strain at ultimate failure is 77% and the average specific compressive stress at ultimate failure is 3.89 x 10(exp 5) N m kg(sup -1). The specific flexural strength is 2.16 x 10(exp 4) N m kg(sup -1). Effects on the compressive behavior of strain rate and low temperature were also evaluated.

Katti, Atul↗
NASA NTRS2012

Discovery of a Plains Caldera Complex and Extinct Lava Lake in Arabia Terra, Mars: Implications for the Discovery of Additional Highland Volcanic Source Regions

Several irregularly shaped topographic depressions occur near the dichotomy boundary in northern Arabia Terra, Mars. The geomorphology of these features suggests that they formed by collapse, opposed to meteor impact. At least one depression (approx.55 by 85 km) displays geologic features indicating a complex, multi-stage collapse history. Features within and around the collapse structure indicate volcanic processes. The complex occurs within Hesperian ridged plains of likely volcanic origin and displays no crater rim or evidence for ejecta. Instead the depression consists of a series of circumferential graben and down-dropped blocks which also display upper surfaces similar to ridged plain lavas. Large blocks within the depression are tilted towards the crater center, and display graben that appear to have originally been linked with circumferential graben outside of the complex related to earlier collapse events. A nearly 700 m high mound exists along a graben within the complex that might be a vent. The deepest depression displays two sets of nearly continuous terraces, which we interpret as high-stands of a drained lava lake. These features appear similar to the black ledge described during the Kilauea Iki eruption in 1959. A lacustrine origin for the terraces seems unlikely because of the paucity of channels found in or around the depression that could be linked to aqueous surface processes. In addition, there is no obvious evidence for lacustrine sediments within the basin. Together with the presence of significant faulting that is indicative of collapse we conclude that this crater complex represents a large caldera formed in the Late Noachian to Early Hesperian. Other linear and irregular depressions in the region also might be linked to ancient volcanism. If that hypothesis is correct, it suggests that northern Arabia Terra could contain a large, previously unrecognized highland igneous province. Evacuation of magma via explosive and effusive activity produced localized collapse, might have contributed to nearby ridged plains, and pyroclastic materials erupted from these vents might have supplied sediments in fretted terrain and other deposits. The recognition of volcanoes within Arabia Terra expands the known extent of Noachian-Hesperian volcanism to cover much of the preserved martian highland crust.

Bleacher, Jacob↗
DOE OSTI2024

Pressure-tuned quantum criticality in the large-D antiferromagnet DTN

Strongly correlated spin systems can be driven to quantum critical points via various routes. In particular, gapped quantum antiferromagnets can undergo phase transitions into a magnetically ordered state with applied pressure or magnetic field, acting as tuning parameters. These transitions are characterized by z = 1 or z = 2 dynamical critical exponents, determined by the linear and quadratic low-energy dispersion of spin excitations, respectively. Employing high-frequency susceptibility and ultrasound techniques, we demonstrate that the tetragonal easy-plane quantum antiferromagnet NiCl 2 · 4SC(NH 2 ) 2 (aka DTN) undergoes a spin-gap closure transition at about 4.2 kbar, resulting in a pressure-induced magnetic ordering. The studies are complemented by high-pressure-electron spin-resonance measurements confirming the proposed scenario. Powder neutron diffraction measurements revealed that no lattice distortion occurs at this pressure and the high spin symmetry is preserved, establishing DTN as a perfect platform to investigate z = 1 quantum critical phenomena. The experimental observations are supported by DMRG calculations, allowing us to quantitatively describe the pressure-driven evolution of critical fields and spin-Hamiltonian parameters in DTN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Planetary-scale waves in the Venus atmosphere

Observed wave-like cloud features on Venus, which at times form a Y-like structure which encircles the planet, were modelled numerically. Linearized primitive equations for a shallow, hydrostatic atmosphere are defined, along with upper boundary layer conditions for the 4-6 day periods detected in the atmospheric response to forcing. The basic state of the Venus atmosphere was considered to display variations in static stability and the mean zonal wind as a function of altitude. Forcing was introduced over a wide range of frequencies in order to produce the target oscillation modes. Combination of a midlatitude Rossby wave and an equatorial Kelvin wave was found to yield the observed Y-shape, which could be preserved with nonlinear coupling.

Covey, C.↗
NASA NTRS1988

On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices

Conjugate gradient type methods are considered for the solution of large linear systems Ax = b with complex coefficient matrices of the type A = T + i(sigma)I where T is Hermitian and sigma, a real scalar. Three different conjugate gradient type approaches with iterates defined by a minimal residual property, a Galerkin type condition, and an Euclidian error minimization, respectively, are investigated. In particular, numerically stable implementations based on the ideas behind Paige and Saunder's SYMMLQ and MINRES for real symmetric matrices are proposed. Error bounds for all three methods are derived. It is shown how the special shift structure of A can be preserved by using polynomial preconditioning. Results on the optimal choice of the polynomial preconditioner are given. Also, some numerical experiments for matrices arising from finite difference approximations to the complex Helmholtz equation are reported.

Freund, Roland↗
DOE OSTI2024

Topological Optimization with Big Steps

Using persistent homology to guide optimization has emerged as a novel application of topological data analysis. Existing methods treat persistence calculation as a black box and backpropagate gradients only onto the simplices involved in particular pairs. We show how the cycles and chains used in the persistence calculation can be used to prescribe gradients to larger subsets of the domain. In particular, we show that in a special case, which serves as a building block for general losses, the problem can be solved exactly in linear time. This relies on another contribution of this paper, which eliminates the need to examine a factorial number of permutations of simplices with the same value. Here, we present empirical experiments that show the practical benefits of our algorithm: the number of steps required for the optimization is reduced by an order of magnitude.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2023

Electronic structure of Co 3 d states in the Kitaev material candidate honeycomb cobaltate Na 3 Co 2 SbO 6 probed with x-ray dichroism

The recent prediction that honeycomb lattices of Co 2+ (3d 7 ) ions could host dominant Kitaev interactions provides an exciting direction for exploration of new routes to stabilizing Kitaev’s quantum spin liquid in real materials. Na 3 Co 2 SbO 6 has been singled out as a potential material candidate provided that spin and orbital moments couple into a J eff = $\frac{1}{2}$ ground state, and that the relative strength of trigonal crystal field and spin-orbit coupling acting on Co ions can be tailored. Using x-ray linear dichroism (XLD) and x-ray magnetic circular dichroism (XMCD) experiments, alongside configuration interaction calculations, we confirm the counterintuitive positive sign of the trigonal crystal field acting on Co 2+ ions and test the validity of the J eff = $\frac{1}{2}$ description of the electronic ground state. Furthermore, the results lend experimental support to recent theoretical predictions that a compression (elongation) of CoO 6 octahedra along (perpendicular to) the trigonal axis would drive this cobaltate toward the Kitaev limit, assuming the J eff = $\frac{1}{2}$ character of the electronic ground state is preserved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗