Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Exact two component”

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 181 records · Page 10

Satellite-Based Investigations of the Transition from an Oceanic to Continental Transform Margin

Detailed characterization of neotectonics evolution of the Valle de San Felipe and Arroyo Grande regions in northern Baja California. Reoccupied GEOMEX GPS sites, and occupied a regional GPS (Global Positioning System) network. The Baja California peninsula in Mexico offers a unique setting for studying the kinematic evolution of a complex, active strike-slip/rift plate boundary. We are currently conducting remote sensing, geologic, and geodetic studies of this boundary. The combined data sets will yield instantaneous and time integrated views of its evolution. This proposal solicits renewed funding from NASA to support remote sensing and geologic studies. During the late Cenozoic, Baja California has been the locus of changing fault geometry that has accommodated components of the relative motion between the North America and Pacific plates. Contemporary slip between the two plates occurs in a broad zone that encompasses much of southern California and the Baja California Peninsula. The transfer of slip across this zone in southern California is relatively well understood. South of the border, the geometry and role of specific faults and structural provinces in transferring plate margin deformation across the peninsula is enigmatic. Results We use Landsat Thematic Mapper imagery of the Baja California Peninsula to identify recent and active faults, and then conduct field studies that characterize the temporal and spatial structural evolution of the plate margin. These data address questions concerning the neotectonic development of the Gulf of California, the Baja California Peninsula, and their role in evolution of the post-Miocene Pacific - North American plate boundary. Moreover, these studies provide constraints on the geometry of active faults, allowing more exact understanding of the results of ongoing NASA-supported geodetic experiments. In addition, anticipated publication of the TM scenes will provide a widely available geological data base for relatively little-known peninsula California. Achievements include development of an ArcInfo data base of Landsat and SPOT imagery, detailed field studies of Neogene structures in northeastern Baja California, and new constraint on Pacific - North America plate motion at Baja California latitudes. These results are reported in maps, manuscripts and data products which are published or near completion.

Miller, M. Meghan↗

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Upper Bound on Orbital Debris Collision Probability When Only One Object has Position Uncertainty Information

Upper bounds on high speed satellite collision probability, P (sub c), have been investigated. Previous methods assume an individual position error covariance matrix is available for each object. The two matrices being combined into a single, relative position error covariance matrix. Components of the combined error covariance are then varied to obtain a maximum P (sub c). If error covariance information for only one of the two objects was available, either some default shape has been used or nothing could be done. An alternative is presented that uses the known covariance information along with a critical value of the missing covariance to obtain an approximate but useful P (sub c) upper bound. There are various avenues along which an upper bound on the high speed satellite collision probability has been pursued. Typically, for the collision plane representation of the high speed collision probability problem, the predicted miss position in the collision plane is assumed fixed. Then the shape (aspect ratio of ellipse), the size (scaling of standard deviations) or the orientation (rotation of ellipse principal axes) of the combined position error ellipse is varied to obtain a maximum P (sub c). Regardless as to the exact details of the approach, previously presented methods all assume that an individual position error covariance matrix is available for each object and the two are combined into a single, relative position error covariance matrix. This combined position error covariance matrix is then modified according to the chosen scheme to arrive at a maximum P (sub c). But what if error covariance information for one of the two objects is not available? When error covariance information for one of the objects is not available the analyst has commonly defaulted to the situation in which only the relative miss position and velocity are known without any corresponding state error covariance information. The various usual methods of finding a maximum P (sub c) do no good because the analyst defaults to no knowledge of the combined, relative position error covariance matrix. It is reasonable to think, given an assumption of no covariance information, an analyst might still attempt to determine the error covariance matrix that results in an upper bound on the P (sub c). Without some guidance on limits to the shape, size and orientation of the unknown covariance matrix, the limiting case is a degenerate ellipse lying along the relative miss vector in the collision plane. Unless the miss position is exceptionally large or the at-risk object is exceptionally small, this method results in a maximum P (sub c) too large to be of practical use. For example, assuming that the miss distance is equal to the current ISS alert volume along-track (+ or -) distance of 25 kilometers and that the at-risk area has a 70 meter radius. The maximum (degenerate ellipse) P (sub c) is about 0.00136. At 40 kilometers, the maximum P (sub c) would be 0.00085 which is still almost an order of magnitude larger than the ISS maneuver threshold of 0.0001. In fact, a miss distance of almost 340 kilometers is necessary to reduce the maximum P (sub c) associated with this degenerate ellipse to the ISS maneuver threshold value. Such a result is frequently of no practical value to the analyst. Some improvement may be made with respect to this problem by realizing that while the position error covariance matrix of one of the objects (usually the debris object) may not be known the position error covariance matrix of the other object (usually the asset) is almost always available. Making use of the position error covariance information for the one object provides an improvement in finding a maximum P (sub c) which, in some cases, may offer real utility. The equations to be used are presented and their use discussed.

Frisbee, Joseph H., Jr.↗

Visualizing Comparisons of Bill of Materials

Protecting critical infrastructure from cyber attacks, natural disasters, and other disruptions is a priority of the U.S. Government. Critical infrastructure includes providing electricity to homes and businesses, supplying natural gas for heating, and producing renewable energy sources. A loss of these services, as seen in the Solarwinds supply chain attack in 2020 , Texas snowstorm of 2021, the Colonial Pipeline cyber incident of 2021, and the Washington power substation attacks in 2022 result in high costs to consumers, disruption of everyday life, and even death. To protect the infrastructure, we first have to know what equipment we are protecting. The complexity of distributed manufacturing and development coupled with the increasing prevalence of cyber and supply chain attacks necessitates a greater understanding of the hardware and software components that comprise equipment in critical infrastructure. When a vulnerability in a single software library can have disastrous consequences, it is vital to understand critical equipment and systems at a granular level. This need has led to increased energy around the development and incorporation of bill-of-materials (BOM) into existing asset management practices to aid in mitigating, and responding to future attacks \cite{noauthor_software_nodate}. While much of the current research is devoted to creating BOMs, it is equally important to develop methodologies for leveraging BOMs to answer questions, such as: How has my software changed? Are two pieces of equipment equivalent? Does this piece of equipment that just arrived match my historical information? In this work, we demonstrate how BOMs can be represented by graph structures. We then describe how these structures can be fed into a graph comparison algorithm to produce a novel interactive visualization that allows us to not only identify differences in BOMs, but show exactly where they are in the product.

Jones, Rebecca D.↗

Comparison of V-4 and V-5 Exercise/Oxygen Prebreathe Protocols to Support Extravehicular Activity in Microgravity

The Prebreathe Reduction Program (PRP) used exercise during oxygen prebreathe to reduce necessary prebreathe time prior to depressurizing to work in a 4.3 psi suit during extravehicular activity (EVA). Initial testing produced a two-hour protocol incorporating ergometry exercise and a 30 min cycle of depress/repress to 10.2 psi where subjects breathed 26.5% oxygen/balance nitrogen (Phase II - 10 min at 75% peak oxygen consumption [VO2 peak] followed by 40 min intermittent light exercise [ILE] [approx. 5.8 mL-per kilogram- per minute], then 50 min of rest). The Phase II protocol (0/45 DCS) was approved for operations and has been used on 40 EVAs, providing significant time savings compared to the standard 4 h resting oxygen prebreathe. The Phase V effort focused on performing all light in-suit exercise. Two oxygen prebreathe protocols were tested sequentially: V-4) 160 min prebreathe with 150 min of continuous ILE. The entire protocol was completed at 14.7 psi. All exercise involved upper body effort. Exercise continued until decompression. V-5) 160 min prebreathe with 140 min of ILE - first 40 min at 14.7 psi, then 30 min at 10.2 psi (breathing 26.5% oxygen) after a 20 min depress, simulating a suit donning period. Subjects were then repressed to 14.7 psi and performed another 50 min of lower body ILE, followed by 50 min rest before decompression. The V-4 protocol was rejected with 3 DCS/6 person-exposures. Initial V-5 testing has produced 0 DCS/11 person-exposures (ongoing trials). The difference in DCS rate was significant (Fisher Exact p=0.029). The observations of DCS were significantly lower in early V-5 trials than in V-4 trials. Additional studies are required to evaluate the relative contribution of the variables in exercise distribution, the 10.2 psi depress/repress component, pre-decompression rest, or possible variation in total oxygen consumption.

Pollock, N. W.↗

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗

Active Shielding and Control of Environmental Noise

In the framework of the research project supported by NASA under grant # NAG-1-01064, we have studied the mathematical aspects of the problem of active control of sound, i.e., time-harmonic acoustic disturbances. The foundations of the methodology are described in our paper [1]. Unlike. many other existing techniques, the approach of [1] provides for the exact volumetric cancellation of the unwanted noise on a given predetermined region airspace, while leaving unaltered those components of the total acoustic field that are deemed as friendly. The key finding of the work is that for eliminating the unwanted component of the acoustic field in a given area, one needs to know relatively little; in particular, neither the locations nor structure nor strength of the exterior noise sources need to be known. Likewise, there is no need to know the volumetric properties of the supporting medium across which the acoustic signals propagate, except, maybe, in a narrow area of space near the perimeter of the protected region. The controls are built based solely on the measurements performed on the perimeter of the domain to be shielded; moreover, the controls themselves (i.e., additional sources) are concentrated also only on or near this perimeter. Perhaps as important, the measured quantities can refer to the total acoustic field rather than to its unwanted component only, and the methodology can automatically distinguish between the two. In [1], we have constructed the general solution for controls. The apparatus used for deriving this general solution is closely connected to the concepts of generalized potentials and boundary projections of Calderon's type. For a given total wave field, the application of a Calderon's projection allows one to definitively tell between its incoming and outgoing components with respect to a particular domain of interest, which may have arbitrary shape. Then, the controls are designed so that they suppress the incoming component for the domain to be shielded or alternatively, the outgoing component for the domain, which is complementary to the one to be shielded.

Tsynkov, S. V.↗

Stability of the particle-hole Pfaffian state and the 5 2 -fractional quantum Hall effect

Here we present a method for the exact construction of the fully particle-hole symmetric Pfaffian (PH-Pfaffian) ground state and its charged excitations on a sphere. We adopt the Moore-Read state, but with a nonholomorphic pairing component as in previous studies, and project it to the lowest Landau level. We study the energetics as well as other properties of these states and find that in a pure system interacting with the Coulomb forces the PH-Pfaffian cannot compete with either the Moore-Read state or its particle-hole conjugate, the anti-Pfaffian state, as an explanation for the 5/2 effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Toward a Unified View of Black-Hole High-Energy States

We present here a review of high-energy (greater than 1 keV) observations of seven black-hole candidates, six of which have estimated masses. In this review we focus on two parameters of interest: the ratio of 'nonthermal' to total luminosity as a function of the total luminosity divided by the Eddington luminosity, and the root-mean-square (rms) variability as a function of the nonthermal-to-total luminosity ratio. Below approx. 10% Eddington luminosity, the sources tend to be strictly nonthermal (the so called 'off' and 'low' states). Above this luminosity the sources become mostly thermal (the 'high' state). with the nonthermal component increasing with luminosity (the 'very high' and 'flare' states). There are important exceptions to this behavior, however, and no steady - as opposed to transient - source has been observed over a wide range of parameter space. In addition, the rms variability is positively correlated with the ratio of nonthermal to total luminosity, although there may be a minimum level of variability associated with 'thermal' states. We discuss these results in light of theoretical models and find that currently no single model describes the full range of black-hole high-energy behavior. In fact, the observations are exactly opposite from what one expects based upon simple notions of accretion disk instabilities.

Nowak, Michael A.↗

Static galactic halo and galactic wind

Although the exact state of the interstellar medium (ISM) in our Galaxy (other galaxies as well) is not clear at all, the 'common consensus' is that a rough pressure balance (or equipartition of energy) exists between different components and phases: cold, warm, hot phases of the ISM, magnetic field, cosmic rays, etc. If the halo of a galaxy is taken to be an extension of the ISM, then its structure is influenced by various ISM components. A 'complete' description of the halo is evidently very complicated. This paper gives a brief account on cosmic ray halo, which emphasizes the role played by cosmic rays. The interaction between cosmic rays and thermal plasma is facilitated by magnetic field. The cosmic rays are scattered by hydromagnetic waves (e.g., Alfven waves) which in turn can be generated by cosmic ray streaming instability. This constitutes a self-consistent picture. Since we are interested in the structure of the halo, we adopted a hydrodynamic model in which the cosmic rays and waves are described by their pressures. In general there are two classes of halos: static and dynamic.

Ko, Chung-Ming↗

Oxidizer heat exchanger component test

The RL10-IIB engine, is capable of multimode thrust operation. The engine operates at two low-thrust levels: tank head idle (THI), approximately 1 to 2 percent of full thrust; and pumped idle, 10 percent of full thrust. Operation at THI provides vehicle propellant settling thrust and efficient thermal conditioning; PI operation provides vehicle tank prepressurization and maneuver thrust for low-g deployment. Stable combustion of the RL10-IIB engine during the low-thrust operating modes can be accomplished by using a heat exchanger to supply gaseous oxygen to the propellant injector. The oxidized heat exchanger (OHE) vaporizes the liquid oxygen using hydrogen as the energy source. This report summarizes the test activity and post-test data analysis for two possible heat exchangers, each of which employs a completely different design philosophy. One design makes use of a low-heat transfer (PHT) approach in combination with a volume to attenuate pressure and flow oscillations. The test data showed that the LHT unit satisfied the oxygen exit quality of 0.95 or greater in both the THI and PI modes while maintaining stability. The HHT unit fulfilled all PI requirements; data for THI satisfactory operation is implied from experimental data that straddle the exact THI operating point.

Kanic, P. G.↗

Adaptive Nulling for Interferometric Detection of Planets

An adaptive-nulling method has been proposed to augment the nulling-optical- interferometry method of detection of Earth-like planets around distant stars. The method is intended to reduce the cost of building and aligning the highly precise optical components and assemblies needed for nulling. Typically, at the mid-infrared wavelengths used for detecting planets orbiting distant stars, a star is millions of times brighter than an Earth-sized planet. In order to directly detect the light from the planet, it is necessary to remove most of the light coming from the star. Nulling interferometry is one way to suppress the light from the star without appreciably suppressing the light from the planet. In nulling interferometry in its simplest form, one uses two nominally identical telescopes aimed in the same direction and separated laterally by a suitable distance. The light collected by the two telescopes is processed through optical trains and combined on a detector. The optical trains are designed such that the electric fields produced by an on-axis source (the star) are in anti-phase at the detector while the electric fields from the planet, which is slightly off-axis, combine in phase, so that the contrast ratio between the star and the planet is greatly decreased. If the electric fields from the star are exactly equal in amplitude and opposite in phase, then the star is effectively nulled out. Nulling is effective only if it is complete in the sense that it occurs simultaneously in both polarization states and at all wavelengths of interest. The need to ensure complete nulling translates to extremely tight demands upon the design and fabrication of the complex optical trains: The two telescopes must be highly symmetric, the reflectivities of the many mirrors in the telescopes and other optics must be carefully tailored, the optical coatings must be extremely uniform, sources of contamination must be minimized, optical surfaces must be nearly ideal, and alignments must be extremely precise. Satisfaction of all of these requirements entails substantial cost.

Lay, Oliver P.↗

Nontrivial fusion of Majorana zero modes in interacting quantum-dot arrays

Motivated by recent experimental reports of Majorana zero modes (MZMs) in quantum-dot systems at the “sweet spot,” where the electronic hopping t ℎ is equal to the superconducting coupling Δ, we study the time-dependent spectroscopy corresponding to the nontrivial fusion of MZMs. The term “nontrivial” refers to the fusion of Majoranas from different original pairs of MZMs, each with well-defined parities. We employ an experimentally accessible time-dependent real-space local density-of-states (LDOS) method to investigate the nontrivial MZM fusion outcomes in canonical chains and in a Y-shaped array of interacting electrons. In the case of quantum-dot chains where two pairs of MZMs are initially disconnected, after fusion we find equal-height peaks in the electron and hole components of the LDOS, signaling nontrivial fusion into both the vacuum I and fermion Ψ channels with equal weight. For π-junction quantum-dot chains, where the superconducting phase has opposite signs on the left and right portions of the chain, after the nontrivial fusion we observed the formation of an exotic two-site MZM near the center of the chain, coexisting with another single-site MZM. Furthermore, we also studied the fusion of three MZMs in the Y-shaped geometry. In this case, after the fusion we observed the novel formation of another exotic multisite MZM, with properties depending on the connection and geometry of the central region of the Y-shaped quantum-dot array.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reionization history constraints from neural network based predictions of high-redshift quasar continua

ABSTRACT Observations of the early Universe suggest that reionization was complete by z ∼ 6, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Lyα damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Lyα, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening intergalactic medium (IGM). In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Lyα. Our Quasar Spectra from Artificial Neural Network based predictive Regression Algorithm(QSANNdRA) improves the error in this reconstruction compared to the state-of-the-art principal component analysis (PCA) based model in the literature by 14.2 per cent on average, and provides an improvement of 6.1 per cent on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1 per cent and 16.8 per cent when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at z = 7.0851 and ULAS J1342+0928 at z = 7.5413, respectively. Using our more accurate reconstructions of these two z > 7 quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{H\, \small{I}} = 0.25^{+0.05}_{-0.05}$ at z = 7.0851 and $\bar{x}_\mathrm{H\, \small{I}} = 0.60^{+0.11}_{-0.11}$ at z = 7.5413. Our results are consistent with the literature and favour a rapid end to reionization.

Ďurovčíková, Dominika↗

Accurate interlaminar stress recovery from finite element analysis

The accuracy and robustness of a two-dimensional smoothing methodology is examined for the problem of recovering accurate interlaminar shear stress distributions in laminated composite and sandwich plates. The smoothing methodology is based on a variational formulation which combines discrete least-squares and penalty-constraint functionals in a single variational form. The smoothing analysis utilizes optimal strains computed at discrete locations in a finite element analysis. These discrete strain data are smoothed with a smoothing element discretization, producing superior accuracy strains and their first gradients. The approach enables the resulting smooth strain field to be practically C1-continuous throughout the domain of smoothing, exhibiting superconvergent properties of the smoothed quantity. The continuous strain gradients are also obtained directly from the solution. The recovered strain gradients are subsequently employed in the integration o equilibrium equations to obtain accurate interlaminar shear stresses. The problem is a simply-supported rectangular plate under a doubly sinusoidal load. The problem has an exact analytic solution which serves as a measure of goodness of the recovered interlaminar shear stresses. The method has the versatility of being applicable to the analysis of rather general and complex structures built of distinct components and materials, such as found in aircraft design. For these types of structures, the smoothing is achieved with 'patches', each patch covering the domain in which the smoothed quantity is physically continuous.

Tessler, Alexander↗

Radiative Impacts of Cloud Heterogeneity and Overlap in an Atmospheric General Circulation Model

The radiative impacts of introducing horizontal heterogeneity of layer cloud condensate, and vertical overlap of condensate and cloud fraction are examined with the aid of a new radiation package operating in the GEOS-5 Atmospheric General Circulation Model. The impacts are examined in terms of diagnostic top-of-the-atmosphere shortwave (SW) and longwave (LW) cloud radiative effect (CRE) calculations for a range of assumptions and parameter specifications about the overlap. The investigation is conducted for two distinct cloud schemes, the one that comes with the standard GEOS-5 distribution, and another which has been recently used experimentally for its enhanced GEOS-5 distribution, and another which has been recently used experimentally for its enhanced cloud microphysical capabilities; both are coupled to a cloud generator allowing arbitrary cloud overlap specification. We find that cloud overlap radiative impacts are significantly stronger for the operational cloud scheme for which a change of cloud fraction overlap from maximum-random to generalized results to global changes of SW and LW CRE of approximately 4 Watts per square meter, and zonal changes of up to approximately 10 Watts per square meter. This is because of fewer occurrences compared to the other scheme of large layer cloud fractions and of multi-layer situations with large numbers of atmospheric being simultaneously cloudy, conditions that make overlap details more important. The impact on CRE of the details of condensate distribution overlap is much weaker. Once generalized overlap is adopted, both cloud schemes are only modestly sensitive to the exact values of the overlap parameters. We also find that if one of the CRE components is overestimated and the other underestimated, both cannot be driven towards observed values by adjustments to cloud condensate heterogeneity and overlap alone.

Oreopoulos, L.↗

In Situ Resource Utilization (ISRU) Experiments for Mars Exploration

In situ resource utilization can best be described as living off the land. In our case the land is the planet Mars. ISRU is based on the idea that some fraction of the consumables, life support and propellant materials do not have to be flown from earth. Rather, they can be manufactured or extracted from resources already present on Mars. The primary resources on Mars are the atmosphere, polar caps and regolith. The atmosphere of Mars is mostly carbon dioxide as shown in the table below. The proportion of oxygen on the other hand is quite small. Still, there is quite a bit of oxygen in the Martian atmosphere, but it is unfortunately tied up with carbon. Thus, one of the goals of ISRU is the separation of breathable oxygen from the carbon dioxide. Several means of separation have been proposed. We have begun experiments on another approach for production of oxygen with carbon monoxide as a useful by product. Our work on a CO2 separator is described later in this report. Regolith melting is another means of obtaining materials. Two materials of interest are iron and silicon. Iron oxide is plentiful on Mars and is of obvious importance for structural components. Silicon is the foundation of solid state devices. Power generation on Mars may be accomplished using silicon solar cells. There is discussion of the feasibility of in situ production of solar cells. This would require a means of extracting silicon from the regolith. We have conducted several experiments concerning melting and glassification of the Mars soil simulant. Other summer faculty fellows have tried various means of processing the stimulant material. These include furnace melting, microwave melting and laser ablation. We have conducted several furnace melting experiments in both air and carbon dioxide environments. We have also carried out experiments to test spark melting in a carbon dioxide atmosphere. These experiments suggest the possibility of using arc melting in a reducing atmosphere. It is important to keep in mind that we are working with a soil stimulant. Any simulant, no matter how chemically similar it is to Martian regolith, may differ in mineralogy. The underlying assumption in this work is that once a glass is formed, any differences between simulant and regolith are unimportant. The exact means of forming the glass do, however, depend on the mineralogy of the regolith. A sample return mission is required to help answer these questions.

Marone, Matt↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗