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 1,621 records · Page 90

Forced MHD turbulence in a uniform external magnetic field

Two-dimensional dissipative MHD turbulence is randomly driven at small spatial scales and is studied by numerical simulation in the presence of a strong uniform external magnetic field. A behavior is observed which is apparently distinct from the inverse cascade which prevails in the absence of an external magnetic field. The magnetic spectrum becomes dominated by the three longest wavelength Alfven waves in the system allowed by the boundary conditions: those which, in a box size of edge 2 pi, have wave numbers (kx' ky) = (1, 1), and (1, -1), where the external magnetic field is in the x direction. At any given instant, one of these three modes dominates the vector potential spectrum, but they do not constitute a resonantly coupled triad. Rather, they are apparently coupled by the smaller-scale turbulence.

Hossain, M.↗

Forced magnetohydrodynamic turbulence in a uniform external magnetic field

Two-dimensional dissipative MHD turbulence is randomly driven at small spatial scales and is studied by numerical simulation in the presence of a strong uniform external magnetic field. A behavior is observed which is apparently distinct from the inverse cascade which prevails in the absence of an external magnetic field. The magnetic spectrum becomes dominated by the three longest wavelength Alfven waves in the system allowed by the boundary conditions: those which, in a box size of edge 2 pi, have wave numbers (kx, ky) = (1, 1), and (1, -1), where the external magnetic field is in the x direction. At any given instant, one of these three modes dominates the vector potential spectrum, but they do not constitute a resonantly coupled triad. Rather, they are apparently coupled by the smaller-scale turbulence.

Hossain, M.↗

Automated remote sensing of sea ice using synthetic aperture radar

Two techniques for automated sea ice tracking: image pyramids (hierarchical correlation) and feature tracking were applied to a pair of SEASAT SAR sea ice images. The results compare well with each other and manually tracked estimates of the ice velocity field. Refinement of the image pyramid technique should include removal of a small number of erroneous velocity vectors using constraints dependent on correlation coefficient, near neighbor velocities, and velocity gradients. Estimates of sea ice velocity are successfully obtained over limited areas using the boundary segment tracking technique. Refinement of the feature tracking technique should include use of pressure ridge features, constrained search, and relaxation techniques.

Vesecky, J. F.↗

Acoustic forcing of a liquid drop

The development of systems such as acoustic levitation chambers will allow for the positioning and manipulation of material samples (drops) in a microgravity environment. This provides the capability for fundamental studies in droplet dynamics as well as containerless processing work. Such systems use acoustic radiation pressure forces to position or to further manipulate (e.g., oscillate) the sample. The primary objective was to determine the effect of a viscous acoustic field/tangential radiation pressure forcing on drop oscillations. To this end, the viscous acoustic field is determined. Modified (forced) hydrodynamic field equations which result from a consistent perturbation expansion scheme are solved. This is done in the separate cases of an unmodulated and a modulated acoustic field. The effect of the tangential radiation stress on the hydrodynamic field (drop oscillations) is found to manifest as a correction to the velocity field in a sublayer region near the drop/host interface. Moreover, the forcing due to the radiation pressure vector at the interface is modified by inclusion of tangential stresses.

Lyell, M. J.↗

Development of a stable electro-optical modulator

The activity of the sun can be studied through the measurement of its magnetic field. The measurement of the polarization state of the sunlight in a narrow spectral bandwidth can be related to the magnetic fields strengths and directions over the active sunspots. Measuring the polarization state can be done using the vector magnetograph. The essential component of the vector magnetograph is the polarization modulators. Electro-optical modulators are used as dynamic retarders. KDP (potassium dideuterium phosphate) crystals provide variable retardation which can be controlled by the applied voltage. KDP's have a high contrast ratio but suffer some problems; they require a high voltage and have a limited modulation lifetime. The KDP modulator forms the heart of the polarimeter assembly which is used for carrying out accurate measurements of linear polarization of the sun. It is therefore important to study the optical and electrical behavior of the KDP crystals so that the ultimate goal of making stable electro-optical modulators in the laboratory is achieved. A number of activities were performed in coating and testing the KDP crystals. A summary of these activities is given.

Abushagur, Mustafa A. G.↗

SSMI Wind Speed Climatology of the Time of Monsoon Wind Offset in the Western Arabian Sea

Forecasting the time of onset of monsoon wind in the western Arabian Sea, which is believed to precede the onset of rainfall along the west coast of India, is an important unsolved problem. Prior to measurements of the surface wind field by satellite, there was an absence of suitable surface wind observations. NASA scatterometer (NSCAT) surface wind vectors revealed that the time of the 1997 onset of 12 m/s southwest monsoon wind speeds in the western Arabian Sea preceded the onset of monsoon rainfall in Goa, India, by 3 - 4 days. Wind speed and direction data were necessary to establish a dynamical mechanism between times of onset of 12 m/s wind speed off Somalia and rainfall in Goa. Except for NSCAT, no satellite scatterometer wind product recorded adequately sampled 2-day 1deg x 1deg averaged wind vectors, which are the required space and time scales, to examine the wind-rain relationship in other years. However, the greater-than-95% steadiness of summer monsoon winds allows an opportunity to use satellite measurements of surface wind speed. The Special Sensor Microwave Imager (SSMI) recorded surface wind speed with adequate sampling to produce a 1-day, 1deg x 1deg data product during 1988 - 1998. SSMI data had been uniformly processed throughout the period. Times of onset of 12 m/s wind speed off Somalia determined with the SSMI data set were 21 May 1988, 24 May 1989, 17 May 1990, 28 May 1991, 8 June 1992, 28 May 1993, 30 May 1994, 7 June 1995, 29 May 1996, 12 June 1997, and 15 May 1998. Uncertainty of the 1992 and 1996 times of onset were increased because of the absence of SSMI data on 6 and 7 June 1992 and on 30 May 1996. Correlations of timing of monsoon wind onset with El Nino will be described. Variability of the time difference between times of onset of 12 m/s wind speed and Goa rainfall will be discussed. At the time of submission of the abstract, the Goa rainfall data have not arrived from the India Meteorological Department.

Halpern, David↗

Trajectory Design Employing Convex Optimization for Landing on Irregularly Shaped Asteroids

Mission proposals that land spacecraft on asteroids are becoming increasingly popular. However, in order to have a successful mission the spacecraft must reliably and softly land at the intended landing site with pinpoint precision. The problem under investigation is how to design a propellant optimal powered descent trajectory that can be quickly computed onboard the spacecraft, without interaction from the ground control. The propellant optimal control problem in this work is to determine the optimal finite thrust vector to land the spacecraft at a specified location, in the presence of a highly nonlinear gravity field, subject to various mission and operational constraints. The proposed solution uses convex optimization, a gravity model with higher fidelity than Newtonian, and an iterative solution process for a fixed final time problem. In addition, a second optimization method is wrapped around the convex optimization problem to determine the optimal flight time that yields the lowest propellant usage over all flight times. Gravity models designed for irregularly shaped asteroids are investigated. Success of the algorithm is demonstrated by designing powered descent trajectories for the elongated binary asteroid Castalia.

Pinson, Robin M.↗

Determination of the viscous acoustic field for liquid drop positioning/forcing in an acoustic levitation chamber in microgravity

The development of acoustic levitation systems has provided a technology with which to undertake droplet studies as well as do containerless processing experiments in a microgravity environment. Acoustic levitation chambers utilize radiation pressure forces to position/manipulate the drop. Oscillations can be induced via frequency modulation of the acoustic wave, with the modulated acoustic radiation vector acting as the driving force. To account for tangential as well as radial forcing, it is necessary that the viscous effects be included in the acoustic field. The method of composite expansions is employed in the determination of the acoustic field with viscous effects.

Lyell, Margaret J.↗

Application of Powder Diffraction Methods to the Analysis of the Atomic Structure of Nanocrystals: The Concept of the Apparent Lattice Parameter (ALP)

The applicability of standard methods of elaboration of powder diffraction data for determination of the structure of nano-size crystallites is analysed. Based on our theoretical calculations of powder diffraction data we show, that the assumption of the infinite crystal lattice for nanocrystals smaller than 20 nm in size is not justified. Application of conventional tools developed for elaboration of powder diffraction data, like the Rietveld method, may lead to erroneous interpretation of the experimental results. An alternate evaluation of diffraction data of nanoparticles, based on the so-called 'apparent lattice parameter' (alp) is introduced. We assume a model of nanocrystal having a grain core with well-defined crystal structure, surrounded by a surface shell with the atomic structure similar to that of the core but being under a strain (compressive or tensile). The two structural components, the core and the shell, form essentially a composite crystal with interfering, inseparable diffraction properties. Because the structure of such a nanocrystal is not uniform, it defies the basic definitions of an unambiguous crystallographic phase. Consequently, a set of lattice parameters used for characterization of simple crystal phases is insufficient for a proper description of the complex structure of nanocrystals. We developed a method of evaluation of powder diffraction data of nanocrystals, which refers to a core-shell model and is based on the 'apparent lattice parameter' methodology. For a given diffraction pattem, the alp values are calculated for every individual Bragg reflection. For nanocrystals the alp values depend on the diffraction vector Q. By modeling different a0tomic structures of nanocrystals and calculating theoretically corresponding diffraction patterns using the Debye functions we showed, that alp-Q plots show characteristic shapes which can be used for evaluation of the atomic structure of the core-shell system. We show, that using a simple model of a nanocrystal with spherical shape and centro-symmetric strain at the surface shell we obtain theoretical alp-Q values which match very well the alp-Q plots determined experimentally for Sic, GaN, and diamond nanopowders. The theoretical models are defined by the lattice parameter of the grain core, thickness of the surface shell, and the magnitude and distribution of the strain field in the surface shell. According to our calculations, the part of the diffraction pattern measured at relatively low diffraction vectors Q (below 10/angstrom) provides information on the surface strain, whle determination of the lattice parameters in the grain core requires measurements at large Q-values (above 15 - 20/angstrom).

Palosz, B.↗

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ~0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

3C 84 (NGC 1275)↗

Toward the large-eddy simulations of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model (in terms of Favre-filtered fields) is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The three dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 to the third power grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. Future applications of these compressible subgrid-scale models to the large-eddy simulation of supersonic aerodynamic flows are discussed briefly.

Erlebacher, G.↗

The wavenumber spectra of scatterometer-derived winds

The accuracy of scatterometer-derived winds using a frequency-domain analysis and simulation is studied. The wavenumber spectra of the Seasat-A-scatterometer (SASS)-derived wind fields have been observed to be accentuated relative to the input wind field. The results of extensive simulations designed to test this observation are reported. Actual SASS measurements of the normalized radar backscatter (NRB) over an orbit (rev) are used as a template to generate simulated NRB measurements. Pointwise estimation of winds from the simulated NRB measurements is accomplished with the ambiguity closest to the true wind selected as the unique wind vector estimate. For comparison, winds are also estimated using a new model-based approach. After wind retrieval, the spectra of the estimated wind fields are computed and compared to the input wind field. The high wavenumber portion of the spectra of the pointwise estimated winds was higher than the spectra of the true winds by an amount which depends on the wind speed variance.

Long, D. G.↗

A Hybrid Symmetric FEM/MOM Formulation Applied to Scattering by Inhomogeneous Bodies of Revolution

A new symmetric formulation of the Hybrid Finite Element Method (HFEM) is described which combines elements of the Electric Field Integral Equation (EFIE) and the Magnetic Field Integral Equation (MFIE) for the exterior region along with the finite element solution for the interior problem. The formulation is applied to scattering by inhomogeneous bodies of revolution. To avoid spurious modes in the interior region a combination of vector and nodal based finite elements are used. Integral equations in the exterior region are used to enforce the Sommerfield radiation condition by matching both the tangential electric and magnetic fields between interior and exterior regions.

Lee, J-F.↗

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗

X-Ray Polarization Measurement of the Gold Standard of Radio-Quiet Active Galactic Nuclei: NGC 1068

Context.NGC 1068 is the most observed radio-quiet active galactic nucleus (AGN) in polarimetry, yet its high-energy polarization has never been probed before due to a lack of dedicated polarimeters.Aims.Using the first X-ray polarimeter sensitive enough to measure the polarization of AGNs, we want to probe the orientation and geometric arrangement of (sub)parsec-scale matter around the X-ray source. Methods. We used the Imaging X-ray Polarimetry Explorer (IXPE) satellite to measure, for the first time, the 2–8 keV polarization of NGC 1068.We pointed IXPE at the target for a net exposure time of 1.15 Ms, in addition to using two Chandra snapshots of∼10 ks each in order to account for the potential impact of several ultraluminous X-ray sources (ULXs) within IXPE’s field of view. Results. We measured a 2–8 keV polarization degree of 12.4%±3.6% and an electric vector polarization angle of 101◦±8◦at a 68% confidence level. If we exclude the spectral region containing bright Fe K lines and other soft X-ray lines where depolarization occurs, the polarization fraction rises to 21.3%±6.7% in the 3.5–6.0 keV band, with a similar polarization angle. The observed polarization angle is found to be perpendicular to the parsec-scale radio jet. Using a combined Chandra and IXPE analysis plus multiwavelength constraints, we estimated that the circum nuclear“torus” may sustain a half-opening angle of 50–55◦(from the vertical axis of the system). Conclusions. Thanks to IXPE, we have measured the X-ray polarization of NGC 1068 and found comparable results, both in terms of the polarization angle orientation with respect to the radio jet and the torus half-opening angle, to the X-ray polarimetric measurement achieved for the other archetypal Compton-thick AGN: the Circinus galaxy. Probing the geometric arrangement of parsec-scale matter in extragalactic objects is now feasible thanks to X-ray polarimetry

active galaxies↗

A Measure from Line-of-Sight Magnetograms for Prediction of Coronal Mass Ejections

From a sample of 17 vector magnetograms of 12 bipolar active regions we have recently found (1) that a measure of the overall nonpotentiality (the overall twist and shear in the magnetic field) of an active region is given by the strong shear length Lss, the length of the portion of the main neutral line on which the observed transverse fields is strong (greater than 150 Guass (G)) and strongly sheared (shear angle greater than 45 degrees), and (2) that L(sub ss) is well correlated with the coronal mass ejection (CME) productivity of the active regions during the plus or minus 2-day time window centered on the day of the magnetogram. In the present paper, from the same sample of 17 vector magnetograms, we show that there is a viable proxy for L(sub ss) that can be measured from a line-of-sight magnetogram. This proxy is the strong gradient length L(sub SG), the length of the portion of the main neutral line on which the potential transverse field is strong (greater than 150 G), and the gradient of the line-of-sight field is sufficiently steep (greater than approximately 50 G/Mm). In our sample of active regions, L(sub SG) is statistically significantly correlated with Lss (correlation confidence level greater than 95%), and L(sub SG) is as strongly correlated with active region CME productivity as is L(sub SS) (correlation confidence level approximately 99.7%). Because L(sub SG) can be measured from line-of-sight magnetograms obtained from conventional magnetographs, such as the magnetograph mode of the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory, it is a dependable substitute for L(sub SS) for use in operational CME forecasting. In addition, via measurement of L(sub SG), the years-long, nearly continuous sequence of 1.5-hour cadence full disk line-of-sight magnetograms from MDI can be used to track the growth and decay of the large-scale nonpotentiality in active regions and to examine the role of this evolution in active region CME productivity.

Falconer, D. A.↗

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗