SEARCH · Search NASA
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.
Various combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are studied by a priori tests and numerical simulations. Consistency between model and filter is found to be essential to ensure accurate results. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.
A series of models of the lunar interior are derived from topographic, gravitational, librational, and seismic data. The librational parameters and low-degree gravity harmonics result primarily from surface height variations and only secondarily from lateral density variations. The moon departs from isostasy, even for the low-degree harmonics, with a maximum superisostatic stress of 200 bars under the major mascon basins. The mean crustal thicknesses under different physiographic regions are: mascons, 30-35 km; irregular maria, 50-60 km; and highlands, 90-110 km. A possible composition consistent with our model is an anorthositic crust, underlain by a predominantly forsterite upper mantle which grades into a refractory rich lower mantle surrounding a pyrrhotite core.
This presentation discusses the importance of consistent evaluations for fission observables. Currently, prompt fission quantities (e.g., neutron multiplicity, neutron energy spectrum, and fission product yields) are evaluated independently of one another, either with separate models or data-only/non-model evaluations (few – if any – shared model inputs) and inconsistencies can arise in evaluated data (e.g. NSE 190, 258 (2018)). They have found that consistent, model-based evaluations lead to more physical constraints on free parameters in the model, more robust predictions for other observables and other isotopes, that there are correlations between the uncertainties on different observables, and that uncertainties can be propagated to unmeasured quantities. They are working on developing models to the point where fission observables are calculated at the quality needed for an evaluation. They are striving to consistently calculate fission observables, such as $\bar{v}$, PFNS, and FPY, to a quality suitable for an evaluation. Current evaluations for fission observables use separate models for each observable (with minimal shared input) or are purely data driven. LANL models, such as the Monte Carlo CGMF and deterministic BeoH, have different strengths; the similar fission fragment initial conditions and decay models leads to consistency between the two codes (FPY can be connected to $\bar{v}$ and PFNS, even when calculated by different codes). Work is underway by Lovell to optimize the models and Neudecker is working to perform detailed uncertainty quantification. The PFNS, in particular, presents a significant challenge.
We examine the possibility of constructing scale-free triaxial logarithmic potentials self-consistently, using Schwarzschild's linear programing method. In particular, we explore the limit of nonaxisymmetric disks. In this case it is possible to reduce the problem to the self-consistent reconstruction of the disk surface density on the unit circle, a considerably simpler problem than the usual 2D or 3D one. Models with surface densities of the form Sigma = (x exp n + (y/q) exp n) exp - 1/n with n = 2 or 4 are investigated. We show that the complicated shapes of the 'boxlet' orbit families (which replace the box orbit family found in potentials with smooth cores) limit the possibility of building self-consistent models, though elliptical disks of axis ratio above 0.7 and a restricted range of boxier models can be constructed. This result relies on using sufficiently fine bins, smaller than the 10 deg bins commonly used in 2D or 3D investigations. It also indicates the need for caution in interpreting N-body models of triaxial halos in which the core of the potential is numerically smoothed.
An approach is disclosed that maintains a consistent view of a virtual address by a local node which writes a first value to the virtual address and, after writing the first value, establishes a snapshot consistency state of the virtual address. The virtual address is shared amongst any number of processes and the processes includes a writing process and other processes that read from the virtual address. After writing the first value, the writing process writes a second value to the virtual address. Even after writing the second value, the first value is still visible to the other processes.
This paper continues presentation and discussion of the results from our new global self-consistent theoretical model of interacting ring current ions and propagating electromagnetic ion cyclotron waves [Khazanov et al., 2006]. To study the effects of electromagnetic ion cyclotron wave propagation and refraction on the wave induced ring current precipitation and heating of the thermal plasmaspheric electrons, we simulate the May 1998 storm. The main findings after a simulation can be summarized as follows. Firstly, the wave induced ring current precipitation exhibits quite a lot of fine structure, and is highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 x 10(exp 6) (cm(raised dot) s(raised dot) sr(raised dot) (sup -1)) are observed during the maill and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not connected to the most intense waves in simple manner. The characteristics of the wave power spectral density distribution over the wave normal angle are extremely crucial for the effectiveness of the ring current ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from RAM [Kozyra et al., 1997a] reveals that although we observe a qualitative agreement between the localizations of the wave induced precipitations in the models, there is no quantitative agreement between the magnitudes of the fluxes. The quantitative differences are mainly due to a qualitative difference between the characteristics of the wave power spectral density distributions over the wave normal angle in RAM and in our model. Thirdly, the heat fluxes to plasmaspheric electrons caused by Landau resonate energy absorption from electromagnetic ion cyclotron waves are observed in the postnoon-premidnight MLT sector, and can reach the magnitude of 10(exp 11) eV/(cm(sup 2)(raised dot)s). The Coulomb energy degradation of the RC H(+) and O(+) ions maximizes at about 10(exp 11) (eV/(cm(sup 2) (raised dot) s), and typically leads to electron energy deposition rates of about 2(raised dot) 10(exp 10) (eV/(cm(sup 2)(raised dot)s) which are observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the ring current with the plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, ring current, and topside ionosphere [Gurgiolo et al., 2005]. Finally, the wave induced intense electron heating has a structure of the spot-like patches along the most enhanced density gradients in the plasmasphere boundary layer and can be a possible driver to the observed but still not explained small-scale structures of enhanced emissions in the stable auroral red arcs.
Electromagnetic ion cyclotron (EMIC) waves are an important magnetospheric emission, which is excited near the magnetic equator with frequencies below the proton gyro-frequency. The source of bee energy for wave growth is provided by temperature anisotropy of ring current (RC) ions, which develops naturally during inward convection from the plasma sheet These waves strongly affect the dynamic s of resonant RC ions, thermal electrons and ions, and the outer radiation belt relativistic electrons, leading to non-adiabatic particle heating and/or pitch-angle scattering and loss to the atmosphere. The rate of ion and electron scattering/heating is strongly controlled by the Wave power spectral and spatial distributions, but unfortunately, the currently available observational information regarding EMIC wave power spectral density is poor. So combinations of reliable data and theoretical models should be utilized in order to obtain the power spectral density of EMIC waves over the entire magnetosphere throughout the different storm phases. In this study, we present the simulation results, which are based on two coupled RC models that our group has developed. The first model deals with the large-scale magnetosphere-ionosphere electrodynamic coupling, and provides a self-consistent description of RC ions/electrons and the magnetospheric electric field. The second model is based on a coupled system of two kinetic equations, one equation describes the RC ion dynamics and another equation describes the power spectral density evolution of EMIC waves, and self-consistently treats a micro-scale electrodynamic coupling of RC and EMIC waves. So far, these two models have been applied independently. However, the large-scale magnetosphere-ionosphere electrodynamics controls the convective patterns of both the RC ions and plasmasphere altering conditions for EMIC wave-particle interaction. In turn, the wave induced RC precipitation Changes the local field-aligned current distributions and the ionospheric conductances, which are crucial for a large-scale electrodynamics. The initial results from this new self-consistent model of the magnetospheric electric field, RC and EMIC waves will be shown in this presentation.
Electromagnetic backscattering from an open ended three dimensional inlet model is analyzed and computed patterns are compared with results of experimental measurements. The model is comprised of two sections. The first section consists of a linearly tapered waveguide with a rectangular opening at one end and the other end is connected to the second section which is a uniform rectangular waveguide with a planar perfectly conducting termination. The model is electrically large so that many propagating modes are excited. The method of analysis contains conventional aperture integration and modal techniques combined with high frequency techniques, which employ concepts such as modal rays, geometrical theory of diffraction and equivalent currents. For the cases considered, it is shown that only a few of the many propagating modes contribute appreciably to the backscattered field. These modes are selected according to their modal ray angle directions.
HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Here using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.
In the rotation-dominated model of the Jovian magnetosphere, proposed in the present paper, the plasma current distribution is dynamically consistent with the magnetic field. The model is constructed under the assumptions of an ionospheric plasma source without pitch angle scattering; loss-free radiation transport through flux tube interchange diffusion; and a static balance between centrifugal force and magnetic stress. Using these assumptions, a dynamic equation is derived for a self-consistent field configuration having a single adjustable parameter that is related to the plasma source strength. The equation is solved by an iterative technique. The self-consistent field resembles a radially distorted dipole field, rather than the flat magnetodisk configuration inferred from Pioneer data. The implication is that such a magnetodisk field would require an equatorial plasma source rather than a source at the feet of the field lines.
A simple, two-dimensional, steady state, viscous model of the dawnside and duskside low-latitude boundary layer is presented. The model described boundary layer currents, magnetic fields, and plasma flow in a self-consistent way. Coupling to the ionosphere via field-aligned currents and associated field-aligned potential drops governed by a simple conductance law are included in the model. The model equations are given and the results from the model are compared with observations.
A short analysis is presented of Earth Radiation Budget Experiment (ERBE) errors inherent in the directional models used for data interpretation. The models were all developed on the basis of experience with the Nimbus-7 ERB experiment, which had a spatial resolution one-third that of ERBE instrumentation. A pseudo-directional model is defined to simulate the ERBE scanner data, using the assumptions that the average radiant exitance for any particular scene is independent of the viewing geometry, geographic location and time the data is collected. The directionality of the view angle and solar zenith angle is accounted for by a method of bins.
Second-order Møller-Plesset perturbation theory is well-known as a computationally inexpensive approach to the electron correlation problem that is size-consistent with a size-consistent reference but fails to be regular. On the other hand, the less well-known many-body version of Brillouin-Wigner perturbation theory has the reverse properties: it is regular but fails to be size-consistent when used with the standard MP partitioning. Consequently, its widespread use remains limited. In this work, we analyze the ways in which it is possible to use alternative non-MP partitions of the Hamiltonian to yield variants of BW2 that are size-consistent as well as regular. We show that there is a vast space of such BW2 theories and also show that it is possible to define a repartitioned BW2 theory from the ground state density alone, which regenerates the exact correlation energy. We also provide a general recipe for deriving regular, size-consistent, and size-extensive partitions from physically meaningful components, and we apply the result to small model systems. The scope of these results appears to further set the stage for a revival of BW2 in quantum chemistry.
Calculations of accreting magnetized neutron star atmospheres heated by the gradual deceleration of protons via Coulomb collisions are presented. Self consistent determinations of the temperature and density structure for different accretion rates are made by assuming hydrostatic equilibrium and energy balance, coupled with radiative transfer. The full radiative transfer in two polarizations, using magnetic cross sections but with cyclotron resonance effects treated approximately, is carried out in the inhomogeneous atmospheres.
Calculations of accreting magnetized neutron star atmospheres heated by the gradual deceleration of Protons via Coulomb collisions are presented. Self consistent determinations of the temperature and density structure for different accretion rates are made by assuming hydrostatic equilibrium and energy balance, coupled with radiative transfer. The full radiative transfer in two polarizations, using magnetic cross sections but with cyclotron resonance effects treated approximately, is carried out in the inhomogeneous atmospheres. Previously announced in STAR as N84-12012
Motivated by the observed correlation between black hole masses M(sub BH) and the velocity dispersion sigma of host galaxies, we develop a theoretical model of black hole formation in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an initial state specified by a uniform rotation rate OMEGA and a density distribution of the form rho = a(sup 2)(sub eff)per2piGR(sup 2)(so that a(sub eff)is an effective transport speed). The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole (for Schwarzschild geometry). This model reproduces the observed correlation between the estimated black hole masses and the velocity dispersions of galactic bulges, i.e., M(sub BH) approximately equal to 10(sup 8) solar mass(sigma per 200 kilometers per second)(sup 4) where sigma = the square root of 2a(sub eff). To obtain this normalization, the rotation rate OMEGA approximately equal to 2 x 10(exp -15) rad per second. The model also defines a bulge mass scale M(sub B). If we identify the scale M(sub B) with the bulge mass, the model determines the ratio mu(sub B) of black hole mass to the host mass: mu(sub B) approximately equal to 0.0024(sigma per 200 kilometer per second), again in reasonable agreement with observed values. In this scenario, supermassive black holes form quickly (in approximately 10(exp 5) yr) and are born rapidly rotating (with a per M approximately 0.9). This paper also shown how these results depend on the assumed initial conditions; the most important quantity is the initial distribution of specific angular momentum in the precollapse state.
Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of ring current ions and ion cyclotron waves in a quasilinear approach. These two equations were solved on a global scale under non steady-state conditions during the May 2-5, 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the wave active zones at three time cuts around initial, main, and late recovery phases of the May 4, 1998 storm phase are presented and discussed in detail. Comparisons of the model wave-ion data with the Polar/HYDRA and Polar/MFE instruments results are presented..
A viscoplastic constitutive model is presented for the estimation of the overall mechanical response of Hastelloy-X polycrystalline metals from a knowledge of single crystal behavior. The behavior of polycrystal is derived from that of single crystals using a self-consistent formulation. The single crystal behavior which has been used was developed by summing postulated slip on crystallographic slip systems. The plasticity and creep are treated coupledly using unified viscoplastic model which includes the interaction effects between rapid and slow deformation at elevated temperature. The validity of the model is directly tested by experiments on Hastelloy-X in both single crystal and polycrystalline versions.