Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector field”

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 145 records · Page 8

Alternative formulations of magnetospheric plasma electrodynamics

The fundamental equations of magnetospheric plasma electrodynamics are considered from a theoretical standpoint that stresses the basic equivalence of various seemingly different formal representations. The mathematical properties of vector fields are reviewed, and their implications in electrodynamics are studied. The irrotational and solenoidal parts of the electric field are associated with two physically distinct types of sources. Relativistic covariance and gauge invariance in electromagnetic theory are reviewed and discussed in the context of an approach in which the mathematical properties of vector fields are taken as primary concepts. Special attention is given to the use and interpretation of the Coulomb gauge potential functions. This choice of gauge is sometimes regarded with undue suspicion, possibly because of a certain paradox concerning causality. The paradox is discussed and resolved. Useful properties of the Coulomb gauge are identified. These need not be limited to the case of slow time variations and can extend beyond the limits of validity of ideal MHD theory.

Cragin, B. L.↗

The seasonal variation of the D region as inferred from propagation characteristics of LF radio waves

The propagation data of JG2AS 40 kHz (Japanese Standard Frequency), Loran C 100 kHz radio waves, and meteorological data were analyzed to study the association of propagation characteristics of LF radio waves with the atmospheric circulation in the mesosphere. The monthly averaged electric fields were depicted on the complex plane for typical summer and winter months, June and November. The locus traced out by the electric field vector during daytime is nearly circular. This is because during daytime the amplitude of the sky wave remains nearly constant while its phase changes in accord with the height change of the reflection layer, and thus the electric field vector traces out a circular locus with its center at the tip of the supposed ground wave vector. The locus has a loop during the sunrise or sunset period, which seems to arise from interference of two waves reflected by two different layers. In June the amplitude of the sky wave decreases rapidly before the dawn or increases after the dusk. In November such rapid change is not observed. During nighttime, the sky wave phase changes in such a way as to suggest that the reflection height moves upwards with time before midnight or lowers after midnight in November. In June it changes similarly before midnight, but after midnight it varies erratically. These characteristics are closely related to the structure of the D region, which is clearly shown by simulating the loci traced out by electric fields.

Ishimine, T.↗

Flare activity, sunspot motions, and the evolution of vector magnetic fields in Hale region 17244

The magnetic and dynamical circumstances leading to the 1B/M4 flare of November 5, 1980 are studied, and a strong association is found between the buildup of magnetic shear and the onset of flare activity within the active region. The development of shear, as observed directly in vector magnetograms, is consistent in detail with the dynamical history of the active region and identifies the precise location of the optical and hard-X-ray kernels of the flare emission.

Neidig, Donald F.↗

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗

Hamiltonian indices and rational spectral densities

Several (global) topological properties of various spaces of linear systems, particularly symmetric, lossless, and Hamiltonian systems, and multivariable spectral densities of fixed McMillan degree are announced. The study is motivated by a result asserting that on a connected but not simply connected manifold, it is not possible to find a vector field having a sink as its only critical point. In the scalar case, this is illustrated by showing that only on the space of McMillan degree = /Cauchy index/ = n, scalar transfer functions can one define a globally convergent vector field. This result holds both in discrete-time and for the nonautonomous case. With these motivations in mind, theorems of Bochner and Fogarty are used in showing that spaces of transfer functions defined by symmetry conditions are, in fact, smooth algebraic manifolds.

Byrnes, C. I.↗

Surface representations of two- and three-dimensional fluid flow topology

We discuss our work using critical point analysis to generate representations of the vector field topology of numerical flow data sets. Critical points are located and characterized in a two-dimensional domain, which may be either a two-dimensional flow field or the tangential velocity field near a three-dimensional body. Tangent curves are then integrated out along the principal directions of certain classes of critical points. The points and curves are linked to form a skeleton representing the two-dimensional vector field topology. When generated from the tangential velocity field near a body in a three-dimensional flow, the skeleton includes the critical points and curves which provide a basis for analyzing the three-dimensional structure of the flow separation. The points along the separation curves in the skeleton are used to start tangent curve integrations to generate surfaces representing the topology of the associated flow separations.

Helman, James L.↗

On the alignment of plasma anisotropies and the magnetic field direction in the solar wind

One year's Imp 6 solar wind plasma and magnetic field data are examined to determine whether anisotropies in particle velocity distributions are aligned with the measured interplanetary magnetic field vector. Alignment of components in the analysis plane was generally found to be excellent whenever plasma parameter magnitudes were larger than determination uncertainties, although some spread exists (typical rms approximately equal to 10 deg). By assuming cylindrical symmetry about the simultaneously measured magnetic field vector during the 1-year interval under study, three-dimensional values of selected solar wind plasma thermal parameters were constructed from the two-dimensional plasma measurements, and the statistical properties of their distributions have been tabulated.

Asbridge, J. R.↗

TROPHY: A Topologically Robust Physics-Informed Tracking Framework for Tropical Cyclones

Tropical cyclones (TCs) are among the most destructive weather systems. Realistically and efficiently detecting and tracking TCs are critical for assessing their impacts and risks. In particular, the eye is a signature feature of a mature TC. Therefore, knowing the eyes’ locations and movements is crucial for both operational weather forecasts and climate risk assessments. Recently, a multilevel robustness framework has been introduced to study the critical points of time-varying vector fields. The framework quantifies the robustness (i.e., structural stability) of critical points across varying neighborhoods. By relating the multilevel robustness with critical point tracking, the framework has demonstrated its potential in cyclone tracking. An advantage is that it identifies cyclonic features using only 2D wind vector fields, which is encouraging as most tracking algorithms require multiple dynamic and thermodynamic variables at different altitudes. A disadvantage is that the framework does not scale well computationally for datasets containing a large number of cyclones. Herein this paper introduces a topologically robust physics-informed tracking framework (TROPHY) for TC tracking. The main idea is to integrate physical knowledge of TC to drastically improve the computational efficiency of multilevel robustness framework for large-scale climate datasets. First, during preprocessing, we propose a physics-informed feature selection strategy to filter 90% of critical points that are short-lived and have low stability, thus preserving good candidates for TC tracking. Second, during in-processing, we impose constraints during the multilevel robustness computation to focus only on physics-informed neighborhoods of TCs. We apply TROPHY to 30 years of 2D wind fields from reanalysis data in ERA5 and generate a number of TC tracks. In comparison with the observed tracks, we demonstrate that TROPHY can capture TC characteristics (e.g., frequency, intensity, duration, latitudes with maximum intensity, and genesis) that are comparable to and sometimes even better than a well-validated TC tracking algorithm that requires multiple dynamic and thermodynamic scalar fields.

97 MATHEMATICS AND COMPUTING↗

Effect of polarized radiative transfer on the Hanle magnetic field determination in prominences: Analysis of hydrogen H alpha line observations at Pic-du-Midi

The linear polarization of the Hydrogen H alpha line of prominences has been computed, taking into account the effect of a magnetic field (Hanle effect), of the radiative transfer in the prominence, and of the depolarization due to collisions with the surrounding electrons and protons. The corresponding formalisms are developed in a forthcoming series of papers. In this paper, the main features of the computation method are summarized. The results of computation have been used for interpretation in terms of magnetic field vector measurements from H alpha polarimetric observations in prominences performed at Pic-du-Midi coronagraph-polarimeter. Simultaneous observations in one optically thin line (He I D(3)) and one optically thick line (H alpha) give an opportunity for solving the ambiguity on the field vector determination.

Bommier, V.↗

Demonstrating the viability of Lagrangian in situ reduction on supercomputers

Performing exploratory analysis and visualization of large-scale time-varying computational science applications is challenging due to inaccuracies that arise from under-resolved data. In recent years, Lagrangian representations of the vector field computed using in situ processing are being increasingly researched and have emerged as a potential solution to enable exploration. However, prior works have offered limited estimates of the encumbrance on the simulation code as they consider “theoretical” in situ environments. Further, the effectiveness of this approach varies based on the nature of the vector field, benefitting from an in-depth investigation for each application area. With this study, an extended version of Sane et al. (2021), we contribute an evaluation of Lagrangian analysis viability and efficacy for simulation codes executing at scale on a supercomputer. We investigated previously unexplored cosmology and seismology applications as well as conducted a performance benchmarking study by using a hydrodynamics mini-application targeting exascale computing. Here, to inform encumbrance, we integrated in situ infrastructure with simulation codes, and evaluated Lagrangian in situ reduction in representative homogeneous and heterogeneous HPC environments. To inform post hoc accuracy, we conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. Additionally, our study contributes cost estimates for distributed-memory post hoc reconstruction. In all, we demonstrate viability for each application — data reduction to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 90% of our experiments.

97 MATHEMATICS AND COMPUTING↗

Multipole moments of stellar oscillation modes

The oscillating mass 2(exp l)-pole moment, M(sub nl), of a star in a given (normalized) oscillation mode determines the energy that can be absorbed by the mode in a tidal interaction and the power radiated by the mode in gravitational waves, both of which are proportional to (absolute value of M(sub nl))(exp 2). The coefficients in the expansion of the vector fields del(r(exp l)Y(sub lm)(theta, phi)) in terms of the displacement fields of modes of given l and m are proportional to M(sub nl). This expansion leads to a sum rule sum over n(absolute value of M(sub nl))(exp 2) = constant. For stars of weak to moderate central condensation (such as neutron stars), the f-mode is well approximated by the vector field being expanded, and therefore it takes the lion's share of the sum. Thus the multipole moments of all other modes must be small. In there numerical evaluation, it is necessary to know the shape of the eigenfunctions quite precisly, since a small f-mode contamination can significantly increase the obtained values. This contamination occurs in some `hybrid' numerical computations of neutron star oscillations with relativistic equilibrium stars and Newtonian dynamics (e.g., McDermott et al. 1988). In this case, it is due to a slight inconsistency in the models and leads to a large overestimate of the power radiated in gravitational waves by modes other than the f-mode, although their oscillation periods are nearly unaffected.

Reisenegger, Andreas↗

Detection of open cluster rotation fields from Gaia EDR3 proper motions

Context: Most stars from in groups which with time disperse, building the field population of their host galaxy. In the Milky Way, open clusters have been continuously forming in the disk up to the present time, providing it with stars spanning a broad range of ages and masses. Observations of the details of cluster dissolution are, however, scarce. One of the main difficulties is obtaining a detailed characterisation of the internal cluster kinematics, which requires very high-quality proper motions. For open clusters, which are typically loose groups with tens to hundreds of members, there is the additional difficulty of inferring kinematic structures from sparse and irregular distributions of stars. Aims: Here, we aim to analyse internal stellar kinematics of open clusters, and identify rotation, expansion, or contraction patterns. Methods: We use Gaia Early Data Release 3 (EDR3) astrometry and integrated nested Laplace approximations to perform vector-field inference and create spatio-kinematic maps of 1237 open clusters. The sample is composed of clusters for which individual stellar memberships were already known, thus minimising contamination from field stars in the velocity maps. Projection effects were corrected using EDR3 data complemented with radial velocities from Gaia Data Release 2 and other surveys. Results: We report the detection of rotation patterns in eight open clusters. Nine additional clusters display possible rotation signs. We also observe 14 expanding clusters, with 15 other objects showing possible expansion patterns. Contraction is evident in two clusters, with one additional cluster presenting a more uncertain detection. In total, 53 clusters are found to display kinematic structures. Within these, elongated spatial distributions suggesting tidal tails are found in five clusters. These results indicate that the approach developed here can recover kinematic patterns from noisy vector fields, as those from astrometric measurements of open clusters or other stellar or galactic populations, thus offering a powerful probe for exploring the internal kinematics and dynamics of these types of objects.

79 ASTRONOMY AND ASTROPHYSICS↗

Scalable In Situ Computation of Lagrangian Representations via Local Flow Maps

In situ computation of Lagrangian flow maps to enable post hoc time-varying vector field analysis has recently become an active area of research. However, the current literature is largely limited to theoretical settings and lacks a solution to address scalability of the technique in distributed memory. To improve scalability, we propose and evaluate the benefits and limitations of a simple, yet novel, performance optimization. Our proposed optimization is a communication-free model resulting in local Lagrangian flow maps, requiring no message passing or synchronization between processes, intrinsically improving scalability, and thereby reducing overall execution time and alleviating the encumbrance placed on simulation codes from communication overheads. To evaluate our approach, we computed Lagrangian flow maps for four time-varying simulation vector fields and investigated how execution time and reconstruction accuracy are impacted by the number of GPUs per compute node, the total number of compute nodes, particles per rank, and storage intervals. Our study consisted of experiments computing Lagrangian flow maps with up to 67M particle trajectories over 500 cycles and used as many as 2048 GPUs across 512 compute nodes. In all, our study contributes an evaluation of a communication-free model as well as a scalability study of computing distributed Lagrangian flow maps at scale using in situ infrastructure on a modern supercomputer.

Sane, Sudhanshu↗

Functional Principal Component Analysis of Radio–Optical Reference Frame Tie

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wassersplines for Neural Vector Field-Controlled Animation

Much of computer-generated animation is created by manipulating meshes with rigs. While this approach works well for animating articulated objects like animals, it has limited flexibility for animating less structured free-form objects. Here we introduce Wassersplines, a novel trajectory inference method for animating unstructured densities based on recent advances in continuous normalizing flows and optimal transport. The key idea is to train a neurally-parameterized velocity field that represents the motion between keyframes. Trajectories are then computed by advecting keyframes through the velocity field. We solve an additional Wasserstein barycenter interpolation problem to guarantee strict adherence to keyframes. Our tool can stylize trajectories through a variety of PDE-based regularizers to create different visual effects. We demonstrate our tool on various keyframe interpolation problems to produce temporally-coherent animations without meshing or rigging.

97 MATHEMATICS AND COMPUTING↗

Vector magnetic field evolution, energy storage, and associated photospheric velocity shear within a flare-productive active region

Sheared photospheric velocity fields inferred from spot motions for April 5-7, 1980, are compared with both transverse magnetic field orientation changes and with the region's flare history. Rapid spot motions and high inferred velocity shear coincide with increased field alignment along the longitudinal neutral line and with increased flare activity, while a later decrease in velocity shear precedes a more relaxed magnetic configuration and decrease in flare activity. It is estimated that magnetic reconfiguration produced by the relative velocities of the spots could cause storage of about 10 to the 32nd erg/day, while flares occurring during this time expended no more than about 10 to the 31st erg/day.

Krall, K. R.↗

The photospheric vector magnetic field of a sunspot and its vertical gradient

The results of direct comparisons of photospheric and transition region line-of-sight field observations of sunspots using the SMM UV spectrometer and polarimeter are reported. The analysis accompanying the data is concentrated on demonstrating that the sunspot concentrated magnetic field extends into the transition region. An observation of a sunspot on Oct. 23, 1980 at the S 18 E 03 location is used as an example. Maximum field strengths ranged from 2030-2240 gauss for large and small umbrae viewed and inclination of the field to the line-of-sight was determined for the photosphere and transition region. The distribution of the magnetic field over the sunspot and variation of the line-of-sight gradient are discussed, as are the magnitudes and gradients of the photospheric field across the penumbral-photospheric boundaries.

Hagyard, M. J.↗

Ion acceleration at shocks in interplanetary space - A brief review of recent observations

An assessment is presented of recent observations providing insights into the collisionless shock ion acceleration process as it is observed near the earth's bow shock and near interplanetary shocks generated by solar activity, with emphasis on the way in which the ion acceleration process appears to be the same for both types of shocks. The observed differences in energetic particle distributions found near these shocks appear to be related to the seed populations available for acceleration, the shock's extent, radius of curvature, and the angle beween the local shock normal and the magnetic field vector, the time of field line connection to the shock, and possibly the shock Mach number.

Gosling, J. T.↗