Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field model”

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 631 records · Page 35

Particle orbits in model current sheet with a nonzero B(y) component

The problem of charged particle motions in magnetotaillike model current sheets is revisited with the inclusion of a nonzero dawn-dusk magnetic field component. Three cases are examined considering both trapped and escaped orbits. The results show that a nonzero B(y) component disturbs the particle orbits by destroying orbit symmetry in the phase space about the z = 0 plane. It also changes the bounce frequency of particle orbits. The presence of B(y) thus modifies the Speiser orbits, particularly near the ejection phase. The process of ejected particle such as ejection direction, ejection velocity, and pitch angles are shown to depend on the sign of the charge.

Zhu, Zhongwei↗

Charged dust in planetary magnetospheres: Hamiltonian dynamics and numerical simulations for highly charged grains

We use a combination of analytical and numerical methods to investigate the dynamics of charged dust grains in planetary magnetospheres. Our emphasis is on obtaining results valid for particles that are not necessarily dominated either by gravitational or electromagnetic forces. A Hamiltonian formulation of the problem yields exact results, for all values of charge-to-mass ratio, when we introduce two constraints: particles remain in the equatorial plane and the magnetic field is taken as axially symmetric. In particular, we obtain locations of equilibrium points, the frequencies of stable periodic orbits, the topology of separatrices in phase space, and the rate of longitudinal drift. These results are significant for specific applications: motion in the nearly aligned dipolar field of Saturn, and the trajectories of arbitrarily charged particles in complex magnetic fields for limited periods of time after ejection from parent bodies. Since the model is restrictive, we also use numerical integrations of the full three-dimensional equations of motion and illustrate under what conditions the constrained problem yields reasonable results. We show that a large fraction of the intermediately charged and highly charged (gyrating) particles will always be lost to a planet's atmosphere within a few hundred hours, for motion through tilted-dipole magnetic fields. We find that grains must have a very high charge-to-mass ratio in order to be mirrored back to the ring plane. Thus, except perhaps at Saturn where the dipole tilt is very small, the likely inhabitants of the dusty ring systems are those particles that are either nearly Keplerian (weakly charged) grains or grains whose charges place them in the lower end of the intermediate charge zone. Fianlly, we demonstrate the effect of plasma drag on the orbits of gyrating particles to be a rapid decrease in gyroradius followed by a slow radial evolution of the guiding center.

Schaffer, L.↗

Topological Evolution of a Fast Magnetic Breakout CME in 3-Dimensions

W present the extension of the magnetic breakout model for CME initiation to a fully 3-dimensional, spherical geometry. Given the increased complexity of the dynamic magnetic field interactions in 3-dimensions, we first present a summary of the well known axisymmetric breakout scenario in terms of the topological evolution associated with the various phases of the eruptive process. In this context, we discuss the completely analogous topological evolution during the magnetic breakout CME initiation process in the simplest 3-dimensional multipolar system. We show that an extended bipolar active region embedded in an oppositely directed background dipole field has all the necessary topological features required for magnetic breakout, i.e. a fan separatrix surface between the two distinct flux systems, a pair of spine fieldlines, and a true 3-dimensional coronal null point at their intersection. We then present the results of a numerical MHD simulation of this 3-dimensional system where boundary shearing flows introduce free magnetic energy, eventually leading to a fast magnetic breakout CME. The eruptive flare reconnection facilitates the rapid conversion of this stored free magnetic energy into kinetic energy and the associated acceleration causes the erupting field and plasma structure to reach an asymptotic eruption velocity of greater than or approx. equal to 1100 km/s over an approx.15 minute time period. The simulation results are discussed using the topological insight developed to interpret the various phases of the eruption and the complex, dynamic, and interacting magnetic field structures.

Lynch, B. J.↗

Simulation Assisted Risk Assessment Applied to Launch Vehicle Conceptual Design

A simulation-based risk assessment approach is presented and is applied to the analysis of abort during the ascent phase of a space exploration mission. The approach utilizes groupings of launch vehicle failures, referred to as failure bins, which are mapped to corresponding failure environments. Physical models are used to characterize the failure environments in terms of the risk due to blast overpressure, resulting debris field, and the thermal radiation due to a fireball. The resulting risk to the crew is dynamically modeled by combining the likelihood of each failure, the severity of the failure environments as a function of initiator and time of the failure, the robustness of the crew module, and the warning time available due to early detection. The approach is shown to support the launch vehicle design process by characterizing the risk drivers and identifying regions where failure detection would significantly reduce the risk to the crew.

Mathias, Donovan L.↗

Multifunctional Collaborative Modeling and Analysis Methods in Engineering Science

Engineers are challenged to produce better designs in less time and for less cost. Hence, to investigate novel and revolutionary design concepts, accurate, high-fidelity results must be assimilated rapidly into the design, analysis, and simulation process. This assimilation should consider diverse mathematical modeling and multi-discipline interactions necessitated by concepts exploiting advanced materials and structures. Integrated high-fidelity methods with diverse engineering applications provide the enabling technologies to assimilate these high-fidelity, multi-disciplinary results rapidly at an early stage in the design. These integrated methods must be multifunctional, collaborative, and applicable to the general field of engineering science and mechanics. Multifunctional methodologies and analysis procedures are formulated for interfacing diverse subdomain idealizations including multi-fidelity modeling methods and multi-discipline analysis methods. These methods, based on the method of weighted residuals, ensure accurate compatibility of primary and secondary variables across the subdomain interfaces. Methods are developed using diverse mathematical modeling (i.e., finite difference and finite element methods) and multi-fidelity modeling among the subdomains. Several benchmark scalar-field and vector-field problems in engineering science are presented with extensions to multidisciplinary problems. Results for all problems presented are in overall good agreement with the exact analytical solution or the reference numerical solution. Based on the results, the integrated modeling approach using the finite element method for multi-fidelity discretization among the subdomains is identified as most robust. The multiple-method approach is advantageous when interfacing diverse disciplines in which each of the method's strengths are utilized. The multifunctional methodology presented provides an effective mechanism by which domains with diverse idealizations are interfaced. This capability rapidly provides the high-fidelity results needed in the early design phase. Moreover, the capability is applicable to the general field of engineering science and mechanics. Hence, it provides a collaborative capability that accounts for interactions among engineering analysis methods.

Ransom, Jonathan B.↗

Pulsating type IV solar radio bursts

Several models for pulsating type IV radio bursts are presented based on the assumption that the pulsations are the result of fluctuations in the synchrotron emission due to small variations in the magnetic field of the source. It is shown that a source that is optically thick at low frequencies due to synchrotron self-absorption exhibits pulsations that occur in two bands situated on either side of the spectral peak. The pulsations in the two bands are 180 deg out of phase and the band of pulsations at the higher frequencies is the more intense. In contrast, a synchrotron source that is optically thin at all frequencies and whose low-frequency emission is suppressed due to the Razin effect develops only a single band of pulsations around the frequency of maximum emission. However, the flux density associated with the latter model would be too small to explain the more intense pulsations that have been observed, unless the source area is considerably larger than presently seems reasonable.

Gotwols, B. L.↗

Formation of ferrobasalt at east Pacific midocean spreading centers

A numerical modeling using a least squares petrologic mixing program demonstrates that ridge ferrobasalts from three east Pacific regions are generated from normal midocean ridge balsalt by crystal fractionation of plagioclase, clinopyroxene, and minor olivine at shallow mantle depths. Clinopyroxene is an important fractionating phase in this model in contrast to most previously proposed fractionation models. Calculations demonstrate that the ferrobasalts represent a much smaller percentage of residual liquid than has previously been proposed. Implications of extensive shallow crystal fractionation for the determination of depths of magma generation are discussed. The different fields on an anorthite-olivine-hypersthene ternary diagram for basalts from east Pacific fast and moderate spreading ridges can be explained by variable degrees of shallow crystal fractionation rather than by different depths of magma generation as suggested by Bass (1971) and Scheidegger (1973).

Clague, D. A.↗

Geological mapping in the Amazon Jungle: A challenge to side-looking radar

A brief outline of the physical environment and vegetation characteristics peculiar to the Amazon Jungle is given to provide the background for the discussion of energy target interaction models. Examples of radar images of the central part of the jungle are used to illustrate a common approach to geologic mapping in this environment. During the initial phase of interpretation, the interaction phenomenon is evaluated and textural and tonal characteristics of the radar image are observed. Field checking of selected areas is the next phase of a geologic mapping program. Reinterpretation of radar images and auxiliary data (aerial photographs, Landsat images, and field data) is the last phase of geologic mapping.

Correa, A. C.↗

Modeling of convection phenomena in Bridgman-Stockbarger crystal growth

Thermal convection phenomena in a vertically oriented Bridgman-Stockbarger apparatus were modeled by computer simulations for different gravity conditions, ranging from earth conditions to extremely low gravity, approximate space conditions. The modeling results were obtained by the application of a state-of-the art, transient, multi-dimensional, completely densimetrically coupled, discrete-element computational model which was specifically developed for the simulation of flow, temperature, and species concentration conditions in two-phase (solid-liquid) systems. The computational model was applied to the simulation of the flow and the thermal conditions associated with the convection phenomena in a modified Germanium-Silicon charge enclosed in a stationary fused-silica ampoule. The results clearly indicated that the gravitational field strength influences the characteristics of the coherent vortical flow patterns, interface shape and position, maximum melt velocity, and interfacial normal temperature gradient.

Carlson, F. M.↗

Mesoscale and large-scale variability of the Antarctic circumpolar current

This investigation of the physical oceanography of the Southern Ocean will carry out two parallel efforts during the years preceding the launch of TOPEX/POSEIDON. First, the Geosat data will be used to develop a preliminary descriptive picture of the mesoscale and large-scale, low-frequency surface circulation of the Southern Ocean. Some of this analysis of Geosat data has already begun. For example, as a measure of the geographical distribution of mesoscale variability, a color-coded map of the standard deviation of sea level from two years of Geosat data is shown. Efforts are presently under way to investigate the seasonal and year-to-year variability of this mesoscale energy. The data are also being used to generate low-pass filtered fields of sea level from which the temporal evolution of large-scale variability in the Southern Ocean may be investigated. The second parallel effort is the development and test of modeling and data assimilation techniques that will later be applied to TOPEX/POSEIDON data during the postlaunch phase. One objective of the modeling and data assimilation is to investigate the relation between mesoscale sea level variations and eddy flux in the Southern Ocean. Uncertainties in present estimates of the various components of meridional oceanic heat transport are large. The evidence presented indicates that the very energetic mesoscale variability in the ACC apparently accounts for much of the estimated 0.45x10(exp 15) watts of poleward heat transport across the ACC required to balance the heat budget. Eddy variability is strongly coherent vertically in the ACC, at least in the vicinity of Drake Passage where nearly all of the historical in situ data have been collected.

Chelton, Dudley B.↗

Partial melting of ordinary chondrites: Lost City (H) and St. Severin (LL)

Eucrites and diogenites are examples of asteroidal basalts and orthopyroxenites, respectively. As they are found intermingled in howardites, which are inferred to be regolith breccias, eucrites and diogenites are thought to be genetically related. But the details of this relationship and of their individual origins remain controversial. Work by Jurewicz et al. showed that 1170-1180 C partial melts of the (anhydrous) Murchison (CM) chondrite have major element compositions extremely similar to primitive eucrites, such as Sioux County. However, the MnO contents of these melts were about half that of Sioux County, a problem for the simple partial melting model. In addition, partial melting of Murchison could not produce diogenites, because residual pyroxenes in the Murchison experiments were too Fe- and Ca-rich and were minor phases at all but the lowest temperatures. A parent magma for diogenites needs an expanded low-calcium pyroxene field. In their partial melting study of an L6 chondrite, Kushiro and Mysen found that ordinary chondrites did have an expanded low-Ca pyroxene field over that of CV chondrites (i.e., Allende), probably because ordinary chondrites have lower Mg/Si ratios. This study expands that of both Kushiro and Mysen and Jurewicz et al. to the Lost City (H) and St. Severin (LL) chondrites at temperatures ranging from 1170 to 1325 C, at an fO2 of one log unit below the iron-wuestite buffer (IW-1).

Jurewicz, Amy J. G.↗

Circinus X-1 revisited: Fast-timing properties in relation to spectral state

We have studied the X-ray spectral and fast-timing variations of Cir X-1 by performing a homogenous analysis of all EXOSAT ME data on this source using X-ray hardness-intensity diagrams (HIDs), color-color diagrams (CDs), and power spectra. Cir X-1 exhibits a wide range of power spectral shapes and a large variety in X-ray spectral shapes. At different epochs the power spectra variously resemble those of an atoll source, a Z source, a black-hole candidate, or are unlike any of these. At some epochs one-dimensional connected-branch patterns are seen in HID and CD, and at other times more complex structures are found. We interpret the complex behavior of Cir X-1 in terms of a model where accretion rate, orbital phase and epoch are the main determinants of the source behavior, and where the unique properties of the source are due to two special circumstances: (1) the source is the only known atoll source (accreting neutron star with a very low magnetic field) that can reach the Eddington critical accretion rate, and (2) it has a unique, highly eccentric and probably precessing orbit. Property (1) makes Cir X-1 a very important source for our understanding of the similarities in the observable properties of neutron stars and black holes as it allows to separate out black hole signatures from properties that are merely due to the presence of accretion compact with a low magnetic field.

Oosterbroek, T.↗

Development and Validation of a Supersonic Helium-Air Coannular Jet Facility

Data are acquired in a simple coannular He/air supersonic jet suitable for validation of CFD (Computational Fluid Dynamics) codes for high speed propulsion. Helium is employed as a non-reacting hydrogen fuel simulant, constituting the core of the coannular flow while the coflow is composed of air. The mixing layer interface between the two flows in the near field and the plume region which develops further downstream constitute the primary regions of interest, similar to those present in all hypersonic air breathing propulsion systems. A computational code has been implemented from the experiment's inception, serving as a tool for model design during the development phase.

Carty, Atherton A.↗

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗

Deep Learning Models for Planetary Seismicity Detection

Research in planetary seismology is fundamentally constrained by a lack of data. Seismo-logical science products of future missions can typically only be informed by theoretical signal/noise characteristics of the environment or likely Earth-analogues. Although objectives can be re-assessed after some initial data-collection upon lander arrival, transfer of high-resolution data back to Earth is costly on lander power usage. Over the last several years, development of GPU computing techniques and open-source high-level APIs have led to rapid advances in deep learning within the fields of computer vision, natural language processing, and collaborative filtering. These techniques are actively being adapted in seismology for a variety of tasks, including: earthquake detection, seismic phase discrimination, and ground-motion prediction. Until the recent detection of mars quakes during the Mars InSight mission, the only other measurements of seismicity recorded outside of Earth was on the Moon during the Apollo missions between 1969 to 1977. These unique data sets have been periodically revisited using new seismological methods, including ambient noise interferometry and Hidden Markov Models. Our objective is to develop a deep learning seismic detector and use it to catalog moonquakes from the Apollo 17 Lunar Seismic Profiling Experiment (LSPE) and compare the results with those obtained by other methods. Additionally, we will assess the accuracy tradeoff between using a training set of lunar data and one composed of Earth seismicity. In this document, we present preliminary results using a prototype classifier trained on a small set of earthquakes that was able to obtain detections for LSPE moonquakes with a greater accuracy than a recent study using Hidden Markov Models.

Civilini, F.↗

Testing and Application of Commercial Thermodynamics Software for Solid-Gas Equilibria in Planetary Science Applications

The ability to calculate chemical equilibria from thermochemical constants (e.g., H, S, Cp) or other experimental data is a tool widely employed by scientists. However, phase equilibria for planetary science problems is limited by easy access to libraries of thermochemical data, as well as the software to invert that data. Many academic software packages are focused on specific types of problems (e.g., aqueous solutions or melting and crystallization of magmas). This focus makes them excellent tools for specific uses, but poorer tools for other scenarios, especially those involving elements or chemical species that are less well-studied in that specific scientific field (e.g. NaCl vapor). Amongst the more general thermodynamic software tools, there are closed academic software packages, open-source options, as well as commercial software. We have explored the extent to which commercial chemistry/chemical engineering software may be able to satisfy the need for thermochemical modeling packages that can accommodate the diverse species and conditions of planetary science problems. Here we compare the output of software package HSC Chemistry™ (Metso Outotec) to published vapor-solid phase diagrams for a variety of major and minor elements. HSC Chemistry™ has been used in several recent geochemical studies of terrestrial volcanic systems [2- 4]. However, it has not to our knowledge been benchmarked against peer-reviewed model results for the conditions and compositions relevant to planetary science. We present one such comparison here, the low- pressure environment of the solar nebula, using literature data. In short, we find excellent agreement between HSC Chemistry™ and the well-regarded academic software package CONDOR. These preliminary results suggest that at least some commercial software packages are capable of robust thermochemical calculations for planetary science.

planetary science↗

Gamma-ray pulsars: Radiation processes in the outer magnetosphere

We describe an emission model for gamma ray pulsars based on curvature radiation-reaction limited charges in the outer magnetosphere. We show how pair production on thermal surface flux can limit the acceleration zones. Estimates for the efficiency of GeV photon production eta gamma and the gamma-ray beaming fraction are derived, including their dependence on pulsar parameters. In general eta gamma increases with pulsar age, but is decreased for low magnetic fields and for small magnetic inclinations. We argue that this produces GeV pulse profiles, curvature spectra and detection statistics consistent with the observations. We also describe the optical through X-ray pulsar synchrotron spectrum and the spectral variations with pulsar phase. A test computation for Vela-like parameters reproduces phase-resolved GeV spectra consistent with those observed by EGRET. Finally we comment on very high energy pulsed emission and particle production and note extensions needed to allow a more complete pulsar model.

Romani, Roger W.↗

Low frequency dielectric relaxation in boracites

In order to elucidate the nature of the imperfections which adversely affect pyroelectric processes in boracites, the thermal and dispersive characteristics of the low frequency dielectric response in iron-iodide boracite (Fe3B7O13I) and copper-chloride boracite (Cu3B7O13Cl) have been investigated. These characteristics were measured as a function of crystallographic orientation and applied field in both the ferroelectric and paraelectric states. The low frequency dielectric relaxation of 100 line oriented multi-domain copper-chloride boracite clearly indicates the dipole nature of the lattice imperfections. The activation energies calculated from a noninteracting Debye model, are 0.53 eV in the ferroelectric phase and 0.10 eV in the paraelectric phase.

Kim, Q.↗