Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector fields”

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 127 records · Page 7

The Topology of Symmetric Tensor Fields

Combinatorial topology, also known as "rubber sheet geometry", has extensive applications in geometry and analysis, many of which result from connections with the theory of differential equations. A link between topology and differential equations is vector fields. Recent developments in scientific visualization have shown that vector fields also play an important role in the analysis of second-order tensor fields. A second-order tensor field can be transformed into its eigensystem, namely, eigenvalues and their associated eigenvectors without loss of information content. Eigenvectors behave in a similar fashion to ordinary vectors with even simpler topological structures due to their sign indeterminacy. Incorporating information about eigenvectors and eigenvalues in a display technique known as hyperstreamlines reveals the structure of a tensor field. The simplify and often complex tensor field and to capture its important features, the tensor is decomposed into an isotopic tensor and a deviator. A tensor field and its deviator share the same set of eigenvectors, and therefore they have a similar topological structure. A a deviator determines the properties of a tensor field, while the isotopic part provides a uniform bias. Degenerate points are basic constituents of tensor fields. In 2-D tensor fields, there are only two types of degenerate points; while in 3-D, the degenerate points can be characterized in a Q'-R' plane. Compressible and incompressible flows share similar topological feature due to the similarity of their deviators. In the case of the deformation tensor, the singularities of its deviator represent the area of vortex core in the field. In turbulent flows, the similarities and differences of the topology of the deformation and the Reynolds stress tensors reveal that the basic addie-viscosity assuptions have their validity in turbulence modeling under certain conditions.

Levin, Yingmei↗

Towards an Optimal Noise Versus Resolution Trade-Off in Wind Scatterometry

A scatterometer is a radar that measures the normalized radar cross section sigma(sup 0) of the Earth's surface. Over the ocean this signal is related to the wind via the geophysical model function (GMF). The objective of wind scatterometry is to estimate the wind vector field from sigma(sup 0) measurements; however, there are many subtleties that complicate this problem-making it difficult to obtain a unique wind field estimate. Conventionally, wind estimation is split into two stages: a wind retrieval stage in which several ambiguous solutions are obtained, and an ambiguity removal stage in which ambiguities are chosen to produce an appropriate wind vector field estimate. The most common approach to wind field estimation is to grid the scatterometer swath into wind vector cells and estimate wind vector ambiguities independently for each cell. Then, field wise structure is imposed on the solution by an ambiguity selection routine. Although this approach is simple and practical, it neglects field wise structure in the retrieval step and does not account for the spatial correlation imposed by the sampling. This makes it difficult to develop a theoretically appropriate noise versus resolution trade-off using pointwise retrieval. Fieldwise structure may be imposed in the retrieval step using a model-based approach. However, this approach is generally only practical if a low order wind field model is applied, which may discard more information than is desired. Furthermore, model-based approaches do not account for the structure imposed by the sampling. A more general fieldwise approach is to estimate all the wind vectors for all the WVCs simultaneously from all the measurements. This approach can account for structure of the wind field as well as structure imposed by the sampling in the wind retrieval step. Williams and Long in 2010 developed a fieldwise retrieval method based on maximum a posteriori estimation (MAP). This MAP approach can be extended to perform a noise versus resolution trade-off, and deal with ambiguity selection. This paper extends the fieldwise MAP estimation approach and investigates both the noise versus resolution trade-off as well as ambiguity removal in the fieldwise wind retrieval step. The method is then applied to the Sea Winds scatterometer and the results are analyzed. This paper extends the fieldwise MAP estimation approach and investigates both the noise versus resolution trade-off as well as ambiguity removal in the fieldwise wind retrieval step. The method is then applied to the Sea Winds scatterometer and the results are analyzed.

soil moisture↗

The effect of line damping, magneto-optics and parasitic light on the derivation of sunspot vector magnetic fields

The least square fitting of Stokes observations of sunspots using a Milne-Eddington-Unno model appears to lead, in many circumstances, to various inconsistencies such as anomalously large doppler widths and, hence, small magnetic fields which are significantly below those inferred solely from the Zeeman splitting in the intensity profile. It is found that the introduction of additional physics into the model such as the inclusion of damping wings and magneto-optic birefrigence significantly improves the fit to Stokes parameters. Model fits excluding the intensity profile, i.e., of both magnitude as well as spectral shape of the polarization parameters alone, suggest that parasitic light in the intensity profile may also be a source of inconsistencies. The consequences of the physical changes on the vector properties of the field derived from the Fe I lambda 6173 line for the 17 November 1975 spot as well as on the thermodynamic state are discussed. A Doppler width delta lambda (D) - 25mA is bound to be consistent with a low spot temperature and microturbulence, and a damping constant of a = 0.2.

Skumanich, A.↗

Direct discretization of planar div-curl problems

A control volume method is proposed for planar div-curl systems. The method is independent of potential and least squares formulations, and works directly with the div-curl system. The novelty of the technique lies in its use of a single local vector field component and two control volumes rather than the other way around. A discrete vector field theory comes quite naturally from this idea and is developed. Error estimates are proved for the method, and other ramifications investigated.

Nicolaides, R. A.↗

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

Rotational polarities of sudden impulses in the magnetotail lobe

A sudden impulse (SI) is a sudden change in the magnetic field strength which is caused by a change in the solar wind pressure and is observed throughout the magnetosphere. In this report we have examined the rotations of the magnetic field vectors at times of SIs in the magnetotail lobe, by using IMP 6, 7, and 8 magnetometer data. The following properties have been found: (1) at the time of SI the arrowhead of the magnetic vector tends to rotate in one plane; (2) the plane of rotation tends to include the unperturbed magnetic field vector; (3) the plane of rotation tends to be aligned with the radial direction from the magnetotail axis; and (4) the magnetic vectors have a particular rotational polarity: when the plane of rotation is viewed so that the Sun is to the right of the viewed plane and the magnetotail axis is to the bottom, the arrowhead of the vector tends to rotate counterclockwise in this plane. These magnetic vector properties are consistent with those expected when part of an increase in solar wind lateral pressure squeezes the magnetotail axisymmetrically while moving tailward.

Kawano, H.↗

Split-Block Waveguide Polarization Twist for 220 to 325 GHz

A split-block waveguide circuit that rotates polarization by 90 has been designed with WR-3 input and output waveguides, which are rectangular waveguides used for a nominal frequency range of 220 to 325 GHz. Heretofore, twisted rectangular waveguides equipped with flanges at the input and output have been the standard means of rotating the polarizations of guided microwave signals. However, the fabrication and assembly of such components become difficult at high frequency due to decreasing wavelength, such that twisted rectangular waveguides become impractical at frequencies above a few hundred gigahertz. Conventional twisted rectangular waveguides are also not amenable to integration into highly miniaturized subassemblies of advanced millimeter- and submillimeter-wave detector arrays now undergoing development. In contrast, the present polarization- rotating waveguide can readily be incorporated into complex integrated waveguide circuits such as miniaturized detector arrays fabricated by either conventional end milling of metal blocks or by deep reactive ion etching of silicon blocks. Moreover, the present split-block design can be scaled up in frequency to at least 5 THz. The main step in fabricating a splitblock polarization-rotating waveguide of the present design is to cut channels having special asymmetrically shaped steps into mating upper and lower blocks (see Figure 1). The dimensions of the steps are chosen to be consistent with the WR-3 waveguide cross section, which is 0.864 by 0.432 mm. The channels are characterized by varying widths with constant depths of 0.432, 0.324, and 0.216 mm and by relatively large corner radii to facilitate fabrication. The steps effect both a geometric transition and the corresponding impedance-matched electromagnetic-polarization transition between (1) a WR-3 rectangular waveguide oriented with the electric field vector normal to the block mating surfaces and (2) a corresponding WR-3 waveguide oriented with its electric field vector parallel to the mating surfaces of the blocks. A prototype has been built and tested. Figure 2 presents test results indicative of good performance over nearly the entire WR-3 waveguide frequency band.

Ward, John↗

An Efficient High-Order Solver for Diffusion Equations with Strong Anisotropy on Non-Anisotropy-Aligned Meshes

This paper concerns numerical solution of the diffusion equation with strong anisotropy on meshes not aligned with the anisotropic vector field. In order to resolve the numerical pollution for simulations on a non-anisotropy-aligned mesh and reduce the associated high computational cost we propose an effective preconditioner, extending our previous work. Similar to the anisotropy-aligned mesh case, we apply the auxiliary space preconditioning framework to design a preconditioner where a continuous finite element space is used as the auxiliary space for the discontinuous finite element space. The key component is an effective line smoother that can mitigate the high-frequency errors perpendicular to the magnetic field. We design a graph-based approach to find such a line smoother that is approximately perpendicular to the vector fields when the mesh does not align with the anisotropy. Finally, numerical experiments for several benchmark problems are presented, demonstrating the effectiveness and robustness of the proposed preconditioner when applied to Krylov iterative methods.

97 MATHEMATICS AND COMPUTING↗

Evaluation of two January - June 1992 ERS-1 AMI wind vector data sets

A Jan. - Jun. 1992 10 m height wind vector field was computed with the CMOD2 model function. Freilich and Dunbar used ESA measurements of AMI (Active Microwave Instrument) sigma-naught to derive a refined model function and to produce an alternate Jan. - Jun. 1992 10 m height wind vector field, named FD data. The ESA and FD model functions and algorithms were substantially different. The FD wind vectors were in better agreement with the climatological mean wind field and with moored buoy wind measurements. The average correlation coefficients between monthly mean and daily averaged FD and moored buoy matchups of the Cartesian wind components were 0.84 and 0.74, respectively. The average root mean square differences between the monthly mean and daily averaged FD and moored buoy matchups of the Cartesian wind components were 1.25 and 3.0 m/s, respectively.

Halpern, D.↗

Magnetic Test Performance Capabilities at the Goddard Space Flight Center as Applied to the Global Geospace Science Initiative

Goddard Space Flight Center's (GSFC) Spacecraft Magnetic Test Facility (SMTF) is a historic test facility that has set the standard for all subsequent magnetic test facilities. The SMTF was constructed in the early 1960's for the purpose of simulating geomagnetic and interplanetary magnetic fields. Additionally, the facility provides the capability for measuring spacecraft generated magnetic fields as well as calibrating magnetic attitude control systems and science magnetometers. The SMTF was designed for large, spacecraft level tests and is currently the second largest spherical coil system in the world. The SMTF is a three-axis Braunbek system composed of four coils on each of three orthogonal axes. The largest coils are 12.7 meters (41.6 feet) in diameter. The three-axis Braunbek configuration provides a highly uniform cancellation of the geomagnetic field over the central 1.8 meter (6 foot) diameter primary test volume. Cancellation of the local geomagnetic field is to within +/-0.2 nanotesla with a uniformity of up to 0.001% within the 1.8 meter (6 foot) diameter primary test volume. Artificial magnetic field vectors from 0-60,000 nanotesla can be generated along any axis with a 0.1 nanotesla resolution. Oscillating or rotating field vectors can also be produced about any axis with a frequency of up to 100 radians/second. Since becoming fully operational in July of 1967, the SMTF has been the site of numerous spacecraft magnetics tests. Spacecraft tested at the SMTF include: the Solar Maximum Mission (SMM), Magsat, LANDSAT-D, the Fast Aurora] Snapshot (FAST) Explorer and the Sub-millimeter-Wave-Astronomy Satellite (SWAS) among others. This paper describes the methodology and sequencing used for the Global Geospace Science (GGS) initiative magnetic testing program in the Goddard Space Flight Center's SMTF. The GGS initiative provides an exemplary model of a strict and comprehensive magnetic control program.

Mitchell, Darryl R.↗

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗

Lie theoretic aspects of the Riccati equation

Various features of the application of Lie theory to matrix Riccati equations, of basic importance in control and system theories, are discussed. Particular consideration is given to centralizer foliation, the Cartan decomposition, matrix Riccati equations as Lie systems on Grassmanians, local analysis near a zero point of a vector field, linearization in homogeneous space, the tangent bundle in terms of partitioned matrices, and stability properties of fixed points of Riccati vector fields.

Hermann, R.↗

Finding a contaminant source in a volume of flowing fluid

A method includes: mapping a boundary of a volume of flowing fluid; partitioning the volume by a computational mesh; finding a contaminant location at a first sensor that is disposed within the volume; obtaining a measured velocity of the flowing fluid within the volume; generating a reversed velocity vector field within the mesh, in response to the measured velocity; time stepping the contaminant location from the first sensor along the reversed velocity vector field until the contaminant location intersects the boundary of the volume; and finding a contaminant source at the intersection of the time stepped contaminant location with the boundary of the volume.

Lopez-Marrero, Vanessa↗

Effective field theory of Stückelberg vector bosons

We explore the effective field theory of a vector field $X^μ$ that has a Stückelberg mass. The absence of a gauge symmetry for $X^μ$ implies Lorentz-invariant operators are constructed directly from Xμ. Beyond the kinetic and mass terms, allowed interactions at the renormalizable level include $X_μX^μH^†H, (X_μX^μ)^2$, and $X_μj^μ$, where $j^μ$ is a global current of the SM or of a hidden sector. We show that all of these interactions lead to scattering amplitudes that grow with powers of $\sqrt{s}/m_X$, except for the case of $X_μj^μ$ where $j^μ$ is a nonanomalous global current. The latter is well-known when $\textit{X}$ is identified as a dark photon coupled to the electromagnetic current, often written equivalently as kinetic mixing between $\textit{X}$ and the photon. The power counting for the energy growth of the scattering amplitudes is facilitated by isolating the longitudinal enhancement. We examine in detail the interaction with an anomalous global vector current $X_μj^μ_{\text{anom}}$, carefully isolating the finite contribution to the fermion triangle diagram. We calculate the longitudinally-enhanced observables $Z → X_γ$ (when $m_X < m_Z), f\bar{f} → X_γ$, and $Z_γ → Z_γ$ when $\textit{X}$ couples to the baryon number current. Introducing a "fake" gauge-invariance by writing $X^μ = A^μ – ∂^μπ/m_X$, the would-be gauge anomaly associated with $A_μj^μ_{\text{anom}}$ is canceled by $j^μ_{\text{anom}}∂_μπ/m_X$; this is the four-dimensional Green-Schwarz anomaly-cancellation mechanism at work. Our analysis demonstrates there is a much larger set of possible interactions that an EFT with a Stückelberg vector field can have, revealing scattering amplitudes that grow with energy. The growth of these amplitudes can be tamed by a dark Higgs sector, but this requires dark Higgs boson interactions (and reintroduces fine-tuning in the dark Higgs sector) that can be separated from X interactions only in the limit $g \ll 1$ .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bifurcations of relative equilibria

The characteristics of equivariant dynamical systems near relative equilibria (RE: group orbits which are invariant in the flow of an equivariant vector field) are investigated analytically, with a focus on the dynamics and bifurcation (B) behavior. The principles of Lie-group theory are reviewed; the decomposition of the vector field is explained; and particular attention is given to the Bs of RE occurring when an eigenvalue passes through zero, Hopf Bs of RE, the classification of generic secondary steady-state and Hopf Bs with symmetry group O(2), Bs of the zero solution of the Kuramoto-Shivashinsky equation, and possible generic steady-state Bs in the two-dimensional Benard problem. In the latter case, it is shown that the primary generic Bs are to two types of equilibria (hexagons and rolls), while the secondary Bs result in trajectories which are either equilibria or rotating waves.

Krupa, Martin↗

Studies of high latitude current systems using MAGSAT vector data

The magnetic disturbance fields caused by global external current systems are considered with particular emphasis on improving the understanding of the physical processes which control high latitude current systems. Following processing the MAGSAT data were routinely plotted in the Universal Time (UT) format as well as in a polar plot format. The H'D'U' coordinate system, was adopted as the standard for representing the MAGSAT residual magnetic field vectors. A data file was generated and the TPOLAR computer code was developed to determine from the orbital elements, the time, latitude, and MLT of the extremum latitude of each transpolar segment of orbit. The precision of the vector data set from MAGSAT prompted an extended exploratory phase for data analysis procedures, modeling techniques and phenomenology.

Burrows, J. R.↗

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.↗