Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Saturation problem”

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 163 records · Page 9

Attitude Determination and Control Subsystem (ADCS) Preparations for the EPOXI Flyby of Comet Haley 2

On November 4, 2010 the already "in-flight" Deep Impact spacecraft flew within 700km of comet 103P/Hartley 2 as part of its extended mission EPOXI, the 5th time to date any spacecraft visited a comet. In 2005, the spacecraft had previously imaged a probe impact comet Tempel 1. The EPOXI flyby marked the first time in history that two comets were explored with the same instruments on a re-used spacecraft-with hardware and software originally designed and optimized for a different mission. This made the function of the attitude determination and control subsystem (ADCS) critical to the successful execution of the EPOXI flyby. As part of the spacecraft team preparations, the ADCS team had to perform thorough sequence reviews, key spacecraft activities and onboard calibrations. These activities included: review of background sequences for the initial conditions vector, sun sensor coefficients, and reaction wheel assembly (RWA) de-saturations; design and execution of 10 trajectory correction maneuvers; science calibration of the two telescope instruments; a flight demonstration of the fastest turns conducted by the spacecraft between Earth and comet point; and assessment of RWA health (given RWA problems on other spacecraft).

Harley↗

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Generation of Bright Soft X-ray Laser Beams (Final Report)

The project goal was to demonstrate new compact soft x-ray lasers emitting high energy ultrashort pulses at shorter wavelengths and increased repetition rates by amplification in atomic transitions in high density laser-created plasmas. To demonstrate these lasers we developed new high power solid state lasers. The project succeeded in demonstrating the shortest wavelength compact soft x-ray laser to date, 6.85 nm wavelength, and the highest repetition rate table-top soft x-ray laser, 400 Hz. We also applied these compact soft x-ray laser sources to determine the composition of materials with nanoscale resolution using a new technique that combines soft x-ray laser ablation with mass spectrometry, and we demonstrated that these lasers are capable of acquiring nanoscale resolution holograms with picosecond time resolution. The combination of an increased pumping efficiency with increased repetition rates can yield soft x-ray laser beams with an unprecedented average power on a table-top for applications demanding a high photon flux. We have demonstrated that the efficiency of plasma-based laser-pumped soft x-ray lasers can be significantly increased by detailed pulse shaping of the pump laser using a new pulse synthesizer. Results include the demonstration of 18.9 nm wavelength laser operation at a record 400 Hz repetition rate, and the extension of repetitive gain-saturated plasma-based table-top x-ray lasers to the shortest wavelength to date, 6.85 nm. Laser amplification was also observed at wavelengths as short at 5.8 nm by isolectronic scaling in Ni-like lanthanide ions. We have also conducted simulations to explore the possibility of demonstrating gain saturated atomic x-ray lasers near the carbon K-edge, in the “carbon window” for the first time. The work conducted under this program also includes the use of compact soft x-ray lasers in applications to problems of chemical interest, such as isotopic composition mapping of micron-size uranium particles with nanoscale resolution, and nanoscale resolution depth profiling analysis of CoNCN-coated electrodes for water oxidation catalysis. We also demonstrated single-shot high resolution soft x-ray Fourier transform holography over a broad 7 µm diameter field of view with ~ 5 ps temporal resolution.

74 ATOMIC AND MOLECULAR PHYSICS↗

The Case for a Heat-Pipe Phase of Planet Evolution on the Moon

The prevalence of anorthosite in the lunar highlands is generally attributed to the flotation of less dense plagioclase in the late stages of the solidification of the lunar magma ocean. It is not clear, however, that these models are capable of producing the extremely high plagioclase contents (near 100%) observed in both Apollo samples and remote sensing data, since a mostly solid lithosphere forms (at 60-70% solidification) before plagioclase feldspar reaches saturation (at approximately 80% solidification). Formation as a floating cumulate is made even more problematic by the near uniformity of the alkali composition of the plagioclase, even as the mafic phases record significant variations in Mg/(Mg+Fe) ratios. These problems can be resolved for the Moon if the plagioclase-rich crust is produced and refined through a widespread episode of heat-pipe magmatism rather than a process dominated by density-driven plagioclase flotation. Heat-pipes are an important feature of terrestrial planets at high heat flow, as illustrated by Io's present activity. Evidence for their operation early in Earth's history suggests that all terrestrial bodies should experience an early episode of heat-pipe cooling. As the Moon likely represents the most wellpreserved example of early planetary thermal evolution in our solar system, studies of the lunar surface and of lunar materials provide useful data to test the idea of a universal model of the way terrestrial bodies transition from a magma ocean state into subsequent single-plate, rigid-lid convection or plate tectonic phases.

Simon, J. I.↗

Hanford Site Composite Analysis: LLBG-200-W B Vadose Zone Model

The objectives of the vadose modeling for the updated Hanford Site Composite Analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates to the CA saturated zone model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP 1 ) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the saturated zone model. The rates will be summed over the 100 by 100 m saturated zone model grid cells that fall within the vadose zone model source domain.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Problems Associated with Statistical Pattern Recognition of Acoustic Emission Signals in a Compact Tension Fatigue Specimen

Acoustic emission (AE) data were acquired during fatigue testing of an aluminum 2024-T4 compact tension specimen using a commercially available AE system. AE signals from crack extension were identified and separated from noise spikes, signals that reflected from the specimen edges, and signals that saturated the instrumentation. A commercially available software package was used to train a statistical pattern recognition system to classify the signals. The software trained a network to recognize signals with a 91-percent accuracy when compared with the researcher's interpretation of the data. Reasons for the discrepancies are examined and it is postulated that additional preprocessing of the AE data to focus on the extensional wave mode and eliminate other effects before training the pattern recognition system will result in increased accuracy.

Hinton, Yolanda L.↗

One-carbon (bio ?) Geochemistry in Subsurface Waters of the Serpentinizing Coast Range Ophiolite

Serpentinization - the aqueous alteration of ultramafic rocks - typically imparts a highly reducing and alkaline character to the reacting fluids. In turn, these can influence the speciation and potential for metabolism of one-carbon compounds in the system. We examined the aqueous geochemistry and assessed the biological potential of one-carbon compounds in the subsurface of the McLaughlin Natural Reserve (Coast Range Ophiolite, California, USA). Fluids from wells sunk at depths of 25-90 meters have pH values ranging from 9.7 to 11.5 and dissolved inorganic carbon (DIC concentrations) generally below 60 micromolar. Methane is present at concentrations up to 1.3 millimolar (approximately one-atmosphere saturation), and hydrogen concentrations are below 15 nanomolar, suggesting active consumption of H2 and production of CH4. However, methane production from CO2 is thermodynamically unfavorable under these conditions. Additionally, the speciation of DIC predominantly into carbonate at these high pH values creates a problem of carbon availability for any organisms that require CO2 (or bicarbonate) for catabolism or anabolism. A potential alternative is carbon monoxide, which is present in these waters at concentrations 2000-fold higher than equilibrium with atmospheric CO. CO is utilized in a variety of metabolisms, including methanogenesis, and bioavailability is not adversely affected by pH-dependent speciation (as for DIC). Methanogenesis from CO under in situ conditions is thermodynamically favorable and would satisfy biological energy requirements with respect to both Gibbs Energy yield and power.

McLaughlin Natural Reserve↗

Correlated Electromagnetic Levitation Actuator: A Reaction Sphere Based Attitude Control System

To address problems experienced by current reaction wheels and control moment gyroscopebased attitude control systems (ACS), researchers at NASA’s Marshall Space Flight Center have begun developing a reaction sphere actuator based on correlated electromagnetic levitation that will be immune to destructive bearing friction, momentum saturation, and gimbal lock. The Correlated Electromagnetic Levitation Actuator (CELA) advances the state of the art of reaction sphere ACSs by employing the concept of correlated magnetics. It is a frictionless, direct-drive reaction sphere that harnesses a unique technology with an array of applications across multiple disciplines. Correlated electromagnets function in a manner that is analogous to a matched filter; the convolution of two signals is peaked at the index representing the greatest match. For CELA, the signals are the patterns of magnetic flux density as a function of position. The magnitude of the convolution equates to an attractive or repulsive force, and these forces can be azimuthal or radial. The development of CELA is based in four distinct disciplines: Advanced Manufacturing, Prototype Development, Electromagnetic Modeling, and Controls. We are developing novel manufacturing techniques required to build arrays of permanent and electromagnet dipoles on curved surfaces. To print the permanent magnetic array, we have developed a probe with pyramidal magnets that will reside on a robotic arm to induce localized magnetic fields on a surface. The probe also includes the ability to erase dipole patterns from a permanent magnet by heating the surface to its Curie temperature. A number of test articles and prototypes have been developed using additive manufacturing methods. These prototypes have included hemispherical motors to test the drive algorithm, and a levitation test bed that demonstrates a magnetic bearing method based on attractive magnetic forces and ratiometric Hall effect sensors. We developed an array of electromagnetic dipoles on a printed circuit board (PCB) with individual H-bridges controlling each coil. This device created various flux density patterns and we measured their magnetic fields using a custom Hall effect 3-D probe and a LabVIEW virtual instrument. These data will serve as a benchmark for characterizing the accuracy of future models. Current work is focused on modeling the magnetic fields of our prototype arrays using COMSOL Finite Element Analysis and verifying the model against our test data. Accurate modeling will allow us to quickly test new patterns of electromagnets and their macro behavior. Eventually, the magnetic field models will be implemented in our controls simulations to facilitate precise control of the reaction sphere. Initial model results agree with field measurements to within 1 G (5% of measured flux density). Currently, we are testing different material properties of the electromagnets and their magnetic fields and thermal effects. These results will be used to refine the design of the electromagnetic dipoles. Our control efforts have centered on developing commutation, levitation, and field pattern shaping hardware in the form of breadboards and PCBs with software running on a local microcontroller. In addition, our partners developed MATLAB Simulink models to demonstrate a PID controller thatmitigates disturbance forces resulting from the interaction of drive and levitation magnetics. Finally, we have designed a three-axis test stand that will be used in future work to demonstrate CELA’s orientation control capability.

controls↗

Acceleration of electrons and ions by strong lower-hybrid turbulence in solar flares

One of the outstanding problems in solar flare theory is how to explain the 10-20 keV and greater hard x-ray emissions by a thick target bremsstrahlung model. The model requires the acceleration mechanism to accelerate approximately 10(exp 35) electrons sec(exp -l) with comparable energies, without producing a large return current which persists for long time scales after the beam ceases to exist due to Lenz's law, thereby, producing a self-magnetic field of order a few mega-Gauss. In this paper, we investigate particle acceleration resulting from the relaxation of unstable ion ring distributions, producing strong wave activity at the lower hybrid frequency. It is shown that strong lower hybrid wave turbulence collapses in configuration space producing density cavities containing intense electrostatic lower hybrid wave activity. The collapse of these intense nonlinear wave packets saturate by particle acceleration producing energetic electron and ion tails. There are several mechanisms whereby unstable ion distributions could be formed in the solar atmosphere, including reflection at perpendicular shocks, tearing modes, and loss cone depletion. Numerical simulations of ion ring relaxation processes, obtained using a 2 1/2-D fully electromagnetic, relativistic particle in cell code are discussed. We apply the results to the problem of explaining energetic particle production in solar flares. The results show the simultaneous acceleration of both electrons and ions to very high energies: electrons are accelerated to energies in the range 10-500 keV, while ions are accelerated to energies of the order of MeVs, giving rise to x-ray emission and gamma-ray emission respectively. Our simulations also show wave generation at the electron cyclotron frequency. We suggest that these waves are the solar millisecond radio spikes. The strong turbulence collapse process leads to a highly filamented plasma producing many localized regions for particle acceleration and resulting in approximately 10(exp 17) electron 'beamlets' of width approximately equal to 10 lambda sub De which eliminates the production of large magnetic fields. In this paper, we demonstrate that the model produces an energetic electron spectrum with the right flux to account for the hard x-ray observations.

Spicer, D. S.↗

A Search for Coronal Emission at the Bottom of the Main-Sequence: Stars and Brown Dwarf Candidates with Spectral Types Later than M7 and the Rotation-Activity Relation

This program intended to test whether the lowest mass stars at the bottom end of the main sequence and the lower mass brown dwarfs have coronae. If they have coronae, what are the coronal characteristics and what drives them? In the classical dynamo picture, the closed magnetic loop structure is generated near the boundary of the convective envelope and the radiative core. Stars with mass below 0.30 Msun however are fully convective, and the nature of the dynamo responsible for the generation of the coronae in this regime is poorly understood. Previous results from the ROSAT mission (e.g., Fleming et al. 1993, 1995; Schmitt et al. 1995) had confirmed three very important characteristics of M-star coronae: (1) a very high percentage of all M dwarfs have coronae (of order 85% in the local 7 pc sample), (2) those M dwarfs showing high chromospheric activity, such as having the Balmer series in emission or large/numerous optical flaring, indeed exhibit the highest coronal activity, and (3) that the maximum saturation boundary in X-ray luminosity, which amounts to 0.0001-0.001 for Lx/Lbol for the dMe stars, extends down to the current detection limit, through spectral types M7. It was likely that the incompleteness noted for result (1) above was simply a detection limit problem; for more distant sources, the X-ray fainter dM stars will drop below detection thresholds before the more X-ray luminous dMe stars. The latest stars for which direct detection of the corona had been successful were of spectral type dM7 (e.g., VB8, LHS 3003). This program proposed to obtain ROSAT HRI observations for a large number of the coolest known (at that time) stars at the bottom of the main-sequence, which had spectral types of M9 or later. Three stars were approved for observations with ROSAT-HRI totaling 180 ksec. The goal was to obtain X-ray detections or low upper limits for the three approved stars.

Stringfellow, Guy↗

Ground-based Photon Path Measurements from Solar Absorption Spectra of the O2 A-band

High-resolution solar absorption spectra obtained from Table Mountain Facility (TMF, 34.38degN, 117.68degW, 2286 m elevation) have been analyzed in the region of the O2 A-band. The photon paths of direct sunlight in clear sky cases are retrieved from the O2 absorption lines and compared with ray-tracing calculations based on the solar zenith angle and surface pressure. At a given zenith angle, the ratios of retrieved to geometrically derived photon paths are highly precise (approx.0.2%), but they vary as the zenith angle changes. This is because current models of the spectral lineshape in this band do not properly account for the significant absorption that exists far from the centers of saturated lines. For example, use of a Voigt function with Lorentzian far wings results in an error in the retrieved photon path of as much as 5%, highly correlated with solar zenith angle. Adopting a super-Lorentz function reduces, but does not completely eliminate this problem. New lab measurements of the lineshape are required to make further progress.

Yang, Z.↗

Fourier-MIONet: Fourier-enhanced multiple-input neural operators for multiphase modeling of geological carbon sequestration

Geologic carbon sequestration (GCS) is a safety-critical technology that aims to reduce the amount of carbon dioxide in the atmosphere, which also places high demands on reliability. Multiphase flow in porous media is essential to understand CO 2 migration and pressure fields in the subsurface associated with GCS. However, numerical simulation for such problems in 4D is computationally challenging and expensive, due to the multiphysics and multiscale nature of the highly nonlinear governing partial differential equations (PDEs). It prevents us from considering multiple subsurface scenarios and conducting real-time optimization. Here, we develop a Fourier-enhanced multiple-input neural operator (Fourier-MIONet) to learn the solution operator of the problem of multiphase flow in porous media. Fourier-MIONet utilizes the recently developed framework of the multiple-input deep neural operators (MIONet) and incorporates the Fourier neural operator (FNO) in the network architecture. Once Fourier-MIONet is trained, it can predict the evolution of saturation and pressure of the multiphase flow under various reservoir conditions, such as permeability and porosity heterogeneity, anisotropy, injection configurations, and multiphase flow properties. Compared to the enhanced FNO (U-FNO), the proposed Fourier-MIONet has 90% fewer unknown parameters, and it can be trained in significantly less time (about 3.5 times faster) with much lower CPU memory (<15%) and GPU memory (<35%) requirements, to achieve similar prediction accuracy. In addition to the lower computational cost, Fourier-MIONet can be trained with only 6 snapshots of time to predict the PDE solutions for 30 years. Furthermore, we observed that Fourier-MIONet can maintain good accuracy when predicting out-of-distribution (OOD) data. The excellent generalizability of Fourier-MIONet is enabled by its adherence to the physical principle that the solution to a PDE is continuous over time. Furthermore, the developed Fourier-MIONet makes it possible to solve the long-time evolution of geological carbon sequestration in a large-scale three-dimensional space accurately and efficiently.

97 MATHEMATICS AND COMPUTING↗

Remarks on the overall distribution of hydrogen in the galactic disk

The degree of saturation characterizing low altitude emission observations of H(I) and the optical depth corrections to the derived column and volume densities are discussed. The amount of fine-scale velocity and spatial structure diluted by the instrumental limitations of the presently available surveys are described. The general problem of detailed mapping of H(I) in the galaxy is explored. Comparison is made between the distribution of H(I) and that of CO and several other galactic tracers. Atomic hydrogen is unique in its distribution, instead of being typical of many Population I constituents. As defined by atomic hydrogen, the galactic disk has a diameter fully twice as large as that defined by the ionized and molecular states of hydrogen, as well as by other molecules, supernova remnants, pulsars, gamma-radiation, synchroton radiation, and the youngest stars. It is also less confined to the galactic equator than most of the other constituents. The degree of small scale structure apparent in the molecular observations is much greater than that in the H(I) observations.

Burton, W. B.↗

Restoration of HST images with missing data

Missing data are a fairly common problem when restoring Hubble Space Telescope observations of extended sources. On Wide Field and Planetary Camera images cosmic ray hits and CCD hot spots are the prevalent causes of data losses, whereas on Faint Object Camera images data are lossed due to reseaux marks, blemishes, areas of saturation and the omnipresent frame edges. This contribution discusses a technique for 'filling in' missing data by statistical inference using information from the surrounding pixels. The major gain consists in minimizing adverse spill-over effects to the restoration in areas neighboring those where data are missing. When the mask delineating the support of 'missing data' is made dynamic, cosmic ray hits, etc. can be detected on the fly during restoration.

Adorf, Hans-Martin↗

Elucidating the Onset of Plasticity in Sliding Contacts Using Differential Computational Orientation Tomography

Abstract Depending on the mechanical and thermal energy introduced to a dry sliding interface, the near-surface regions of the mated bodies may undergo plastic deformation. In this work, we use large-scale molecular dynamics simulations to generate “differential computational orientation tomographs” (dCOT) and thus highlight changes to the microstructure near tribological FCC alloy surfaces, allowing us to detect subtle differences in lattice orientation and small distances in grain boundary migration. The analysis approach compares computationally generated orientation tomographs with their undeformed counterparts via a simple image analysis filter. We use our visualization method to discuss the acting microstructural mechanisms in a load- and time-resolved fashion, focusing on sliding conditions that lead to twinning, partial lattice rotation, and grain boundary-dominated processes. Extracting and laterally averaging the color saturation value of the generated tomographs allows us to produce quantitative time- and depth-resolved maps that give a good overview of the progress and severity of near-surface deformation. Corresponding maps of the lateral standard deviation in the color saturation show evidence of homogenization processes occurring in the tribologically loaded microstructure, frequently leading to the formation of a well-defined separation between deformed and undeformed regions. When integrated into a computational materials engineering framework, our approach could help optimize material design for tribological and other deformation problems. Graphic Abstract .

36 MATERIALS SCIENCE↗

Magnetized thermal conduction fronts

The evolution of planar thermal conduction fronts in the presence of a dynamically weak, but otherwise self-consistent, magnetic field is considered. The field is assumed to be connected and untangled. In the diffusion limit for the thermal conductivity, these fronts exhibit self-similar behavior, even in the presence of a field. The role of the field is restricted to channeling the heat flux along its lines of force, and it enters into the problem as a dimensionless angle variable. 'Combing' (or opening) of insulating field lines by the evaporative flow is explicitly demonstrated. Unless the field is nearly perpendicular to the front normal in the hot gas, insulating effects are not profound. Self-similarity breaks down if the front becomes saturated, and under certain conditions magnetized saturated conduction fronts cannot propagate: the solution characteristics of the wave equation form caustics. The physical resolution is the advent of two-fluid (nonlocal) heating. Such Coulomb-heated fronts are expected to be relatively rare in typical astrophysical systems. The large-scale effects of a magnetic field on cloud evaporation in the interstellar medium are briefly discussed, and it is suggested that these fields preclude the presence of time-independent evaporative solutions. Thermal interfaces may then continue to evolve until radiative cooling halts their development; large tracts of warm 10,000 K gas may result.

Balbus, S. A.↗

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.↗

A Path to Improving Simulated Properties of Low Clouds over the Beaufort Sea using Airborne In Situ Observations of Subgrid-Scale Variability

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). Multiple models/reanalyses use SGS supersaturated total water (relative to GS saturation) as a threshold for partitioning available water for condensation, and indeed we find significant correlation between GS cloud water and SGS supersaturation. However, we also find that the assumption of a static threshold of 100% saturation to be unrealistic. Empirical calculations from the ARISE data show a large sensitivity of this threshold to GS relative humidity, and so a microphysical parameterization allowing the threshold to vary according to GS thermodynamic properties may result in more realistic GS cloud water values. Finally, to determine how sensitive the ARISE-derived results are to that particular campaign, we include additional data from the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) conducted in 1998. The inclusion of the second dataset will help with demonstrating the robustness of the results and their utility in improving the representation of Arctic clouds in models and reanalyses.

J Brant Dodson↗