SEARCH · Engineering Papers
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.
The motion of charged particles is studied numerically and analytically using a test particle approach. Acceleration near the X-type neutral line is described in terms of adiabatic invariants and separatrix crossings. The relation between regular and chaotic properties of the process is discussed for a realistic source in the plasma mantle.
A thick, quasi-parallel bow shock structure was observed with field and particle detectors of both HEOS 1 and OGO 5. The typical magnetic pulsation structure was at least 1 to 2 earth radii thick radially and was accompanied by irregular but distinct plasma distributions characteristic of neither the solar wind nor the magnetosheath. Waves constituting the large pulsations were polarized principally in the plane of the nominal shock, therefore also in the plane perpendicular to the average interplanetary field. A separate interpulsation regime detected between bursts of large amplitude oscillations was similar to the upstream wave region magnetically, but was characterized by disturbed plasma flux and enhanced noise around the ion plasma frequency. The shock structure appeared to be largely of an oblique, whistler type, probably complicated by counterstreaming high energy protons. Evidence for firehose instability-based structure was weak at best and probably negative.
The field of elementary particle physics explores the universe at the smallest and largest scales to study the fundamental constituents of matter and energy. The United States plays the leading role in this global effort, including some of the world's most extensive scientific collaborations and most complex machines. Current U.S. investments in particle physics are guided by the strategy and priorities identified in the 2014 report of the Particle Physics Project Prioritization Panel (P5), the culmination of a multi-year community-driven process.
The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.
In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).
The fluctuations in magnetic field and plasma velocity in solar wind, which possess many features of fully developed magnetohydrodynamic (MHD) turbulence, are discussed. Direct spacecraft observations from 0.3 to over 20 AU, remote sensing radio scintillation observations, numerical simulations, and various models provide complementary methods that show that the fluctuations in the wind parameters undergo significant dynamical evolution independent of whatever turbulence might exist in the solar photosphere and corona. The Cluster mission, with high time resolution particle and field measurements and its variable separation strategies, should be able to provide data for answering many questions on MHD turbulence.
Despite being a quantum two-fluid system, superfluid helium-4 (He II) is observed to behave similarly to classical fluids when a flow is generated by mechanical forcing. This similarity has brought up the feasibility of utilizing He II for high Reynolds number classical turbulence research, considering the small kinematic viscosity of He II. However, it has been suggested that the nonclassical dissipation mechanism in He II at small scales may alter its turbulent statistics and intermittency. In this work, we report our study of a nearly homogeneous and isotropic turbulence (HIT) generated by a towed grid in He II. We measure the velocity field using particle tracking velocimetry with solidified deuterium particles as the tracers. By correlating the velocities measured simultaneously on different particle trajectories or at different times along the same particle trajectory, we are able to conduct both Eulerian and Lagrangian flow analyses. Spatial velocity structure functions obtained through the Eulerian analysis show scaling behaviors in the inertial subrange similar to that for classical HIT but with enhanced intermittency. The Lagrangian analysis allows us to examine the flow statistics down to below the dissipation length scale. Interestingly, abnormal deviations from the classical scaling behaviors are observed in this regime. Lastly, we discuss how these deviations may relate to the motion of quantized vortices in the superfluid component in He II.
The Alfven critical velocity mechanism for ionization of a neutral gas streaming across the magnetic field has been demonstrated in laboratory experiments. In March 1983, two rocket-borne experiments with Ba and Sr tested the effect in the wall-less laboratory of space from Punto Lobos, Peru, near 430 km altitude. 'Star of Lima' used a conical Ba shaped charge aimed at an instrument payload about 2 km away. Because of rocket overperformance the detonation occurred in partial sunlight, so that less than 21.6 percent of the ionizing UV was present. Particle and field measurements indicate the production of hot electrons and waves in the energy and frequency range that are respectively predicted to produce a cascade of ionization by the Alfven mechanism. However, the ionization fluxes and wave energy density did not reach cascade levels, and optical observations indicate that only 2.5 to 5 x 10 to the 20th Ba ions were produced. A substantial portion and perhaps all of the ionization could have been produced by solar UV. The failure of the Alfven process in this experiment is not well understood.
Alpha particles are known to be a major source of particles creating soft errors in semiconductor devices, such as content flipping in Static Random-Access Memory (SRAM). Recent advancements in transistor nodes have led to the introduction of Gate-All-Around Field Effect Transistors (GAA-FETs), which have better gate control, thus better electrostatics. Moreover, the introduction of bottom dielectric isolation (BDI) eliminates substrate leakage and thus is expected to enhance its radiation hardness. It is thus important to explore if one can design an SRAM that is completely radiation-hard to alpha particles. In this paper, using 3D Technology Computer-Aided-Design (TCAD) simulations, we show that it is possible to design an SRAM using GAA-FET technology so that it is immune to single alpha particle radiation error. In other words, with the design, there will be no single-event upset (SEU) due to alpha particles. We first use ab initio calculations in PHITS to show that there is a maximum linear energy transfer (LET), LET max , for the alpha particle in Si and Si x Ge 1-x . Based on that, by de signing a sub-7nm GAA-FET-based SRAM with BDI, we show that the SRAM does not flip even if the particle strike is in the worst-case scenario for LET > LET max .
The nonlinear coalescence instability of current carrying solar loops can explain many of the characteristics of the solar flares such as their impulsive nature, heating and high energy particle acceleration, amplitude oscillations of electromagnetic emission as well as the characteristics of 2-D microwave images obtained during a solar flare. The physical characteristics of the explosive coalescence of currents are presented in detail through computer simulation and theory. Canonical characteristics of the explosive coalescence are: (1) a large amount of impulsive increase of kinetic energies of electrons and ions; (2) simultaneous heating and acceleration of electrons and ions in high and low energy spectra; (3) ensuing quasi-periodic amplitude oscillations in fields and particle quantities; and (4) the double peak (or triple peak) structure in these profiles, participate in the coalescence process, yielding varieties of phenomena.
The study of relativistic electron–positron pair plasmas is both of fundamental physics interest and important to understand the processes that shape the magnetic field dynamics, particle acceleration, and radiation emission in high-energy astrophysical environments. Although it is highly desirable to study relativistic pair plasmas in the laboratory, their generation and control constitutes a critical challenge. Significant experimental and theoretical progress has been made over recent years to explore the use of intense lasers to produce dense relativistic pair plasma in the laboratory and study the basic collective plasma processes associated with these systems. Important challenges remain in terms of improving the number of pairs, system size, and control over the charge neutrality required to establish laboratory platforms that can expand our understanding of relativistic pair plasma and help validate underlying models in conditions relevant to high-energy astrophysical phenomena. Furthermore, we highlight recent progress in this field, discuss the main challenges, and the exciting prospects for studying relativistic pair plasmas and astrophysics relevant instabilities in the laboratory in the near future.
Abstract Particle‐resolved direct numerical simulations (PR‐DNS) play an increasing role in investigating aerosol‐cloud‐turbulence interactions at the most fundamental level of processes. However, the high computational cost associated with high resolution simulations poses considerable challenges for large domain or long duration simulation using PR‐DNS. To address these issues, here we present an emulator of the complex physics‐based PR‐DNS developed by use of the data‐driven Fourier Neural Operator (FNO) method. The effectiveness of the method is showcased by presenting turbulence and temperature fields in a two‐dimensional space. The results demonstrate high accuracy at various resolutions and the emulator is two orders of magnitude cheaper in terms of computational demand compared to the physics‐based PR‐DNS model. Furthermore, the FNO emulator exhibits strong generalization capabilities for different initial conditions and ultra‐high‐resolution without the need to retrain models. These findings highlight the potential of the FNO method as a promising tool to simulate complex fluid dynamics problems with high accuracy, computational efficiency, and generalization capabilities, enhancing our understanding of the aerosol‐cloud‐precipitation system.
Solar-terrestrial physics concerns the study of the generation, flow, and dissipation of mass, momentum, and energy between the Sun and the Earth. Mass, momentum, and energy are carried by charged particles that compose the solar wind. When the solar wind reaches the Earth, some solar-wind particles enter the magnetosphere; this coupling between the solar wind and the Earth means that the solar wind can influence the Earth's upper atmosphere. As the first step in addressing the behavior of this solar-terrestrial system, the Global Geospace Science (GGS) Initiative will use the Wind and Polar satellites, provided by NASA, and the Geotail satellite provided by the Japanese Institute of Space and Astronautical Science (ISAS), to perform simultaneous and closely coordinated measurements of the key geospace regions and will add data from equatorial missions. Magnetic field and particle changes that occur when particles are energized during auroral events will be monitored. The intention behind the GGS Initiative is to understand the physical mechanisms and various regions controlling the transport of mass, momentum, and energy in geospace. A summary of the GGS Initiative is presented.
The Martian magnetotail is a complex regime through which atmospheric particles are lost to space. Our current understanding of Mars' tail continues to develop with the comprehensive particle and field data collected by Mars Atmosphere and Volatile EvolutioN (MAVEN). In this work, we identify periods when MAVEN encounters multiple current sheet crossings through a single tail traversal in order to understand tail dynamics. We apply an analysis technique that has been developed and validated by using multipoint measurements in order to separate the spatial and temporal properties associated with current sheet flapping. Events are classified into periods of steady flapping, due to a global motion of the current sheet, and kink-like flapping, resulting from localized wave propagation along the tail current sheet. Out of 106 periods during which multiple current sheet crossings were observed, 20 were due to steady flapping and 10 from kink-like flapping. A majority of the kink-like events resulted from waves propagating in the opposite direction of the solar wind convection electric field, regardless of their location in the tail, unlike at Earth and Venus. This finding suggests that possible magnetosphere energy sources, whereby plasma is accelerated and removed from the Martian environment, are not located in the central magnetotail; rather, these waves may be driven by a source located at the tail flank based on the direction of the solar wind electric field. Therefore, by identifying potential sources of impulsive energy release in the tail, we may better understand mechanisms that drive atmospheric loss at Mars.
The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.
Scheduled to launch in May 2012, NASA's dual spacecraft Living With a Star Radiation Belt Storm Probe mission carries the field and particle instrumentation needed to determine the processes that produce enhancements in radiation belt ion and electron fluxes, the dominant mechanisms that cause the loss of relativistic electrons, and the manner by which the ring current and other geomagnetic phenomena affect radiation belt behavior. The two spacecraft will operate in low-inclination elliptical lapping orbits around the Earth, within and immediately exterior to the Van Allen radiation belts. During course of their two year primary mission, they will cover the full range of local times, measuring both AC and DC electric and magnetic fields to 10kHz, as well as ions from 50 eV to 1 GeV and electrons with energies ranging from 50 eV to 10 MeV.
Non-invertible global symmetry often predicts degeneracy in axion potentials and carries important information about the global form of the gauge group. When these symmetries are spontaneously broken they can lead to the formation of stable axion domain wall networks which support topological degrees of freedom on their worldvolume. Such non-invertible symmetries can be broken by embedding into appropriate larger UV gauge groups where small instanton contributions lift the vacuum degeneracy, and provide a possible solution to the domain wall problem. We explain these ideas in simple illustrative examples and then apply them to the Standard Model, whose gauge algebra and matter content are consistent with several possible global structures. Each possible global structure leads to different selection rules on the axion couplings, and various UV completions of the Standard Model lead to more specific relations. As a proof of principle, we also present an example of a UV embedding of the Standard Model which can solve the axion domain wall problem. The formation and annihilation of the long-lived axion domain walls can lead to observables, such as gravitational wave signals. Observing such signals, in combination with the axion coupling measurements, can provide valuable insight into the global structure of the Standard Model, as well as its UV completion.
The shape of the dayside Venus ionopause, and its dependence on solar wind parameters, is examined using Pioneer Venus Orbiter field and particle data. The ionopause is defined here as the altitude of pressure equality between magnetosheath pressure and ionospheric thermal pressure; its typical altitudes range from about 300 km near the subsolar point to about 900 km near the terminator. A strong correlation between ionopause altitude and magnetosheath magnetic pressure is demonstrated; correlation between magnetic pressure and the normally incident component of solar wind dynamic pressure is also evident. The data support the hypothesis of control of the ionopause altitude by solar wind dynamic pressure, manifested in the sheath as magnetic pressure. The presence of large scale magnetic fields in the ionosphere is observed primarily when dynamic pressure is high and the ionopause is low.