Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable physics”

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

Turbulence and the Stabilization Principle

Further results of research, reported in several previous NASA Tech Briefs articles, were obtained on a mathematical formalism for postinstability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence). To recapitulate: Fictitious control forces are introduced to couple the dynamical equations with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in ordinary perceived three-dimensional space is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. Con sequently, the postinstability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable. The previously reported findings are analyzed from the perspective of the authors Stabilization Principle, according to which (1) stability is recognized as an attribute of mathematical formalism rather than of underlying physics and (2) a dynamical system that appears unstable when modeled by differentiable functions only can be rendered stable by modifying the dynamical equations to incorporate intrinsic stochasticity.

Zak, Michail↗

Investigation of Global Particulate Nitrate from the AeroCom Phase III Experiment

An assessment of global particulate nitrate and ammonium aerosol based on simulations from nine models participating in the Aerosol Comparisons between Observations and Models (AeroCom) phase III study is presented. A budget analysis was conducted to understand the typical magnitude, distribution, and diversity of the aerosols and their precursors among the models. To gain confidence regarding model performance, the model results were evaluated with various observations globally, including ground station measurements over North America, Europe, and east Asia for tracer concentrations and dry and wet depositions, as well as with aircraft measurements in the Northern Hemisphere mid-to-high latitudes for tracer vertical distributions. Given the unique chemical and physical features of the nitrate occurrence, we further investigated the similarity and differentiation among the models by examining (1) the pH-dependent NH3 wet deposition; (2) the nitrate formation via heterogeneous chemistry on the surface of dust and sea salt particles or thermodynamic equilibrium calculation including dust and sea salt ions; and (3) the nitrate coarse-mode fraction (i.e., coarse/total). It is found that HNO3, which is simulated explicitly based on full O3-HOx-NOx-aerosol chemistry by all models, differs by up to a factor of 9 among the models in its global tropospheric burden. This partially contributes to a large difference in NO3(-), whose atmospheric burden differs by up to a factor of 13. The atmospheric burdens of NH3 and NHC 4 differ by 17 and 4, respectively. Analyses at the process level show that the large diversity in atmospheric burdens of NO3(-), NH3, and NHC4(+) is also related to deposition processes. Wet deposition seems to be the dominant process in determining the diversity in NH3 and NHC 4 lifetimes. It is critical to correctly account for contributions of heterogeneous chemical production of nitrate on dust and sea salt, because this process overwhelmingly controls atmospheric nitrate production (typically greater than 80 %) and determines the coarse- and fine-mode distribution of nitrate aerosol.

Ammonium aerosol based on simulations↗

An analytical transformation technique for generating uniformly spaced computational mesh

An analytical transformation method which can map arbitrary physical coordinate grid distribution into desired computational coordinate with uniform grid distribution is derived. The transformation function and its higher derivatives are differentiable. Salient features include; (1) precise control of grid sizes; (2) more than one location of clustered grids; (3) exact positioning of particular computational nodes in the physical plane; and (4) ensuring several patches of uniformly spaced grids in the physical plane for the higher accuracies (such as at the boundaries), etc., while keeping the variation of grid spacing continuous to avoid numerical instability.

Oh, Y. H.↗

On determining important aspects of mathematical models: Application to problems in physics and chemistry

The use of parametric and functional gradient sensitivity analysis techniques is considered for models described by partial differential equations. By interchanging appropriate dependent and independent variables, questions of inverse sensitivity may be addressed to gain insight into the inversion of observational data for parameter and function identification in mathematical models. It may be argued that the presence of a subset of dominantly strong coupled dependent variables will result in the overall system sensitivity behavior collapsing into a simple set of scaling and self similarity relations amongst elements of the entire matrix of sensitivity coefficients. These general tools are generic in nature, but herein their application to problems arising in selected areas of physics and chemistry is presented.

Rabitz, Herschel↗

Simulation of Ultra-Small MOSFETs Using a 2-D Quantum-Corrected Drift-Diffusion Model

The continued down-scaling of electronic devices, in particular the commercially dominant MOSFET, will force a fundamental change in the process of new electronics technology development in the next five to ten years. The cost of developing new technology generations is soaring along with the price of new fabrication facilities, even as competitive pressure intensifies to bring this new technology to market faster than ever before. To reduce cost and time to market, device simulation must become a more fundamental, indeed dominant, part of the technology development cycle. In order to produce these benefits, simulation accuracy must improve markedly. At the same time, device physics will become more complex, with the rapid increase in various small-geometry and quantum effects. This work describes both an approach to device simulator development and a physical model which advance the effort to meet the tremendous electronic device simulation challenge described above. The device simulation approach is to specify the physical model at a high level to a general-purpose (but highly efficient) partial differential equation solver (in this case PROPHET, developed by Lucent Technologies), which then simulates the model in 1-D, 2-D, or 3-D for a specified device and test regime. This approach allows for the rapid investigation of a wide range of device models and effects, which is certainly essential for device simulation to catch up with, and then stay ahead of, electronic device technology of the present and future. The physical device model used in this work is the density-gradient (DG) quantum correction to the drift-diffusion model [Ancona, Phys. Rev. B 35(5), 7959 (1987)]. This model adds tunneling and quantum smoothing of carrier density profiles to the drift-diffusion model. We used the DG model in 1-D and 2-D (for the first time) to simulate both bipolar and unipolar devices. Simulations of heavily-doped, short-base diodes indicated that the DG quantum corrections do not have a large effect on the IN characteristics of electronic devices without heteroj unction s. On the other hand, ultra-small MOSFETs certainly exhibit important quantum effects that the DG model will include: quantum repulsion of the inversion and gate charges from the oxide interfaces, and quantum tunneling through thin gate oxides. We present initial results of 2-D DG simulations of ultra-small MOSFETs. Subtle but important issues involving the specification of the model, boundary conditions, and interface constraints for DG simulation of MOSFETs will also be illuminated.

Biegal, Bryan A.↗

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

Modeling of a Stewart Platform for Analyzing One Directional Dynamics for Spacecraft Docking Operations

A one-directional dynamic model of a Stewart Platform was developed to assist NASA in analyzing the dynamic response in spacecraft docking operations. A simplified mechanical drawing was created, capturing the physical structure's main features. A simplified schematic diagram was developed in a lumped mass model from the mechanical drawing. Three differential equations were derived according to the schematic diagram. A Simulink diagram was created using MATLAB to represent the three equations. System parameters, including spring constants and masses, are derived in detail from the physical system. The model can be used for further analysis via computer simulation in predicting dynamic response in its main docking direction, i.e., up-and-down motion.

Stewart Platform↗

Galactic civilizations: Population dynamics and interstellar diffusion

The interstellar diffusion of galactic civilizations is reexamined by potential theory; both numerical and analytical solutions are derived for the nonlinear partial differential equations which specify a range of relevant models, drawn from blast wave physics, soil science, and, especially, population biology. An essential feature of these models is that, for all civilizations, population growth must be limited by the carrying capacity of the environment. Dispersal is fundamentally a diffusion process; a density-dependent diffusivity describes interstellar emigration. Two models are considered: the first describing zero population growth (ZPG), and the second which also includes local growth and saturation of a planetary population, and for which an asymptotic traveling wave solution is found.

Newman, W. I.↗

Medical support and technology for long-duration space missions

The current philosophy and development directions being taken towards realization of medical systems for use on board space stations are discussed. Data was gained on the performance of physical examinations, venipuncture and blood flow, blood smear and staining, white blood cell differential count, throat culture swab and colony count, and microscopy techniques during a 28-day period of the Skylab mission. It is expected that the advent of Shuttle flights will rapidly increase the number of persons in space, create a demand for in-space rather than on-earth medical procedures, and necessitate treatments for disorders without the provision for an early return to earth. Attention is being given to pressurized environment and extravehicular conditions of treatment, the possibilities of the use of the OTV for moving injured or ill crewmembers to other space stations, and to isolation of persons with communicable diseases from station crews.

Furukawa, S.↗

Elliptic generation of composite three-dimensional grids about realistic aircraft

An elliptic method for generating composite grids about realistic aircraft is presented. A body-conforming grid is first generated about the entire aircraft by the solution of Poisson's differential equation. This grid has relatively coarse spacing, and it covers the entire physical domain. At boundary surfaces, cell size is controlled and cell skewness is nearly eliminated by inhomogeneous terms, which are found automatically by the program. Certain regions of the grid in which high gradients are expected, and which map into rectangular solids in the computational domain, are then designated for zonal refinement. Spacing in the zonal grids is reduced by adding points with a simple, algebraic scheme. Details of the grid generation method are presented along with results of the present application, a wing-body configuration based on the F-16 fighter aircraft.

Sorenson, R. L.↗

Elliptic generation of composite three-dimensional grids about realistic aircraft

An elliptic method for generating composite grids about realistic aircraft is presented. A body-conforming grid is first generated about the entire aircraft by the solution of Poisson's differential equation. This grid has relatively coarse spacing, and it covers the entire physical domain. At boundary surfaces, cell size is controlled and cell skewness is nearly eliminated by inhomogeneous terms, which are found automatically by the program. Certain regions of the grid in which high gradients are expected, and which map into rectangular solids in the computational domain, are then designated for zonal refinement. Spacing in the zonal grids is reduced by adding points with a simple, algebraic scheme. Details of the grid-generation method are presented along with results of the present application, a wing/body configuration based on the F-16 fighter aircraft.

Sorenson, Reese L.↗

Measurements of fugitive hydrocarbon emissions with a tunable infrared DIAL

A tunable infrared differential absorption lidar (DIAL) system has been designed and developed at the National Physics Lab (NPL) which is capable of making measurements throughout the spectral region 3.0 to 4.2 micro-m. It is ideally suited to measuring a range of organic and inorganic species including methane, propane, and butane. The system also has an ultraviolet channel that is capable of making simultaneous measurements of aromatic hydrocarbons such as Toluene and benzene. This paper describes the source and detection system, together with some measurements of fugitive hydrocarbon emissions performed at various petrochemical plants.

Milton, M. J. T.↗

Sparking young minds with Moon rocks and meteorites

What could be more exciting than seeing pieces of other worlds? The Apollo program left a legacy of astounding accomplishments and precious samples. Part of the thrill of those lunar missions is brought to schools by the lunar sample educational disks, which contain artifacts of six piloted trips to the Moon. Johnson Space Center (JSC) is preparing 100 new educational disks containing pieces of meteorites collected in Antarctica. These represent chunks of several different asteroids, that were collected in one of the most remote, forbidding environments on Earth. These pieces of the Moon and asteroids represent the products of basic planetary processes (solar nebular processes, initial differentiation, volcanism, and impact), and, in turn, these processes are controlled by basic physical and chemical processes (energy, energy transfer, melting, buoyancy, etc.). Thus, the lunar and meteorite sample disks have enormous educational potential. New educational materials are being developed to accompany the disks. Present materials are not as effective as they could be, especially in relating samples to processes and to other types of data such as spectral studies and photogeology. Furthermore, the materials are out of date. New background materials will be produced for teachers, assembling slide sets with extensive captions, and devising numerous hands-on classroom activities to do while the disks are at a school and before and after they arrive. The classroom activities will be developed by teams of experienced teachers working with lunar and meteorite experts.

Taylor, G. Jeffrey↗

Shape-induced Gravitational Sorting of Saharan Dust During Transatlantic Voyage: Evidence from CALIOP Lidar Depolarization Measurements

Motivated by the physical picture of shape-dependent air resistance and, consequently, shape-induced differential sedimentation of dust particles, we searched for and found evidence of dust particle asphericity affecting the evolution and distribution of dust-scattered light depolarization ratio (delta). Specifically, we examined a large data set of Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) observations of Saharan dust from June to August 2007. Observing along a typical transatlantic dust track, we find that (1) median delta is uniformly distributed between 2 and 5 km altitudes as the elevated dust leaves the west coast of Africa, thereby indicating uniformly random mixing of particle shapes with height; (2) vertical homogeneity of median delta breaks down during the westward transport: between 2 and 5 km delta increases with altitude and this increase becomes more pronounced with westward progress; (3) delta tends to increase at higher altitude (greater than 4 km) and decrease at lower altitude (less than 4 km) during the westward transport. All these features are captured qualitatively by a minimal model (two shapes only), suggesting that shape-induced differential settling and consequent sorting indeed contribute significantly to the observed temporal evolution and vertical stratification of dust properties. By implicating particle shape as a likely cause of gravitational sorting, these results will affect the estimates of radiative transfer through Saharan dust layers.

CALIPSO↗

Responses of Microbes to Modeled Space Radiation

The built environment of spaceships is host to a microbial community that affects crew and craft alike. While the static composition of this community has been characterized and its temporal dynamics examined, the mechanisms controlling its make-up and evolutionary trajectory are not understood. Systematic analyses of microbial diversity show consistent patterns in community composition and function. Understanding these patterns' ecological origins remains a significant challenge, as it requires connecting processes at varying temporal and spatial scales. However, it is clear that the state and trajectories of microbial communities are in-part determined by their physical environment. In this regard, the spaceflight environment includes numerous interacting factors that differentiate it from Earth environments, including an altered atmospheric composition, reduced gravity (and thus altered fluid dynamics), and increased ionizing radiation. These factors impart selective pressures on microbial communities that affect their evolutionary trajectories and thus the risks and benefits these communities represent to crew and craft. The radiation environment of space leads to chronic exposure to low doses and is difficult to mimic on Earth. Thus, little is known about how microbial communities in spacecraft will respond and evolve. Therefore, given the limitations of existing studies, we aim to empirically determine how exposure to low doses of ionizing radiation for thousands of cell divisions affects rates of mutation accumulation in bacteria and the trajectory of their evolution. In this way, we will provide a critical set of data for designing safe and robust space missions. Here we discuss our progress towards this aim, including the construction of exposure facilities, our culturing and analysis approach, and preliminary data.

radiation↗

Ejecta patterns of Meteor Crater, Arizona derived from the linear un-mixing of TIMS data and laboratory thermal emission spectra

Accurate interpretation of thermal infrared data depends upon the understanding and removal of complicating effects. These effects may include physical mixing of various mineralogies and particle sizes, atmospheric absorption and emission, surficial coatings, geometry effects, and differential surface temperatures. The focus is the examination of the linear spectral mixing of individual mineral or endmember spectra. Linear addition of spectra, for particles larger than the wavelength, allows for a straight-forward method of deconvolving the observed spectra, predicting a volume percent of each endmember. The 'forward analysis' of linear mixing (comparing the spectra of physical mixtures to numerical mixtures) has received much attention. The reverse approach of un-mixing thermal emission spectra was examined with remotely sensed data, but no laboratory verification exists. Understanding of the effects of spectral mixing on high resolution laboratory spectra allows for the extrapolation to lower resolution, and often more complicated, remotely gathered data. Thermal Infrared Multispectral Scanner (TIMS) data for Meteor Crater, Arizona were acquired in Sep. 1987. The spectral un-mixing of these data gives a unique test of the laboratory results. Meteor Crater (1.2 km in diameter and 180 m deep) is located in north-central Arizona, west of Canyon Diablo. The arid environment, paucity of vegetation, and low relief make the region ideal for remote data acquisition. Within the horizontal sedimentary sequence that forms the upper Colorado Plateau, the oldest unit sampled by the impact crater was the Permian Coconino Sandstone. A thin bed of the Toroweap Formation, also of Permian age, conformably overlays the Coconino. Above the Toroweap lies the Permian Kiabab Limestone which, in turn, is covered by a thin veneer of the Moenkopi Formation. The Moenkopi is Triassic in age and has two distinct sub-units in the vicinity of the crater. The lower Wupatki member is a fine-grained sandstone, while the upper Moqui member is a fissile siltstone. Ejecta from these units are preserved as inverted stratigraphy up to 2 crater radii from the rim. The mineralogical contrast between the units, relative lack of post-emplacement erosion and ejecta mixing provide a unique site to apply the un-mixing model. Selection of the aforementioned units as endmembers reveals distinct patterns in the ejecta of the crater.

Ramsey, Michael S.↗

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

Test of the Equivalence Principle in an Einstein Elevator

The scientific goal of the experiment is to test the equality of gravitational and inertial mass (i.e., to test the Principle of Equivalence) by measuring the independence of the rate of fall of bodies from their compositions. The measurement is accomplished by measuring the relative displacement (or equivalently acceleration) of two falling bodies of different materials which are the proof masses of a differential accelerometer spinning about an horizontal axis to modulate a possible violation signal. A non-zero differential acceleration appearing at the signal frequency will indicate a violation of the Equivalence Principle. The goal of the experiment is to measure the Eotvos ratio og/g (differential acceleration/common acceleration) with a targeted accuracy that is about two orders of magnitude better than the state of the art (presently at several parts in 10(exp 13). The analyses carried out during this first grant year have focused on: (1) evaluation of possible shapes for the proof masses to meet the requirements on the higher-order mass moment disturbances generated by the falling capsule; (2) dynamics of the instrument package and differential acceleration measurement in the presence of errors and imperfections; (3) computation of the inertia characteristic of the instrument package that enable a separation of the signal from the dynamics-related noise; (4) a revised thermal analysis of the instrument package in light of the new conceptual design of the cryostat; (5) the development of a dynamic and control model of the capsule attached to the gondola and balloon to define the requirements for the leveling mechanism (6) a conceptual design of the leveling mechanism that keeps the capsule aligned before release from the balloon; and (7) a new conceptual design of the customized cryostat and a preliminary valuation of its cost. The project also involves an international cooperation with the Institute of Space Physics (IFSI) in Rome, Italy. The group at IFSI is in charge of prototyping the differential accelerometer and carrying out precursor laboratory measurements. During this grant year, our partners analyzed and then designed a new prototype of differential accelerometer that has several characteristics in common with the flight accelerometer at this point of the instrument development. The highlights of these activities are documented in a section of this report.

Shapiro, Irwin I.↗