SEARCH · Search NASA

Results for “fields simulating electromagnetic 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

341 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Optimal Transport Based Control of Granular Imaging System in Space

"Orbiting Rainbows” is a paradigm for creating a space-based observatory (telescope) from confined and aligned granular media. The overall objective is to construct a very large and lightweight aperture in space using a cloud of micron-sized particles held in position and aligned by the application of external electromagnetic fields, thus alleviating the extreme expense of deploying and stabilizing large monolithic telescopes by reducing the mass of the primary aperture. Adaptive optics techniques for wavefront stabilization and computational imaging techniques are then used to remove the optical noise from the scattering reflections of the granular medium, and finally producing clean images within the diffraction limit. In this paper, we address one of the key challenges of this new paradigm, namely controlling the position and orientation of the granules or particles so that the cloud of particles achieves the desired shape of an ultra-lightweight telescope and points into the right orientation. A further complication arises due to the fact that the particles do not have any onboard actuators. A number of external electric-field-based actuators are used to control all the particles simultaneously. Moreover, the number of actuators are orders-of-magnitude smaller than the number of particles, making this system largely under-actuated. Since the cloud of particles in any plausible configuration can be modeled as a probability distribution over the workspace, we formulate this control challenge as an optimal transport problem and then show that our optimal transport based control law drives the cloud of particles from any initial distribution to the desired position distribution (optical surface) and angular distribution (polarization state). Simulation results demonstrate the effectiveness of our approach.

Bandyopadhyay, Saptarshi↗

NASA NTRSDate not supplied

Suborbital Lunar Gravity Experiment of an Electrodynamic Regolith Conveyor

NASA’s Swamp Works Electrostatics and Surface Physics Laboratory (ESPL) is planning to test a 4-phase electrodynamic regolith conveyor in simulated lunar gravity on a suborbital flight in 2023. A 4-phase electrodynamic regolith conveyor could convey regolith without the risk of rotating or vibratory actuation, which could jam or require regular maintenance due to the abrasive nature of lunar regolith. Another advantage of electrodynamic conveying is the reduction of conveying power, which is important considering the limited scale of available early-stage lunar power systems. The current state of the art(SOA) for lunar regolith conveying is based on recent NASA system studies for oxygen production plants. These plant designs require conveying rates around 100 kg/hr, to produce 10 mT/yr of oxygen from the regolith. To accomplish this, conventional augers and vibratory spiral conveyors have been identified as the SOA or the leading candidates due to their extensive use in the terrestrial material handling industry. At NASA KSC, the use of dynamic electric fields, generated by alternating high voltage on electrodes, has been developed as a dust mitigation solution known as the Electrodynamic Dust Shield (EDS). The EDS is being developed for lenses, solar panels, radiators, fabric and seals. ESPL researchers have shown the ability to move thin layers (~ a few mm) of dust with mW of power. In Academia, researchers have shown the ability to electrodynamically convey regolith up to ~1 kg/hr with a4-phase EDS and electromagnetically convey particles that would suggest a lunar regolith conveying rate of ~25 kg/hr. This paper will describe the design and expected operation of a 4-phase electrodynamic regolith conveyor in lunar gravity. This paper will also describe how an electrodynamic conveyor could be scaled for transporting ISRU relevant flow rates.

Aaron D S Olson↗

NASA NTRS1996

Evaluation of Two Computational Techniques of Calculating Multipath Using Global Positioning System Carrier Phase Measurements

Two computational techniques are used to calculate differential phase errors on Global Positioning System (GPS) carrier war phase measurements due to certain multipath-producing objects. The two computational techniques are a rigorous computati electromagnetics technique called Geometric Theory of Diffraction (GTD) and the other is a simple ray tracing method. The GTD technique has been used successfully to predict microwave propagation characteristics by taking into account the dominant multipath components due to reflections and diffractions from scattering structures. The ray tracing technique only solves for reflected signals. The results from the two techniques are compared to GPS differential carrier phase ns taken on the ground using a GPS receiver in the presence of typical International Space Station (ISS) interference structures. The calculations produced using the GTD code compared to the measured results better than the ray tracing technique. The agreement was good, demonstrating that the phase errors due to multipath can be modeled and characterized using the GTD technique and characterized to a lesser fidelity using the DECAT technique. However, some discrepancies were observed. Most of the discrepancies occurred at lower devations and were either due to phase center deviations of the antenna, the background multipath environment, or the receiver itself. Selected measured and predicted differential carrier phase error results are presented and compared. Results indicate that reflections and diffractions caused by the multipath producers, located near the GPS antennas, can produce phase shifts of greater than 10 mm, and as high as 95 mm. It should be noted tl the field test configuration was meant to simulate typical ISS structures, but the two environments are not identical. The GZ and DECAT techniques have been used to calculate phase errors due to multipath o the ISS configuration to quantify the expected attitude determination errors.

Gomez, Susan F.↗

NASA NTRS2014

Open Circuit Resonant (SansEC) Sensor Technology for Lightning Mitigation and Damage Detection and Diagnosis for Composite Aircraft Applications

Traditional methods to protect composite aircraft from lightning strike damage rely on a conductive layer embedded on or within the surface of the aircraft composite skin. This method is effective at preventing major direct effect damage and minimizes indirect effects to aircraft systems from lightning strike attachment, but provides no additional benefit for the added parasitic weight from the conductive layer. When a known lightning strike occurs, the points of attachment and detachment on the aircraft surface are visually inspected and checked for damage by maintenance personnel to ensure continued safe flight operations. A new multi-functional lightning strike protection (LSP) method has been developed to provide aircraft lightning strike protection, damage detection and diagnosis for composite aircraft surfaces. The method incorporates a SansEC sensor array on the aircraft exterior surfaces forming a "Smart skin" surface for aircraft lightning zones certified to withstand strikes up to 100 kiloamperes peak current. SansEC sensors are open-circuit devices comprised of conductive trace spiral patterns sans (without) electrical connections. The SansEC sensor is an electromagnetic resonator having specific resonant parameters (frequency, amplitude, bandwidth & phase) which when electromagnetically coupled with a composite substrate will indicate the electrical impedance of the composite through a change in its resonant response. Any measureable shift in the resonant characteristics can be an indication of damage to the composite caused by a lightning strike or from other means. The SansEC sensor method is intended to diagnose damage for both in-situ health monitoring or ground inspections. In this paper, the theoretical mathematical framework is established for the use of open circuit sensors to perform damage detection and diagnosis on carbon fiber composites. Both computational and experimental analyses were conducted to validate this new method and system for aircraft composite damage detection and diagnosis. Experimental test results on seeded fault damage coupons and computational modeling simulation results are presented. This paper also presents the shielding effectiveness along with the lightning direct effect test results from several different SansEC LSP and baseline protected and unprotected carbon fiber reinforced polymer (CFRP) test panels struck at 40 and 100 kiloamperes following a universal common practice test procedure to enable damage comparisons between SansEC LSP configurations and common practice copper mesh LSP approaches. The SansEC test panels were mounted in a LSP test bed during the lightning test. Electrical, mechanical and thermal parameters were measured during lightning attachment and are presented with post test nondestructive inspection comparisons. The paper provides correlational results between the SansEC sensors computed electric field distribution and the location of the lightning attachment on the sensor trace and visual observations showing the SansEC sensor's affinity for dispersing the lightning attachment.

Szatkowski, George N.↗

NASA NTRS2004

Modeling of a Variable Focal Length Flat Lens Using Left Handed Metamaterials

Left Handed Metamaterials (LHM) were originally purposed by Victor Veselago in1968. These substances would allow a flat structure to focus electromagnetic (EM) waves because they have a negative index of refraction. A similar structure made from conventional materials, those with a positive index of refraction, would disperse the waves. But until recently, these structures have been purely theoretical because substances with both a negative permittivity and negative permeability, material properties necessary for a negative index of refraction, do not naturally exist, Recent developments have produced a structure composed of an array of thin wires and split ring resonators that shows a negative index of refraction. area smaller than a square wavelength. How small the area is can be determined by how perfectly the lens is polished and how pure the substance is that composes the lens. These lenses must also be curved for focusing to occur. The focal length is determined by the curvature of the lens and the material. On the other hand, a flat structure made from LHM would focus light because of the effect of a negative index of refraction in Snell s law. The focal length could also be varied by simply adjusting the distance of the lens from the source of radiation. This could create many devices that are adjustable to different situations in fields such as biomedical imaging and communication. the software package XFDTD which solves Maxwell s equations in the frequency domain as well as the time domain. The program used Drude models of materials to simulate the effect of negative permittivity and negative permeability. Because of this, a LHM can be simulated as a solid block of material instead of an array of wires and split ring resonators. After a flat lens is formed, I am to examine the focusing effect of the lens and determine if a higher resolution flat lens can be developed. Traditional lenses made from conventional materials cannot focus an EM wave onto an My goal was to model LHMs and create a flat lens from them. This was to be done using

Reinert, Jason↗

NASA NTRS2012

Implementation and Testing of the First Prompt Search for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (OW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient OW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize OW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 172009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify OW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency OW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with approx.50% or better probability with a few pointings of wide-field telescopes. Key words. gravitational waves - methods: observational

Camp, Jordan↗

NASA NTRS2012

Antenna Measurements: Test & Analysis of the Radiated Emissions/Immunity of the NASA/Orion Spacecraft Dart Parachute Simulator & Prototype Capsule - The Crew Exploration Vehicle

For future NASA Manned Space Exploration of the Moon and Mars, a blunt body capsule, called the Orion Crew Exploration Vehicle (CEV), composed of a Crew Module (CM) and a Service Module (SM), with a parachute decent assembly is planned for reentry back to Earth. A Capsule Parachute Assembly System (CPAS) is being developed for preliminary prototype parachute drop tests at the Yuma Proving Ground (YPG) to simulate high-speed reentry to Earth from beyond Low-Earth-Orbit (LEO) and to provide measurements of position, velocity, acceleration, attitude, temperature, pressure, humidity, and parachute loads. The primary and secondary (backup) avionics systems on CPAS also provide mission critical firing events to deploy, reef, and release the parachutes in three stages (extraction, drogues, mains) using mortars and pressure cartridge assemblies. In addition, a Mid-Air Delivery System (MDS) is used to separate the capsule from the sled that is used to eject the capsule from the back of the drop plane. Also, high-speed and high-definition cameras in a Video Camera System (VCS) are used to film the drop plane extraction and parachute landing events. Intentional and unintentional radiation emitted from and received by antennas and electronic devices on/in the CEV capsule, the MDS sled, and the VCS system are being tested for radiated emissions/immunity (susceptibility) ~ (RE/RS). To verify Electromagnetic Compatibility (EMC) of the Orion capsule, Electromagnetic Interference (EMI) measurements are being made inside a semi-anechoic chamber at NASA/JSC on the components of the CPAS system. Measurements are made at 1m from the components-under-test (CUT). In addition, EMI measurements of the integrated CEV system are being made inside a hanger at YPG. These measurements are made in a complete circle, at 30◦ angles or less, around the Orion Capsule, the spacecraft system under-test (SUT). Near-field B-Dot probe measurements on the surface of the Orion capsule are being extrapolated outward to the 1m standard distance for comparison to the MIL-STD radiated emissions limit, and far-field hybrid antenna measurements at 3m are being extrapolated inward to the 1m distance for similar comparisons.

Norgard, John D.↗

NASA NTRS2011

Implementation and Testing of the First Prompt for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW-event candidates and to reconstruct-maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline s ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with 50% or better probability with a few pointings of wide-field telescopes.

Abadie, J.↗

NASA NTRS2003

On the Determination of Local Sites for Shallow Ground Water Prospection Using Low Frequency Sounding Radars on Mars

Low frequency sounding radars should be able to probe the Martian subsurface layers down to varying depths, depending on the geo-electrical properties of the sounded sites. We present in this work four frequency dependent geo-electrical models of the Martian subsurface in the 1-20 MHz frequency band, based on laboratory electromagnetic characterization of Martian soil analogues. Those models correspond to local Martian sites, where we considered particular interest for the search of water using mainly the Ground Penetrating Radar (GPR) instrument of the Netlander mission. Results and discussion are also valid for both sounding experiments MARSIS and SHARAD. The four models of the Martian subsurface are designed to represent terrains where recent fluvial like features suggest the presence of near subsurface ground ice and probably liquid water. We performed measurements on volcanic and sedimentary materials that may be present on these sites under the appropriate geophysical conditions that may exist in those terrains. We then simulated the backscattered radar echo arising from each site in the 2 MHz frequency band, using the Finite Difference Time Domain (FDTD) algorithm, in order to evaluate the instrument performances to probe the subsurface stratigraphy of each site. Our results confirm that the near subsurface rich iron oxide mineralogy controls the instrument performances in terms of penetration depth and signal to noise ratio in the 2 MHz frequency band. We also discuss the geophysical and geo-electrical sounding conditions that could lead to an ambiguous detection of shallow subsurface water on Mars for the Netlander GPR. Finally we hope to present by the conference time a terrestrial test analogue site selected according to those creteria where deep sub-surface water have been detected 600 m deep using a prototype of the 2 MHz Netlander GPR in a field survey performed last febrary in the eastern part of the Egyptian desert.

Heggy, E.↗

NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

NASA NTRS1991

Remote sensing of earth terrain

In remote sensing, the encountered geophysical media such as agricultural canopy, forest, snow, or ice are inhomogeneous and contain scatters in a random manner. Furthermore, weather conditions such as fog, mist, or snow cover can intervene the electromagnetic observation of the remotely sensed media. In the modelling of such media accounting for the weather effects, a multi-layer random medium model has been developed. The scattering effects of the random media are described by three-dimensional correlation functions with variances and correlation lengths corresponding to the fluctuation strengths and the physical geometry of the inhomogeneities, respectively. With proper consideration of the dyadic Green's function and its singularities, the strong fluctuation theory is used to calculate the effective permittivities which account for the modification of the wave speed and attenuation in the presence of the scatters. The distorted Born approximation is then applied to obtain the correlations of the scattered fields. From the correlation of the scattered field, calculated is the complete set of scattering coefficients for polarimetric radar observation or brightness temperature in passive radiometer applications. In the remote sensing of terrestrial ecosystems, the development of microwave remote sensing technology and the potential of SAR to measure vegetation structure and biomass have increased effort to conduct experimental and theoretical researches on the interactions between microwave and vegetation canopies. The overall objective is to develop inversion algorithms to retrieve biophysical parameters from radar data. In this perspective, theoretical models and experimental data are methodically interconnected in the following manner: Due to the complexity of the interactions involved, all theoretical models have limited domains of validity; the proposed solution is to use theoretical models, which is validated by experiments, to establish the region in which the radar response is most sensitive to the parameters of interest; theoretically simulated data will be used to generate simple invertible models over the region. For applications to the remote sensing of sea ice, the developed theoretical models need to be tested with experimental measurements. With measured ground truth such as ice thickness, temperature, salinity, and structure, input parameters to the theoretical models can be obtained to calculate the polarimetric scattering coefficients for radars or brightness temperature for radiometers and then compare theoretical results with experimental data. Validated models will play an important role in the interpretation and classification of ice in monitoring global ice cover from space borne remote sensors in the future. We present an inversion algorithm based on a recently developed inversion method referred to as the Renormalized Source-Type Integral Equation approach. The objective of this method is to overcome some of the limitations and difficulties of the iterative Born technique. It recasts the inversion, which is nonlinear in nature, in terms of the solution of a set of linear equations; however, the final inversion equation is still nonlinear. The derived inversion equation is an exact equation which sums up the iterative Neuman (or Born) series in a closed form and, thus, is a valid representation even in the case when the Born series diverges; hence, the name Renormalized Source-Type Integral Equation Approach.

Yueh, Herng-Aung↗

NASA NTRSDate not supplied

A Dual Linear/Circular Polarized Wideband Flat Panel Phased Array Antenna for K/Ka-band 6U CubeSat Communications Applications

Circular polarized (CP) electromagnetic waves are preferred for space-to-space and space-to-ground satellite communications (SATCOM). CubeSats with 6U sizes (1U represents a 10 cm x 10 cm x 10 cm), that will rotate together with the Moon as it orbits Earth require high gain CP antenna systems to maintain reliable communication links. Here, we have designed an 8x8 dual linear and dual circularly polarized phased array antenna (PAA) aperture that offers 22-28 GHz bandwidth, with maximum gain of 23 dBic, and beam scanning capability of ±50o with acceptable axial ratio over all scan angles. A stacked patch topology is selected to achieve wide impedance bandwidth. Each driven patch has two via feed point locations that are connected to grounded coplanar waveguide (CPW) feed structures. By utilizing phase shifters from the Anokiwave RFICs (AWMF-0165) employed for the beamforming network (BFN), we have generated both dual linear polarization and dual circular polarization in each patch. Each driven and parasitic patch is innovatively shaped to achieve excellent isolation between two feed points which, in turn, leads to low cross polarization, high efficiency and good axial ratio over the entire bandwidth. Each chip has eight outputs which feed one 2x2 sub-array consisting of four dual linear or dual circularly polarized driven patches, simultaneously. Each 2x2 sub-array of the 8x8 array are sequentially rotated which improves axial ratio performance when the antenna operates in the dual circularly polarized mode. By utilizing a mirror effect, cross-polarized fields are also improved when operated in a dual linear polarized mode. The 8x8 array aperture is simulated taking into account the effect of the multilayer printed circuit board (PCB) stack-up of the BFN and co-simulated with the power splitter outputs and chip outputs and results are compared with only passive array performance. The performance of the antenna is still excellent when non-ideal chip performances are included in the simulation analysis. The fabricated array will be tested for radiation patterns in the Antenna and Microwave Laboratory (AML) of the San Diego State University and the NASA Glenn Research Center (GRC) spherical near field chambers. The design of the 8x8 antenna array with BFN layout is shown in Fig. 1.

5G↗

NASA NTRS1997

Retrieval of Soil Moisture and Roughness from the Polarimetric Radar Response

The main objective of this investigation was the characterization of soil moisture using imaging radars. In order to accomplish this task, a number of intermediate steps had to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces was considered and a hybrid model that relates the radar backscattering coefficient to soil moisture and surface roughness was developed. This model is based on extensive experimental measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave frequencies over a wide range of moisture conditions and roughness scales in conjunction with existing theoretical surface scattering models in limiting cases (small perturbation, physical optics, and geometrical optics models). Also a simple inversion algorithm capable of providing accurate estimates of soil moisture content and surface rms height from single-frequency multi-polarization radar observations was developed. The accuracy of the model and its inversion algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized phase difference into the model. Experimental data in conjunction with numerical simulations are used to relate the soil moisture content and surface roughness to the phase difference statistics. For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random surfaces was developed. Finally the scattering problem of short vegetation cover above a rough soil surface was considered. A general scattering model for grass-blades of arbitrary cross section was developed and incorporated in a first order random media model. The vegetation model and the bare-soil model are combined and the accuracy of the combined model is evaluated against experimental observations from a wheat field over the entire growing season. A complete set of ground-truth data and polarimetric backscatter data were collected. Also an inversion algorithm for estimating soil moisture and surface roughness from multi-polarized multi-frequency observations of vegetation-covered ground is developed.

Sarabandi, Kamal↗

NASA NTRS2001

Effects of Traveling Magnetic Field on Dynamics of Solidification

TMF is based on imposing a controlled phase-shift in a train of electromagnets, forming a stack. Thus, the induced magnetic field can be considered to be travelling along the axis of the stack. The coupling of this traveling wave with an electrically conducting fluid results in a basic flow in a form of a single axisymmetric roll. The magnitude and direction of this flow can be remotely controlled. Furthermore, it is possible to localize the effect of this force field though activating only a number of the magnets. This force field generated in the fluid can, in principle, be used to control and modify convection in the molten material. For example, it can be used to enhance convective mixing in the melt, and thereby modify the interface shape, and macrosegregation. Alternatively, it can be used to counteract thermal and/or solutal buoyancy forces. High frequency TMF can be used in containerless processing techniques, such as float zoning, to affect the very edge of the fluid so that Marangoni flow can be counter balanced. The proposed program consists of basic fundamentals and applications. Our goal in conducting the following experiments and analyses is to establish the validity of TMF as a new tool for solidification processes. Due to its low power consumption and simplicity of design, this tool may find wide spread use in a variety of space experiments. The proposed ground based experiments are intended to establish the advantages and limitations of employing this technique. In the fundamentals component of the proposed program, we will use theoretical tools and experiments with mercury to establish the fundamental aspects of TMF-induced convection through a detailed comparison of theoretical predictions and experimental measurements of flow field. In this work, we will conduct a detailed parametric study involving the effects of magnetic field strength, frequency, wave vector, and the fluid geometry. The applications component of this work will be focused on investigating the effect of TMF on the following solidification and pre-directional solidification processes: (1) Bridgman growth of Ga:Ge with the goal of counteracting the buoyancy-driven convection; (2) Mixing of Pb-Ga and Pb-Sn alloys with the aim of initiating and maintaining a uniform melt prior to solidification processing; and (3) Float Zone growth with the aim of identifying, through simulations and model experiments, conditions needed to counteract Marangoni flow in a microgravity environment. The proposed research has strong relevance to microgravity research and the objectives of the NRA. TMF can provide a unique and accurate mechanism for generation and control of desirable flow patterns for microgravity research. These attributes have significant relevance to 1) Alloy mixing prior to solidification in a microgravity environment. TMF can provide this mixing with a low level of power consumption; (2) TMF can offset the deleterious effects of Marangoni convection in microgravity containerless processing. Thus, TMF can be instrumental in further understanding this phenomena; (3) Generation of controlled flows will allow the investigation of the effect of these flows on growth morphology and growth kinetics; and (4) On Earth, TMF has the potential to significantly counter-balance thermosolutal convection, thereby creating conditions similar to those obtained in microgravity. Once demonstrated, this new tool for use in solidification has the strong potential to find applications in a host of microgravity material research projects.

Mazuruk, Konstantin↗

NASA NTRS2016

Photospheric Current Spikes and Their Possible Association with Flares - Results from an HMI Data Driven Model

A data driven, near photospheric magnetohydrodynamic model predicts spikes in the horizontal current density, and associated resistive heating rate per unit volume Q. The spikes appear as increases by orders of magnitude above background values in neutral line regions (NLRs) of active regions (ARs). The largest spikes typically occur a few hours to a few days prior to M or X flares. The spikes correspond to large vertical derivatives of the horizontal magnetic field. The model takes as input the photospheric magnetic field observed by the Helioseismic & Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) satellite. This 2.5 D field is used to determine an analytic expression for a 3 D magnetic field, from which the current density, vector potential, and electric field are computed in every AR pixel for 14 ARs. The field is not assumed to be force-free. The spurious 6, 12, and 24 hour Doppler periods due to SDO orbital motion are filtered out of the time series of the HMI magnetic field for each pixel using a band pass filter. The subset of spikes analyzed at the pixel level are found to occur on HMI and granulation scales of 1 arcsec and 12 minutes. Spikes are found in ARs with and without M or X flares, and outside as well as inside NLRs, but the largest spikes are localized in the NLRs of ARs with M or X flares. The energy to drive the heating associated with the largest current spikes comes from bulk flow kinetic energy, not the electromagnetic field, and the current density is highly non-force free. The results suggest that, in combination with the model, HMI is revealing strong, convection driven, non-force free heating events on granulation scales, and that it is plausible these events are correlated with subsequent M or X flares. More and longer time series need to be analyzed to determine if such a correlation exists. Above an AR dependent threshold value of Q, the number of events N(Q) with heating rates greater than or equal to Q obeys a scale invariant power law distribution for each AR given by N(Q) varies Q(sup -s), where 0.40 less than or equal to S less than or equal to 0.53, with a mean and standard deviation across the 14 ARs of 0.47 and 0.045, showing there is little variation of S from one AR to another. These properties of N(Q) are in close agreement with those of the distribution N(E) for the total energy E of solar flares, determined from observations to be N(E) = constant x E(sup -alpha). From observations of nanoflares in the 0.7 to 4 MK range, and from observations of flares in hard X-rays, it is found that 0.51 less than or equal to alpha less than or equal to 0.57, and 0.4 less than or equal to alpha less than or equal to 0.6, respectively (Crosby et al. 1993, Sol. Phys., 143, 275; Aschwanden & Parnell 2002, ApJ, 572, 1048). Observations also show that, as is found here for the exponent S, there is little variation of alpha with AR (Wheatland 2000, ApJ, 532, 1209), indicating N(E) and N(Q) are largely independent of individual properties of ARs such as area, total magnetic flux, and distribution of current density (i.e. non-potentiality). Therefore the power law scaling of the photospheric heating rate Q computed here on granulation scales is essentially identical to that found for coronal observations of flare energies on scales 1-2 orders of magnitude larger. This suggests the physical mechanisms that cause Q and coronal flares are closely related. It seems likely that Q is the signature of a magnetic reconnection process in an energy range and volume orders of magnitude smaller than those of flares. In this context, at least the larger spikes in Q might be signatures of UV photospheric or lower chromospheric bombs in which plasma is heated to temperatures approximately 10(exp -5) K (Peter et al. 2014, Science 346, 1255726; Judge 2015, ApJ, 808, 116). In addition, lattice based avalanche simulations of flare energy release predict 0.4 less than or equal to alpha less than or equal to 0.5, while analytic, fractal-diffusive self-organized criticality models predict 0.4 less than or equal to alpha less than or equal to 0.67, in excellent agreement with observations, and the results presented here (Aschwanden & Parnell 2002, ApJ, 572, 1048; Aschwanden 2012, A&A, 539, A2; Aschwanden 2013, in "Self Organized Criticality Systems"; Aschwanden et al. 2016, SSR, 198, 47).

Goodman, Michael↗

NASA NTRS2016

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

The numerical simulation of many aerodynamic non-periodic flows of practical interest involves discretized computational domains that often must be artificially truncated. Appropriate boundary conditions are required at these truncated domain boundaries, and ideally, these boundary conditions should be perfectly "absorbing" or "nonreflecting" so that they do not contaminate the flow field in the interior of the domain. The proper specification of these boundaries is critical to the stability, accuracy, convergence, and quality of the numerical solution, and has been the topic of considerable research. The need for accurate boundary specification has been underscored in recent years with efforts to apply higher-fidelity methods (DNS, LES) in conjunction with high-order low-dissipation numerical schemes to realistic flow configurations. One of the most popular choices for specifying these boundaries is the characteristics-based boundary condition where the linearized flow field at the boundaries are decomposed into characteristic waves using either one-dimensional Riemann or other multi-dimensional Riemann approximations. The values of incoming characteristics are then suitably modified. The incoming characteristics are specified at the in flow boundaries, and at the out flow boundaries the variation of the incoming characteristic is zeroed out to ensure no reflection. This, however, makes the problem ill-posed requiring the use of an ad-hoc parameter to allow small reflections that make the solution stable. Generally speaking, such boundary conditions work reasonably well when the characteristic flow direction is normal to the boundary, but reflects spurious energy otherwise. An alternative to the characteristic-based boundary condition is to add additional "buffer" regions to the main computational domain near the artificial boundaries, and solve a different set of equations in the buffer region in order to minimize acoustic reflections. One approach that has been used involves modeling the pressure fluctuations as acoustic waves propagating in the far-field relative to a single noise-source inside the buffer region. This approach treats vorticity-induced pressure fluctuations the same as acoustic waves. Another popular approach, often referred to as the "sponge layer," attempts to dampen the flow perturbations by introducing artificial dissipation in the buffer region. Although the artificial dissipation removes all perturbations inside the sponge layer, incoming waves are still reflected from the interface boundary between the computational domain and the sponge layer. The effect of these refkections can be somewhat mitigated by appropriately selecting the artificial dissipation strength and the extent of the sponge layer. One of the most promising variants on the buffer region approach is the Perfectly Matched Layer (PML) technique. The PML technique mitigates spurious reflections from boundaries and interfaces by dampening the perturbation modes inside the buffer region such that their eigenfunctions remain unchanged. The technique was first developed by Berenger for application to problems involving electromagnetic wave propagation. It was later extended to the linearized Euler, Euler and Navier-Stokes equations by Hu and his coauthors. The PML technique ensures the no-reflection property for all waves, irrespective of incidence angle, wavelength, and propagation direction. Although the technique requires the solution of a set of auxiliary equations, the computational overhead is easily justified since it allows smaller domain sizes and can provide better accuracy, stability, and convergence of the numerical solution. In this paper, the PML technique is developed in the context of a high-order spectral-element Discontinuous Galerkin (DG) method. The technique is compared to other approaches to treating the in flow and out flow boundary, such as those based on using characteristic boundary conditions and sponge layers. The superiority of the current PML technique over other approaches is demonstrated for a range of test cases, viz., acoustic pulse propagation, convective vortex, shear layer flow, and low-pressure turbine cascade flow. The paper is structured as follows. We first derive the PML equations from the non{linear Euler equations. A short description of the higher-order DG method used is then described. Preliminary results for the four test cases considered are then presented and discussed. Details regarding current work that will be included in the final paper are also provided.

Spectral↗

NASA NTRSDate not supplied

Hall Magnetohydrodynamic Power Generation and Drag Augmentation Using A Coaxial-Electrode Configuration During Hypersonic Entry

For interplanetary missions, a spacecraft must reduce its relative velocity from orbital speeds of multiple kilometers per second to zero in order to safely land on the surface of a planetary body. During this phase known as planetary reentry, the spacecraft undergoes high speed above Mach 5 up to Mach 30. This hyper-sonic flight regime imposes a strong bow shock in front of the vehicle that imparts extreme aerodynamic and thermal effects. In addition, during portions of this regime, a plasma is formed in the post-shock flow-field around the vehicle due to these high temperatures. Interaction between an applied magnetic field with this electrically conductive fluid, or magnetohydrodynamic (MHD) interaction, can be utilized to convert kinetic energy from the flow into storable electrical energy. Onboard MHD energy generation can benefit the spacecraft by increasing the power of active thermal control components and vehicle control systems during reentry. MHD interaction also exerts a Lorentz body force on the plasma flow that can be manipulated to act as an additional drag force through the vehicle body. MHD drag augmentation can benefit the spacecraft with an additional control mechanism that operates without moving parts. This can be utilized at higher altitudes and subsequently reduce convective heat transfer to the vehicle. MHD interaction with plasma flow suggests that both power generation and drag augmentation can be integrated together in one design. Modern numerical studies have only demonstrated the viability and potential of MHD energy generation and drag modulation separately for implementation into future spacecraft designs [1] [2] [3]. Steeves et. al. presented one possible generator design specifically for reentry vehicles that utilized a modular panel with two embedded electromagnets and two extruding electrodes for air plasma to flow between. Fujino et. al. presented results that indicated increased total drag with a permanent magnet embedded in the blunt body of the Orbital Reentry Experiment (OREX) trajectory. Furthermore, experimental investigations of these two MHD applications have been limited and separate due to the difficult technical nature of experimental de-sign and analysis. For MHD energy generation, there have been two designs which used a cylindrical fore-body with two embedded permanent magnets and two outward facing curved electrodes for tangency to an artificially ionized plasma that flowed along the body length [4] [5]. The first design was studied using an air microwave (MW) ionized supersonic plasma wind tunnel and the second design used radio frequency (RF) ionization of Argon, but both physically demonstrated the feasibility of MHD energy generation in reentry plasma conditions. For MHD drag modulation, one study involved a small-scale cylindrical blunt body embedded with permanent magnets immersed in plasma flow in an arcjet tunnel [6], and another utilized a spherical permanent magnet enclosed with a spherical forebody that was evaluated in a shock expansion tunnel [7]. Both experimental studies provided physical evidence of drag augmentation due to MHD interactions. Therefore, because of the multi-faceted effects of MHD along with both applications having been deemed feasible and beneficial during reentry, this indicates that a coupled power generation and drag modulation design could be achievable. The axi-symmetric nature of blunt bodies for reentry suggests that a spherical set of two ring electrodes near the nose along the body axis. Based on first principles, flow along the walls of the spherical body can then produce MHD interactions between the electrodes. As a result, two currents are generated: a Hall current that produces power and a primary inductive current that augments the drag force. The goal of this work is to provide a simulated design of a coaxial-electrode con-figuration for power generation and drag modulation for experimental testing.

E. Leong↗