Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear 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 253 records · Page 14

Massive graviton on arbitrary background: derivation, syzygies, applications

We give the detailed derivation of the fully covariant form of the quadratic action and the derived linear equations of motion for a massive graviton in an arbitrary background metric (which were presented in arXiv:1410.8302 [hep-th]). Our starting point is the de Rham-Gabadadze-Tolley (dRGT) family of ghost free massive gravities and using a simple model of this family, we are able to express this action and these equations of motion in terms of a single metric in which the graviton propagates, hence removing in particular the need for a “reference metric' which is present in the non perturbative formulation. We show further how 5 covariant constraints can be obtained including one which leads to the tracelessness of the graviton on flat space-time and removes the Boulware-Deser ghost. This last constraint involves powers and combinations of the curvature of the background metric. The 5 constraints are obtained for a background metric which is unconstrained, i.e. which does not have to obey the background field equations. We then apply these results to the case of Einstein space-times, where we show that the 5 constraints become trivial, and Friedmann-Lemaître-Robertson-Walker space-times, for which we correct in particular some results that appeared elsewhere. To reach our results, we derive several non trivial identities, syzygies, involving the graviton fields, its derivatives and the background metric curvature. These identities have their own interest. We also discover that there exist backgrounds for which the dRGT equations cannot be unambiguously linearized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deconstructing supergravity: massive supermultiplets

Given the success of the deconstruction program in obtaining ghost-free massive gravity from 5-D Einstein gravity, we propose a modification of the deconstruction procedure that incorporates supersymmetry at the linear level. We discuss the relevant limits of a conjectured interacting theory of a massive spin 2 supermultiplet, and determine the linear theory to be the N=1 Zinoviev theory, a supersymmetric extension of Fierz-Pauli theory. We develop a family of 1-site deconstruction procedures for fermionic fields (yielding Dirac and Majorana mass terms). The deconstruction procedure appropriate for giving fermions a Dirac mass is found to preserve half of the supersymmetry of the 5-D theory. We explicitly check this by deconstructing 5-D N=2 super-Maxwell theory down to 4-D N=1 super-Proca theory, and deconstructing linear 5-D N=2 supergravity down to 4-D N=1 Zinoviev theory, and derive the full 4-D supersymmetry algebras and Stückelberg symmetries from the 5-D superalgebras and gauge symmetries, respectively. We conjecture that this procedure should admit a generalization to fully non-linear theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Technical Evaluation of the NASA Model for Cancer Risk to Astronauts Due to Space Radiation

At the request of NASA, the National Research Council's (NRC's) Committee for Evaluation of Space Radiation Cancer Risk Model reviewed a number of changes that NASA proposes to make to its model for estimating the risk of radiation-induced cancer in astronauts. The NASA model in current use was last updated in 2005, and the proposed model would incorporate recent research directed at improving the quantification and understanding of the health risks posed by the space radiation environment. NASA's proposed model is defined by the 2011 NASA report Space Radiation Cancer Risk Projections and Uncertainties 2010 (Cucinotta et al., 2011). The committee's evaluation is based primarily on this source, which is referred to hereafter as the 2011 NASA report, with mention of specific sections or tables cited more formally as Cucinotta et al. (2011). The overall process for estimating cancer risks due to low linear energy transfer (LET) radiation exposure has been fully described in reports by a number of organizations. They include, more recently: (1) The "BEIR VII Phase 2" report from the NRC's Committee on Biological Effects of Ionizing Radiation (BEIR) (NRC, 2006); (2) Studies of Radiation and Cancer from the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR, 2006), (3) The 2007 Recommendations of the International Commission on Radiological Protection (ICRP), ICRP Publication 103 (ICRP, 2007); and (4) The Environmental Protection Agency s (EPA s) report EPA Radiogenic Cancer Risk Models and Projections for the U.S. Population (EPA, 2011). The approaches described in the reports from all of these expert groups are quite similar. NASA's proposed space radiation cancer risk assessment model calculates, as its main output, age- and gender-specific risk of exposure-induced death (REID) for use in the estimation of mission and astronaut-specific cancer risk. The model also calculates the associated uncertainties in REID. The general approach for estimating risk and uncertainty in the proposed model is broadly similar to that used for the current (2005) NASA model and is based on recommendations by the National Council on Radiation Protection and Measurements (NCRP, 2000, 2006). However, NASA's proposed model has significant changes with respect to the following: the integration of new findings and methods into its components by taking into account newer epidemiological data and analyses, new radiobiological data indicating that quality factors differ for leukemia and solid cancers, an improved method for specifying quality factors in terms of radiation track structure concepts as opposed to the previous approach based on linear energy transfer, the development of a new solar particle event (SPE) model, and the updates to galactic cosmic ray (GCR) and shielding transport models. The newer epidemiological information includes updates to the cancer incidence rates from the life span study (LSS) of the Japanese atomic bomb survivors (Preston et al., 2007), transferred to the U.S. population and converted to cancer mortality rates from U.S. population statistics. In addition, the proposed model provides an alternative analysis applicable to lifetime never-smokers (NSs). Details of the uncertainty analysis in the model have also been updated and revised. NASA's proposed model and associated uncertainties are complex in their formulation and as such require a very clear and precise set of descriptions. The committee found the 2011 NASA report challenging to review largely because of the lack of clarity in the model descriptions and derivation of the various parameters used. The committee requested some clarifications from NASA throughout its review and was able to resolve many, but not all, of the ambiguities in the written description.

Source record↗

Multidecadal Changes in Lower Stratospheric Ozone: Variability Vs. Trends

As upper stratospheric ozone appears to be recovering as a result of decreasing chlorine loading following the implementation of the Montreal Protocol and its amendments and in agreement with model projections, several recent studies report an apparent decline of ozone concentrations in the lower stratosphere in the last two decades, particularly in the extratropics. Our previous work as well as at least two other studies provide evidence that this decline results from transport changes rather than an intensification of chemical depletion. It remains unclear whether these changes represent long-term internal variability or are a consequence of a climate forcing. Here we perform free-running ensembles of the recent past (1980-2016) using the Goddard Earth Observing System Model (GEOS) at the cubed sphere C180 (approximately half degree) resolution. Two suites of 10-member ensembles are performed, one in which observed sea surface temperature (SSTs) are fully prescribed, and the other in which the linear SST trend over the recent past is removed so as to only retain internal variability. We evaluate the trends in both ozone as well as two idealized tracers with prescribed uniform loss that are used to isolate the role of transport from chemistry and emissions. Probability-distribution-functions of the trends in both ozone and idealized tracers are compared among ensemble members and with observed trends in order to evaluate the likelihood of recent observed declines in lower stratospheric ozone, relative to large internal variability. Moreover, comparisons among simulations with and without imposed SST trends indicate the extent to which dynamically-driven ozone trends reflect forced trends or internal variability in lower stratospheric dynamics.

Wargan, Krzysztof↗

Direct simulation of proton-coupled electron transfer across multiple regimes

The coupled transfer of electrons and protons is a central feature of biological and molecular catalysis, yet fundamental aspects of these reactions remain poorly understood. In this study, we extend the ring polymer molecular dynamics (RPMD) method to enable direct simulation of proton-coupled electron transfer (PCET) reactions across a wide range of physically relevant regimes. In a system-bath model for symmetric, co-linear PCET in the condensed phase, RPMD trajectories reveal distinct kinetic pathways associated with sequential and concerted PCET reaction mechanisms, and it is demonstrated that concerted PCET proceeds by a solvent-gating mechanism in which the reorganization energy is mitigated by charge cancellation among the transferring particles. We further employ RPMD to study the kinetics and mechanistic features of concerted PCET reactions across multiple coupling regimes, including the fully non-adiabatic (both electronically and vibrationally non-adiabatic), partially adiabatic (electronically adiabatic, but vibrationally non-adiabatic), and fully adiabatic (both electronically and vibrationally adiabatic) limits. Comparison of RPMD with the results of PCET rate theories demonstrates the applicability of the direct simulation method over a broad range of conditions; it is particularly notable that RPMD accurately predicts the crossover in the thermal reaction rates between different coupling regimes while avoiding a priori assumptions about the PCET reaction mechanism. Finally, by utilizing the connections between RPMD rate theory and semiclassical instanton theory, we show that analysis of ring-polymer configurations in the RPMD transition path ensemble enables the a posteriori determination of the coupling regime for the PCET reaction. This analysis reveals an intriguing and distinct 'transient-proton-bridge' mechanism for concerted PCET that emerges in the transition between the proton-mediated electron superexchange mechanism for fully non-adiabatic PCET and the hydrogen atom transfer mechanism for partially adiabatic PCET. Taken together, these results provide a unifying picture of the mechanisms and physical driving forces that govern PCET across a wide range of physical regimes, and they raise the possibility for PCET mechanisms that have not been previously reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

What is the distance to the CMB?

The success of precision cosmology depends not only on accurate observations, but also on the theoretical model --- which must be understood to at least the same level of precision. Subtle relativistic effects can lead to biased measurements if they are neglected. One such effect gives a systematic shift in the distance-redshift relation away from its background value, due to the non-linear relativistic conservation of total photon flux. We also show directly how this shift follows from a fully relativistic analysis of the geodesic deviation equation. We derive the expectation value of the shift using second-order perturbations about a concordance background, and show that the distance to last scattering is increased by 1%. We argue that neglecting this shift could lead to a significant bias in the background cosmological parameters, because it alters the meaning of the background model. A naive adjustment of CMB parameter estimation if this shift is really a correction to the background would raise the H{sub 0} value inferred from the CMB by 5%, potentially removing the tension with local measurements of H{sub 0}. Other CMB parameters which depend on the distance would also be shifted by ∼ 1σ when combined with local H{sub 0} data. While our estimations rely on a simplistic analysis, they nevertheless illustrate that accurately defining the background model in terms of the expectation values of observables is critical when we aim to determine the model parameters at the sub-percent level.

79 ASTRONOMY AND ASTROPHYSICS↗

A Full Demand Response Model in Co-Optimized Energy and

It has been widely accepted that demand response will play an important role in reliable and economic operation of future power systems and electricity markets. Demand response can not only influence the prices in the energy market by demand shifting, but also participate in the reserve market. In this paper, we propose a full model of demand response in which demand flexibility is fully utilized by price responsive shiftable demand bids in energy market as well as spinning reserve bids in reserve market. A co-optimized day-ahead energy and spinning reserve market is proposed to minimize the expected net cost under all credible system states, i.e., expected total cost of operation minus total benefit of demand, and solved by mixed integer linear programming. Numerical simulation results on the IEEE Reliability Test System show effectiveness of this model. Compared to conventional demand shifting bids, the proposed full demand response model can further reduce committed capacity from generators, starting up and shutting down of units and the overall system operating costs.

bidding mechanism↗

Groundwater flow and heat transport for systems undergoing freeze-thaw: Intercomparison of numerical simulators for 2D test cases

In high-elevation, boreal and arctic regions, hydrological processes and associated water bodies can be strongly influenced by the distribution of permafrost. Recent field and modeling studies indicate that a fully-coupled multidimensional thermo-hydraulic approach is required to accurately model the evolution of these permafrost-impacted landscapes and groundwater systems. However, the relatively new and complex numerical codes being developed for coupled non-linear freeze-thaw systems require verification. Here in this paper, this issue is addressed by means of an intercomparison of thirteen numerical codes for two-dimensional test cases with several performance metrics (PMs). These codes comprise a wide range of numerical approaches, spatial and temporal discretization strategies, and computational efficiencies. Results suggest that the codes provide robust results for the test cases considered and that minor discrepancies are explained by computational precision. However, larger discrepancies are observed for some PMs resulting from differences in the governing equations, discretization issues, or in the freezing curve used by some codes.

54 ENVIRONMENTAL SCIENCES↗

A comparison of neural network architectures for data-driven reduced-order modeling

The popularity of deep convolutional autoencoders (CAEs) has engendered new and effective reduced-order models (ROMs) for the simulation of large-scale dynamical systems. Despite this, it is still unknown whether deep CAEs provide superior performance over established linear techniques or other network-based methods in all modeling scenarios. To elucidate this, the effect of autoencoder architecture on its associated ROM is studied through the comparison of deep CAEs against two alternatives: a simple fully connected autoencoder, and a novel graph convolutional autoencoder. Through benchmark experiments, it is shown that the superior autoencoder architecture for a given ROM application is highly dependent on the size of the latent space and the structure of the snapshot data, with the proposed architecture demonstrating benefits on data with irregular connectivity when the latent space is sufficiently large.

42 ENGINEERING↗

Natural Covariant Planck Scale Cutoffs and the Cosmic Microwave Background Spectrum

We calculate here the impact of quantum gravity–motivated ultraviolet cutoffs on inflationary predictions for the cosmic microwave background spectrum. We model the ultraviolet cutoffs fully covariantly to avoid possible artifacts of covariance breaking. Imposing these covariant cutoffs results in the production of small, characteristically k-dependent oscillations in the spectrum. The size of the effect scales linearly with the ratio of the Planck to Hubble lengths during inflation. As a result, the relative size of the effect could be as large as one part in 10 5 ; i.e., eventual observability may not be ruled out.

79 ASTRONOMY AND ASTROPHYSICS↗

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

Perturbational treatment of spin-orbit coupling for generally applicable high-level multi-reference methods

An efficient perturbational treatment of spin-orbit coupling within the framework of high-level multi-reference techniques has been implemented in the most recent version of the COLUMBUS quantum chemistry package, extending the existing fully variational two-component (2c) multi-reference configuration interaction singles and doubles (MRCISD) method. The proposed scheme follows related implementations of quasi-degenerate perturbation theory (QDPT) model space techniques. Our model space is built either from uncontracted, large-scale scalar relativistic MRCISD wavefunctions or based on the scalar-relativistic solutions of the linear-response-theory-based multi-configurational averaged quadratic coupled cluster method (LRT-MRAQCC). The latter approach allows for a consistent, approximatively size-consistent and size-extensive treatment of spin-orbit coupling. The approach is described in detail and compared to a number of related techniques. The inherent accuracy of the QDPT approach is validated by comparing cuts of the potential energy surfaces of acrolein and its S, Se, and Te analoga with the corresponding data obtained from matching fully variational spin-orbit MRCISD calculations. The conceptual availability of approximate analytic gradients with respect to geometrical displacements is an attractive feature of the 2c-QDPT-MRCISD and 2c-QDPT-LRT-MRAQCC methods for structure optimization and ab inito molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Do responses to different anthropogenic forcings add linearly in climate models?

Many detection and attribution and pattern scaling studies assume that the global climate response to multiple forcings is additive: that the response over the historical period is statistically indistinguishable from the sum of the responses to individual forcings. Here, we use the NASA Goddard Institute for Space Studies (GISS) and National Center for Atmospheric Research Community Climate System Model (CCSM4) simulations from the CMIP5 archive to test this assumption for multi-year trends in global-average, annual-average temperature and precipitation at multiple timescales. We find that responses in models forced by pre-computed aerosol and ozone concentrations are generally additive across forcings. However, we demonstrate that there are significant nonlinearities in precipitation responses to different forcings in a configuration of the GISS model that interactively computes these concentrations from precursor emissions. We attribute these to differences in ozone forcing arising from interactions between forcing agents. Lastly, our results suggest that attribution to specific forcings may be complicated in a model with fully interactive chemistry and may provide motivation for other modeling groups to conduct further single-forcing experiments.

54 ENVIRONMENTAL SCIENCES↗

Drift-ideal magnetohydrodynamic simulations of m = 0 modes in Z-pinch plasmas

In this paper, the effects of m = 0 modes on equilibrium Z-pinch plasmas are studied using a drift-ideal magnetohydrodynamic (MHD) model. The model equations are an extension of ideal MHD to include finite-ion-inertial-length/cyclotron-frequency (Ω i ) effects in Ohm's law and in the electron and ion heat transport equations. The linear modes contained in this model include the ideal interchange (sausage) mode and in the magnetized limit, Ω i τ i »1 with τ i the ion collision time, nonideal entropy modes. It is well known that these two modes are decoupled in the kρ s « 1 limit, where k is the axial mode number and ρ s = c s /Ω i is the gyro-Bohm scale with c s the sound speed [B. Kadomtsev, Sov. Phys. JETP-USSR 10, 780 (1960)]. For Bennett equilibrium profiles, it is shown that the regions of stability for both modes are completely governed by the adiabatic coefficient γ in these limits. Equilibria with Bennett profiles are stable to entropy modes for γ < 2 but unstable to ideal modes and vice versa for γ > 2. However, these modes are no longer decoupled when kρ s ≳1. The simulation results of the fully nonlinear set of equations in the magnetized limit show that seeded modes with kρ s ≳1 and γ = 5/3 display the characteristics of both ideal and entropy modes. The general heat flux for both ions and electrons as a function of the species magnetization is retained in the model. Both the linear and nonlinear behaviors of seeded modes for kρ s ≳1 display a strong dependence on the magnetization of the ions. Lastly, the growth rate increases linearly with k at large kρ s when the ions are magnetized but decreases with increasing k when Ω i τ i ≲1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Satellite Constellation Orbit Design to Enable a Space-Based Radio Interferometer

Two different design methods for a networked constellation of N small satellites are presented. In the first method, the (linear) Clohessy-Wilthsire equations are used as an initial design tool, followed by conversion to a two-body model. Discrepancies between the linear and nonlinear solutions are minimized in the conversion process. The second method utilizes invariant manifold theory, by perturbing a reference trajectory in different directions along the center eigenvectors. Both methods require 5N parameters to fully define a constellation. In a relative, rotating frame of the reference path, the spacecraft appear as periodic ellipses of varying sizes. Deployment, reconfiguration using propulsive maneuvers, and station keeping costs for an example mission scenario are addressed.

Chien, Steve A.↗

Parallelized POD-based suboptimal economic model predictive control of a state-constrained Boussinesq approximation

Motivated by an energy efficient building application, we want to optimize a quadratic cost functional subject to the Boussinesq approximation of the Navier-Stokes equations and to bilateral state and control constraints. Since the computation of such an optimal solution is numerically costly, we design an efficient strategy to compute a sub-optimal (but applicationally acceptable) solution with significantly reduced computational effort. We employ an economic Model Predictive Control (MPC) strategy to obtain a feedback control. The MPC sub-problems are based on a linear-quadratic optimal control problem subjected to mixed control and state constraints and a convection-diffusion equation, reduced with proper orthogonal decomposition. Finally, to solve each sub-problem, we apply a primal-dual active set strategy. The method can be fully parallelized, which enables the solution of large problems with real-world parameters.

97 MATHEMATICS AND COMPUTING↗

Tidal alignment of galaxies

We develop an analytic model for galaxy intrinsic alignments (IA) based on the theory of tidal alignment. We calculate all relevant nonlinear corrections at one-loop order, including effects from nonlinear density evolution, galaxy biasing, and source density weighting. Contributions from density weighting are found to be particularly important and lead to bias dependence of the IA amplitude, even on large scales. This effect may be responsible for much of the luminosity dependence in IA observations. The increase in IA amplitude for more highly biased galaxies reflects their locations in regions with large tidal fields. We also consider the impact of smoothing the tidal field on halo scales. We compare the performance of this consistent nonlinear model in describing the observed alignment of luminous red galaxies with the linear model as well as the frequently used 'nonlinear alignment model,' finding a significant improvement on small and intermediate scales. We also show that the cross-correlation between density and IA (the 'GI' term) can be effectively separated into source alignment and source clustering, and we accurately model the observed alignment down to the one-halo regime using the tidal field from the fully nonlinear halo-matter cross correlation. Inside the one-halo regime, the average alignment of galaxies with density tracers no longer follows the tidal alignment prediction, likely reflecting nonlinear processes that must be considered when modeling IA on these scales. Finally, we discuss tidal alignment in the context of cosmic shear measurements.

79 ASTRONOMY AND ASTROPHYSICS↗