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 1,513 records · Page 84

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗

Measurement and Modeling of Steep Ocean Wave Slopes

Our study emphasizes the importance of identifying and quantifying the distribution variance, skewness and kurtosis from optical and microwave scattering observations. Recent field measurements of the sea slope distribution for intermediate-to-long scale gravity waves will be presented. These data were collected using an airborne laser range system designed to estimate the surface slope vector at horizontal scales of 1-2 m. The observed slope distribution tail indicates that the occurrence of steep waves substantially exceeds a Gaussian prediction. This measured peakedness is present over the wide range of sea state and wind speed conditions encountered. Data are further evaluated within the context of Cox and Munk's well-known sea slope investigations. Based on a re-evaluation of the Cox and Munk's reported parameters, we find a consistent picture develops wherein data are shown to consistently indicate non-Gaussian statistics. One fundamental application of such a non-Gaussian slope observation is its place in modifying predicted wave breaking probability to help to better quantify gas transfer processes at the sea surface.

Chapron, B.↗

Helioseismic Observations of the Structure and Dynamics of a Rotating Sunspot Beneath the Solar Surface

Time-distance helioseismology is applied to study the subphotospheric structures and dynamics of an unusually fast-rotating sunspot observed by the Michelson Doppler Imager on bead SOH0 in 2000 August. The subsurface sound speed structures and velocity fields are obtained for the sunspot region at different depths from 0 to 12 Mm. By comparing the subsurface sound speed variations with the surface magnetic field, we find evidence for structural twists beneath the visible surface of this active region, which may indicate that magnetic twists often seen at the photosphere also exist beneath the photosphere. We also report on the observation of subsurface horizontal vortical flows that extend to a depth of 5 Mm around this rotating sunspot and present evidence that opposite vortical flows may exist below 9 Mm. It is suggested that the vortical flows around this active region may build up a significant amount of magnetic helicity and energy to power solar eruptions. Monte Carlo simulation has been performed to estimate the error propagation, and in addition the sunspot umbra is masked to test the reliability of our inversion results. On the basis of the three-dimensional velocity fields obtained from the time-distance helioseismology inversions, we estimate the subsurface kinetic helicity at different depths for the first time and conclude that it is comparable to the current helicity estimated from vector magnetograms.

Zhao, Junwei↗

Estimating Total Open Heliospheric Magnetic Flux

Over the solar-activity cycle, there are extended periods where significant discrepancies occur between the spacecraft-observed total (unsigned) open magnetic flux and that determined from coronal models. In this article, the total open heliospheric magnetic flux is computed using two different methods and then compared with results obtained from insitu interplanetary magnetic-field observations. The first method uses two different types of photospheric magnetic-field maps as input to the Wang–Sheeley–Arge (WSA) model: i) traditional Carrington or diachronic maps, and ii) Air Force Data Assimilative Photospheric Flux Transport model synchronic maps. The second method uses observationally derived helium and extreme-ultraviolet coronal-hole maps overlaid on the same magnetic-field maps in order to compute total open magnetic flux. The diachronic and synchronic maps are both constructed using magnetograms from the same source, namely the National Solar Observatory Kitt Peak Vacuum Telescope and Vector Spectromagnetograph. The results of this work show that the total open flux obtained from observationally derived coronal holes agrees remarkably well with that derived from WSA, especially near solar minimum. This suggests that, on average, coronal models capture well the observed large-scale coronal-hole structure over most of the solar cycle. Both methods show considerable deviations from total open flux deduced from spacecraft data, especially near solar maximum, pointing to something other than poorly determined coronal-hole area specification as the source of these discrepancies.

Wallace, S.↗

Linear stability of the H(+) - O(+) two-stream interaction in a magnetized plasma

The Bergman and Lotko (1986) analysis of the linear stability of the H(+) - O(+) two-stream interactions for parallel propagating wave modes is extended to an analysis of the stability of wave modes that propagate at an angle to the terrestrial magnetic field. The analaysis provides an insight into the manner in which the oblique modes linearly couple and become unstable; the conditions under which the particular couplings are important; the qualitative behavior of the frequency, wave vector, and growth rate of the most unstable mode; and the changes to these introduced by finite ion temperatures. Consideration is also given to the relative importance of obliquely propagating and parallel propagating waves for the condition when the velocity difference between two ion beams is below the upper limit.

Bergmann, R.↗

Subphotospheric current systems and flares

Subphotospheric current systems inferred from recent vector magnetograph observations imply the existence of electric currents penetrating the photosphere and thus flowing deep in the solar convection zone. These currents presumably originate in an internal dynamo that supplies the observed photospheric magnetic fields through the buoyant motions of the initially deeply-buried flux tubes. The coronal fields resulting from this process therefore must carry slowly-varying currents driven by emfs remote from the surface. These currents may then drive solar-flare energy release. This paper discusses the consequences of such a deep origin of the coronal parallel currents. Simple estimates for a large active region suggest a mean current-closure depth of not less than 10,000 km, with a subphotospheric inductance of not less than 100 H and a subphotospheric stored energy of not less than 10 to the 33rd ergs.

Hudson, Hugh S.↗

Stress and vibraton analyses of anisotropic shells of revolution

An efficient computational strategy is presented for reducing the cost of the stress and free vibration analyses of laminated anisotropic shells of revolution. The analytical formulation is based on a form of the Sanders-Budiansky shell theory including the effects of both the transverse shear deformation and the laminated anisotropic material response. The fundamental unknowns consist of the eight strain components, the eight stress resultants and the five generalized displacements of the shell. Each of the shell variables is expressed in terms of trigonometric functions (Fourier series) in the circumferential co-ordinate, and a three-field mixed finite element model is used for the discretization in the meridional direction. The shell response associated with a range of Fourier harmonics is approximated by a linear combination of a few global approximation vectors, which are generated at a particular value of the Fourier harmonic, within that range. The full equations of the finite element model are solved for only a single Fourier harmonic, and the response corresponding to the other Fourier harmonics is generated using a reduced system of equations with considerably fewer degrees of freedom.

Noor, Ahmed K.↗

Convection phenomena in low-gravity processing - The GTE GaAs space experiment

A numerical model for the proposed GTE GaAs space experiment is presented. The model includes the actual ampule design and the experimentally obtained furnace thermal profiles. Results are presented based on both the traditional steady-state models and on fully time dependent simulations. It is demonstrated that the commonly presumed axisymmetric flows in numerical studies are physically accurate only when the local gravitational acceleration vector is parallel to the ampules's axis. Other gravitational orientations result in nontrivial three-dimensional flows which must be resolved using full three-dimensional models. Implications of the resulting flow and thermal fields on the growing crystal at conditions available during space processing are also discussed.

Arnold, W. A.↗

Microwave absorption by a lossy dielectric sphere in a rectangular cavity

A new theory of absorption of microwave power by a lossy dielectric sphere in a resonant rectangular chamber is developed here. In this treatment, resonant modes of the cavity are resolved into traveling plane-wave fields, and Mie's theory of scattering and absorption of a plane electromagnetic wave by a sphere is utilized. The fields in the cavity containing the sphere are found by superposition of results for the individual scattering problem for each traveling plane wave. Absorbed power is then calculated by evaluating the time-averaged Poynting vector and integrating over the surface of the sphere. Applications of microwave power absorption are discussed, with particular emphasis on processing of materials. Numerical studies based on the newly derived formulas are presented. The results exhibit the manner in which microwave absorption is influenced by a number of different experimentally controllable parameters. The examples considered were chosen with a view toward designing experiments and apparatus for containerless processing of materials in the microgravity environment of space.

Jackson, H. W.↗

Data Simulation Using VLIDORT and PCRTM and its Use for Radiometric Calibration of Hyperspectral UV-NIR-SWIR Sensors

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer discrete ordinate scattering codes, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate the radiation fields for hyperspectral UV-NIR-SWIR sensors from 250nm to 2500nm with a very high speed. Since the accurate gas absorptions, molecules and aerosol/cloud particle scatterings, and the surface reflectance have been considered in these simulations, these model data can be used for satellite or airborne sensor design, characterization, or end-to-end system level uncertainty assessment. Comparison of the modelled data with measured results has been a general practice for operational satellite sensor calibrations, such as the Ozone Mapping and Profiler Suite (OMPS) onboard on SNPP and NOAA-20. This presentation will give more details of this simulation system and models. Some examples of the generation of simulations data for CLARREO PathFinder (CPF) and the Surface Biology and Geology (SBG), as well as some examples to use these data for CPF and OMPS radiometric calibration will be presented.

Xiaozhen Xiong↗

Characterizing Spatiotemporal Uncertainty in Interpolated Meteorological Data

Interpolated meteorological data invariably contain errors. These errors have structure in time and space, particularly autocorrelation, which can cause the effects of errors to compound when model outputs are aggregated temporally or spatially. One way to account for this uncertainty is with a probabilistic model from which samples can be drawn that are coherent with respect to underlying spatial and temporal covariance structure. This work describes a probabilistic method for spatial interpolation of point-wise meteorological time series. Observational data from weather stations are generally sparse in space and dense in time (but sometimes missing). The method works by projecting time series onto orthogonal basis vectors and spatially interpolating each resulting component independently. Under suitable assumptions, and data transformations to better satisfy those assumptions, Gaussian process regression provides a complete description of the joint predictive distribution over a Gaussian random field. Spatiotemporally coherent realizations are generated as the sum of conditional (spatial) simulations of each orthogonal (temporal) component. Data-derived and generic orthogonal bases are considered. In addition to spatial interpolation, imputation of missing observational data is examined. The method is applied using near-surface air temperature over the Western United States and validated by comparing theoretical versus actual coverage of predictive distributions and analyzing the degree to which spatial and temporal covariance structure is reproduced. Computational considerations, relating to conditional simulation of random fields, are also addressed.

Conor T Doherty↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

Quantifying the Complexity of Flaring Active Regions

While solar physicists have a better understanding of the importance magnetic fields play in the solar heating mechanism, it is still not possible to predict whether or when an active region will flare. In recent decades, qualitative studies of the changes in active region morphology have shown that there is generally an increase in the complexity of the spatial configuration of a solar active region leading up to a flare event. In this study, we quantify the spatial structure of the region using the differential Box-Counting Method (DBC) of fractal analysis. We analyze data from NASA/Marshall Space Flight Centr's vector magnetograph from two flaring active regions: AR 6089 from June 10, 1990, which produced one M1.7 flare, and AR 6659 from June 8, 9 and 10, 1991, this data set including one C5.7 and two M(6.4 and 3.2) flare. (AR 6659 produced several other flares). Several magnetic parameters are studied, including the transverse and longitudinal magnetic field components (Bt and B1), the total field (Bmag), and the magnetic shear, which describes the non-potentiality of the field. Results are presented for the time series of magnetograms in relation to the timing of flare events.

Stark, B.↗

A new concept of plasma motion and planetary magenetic field for Venus

It is shown that the magnetohydrodynamic conditions of the Venus ionosphere near the terminator favor convection of a magnetic field rather than diffusion. Consequently, any planetary magnetic field which Venus may possess will be strongly affected by the global antisunward flow of the ionosphere which has been revealed by the Pioneer-Venus retarding potential analyzer. The magnetic flux from an internal magnetic field will accumulate in the night hemisphere. Details of the structure and dynamics of such accumulations depend on particular details of the magnetic field source and the time-dependent plasma flow pattern, but a simple interpretation of observational data yields a magnetic dipole moment of 7 x 10 to the 20th cu cm directed along the planet spin vector.

Knudsen, W. C.↗

A vectorized solution for incompressible flow

An algorithm is developed to obtain solutions to the unsteady Reynolds-averaged incompressible Navier-Stokes equations in general curvilinear coordinates on a vector processor. The governing equations are in nonconservative form with the velocity and pressure as dependent variables. Two momentum equations and the Poisson equation for pressure form a set of three governing equations for three flow field unknowns: u, v, and p. The governing equations and boundary conditions are expressed in terms of boundary-conforming curvilinear coordinates, and a checkerboard SOR iteration is used to solve the governing equations. Several possible sequences for a checkerboard SOR iteration are investigated for finding the best overall convergence rate. The efficiency and capability of the present algorithm was assessed using the example of an 18 percent thick NACA 66(3)018 airfoil at zero degree angle of attack for chord Reynolds number range 1000-40,000.

Patel, N. R.↗

A large-amplitude rotational wave in the Venusian ionosheath

Results are reported of a preliminary analysis of a large-amplitude rotational wave in the Venusian ionosheath. The results are based on an analysis of the Pioneer Venus Orbiter plasma and magnetic field observations from six orbits in the first (1979) tail season. This wave appears to be the standing super-Alfvenic wing. It is located within the shocked plasma flow outside the boundary of the tail. The rotation of the magnetic field by about 90 deg across the wave occurs through several successive cycles making the wave similar to the group velocity wing consisting of phase velocity waves. The transition through the wave is accompanied by the vector change of the plasma velocity with the magnitude of the plasma velocity jump comparable to the vector jump of the Alfven velocity. The observed super-Alfvenic wing appears to originate upstream and closer to the planet, possibly near the upper boundary of the magnetic barrier on the dayside.

Fedorov, A. O.↗

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗