Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Inverse dynamics of a 3 degree of freedom spatial flexible manipulator

A technique is presented for solving the inverse dynamics and kinematics of 3 degree of freedom spatial flexible manipulator. The proposed method finds the joint torques necessary to produce a specified end effector motion. Since the inverse dynamic problem in elastic manipulators is closely coupled to the inverse kinematic problem, the solution of the first also renders the displacements and rotations at any point of the manipulator, including the joints. Furthermore the formulation is complete in the sense that it includes all the nonlinear terms due to the large rotation of the links. The Timoshenko beam theory is used to model the elastic characteristics, and the resulting equations of motion are discretized using the finite element method. An iterative solution scheme is proposed that relies on local linearization of the problem. The solution of each linearization is carried out in the frequency domain. The performance and capabilities of this technique are tested through simulation analysis. Results show the potential use of this method for the smooth motion control of space telerobots.

Bayo, Eduardo↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Streamlined spatial and environmental expression signatures characterize the minimalist duckweed Wolffia australiana

Single-cell genomics permits a new resolution in the examination of molecular and cellular dynamics, allowing global, parallel assessments of cell types and cellular behaviors through development and in response to environmental circumstances, such as interaction with water and the light–dark cycle of the Earth. Here, we leverage the smallest, and possibly most structurally reduced, plant, the semiaquaticWolffia australiana, to understand dynamics of cell expression in these contexts at the whole-plant level. We examined single-cell-resolution RNA-sequencing data and foundWolffiacells divide into four principal clusters representing the above- and below-water-situated parenchyma and epidermis. Although these tissues share transcriptomic similarity with model plants, they display distinct adaptations thatWolffiahas made for the aquatic environment. Within this broad classification, discrete subspecializations are evident, with select cells showing unique transcriptomic signatures associated with developmental maturation and specialized physiologies. Assessing this simplified biological system temporally at two key time-of-day (TOD) transitions, we identify additional TOD-responsive genes previously overlooked in whole-plant transcriptomic approaches and demonstrate that the core circadian clock machinery and its downstream responses can vary in cell-specific manners, even in this simplified system. Distinctions between cell types and their responses to submergence and/or TOD are driven by expression changes of unexpectedly few genes, characterizingWolffiaas a highly streamlined organism with the majority of genes dedicated to fundamental cellular processes.Wolffiaprovides a unique opportunity to apply reductionist biology to elucidate signaling functions at the organismal level, for which this work provides a powerful resource.

Biochemistry & Molecular Biology↗
NASA NTRS1990

Low-frequency ionization-driven instability of a discrete auroral arc model

The low-frequency (time scales of tens of seconds) dynamics of the auroral ionosphere differs from that of the nonauroral ionosphere by the presence of strong, unstable space- and time-dependent ionospheric ionization produced by precipitating auroral electrons. If recombination is relatively unimportant (as at high ionospheric heights), it is shown that, in general, transport processes cannot remove this ionization as fast as it is created, and no equilibrium is possible. These nonequilibrium phenomena are investigated in the context of a nonlinear adiabatic auroral model, which has previously been studied in static situations. A linearized local perturbation analysis is given of what amounts to a current-driven E x B gradient-drift instability with an ionization source, as well as some exact nonlinear solutions valid in a finite but limited spatial region. These solutions show continuing motion of auroral potential and plasma density, as the aurora tries to shift its ionization problems from one place to another. The analysis gives clues to the possibility of generation of chaos and of fine-scale spatial structure.

Cornwall, John M.↗
NASA NTRS2005

Discrete element simulations of gravitational volcanic deformation: Deformation structures and geometries - 1

We have carried out two-dimensional particle dynamics simulations of granular piles subject to frictional Coulomb failure criteria to gain a first-order understanding of different modes of gravitational deformation within volcanoes. Under uniform basal and internal strength conditions, granular piles grow self-similarly, developing distinctive stratigraphies, morphologies, and structures. Piles constructed upon cohesive substrates exhibit particle avalanching, forming outward dipping strata and angle of repose slopes. Systematic decreases in basal strength lead to progressively deeper and steeper internal detachment faults and slip along a basal decollement; landslide forms grade from shallow slumps to deep-seated landslide and, finally, to axial subsidence and outward flank displacements, or volcanic spreading. Surface slopes decrease and develop concave up morphologies with decreasing decollement strength; depositional layers tilt progressively inward. Spatial variations in basal strength cause lateral transitions in pile structure, stratigraphy, and morphology. This approximation of volcanoes as Coulomb granular piles reproduces the richness of deformational structures and surface morphologies in many volcanic settings. The gentle slopes of Hawaiian volcanoes and Olympus Mons on Mars suggest weak basal decollements that enable volcanic spreading. High-angle normal faults, favored above weak decollements, are interpreted in both settings and may explain catastrophic sector collapse in Hawaii and broad aureole deposits surrounding Olympus Mons. In contrast, steeper slopes and shallow detachment faults predominate in the Canary Islands, thought to lack a weak decollement, favoring smaller, more frequent slope failures than predicted for Hawaii. The numerical results provide a useful predictive tool for interpreting dynamic behavior and associated geologic hazards of active volcanoes.

Morgan, Julia K.↗
DOE OSTI2024

Probabilistic Context Neighborhood model for lattices

Here we present the Probabilistic Context Neighborhood model designed for two-dimensional lattices as a variation of a Markov random field assuming discrete values. In this model, the neighborhood structure has a fixed geometry but a variable order, depending on the neighbors’ values. Our model extends the Probabilistic Context Tree model, originally applicable to one-dimensional space. It retains advantageous properties, such as representing the dependence neighborhood structure as a graph in a tree format, facilitating an understanding of model complexity. Furthermore, we adapt the algorithm used to estimate the Probabilistic Context Tree to estimate the parameters of the proposed model. We illustrate the accuracy of our estimation methodology through simulation studies. Additionally, we apply the Probabilistic Context Neighborhood model to spatial real-world data, showcasing its practical utility.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

How well is Martian ozone inferred with reflectance spectroscopy?

The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms in one of its two spectral channels every 3 seconds with a spectral resolution of 15 Angstroms and an effective field-of-view of approximately 300 km(exp 2). The only gaseous absorption in the 2000 to 3000 Angstrom region was assumed to come from the Hartley band system of ozone, and therefore the amount of ozone was inferred by fitting this absorption feature with laboratory data of ozone absorption. Mars O3 as inferred from these spectra shows strong seasonal and latitudinal variation, with column abundances ranging from 0.2 microns at equatorial latitudes to 60 microns over the northern winter polar latitudes. The detectability limit of the spectrometer was approximately 3 microns. I use a radiative transfer model based on the discrete ordinate method to calculate synthetic radiance spectra. When typical amounts of dust and cloud are present, significant underestimation of O3 occurs. A factor of 3 times as much O3 is needed to generate the same spectrum for cloudy, dusty atmospheres as for a clear atmosphere. If the scattering properties of Martian clouds and dust were well known, then their appearance would not be a problem, as a model would be capable of retrieving the O3 abundance. However, these properties are not well known, which raises doubts about the effectiveness of the current UV spectroscopy technique used to measure O3. Spatial and temporal variability in temperature and water vapor have been claimed to account for the scatter of the data points. However, water vapor is a small source of odd hydrogen in the winter polar atmosphere, and may not account for most of the variability. Masking by clouds and dust may also account for some of the observed O3 variability, because the nature and opacity of the clouds and dust in the polar hood change dramatically in latitude and even on a day-to-day basis. As the maximum O3 abundance resides near the surface, spacecraft must be able to observe through the entire cloud and dust abundance in order to actually see the total O3 column abundance. If reflectance spectroscopy is used, as on Mariner 9, then the cloud and the airborne dust must be traversed twice; first by the incoming solar flux down to the surface, and then once again upon reflection from the surface out to the spacecraft.

Lindner, Bernhard Lee↗
NASA NTRS2013

Gesture-Based Robot Control with Variable Autonomy from the JPL Biosleeve

This paper presents a new gesture-based human interface for natural robot control. Detailed activity of the user's hand and arm is acquired via a novel device, called the BioSleeve, which packages dry-contact surface electromyography (EMG) and an inertial measurement unit (IMU) into a sleeve worn on the forearm. The BioSleeve's accompanying algorithms can reliably decode as many as sixteen discrete hand gestures and estimate the continuous orientation of the forearm. These gestures and positions are mapped to robot commands that, to varying degrees, integrate with the robot's perception of its environment and its ability to complete tasks autonomously. This flexible approach enables, for example, supervisory point-to-goal commands, virtual joystick for guarded teleoperation, and high degree of freedom mimicked manipulation, all from a single device. The BioSleeve is meant for portable field use; unlike other gesture recognition systems, use of the BioSleeve for robot control is invariant to lighting conditions, occlusions, and the human-robot spatial relationship and does not encumber the user's hands. The BioSleeve control approach has been implemented on three robot types, and we present proof-of-principle demonstrations with mobile ground robots, manipulation robots, and prosthetic hands.

BioSleeve↗
DOE OSTI2021

Characterization of a collimated neutron imager for low-rate fast neutron imaging

Spatial localization of special nuclear materials (SNM) via their neutron signatures amidst background requires knowledge of the background neutron environment or a means of separating a source from background based on low amounts of information. This requirement has created the need for characterizing the spatial distribution of the cosmogenic neutron background. Neutron scatter cameras have been developed and optimized for rapid detection of high activity sources, but have low imaging efficiency, making it difficult to use them to characterize low rate diffuse sources, such as the neutron background. The Low Intensity Neutron Imaging System (LINIS) is a collimated neutron imager that has been designed and optimized for imaging diffuse cosmogenic neutron background in the energy range of 0.5–15 MeV. LINIS operates using 16 liquid scintillation detectors shielded by ultra-high molecular weight polyethylene cylindrical collimators in a staggered orientation and rotates to 7 discrete positions, giving it roughly 2π sensitivity. Finally, LINIS has been characterized using (α, n) and fission neutron sources using two imaging techniques for neutron source localization, simple backprojection and Maximum Likelihood Expectation Maximization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Topological Circular Dichroism in Chiral Multifold Semimetals

Uncovering the physical contents of the nontrivial topology of quantum states is a critical problem in condensed matter physics. Here, we study the topological circular dichroism in chiral semimetals using linear response theory and first-principles calculations. We show that, when the low-energy spectrum respects emergent SO(3) rotational symmetry, topological circular dichroism is forbidden for Weyl fermions, and thus is unique to chiral multifold fermions. This is a result of the selection rule that is imposed by the emergent symmetry under the combination of particle-hole conjugation and spatial inversion. Using first-principles calculations, we predict that topological circular dichroism occurs in CoSi for photon energy below about 0.2 eV. In conclusion, our Letter demonstrates the existence of a response property of unconventional fermions that is fundamentally different from the response of Dirac and Weyl fermions, motivating further study to uncover other unique responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Perturbative Stability and Error-Correction Thresholds of Quantum Codes

Topologically ordered phases are stable to local perturbations, and topological quantum error-correcting codes enjoy thresholds to local errors. We connect the two notions of stability by constructing classical statistical mechanics models for decoding general Calderbank-Shor-Steane codes and classical linear codes. Our construction encodes correction success probabilities under uncorrelated bit-flip and phase-flip errors, and simultaneously describes a generalized ℤ 2 lattice-gauge theory with quenched disorder. We observe that the clean limit of the latter is precisely the discretized imaginary-time path integral of the corresponding quantum code Hamiltonian when the errors are turned into a perturbative 𝑋 or 𝑍 magnetic field. Motivated by error-correction considerations, we define general order parameters for all such generalized ℤ 2 lattice-gauge theories, and show that they are generally lower bounded by success probabilities of error correction. For CSS codes satisfying the low-density parity-check condition and with a sufficiently large code distance, we prove the existence of a low-temperature ordered phase of the corresponding lattice-gauge theories, particularly for those lacking Euclidean spatial locality and/or when there is a nonzero code rate. We further argue that these results provide evidence for stable phases in the corresponding perturbed quantum Hamiltonians, obtained in the limit of continuous imaginary time. To do so, we distinguish space- and timelike defects in the lattice-gauge theory. A high free-energy cost of spacelike defects corresponds to a successful “memory experiment” and suppresses the energy splitting among the ground states, while a high free-energy cost of timelike defects corresponds to a successful “stability experiment” and points to a nonzero gap to local excitations.

quantum error correction↗
NASA NTRS1987

Particle simulation of plasmas on the massively parallel processor

Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.

Gledhill, I. M. A.↗
NASA NTRS1985

Multidimensional Modeling of Atmospheric Effects and Surface Heterogeneities on Remote Sensing

The overall goal of this project is to establish a modeling capability that allows a quantitative determination of atmospheric effects on remote sensing including the effects of surface heterogeneities. This includes an improved understanding of aerosol and haze effects in connection with structural, angular, and spatial surface heterogeneities. One important objective of the research is the possible identification of intrinsic surface or canopy characteristics that might be invariant to atmospheric perturbations so that they could be used for scene identification. Conversely, an equally important objective is to find a correction algorithm for atmospheric effects in satellite-sensed surface reflectances. The technical approach is centered around a systematic model and code development effort based on existing, highly advanced computer codes that were originally developed for nuclear radiation shielding applications. Computational techniques for the numerical solution of the radiative transfer equation are adapted on the basis of the discrete-ordinates finite-element method which proved highly successful for one and two-dimensional radiative transfer problems with fully resolved angular representation of the radiation field.

Gerstl, S. A. W.↗
DOE OSTI2022

Circular Economy Modeling Efforts at NREL: Circular Economy Lifecycle Analysis and Visualization (CELAVI) and the Circular Economy Agent-Based Model (ABM)

Circular Economy (CE) aims at decoupling human activities from economic growth and resource use. While CE economic and environmental benefits are often touted by its proponents, increased circularity and improved sustainability do not necessarily go hand-in-hand. The Circular Economy Lifecycle Assessment and Visualization (CELAVI) framework simulates the transition of renewable energy technology supply chains towards circularity and quantifies spatially explicit environmental impacts and supply chain costs. Using wind blade end-of-life (EOL) management as a case study, this presentation will highlight CELAVI's features. Furthermore, envisioned extension to the framework could improve the stakeholder decision model by accounting for logistical and other factors. The CELAVI framework could be leveraged to answer crucial questions regarding renewables EOL management.

circular economy↗
NASA NTRS1994

The effect of parallel currents on auroral micropulsations

Field aligned currents play an important role in the global coupling between the magnetosphere and the ionosphere and in their relationship to the auroral phenomena. Moreover, there exists evidence that resonant oscillations are related to large-scale Birkeland currents. The spatial confinement of the field-aligned currents forms an inhomogeneous system susceptible to low-frequency oscillations, which can be excited due to periodic variations in the solar wind pressure or to the Kelvin-Helmholtz (KH) instability. In this paper we present a study of ultralow-frequency (ULF) oscillations in an inhomogeneous magnetic field formed by a large-scale current. We investigate the effects of the field-aligned currents on the generation of localized Alfven waves. The field oscillations are described by an eigenvalue wave equation which includes the effects of the field aligned currents and which produces a discrete spectrum of Alfven waves. These waves are observed mainly in three regions of the magnetosphere: in the magnetosheath, in the polar cusp, and in the plasmasphere. In the present study we limit our investigation to the auroral region.

Tavares, M.↗
NASA NTRS1999

Directional Solidification of a Binary Alloy into a Cellular Convective Flow: Localized Morphologies

A steady, two dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behavior of the moving front is described by a slow, spatial-temporal dynamics obtained through a multiple-scale analysis. The resulting system has a "parametric-excitation" structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection stabilizes two dimensional disturbances oriented with the flow, but destabilizes three dimensional disturbances in general. When the flow is weak, the morphological instability behaves incommensurably to the flow wavelength, but becomes quantized and forced to fit into the flow-box as the flow gets stronger. At large flow magnitudes the instability is localized, confined in narrow envelopes with cells traveling with the flow. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundary and asymptotics are obtained by the WKB analysis.

Chen, Y.- J.↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗