Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential geometry”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pressure–strain interaction. II. Decomposition in magnetic field-aligned coordinates

In weakly collisional and collisionless magnetized plasmas, the pressure–strain interaction describes the rate of conversion between bulk flow and thermal energy density. In this study, we derive an analytical expression for the pressure–strain interaction in a coordinate system with an axis aligned with the local magnetic field. The result is eight groups of terms corresponding to different physical mechanisms that can contribute to the pressure–strain interaction. We provide a physical description of each term. Here, the results are immediately of interest to weakly collisional and collisionless magnetized plasmas and the fundamental processes that happen therein, including magnetic reconnection, magnetized plasma turbulence, and collisionless shocks. The terms in the field-aligned coordinate decomposition are likely accessible to measurement with satellite observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multiscale modeling of solute diffusion in triblock copolymer membranes

We develop a multiscale simulation model for diffusion of solutes through porous triblock copolymer membranes. The approach combines two techniques: self-consistent field theory (SCFT) to predict the structure of the self-assembled, solvated membrane and on-lattice kinetic Monte Carlo (kMC) simulations to model diffusion of solutes. Solvation is simulated in SCFT by constraining the glassy membrane matrix while relaxing the brush-like membrane pore coating against the solvent. The kMC simulations capture the resulting solute spatial distribution and concentration-dependent local diffusivity in the polymer-coated pores; we parameterize the latter using particle-based simulations. We apply our approach to simulate solute diffusion through nonequilibrium morphologies of a model triblock copolymer, and we correlate diffusivity with structural descriptors of the morphologies. We also compare the model’s predictions to alternative approaches based on simple lattice random walks and find our multiscale model to be more robust and systematic to parameterize. Furthermore, our multiscale modeling approach is general and can be readily extended in the future to other chemistries, morphologies, and models for the local solute diffusivity and interactions with the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Periodic Korteweg–de Vries soliton potentials generate quasisymmetric magnetic field strength in a finite plasma- β equilibrium

Quasisymmetry (QS) is a hidden symmetry of the magnetic field strength, B , that enables effective confinement of charged particles in a fully three-dimensional (3D) toroidal plasma equilibrium. Such equilibria are typically modeled by the ideal magnetohydrostatic (MHS) equations. The nonlinear, overdetermined nature of the QS MHS equations severely complicates our understanding of the interplay between 3D shaping, equilibrium properties such as pressure and rotational transform, and B . Progress has been made through expansions near the magnetic axis; however, a more comprehensive theory is desirable. Using a combination of analysis and regression on a large dataset of numerically optimized quasisymmetric stellarators, we demonstrate that there is a hidden lower dimensionality of B on a magnetic flux surface with connections to the theory of periodic solitons. We show that B on a flux surface is determined by three or at most four flux functions, each of which determines a critical value of the derivative of B along the field line. While being consistent with the near-axis models, our results are global and hold even on the last closed flux surface.

Differential geometry↗
DOE OSTI2025

The absence of ray-effects in the discrete ordinate solution to the transport equation in spherical coordinates in multi-dimensions

The streaming operator of the transport equation is derived for spherical coordinates by starting from Newton’s second law for a free particle expressed in spherical coordinates. We shall show that the partial derivatives with respect to the velocity variables of the particle, which are absent in the Cartesian coordinate formulation of the transport equation, arise in the spherical coordinate formulation of the transport equation in response to the centrifugal force which prevents a free particle from ‘falling into the origin’ of the coordinate system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

SCExAO/CHARIS Near-infrared Integral Field Spectroscopy of the HD 15115 Debris Disk

We present new, near-infrared (1.1–2.4 μm) high-contrast imaging of the debris disk around HD 15115 with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS). The SCExAO/CHARIS resolves the disk down to ρ ∼ 0.″2 (r {sub proj} ∼ 10 au), a factor of ∼3–5 smaller than previous recent studies. We derive a disk position angle of PA ∼ 279.°4–280.°5 and an inclination of i ∼ 85.°3–86.2.°. While recent SPHERE/IRDIS imagery of the system could suggest a significantly misaligned two-ring disk geometry, CHARIS imagery does not reveal conclusive evidence for this hypothesis. Moreover, optimizing models of both one- and two-ring geometries using differential evolution, we find that a single ring having a Hong-like scattering phase function matches the data equally well within the CHARIS field of view (ρ ≲ 1″). The disk’s asymmetry, well evidenced at larger separations, is also recovered; the west side of the disk appears, on average, around 0.4 mag brighter across the CHARIS bandpass between 0.″25 and 1″. Comparing Space Telescope Imaging Spectrograph (STIS) 50CCD optical photometry (2000–10500 Å) with CHARIS near-infrared photometry, we find a red (STIS/50CCD−CHARIS broadband) color for both sides of the disk throughout the 0.″4–1″ region of overlap, in contrast to the blue color reported at similar wavelengths for regions exterior to ∼2″. Further, this color may suggest a smaller minimum grain size than previously estimated at larger separations. Finally, we provide constraints on planetary companions and discuss possible mechanisms for the observed inner disk flux asymmetry and color.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Geometric differences of the mitral apparatus between ischemic and dilated cardiomyopathy with significant mitral regurgitation: real-time three-dimensional echocardiography study

BACKGROUND: This study was conducted to elucidate the geometric differences of the mitral apparatus in patients with significant mitral regurgitation caused by ischemic cardiomyopathy (ICM-MR) and by idiopathic dilated cardiomyopathy (DCM-MR) by use of real-time 3D echocardiography (RT3DE). METHODS AND RESULTS: Twenty-six patients with ICM-MR caused by posterior infarction, 18 patients with DCM-MR, and 8 control subjects were studied. With the 3D software, commissure-commissure plane and 3 perpendicular anteroposterior (AP) planes were generated for imaging the medial, central, and lateral sides of the mitral valve (MV) during mid systole. In 3 AP planes, the angles between the annular plane and each leaflet (anterior, Aalpha; posterior, Palpha) were measured. In ICM-MR, Aalpha measured in the medial and central planes was significantly larger than that in the lateral plane (39+/-5 degrees, 34+/-6 degrees, and 27+/-5 degrees, respectively; P<0.01), whereas Palpha showed no significant difference in any of the 3 AP planes (61+/-7 degrees, 57+/-7 degrees, and 56+/-7 degrees, P>0.05). In DCM-MR, both Aalpha (38+/-8 degrees, 37+/-9 degrees, and 36+/-7 degrees, P>0.05) and Palpha (59+/-6 degrees, 58+/-5 degrees, and 57+/-6 degrees, P>0.05) revealed no significant differences in the 3 planes. CONCLUSIONS: The pattern of MV deformation from the medial to the lateral side was asymmetrical in ICM-MR, whereas it was symmetrical in DCM-MR. RT3DE is a helpful tool for differentiating the geometry of the mitral apparatus between these 2 different types of functional mitral regurgitation.

Evaluation Studies↗
NASA NTRS2020

SCExAO/CHARIS Near-infrared Integral Field Spectroscopy of the HD 15115 Debris Disk

We present new, near-infrared (1.1–2.4 μm) high-contrast imaging of the debris disk around HD 15115 with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS). The SCExAO/CHARIS resolves the disk down to ρ ∼ 0."2 (r(proj) ∼ 10 au), a factor of ∼3–5 smaller than previous recent studies. We derive a disk position angle of PA ∼ 279.°4–280.°5 and an inclination of i ∼ 85.°3–86.2.°. While recent SPHERE/IRDIS imagery of the system could suggest a significantly misaligned two-ring disk geometry, CHARIS imagery does not reveal conclusive evidence for this hypothesis. Moreover, optimizing models of both one- and two-ring geometries using differential evolution, we find that a single ring having a Hong-like scattering phase function matches the data equally well within the CHARIS field of view (ρ ≲ 1''). The disk's asymmetry, well evidenced at larger separations, is also recovered; the west side of the disk appears, on average, around 0.4 mag brighter across the CHARIS bandpass between 0farcs25 and 1''. Comparing Space Telescope Imaging Spectrograph (STIS) 50CCD optical photometry (2000–10500 Å) with CHARIS near-infrared photometry, we find a red (STIS/50CCD−CHARIS broadband) color for both sides of the disk throughout the 0farcs4–1'' region of overlap, in contrast to the blue color reported at similar wavelengths for regions exterior to ∼2''. Further, this color may suggest a smaller minimum grain size than previously estimated at larger separations. Finally, we provide constraints on planetary companions and discuss possible mechanisms for the observed inner disk flux asymmetry and color.

Kellen Lawson↗
DOE OSTI2024

Loops of loops expansion in the amplituhedron

We study a novel geometric expansion for scattering amplitudes in the planar sector of $\mathcal{N}$ = 4 super Yang-Mills theory, in the context of the Amplituhedron which reproduces the all-loop integrand as a canonical differential form on the positive geometry. In a paper by Arkani-Hamed, Henn and one of the authors, it was shown that this result can be recast in terms of negative geometries with a certain hierarchy of loops (closed cycles) in the space of loop momenta, represented by lines in momentum twistor space. One can then calculate an all-loop order result in the approximation where only tree graphs in the space of all loops are considered. Furthermore, using differential equation methods, it is possible to calculate and resum integrated expressions and obtain strong coupling results. In this paper, we provide a more general framework for the ‘loops of loops’ expansion and outline a powerful method for the determination of differential forms for higher-order geometries. We solve the problem completely for graphs with one internal cycle, but the method can be used more generally for other geometries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1983

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗
DOE OSTI2022

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗
DOE OSTI2023

Data-driven Whitney forms for structure-preserving control volume analysis

Control volume analysis models physics via the exchange of generalized fluxes between subdomains. Here, we introduce a scientific machine learning framework adopting a partition of unity architecture to identify physically-relevant control volumes, with generalized fluxes between subdomains encoded via Whitney forms. The approach provides a differentiable parameterization of geometry which may be trained in an end-to-end fashion to extract reduced models from full field data while exactly preserving physics. The architecture admits a data-driven finite element exterior calculus allowing discovery of mixed finite element spaces with closed form quadrature rules. An equivalence between Whitney forms and graph networks reveals that the geometric problem of control volume learning is equivalent to an unsupervised graph discovery problem. The framework is developed for manifolds in arbitrary dimension, with examples provided for H(div) problems in $\mathbb{R}$ establishing convergence and structure preservation properties. Finally, we consider a lithium-ion battery problem where we discover a reduced finite element space encoding transport pathways from high-fidelity microstructure resolved simulations. The approach reduces the 5.89M finite element simulation to 136 elements while reproducing pressure to under 0.1% error and preserving conservation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Kinetics of particles with short-range interactions

Self-assembly is one of the grand challenges of the 21st century – as the devices and materials we would like to build become too complex or small-scale for top-down manufacturing to be efficient, it is increasingly important to find ways to create these through bottom-up, dynamical approaches. Many particles used in self-assembly have very short-ranged attractive interactions, making simulations expensive or impossible. This proposal develops a set of conceptual and computational tools to study the dynamics of self-assembly for particles with short-ranged interactions, harnessing ideas in differential and computational geometry, and stochastic analysis, to accelerate simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1992

Ejecta patterns of Meteor Crater, Arizona derived from the linear un-mixing of TIMS data and laboratory thermal emission spectra

Accurate interpretation of thermal infrared data depends upon the understanding and removal of complicating effects. These effects may include physical mixing of various mineralogies and particle sizes, atmospheric absorption and emission, surficial coatings, geometry effects, and differential surface temperatures. The focus is the examination of the linear spectral mixing of individual mineral or endmember spectra. Linear addition of spectra, for particles larger than the wavelength, allows for a straight-forward method of deconvolving the observed spectra, predicting a volume percent of each endmember. The 'forward analysis' of linear mixing (comparing the spectra of physical mixtures to numerical mixtures) has received much attention. The reverse approach of un-mixing thermal emission spectra was examined with remotely sensed data, but no laboratory verification exists. Understanding of the effects of spectral mixing on high resolution laboratory spectra allows for the extrapolation to lower resolution, and often more complicated, remotely gathered data. Thermal Infrared Multispectral Scanner (TIMS) data for Meteor Crater, Arizona were acquired in Sep. 1987. The spectral un-mixing of these data gives a unique test of the laboratory results. Meteor Crater (1.2 km in diameter and 180 m deep) is located in north-central Arizona, west of Canyon Diablo. The arid environment, paucity of vegetation, and low relief make the region ideal for remote data acquisition. Within the horizontal sedimentary sequence that forms the upper Colorado Plateau, the oldest unit sampled by the impact crater was the Permian Coconino Sandstone. A thin bed of the Toroweap Formation, also of Permian age, conformably overlays the Coconino. Above the Toroweap lies the Permian Kiabab Limestone which, in turn, is covered by a thin veneer of the Moenkopi Formation. The Moenkopi is Triassic in age and has two distinct sub-units in the vicinity of the crater. The lower Wupatki member is a fine-grained sandstone, while the upper Moqui member is a fissile siltstone. Ejecta from these units are preserved as inverted stratigraphy up to 2 crater radii from the rim. The mineralogical contrast between the units, relative lack of post-emplacement erosion and ejecta mixing provide a unique site to apply the un-mixing model. Selection of the aforementioned units as endmembers reveals distinct patterns in the ejecta of the crater.

Ramsey, Michael S.↗
NASA NTRS2006

Dynamic Features in and nearby a Prominence

We will discuss the thermal and dynamic properties of dynamic structures in and around a prominence channel observed on the limb on 17 April 2003. Observations were taken with the Solar and Heliospheric Observatory's Solar Ultraviolet Measurements of Emitted Radiation (SOHO/SUMER) in lines formed at temperatures from 80,000 to 1.6 million K. The instrument was pointed to a single location and took a series of 90 s exposures. Two-dimensional context was provided by the Transition Region and Coronal Explorer (TRACE) in the UV and EUV and by the Kanzelhohe Solar Observatory in H-alpha. Two dynamic features were studied in depth: a prominence activation and repeated motions in a loop nearby the prominence. We were able to calculate three-dimensional geometries and trajectories, differential emission measure, and limits on the average density, kinetic and thermal energies. These observations provide important tests for models of dynamics in prominences and cool (approx.10(exp 5) K) loops.

Kucera, T.↗