Search NASA⌕ Search

SEARCH · Search NASA

Results for “Free energy representation”

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.

83 records · Page 5

Feature extraction of event-related potentials using wavelets: an application to human performance monitoring

This report describes the development and evaluation of mathematical models for predicting human performance from discrete wavelet transforms (DWT) of event-related potentials (ERP) elicited by task-relevant stimuli. The DWT was compared to principal components analysis (PCA) for representation of ERPs in linear regression and neural network models developed to predict a composite measure of human signal detection performance. Linear regression models based on coefficients of the decimated DWT predicted signal detection performance with half as many free parameters as comparable models based on PCA scores. In addition, the DWT-based models were more resistant to model degradation due to over-fitting than PCA-based models. Feed-forward neural networks were trained using the backpropagation algorithm to predict signal detection performance based on raw ERPs, PCA scores, or high-power coefficients of the DWT. Neural networks based on high-power DWT coefficients trained with fewer iterations, generalized to new data better, and were more resistant to overfitting than networks based on raw ERPs. Networks based on PCA scores did not generalize to new data as well as either the DWT network or the raw ERP network. The results show that wavelet expansions represent the ERP efficiently and extract behaviorally important features for use in linear regression or neural network models of human performance. The efficiency of the DWT is discussed in terms of its decorrelation and energy compaction properties. In addition, the DWT models provided evidence that a pattern of low-frequency activity (1 to 3.5 Hz) occurring at specific times and scalp locations is a reliable correlate of human signal detection performance. Copyright 1999 Academic Press.

Evoked Potentials↗

Dynamics and X-ray emission of a galactic superwind interacting with disk and halo gas

There is a general agreement that the conspicuous extranuclear X-ray, optical-line, and radio-contiuum emission of starbursts is associated with powerful galactic superwinds blowing from their centers. However, despite the significant advances in observational studies of superwinds, there is no consensus on the nature of the emitting material and even on the emission mechanisms themselves. This is to a great extent a consequence of a poor understanding of dynamical processes in the starburst superwind regions. To address this issue, we have conducted two-dimensional hydrodynamical simulations of galactic superwinds. While previous similar studies have used a single (disk) component to represent the ISM of the starburst galaxy, we analyze the interaction of the wind with a two-component disk-halo ambient interstellar medium and argue that this two-component representation is crucial for adequate modeling of starbursts. The emphasis of this study is on the geometry and structure of the wind region and the X-ray emission arising in the wind material and the shocked gas in the disk and the halo of the galaxy. The simulation results have shown that a clear-cut bipolar wind can easily develop under a range of very different conditions. On the other hand, a complex 'filamentary' structure associated with the entrained dense disk material is found to arise within the hot bubble blown out by the wind. The flow pattern within the bubble is dominated equally by the central biconic outflow and a system of whirling motions r elated to the origin and development of the 'filaments'. The filament parameters make them a good candidate for optical-emission-line filamentary gas observed in starburst halos. We find that the history of mass and energy deposition in the starburst region of the galaxy is crucial for wind dynamics. A 'mild' early wind, which arises as a result of the cumulative effect of stellar winds from massive stars, produces a bipolar vertical cavity in the disk and halo gas without strongly affecting the gaseous disk, thus creating conditions for virtually free vertical escape of the hot gas at the later, much more violent supernova-dominated phases of the starburst. We calculate the luminosity, mass, and effective temperature of the X-ray emitting gas in the 'soft' (0.1 to 0.7 keV, 0.7 to 2.2 keV, and 0.1 to 2.2 keV) and 'hard' (1.6 to 8.3 keV) energy bands and estimate the contribution of different gaseous components to the X-ray flux in these bands. Analysis of these parameters enables us to make conclusions regarding the nature of the X-ray-emitting material. We have inferred that the bulk of the soft thermal X-ray emission from starbursts arises in the wind-shocked material of the disk and halo gas rather than in the wind material itself. This enables us to predict that the integrated soft X-ray spectra of starbursts need not show an overabundance of heavy elements which are believed to be produced copiously in the centers of starbursts. Unlike soft X-ray emission, the hard component of thermal X-ray emission is found to originate in the wind material ejected from the starburst region. However, the derived ratio of hard-to-soft X-ray luminosities is too small compared to that observed in starbursts. We conclude therefore that the observed hard X-ray emission of starbursts is probably not associated with the thermal emission of hot wind or ambient shocked gas. Typical temperatures of the bulk of the soft X-ray-emitting material in our very different models have been found to agree well with the ones estimated on the basis of the ROSAT data for the soft component of X-ray emission of nearby starbursts. We predict that temperatures of the extranuclear soft X-ray-emitting gas in starburst galaxies with heavy element abundances near solar should be close to T(sub Xs = 2 to 5 x 10(exp 6)K.

Suchkov, Anatoly A.↗

Development of a finite element based delamination analysis for laminates subject to extension, bending, and torsion

Delamination is a common failure mode of laminated composite materials. This type of failure frequently occurs at the free edges of laminates where singular interlaminar stresses are developed due to the difference in Poisson's ratios between adjacent plies. Typically the delaminations develop between 90 degree plies and adjacent angle plies. Edge delamination has been studied by several investigators using a variety of techniques. Recently, Chan and Ochoa applied the quasi-three-dimensional finite element model to the analysis of a laminate subject to bending, extension, and torsion. This problem is of particular significance relative to the structural integrity of composite helicopter rotors. The task undertaken was to incorporate Chan and Ochoa's formulation into a Raju Q3DG program. The resulting program is capable of modeling extension, bending, and torsional mechanical loadings as well as thermal and hygroscopic loadings. The addition of the torsional and bending loading capability will provide the capability to perform a delamination analysis of a general unsymmetric laminate containing four cracks, each of a different length. The solutions obtained using this program are evaluated by comparing them with solutions from a full three-dimensional finite element solution. This comparison facilitates the assessment of three dimensional affects such as the warping constraint imposed by the load frame grips. It wlso facilitates the evaluation of the external load representation employed in the Q3D formulation. Finally, strain energy release rates computed from the three-dimensional results are compared with those predicted using the quasi-three-dimensional formulation.

Hooper, Steven J.↗

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Implications of Lag-Luminosity Relationship for Unified GRB Paradigms

Spectral lags (tau(sub lag)) are deduced for 1437 long (T(sub 90) greater than 2 s) BATSE gamma-ray bursts (GRBs) with peak flux F(sub p) greater than 0.25 photons cm(sup -2)/s, near to the BATSE trigger threshold. The lags are modeled to approximate the observed distribution in the F(sub p)-T(sub lag) plane, realizing a noise-free representation. Assuming a two-branch lag-luminosity relationship, the lags are self-consistently corrected for cosmological effects to yield distributions in luminosity, distance, and redshift. The results have several consequences for GRB populations and for unified gamma-ray/afterglow scenarios which would account for afterglow break times and gamma-ray spectral evolution in terms of jet opening angle, viewing angle, or a profiled jet with variable Lorentz factor: A component of the burst sample is identified - those with few, wide pulses, lags of a few tenths to several seconds, and soft spectra - whose Log[N]-Log[F(sub p)] distribution approximates a -3/2 power-law, suggesting homogeneity and thus relatively nearby sources. The proportion of these long-lag bursts increases from negligible among bright BATSE bursts to approx. 50% at trigger threshold. Bursts with very long lags, approx. 1-2 less than tau(sub lag) (S) less than 10, show a tendency to concentrate near the Supergalactic Plane with a quadrupole moment of approx. -0.10 +/- 0.04. GRB 980425 (SN 1998bw) is a member of this subsample of approx. 90 bursts with estimated distances less than 100 Mpc. The frequency of the observed ultra-low luminosity bursts is approx. 1/4 that of SNe Ib/c within the same volume. If truly nearby, the core-collapse events associated with these GRBs might produce gravitational radiation detectable by LIGO-II. Such nearby bursts might also help explain flattening of the cosmic ray spectrum at ultra-high energies, as observed by AGASA.

Norris, J. P.↗

Labyrinth Seal Flutter Analysis and Test Validation in Support of Robust Rocket Engine Design

High energy-density turbomachines, like the SSME turbopumps, utilize labyrinth seals, also referred to as knife-edge seals, to control leakage flow. The pressure drop for such seals is order of magnitude higher than comparable jet engine seals. This is aggravated by the requirement of tight clearances resulting in possible unfavorable fluid-structure interaction of the seal system (seal flutter). To demonstrate these characteristics, a benchmark case of a High Pressure Oxygen Turbopump (HPOTP) outlet Labyrinth seal was studied in detail. First, an analytical assessment of the seal stability was conducted using a Pratt & Whitney legacy seal flutter code. Sensitivity parameters including pressure drop, rotor-to-stator running clearances and cavity volumes were examined and modeling strategies established. Second, a concurrent experimental investigation was undertaken to validate the stability of the seal at the equivalent operating conditions of the pump. Actual pump hardware was used to construct the test rig, also referred to as the (Flutter Rig). The flutter rig did not include rotational effects or temperature. However, the use of Hydrogen gas at high inlet pressure provided good representation of the critical parameters affecting flutter especially the speed of sound. The flutter code predictions showed consistent trends in good agreement with the experimental data. The rig test program produced a stability threshold empirical parameter that separated operation with and without flutter. This empirical parameter was used to establish the seal build clearances to avoid flutter while providing the required cooling flow metering. The calibrated flutter code along with the empirical flutter parameter was used to redesign the baseline seal resulting in a flutter-free robust configuration. Provisions for incorporation of mechanical damping devices were introduced in the redesigned seal to ensure added robustness

El-Aini, Yehia↗

Efficient Four-Component Dirac–Coulomb–Gaunt Hartree–Fock in the Pauli Spinor Representation

Four-component Dirac Hartree--Fock is an accurate mean-field method for treating molecular systems where relativistic effects are important. However, the computational cost and complexity of the two-electron interaction makes this method less common, even though we can consider the Dirac Hartree--Fock Hamiltonian the ``ground truth" of electronic structure, barring explicit quantum-electrodynamical effects. Being able to calculate these effects is then vital to the design of lower scaling methods for accurate predictions in computational spectroscopy and properties of heavy element complexes that must include relativistic effects for even qualitative accuracy. In this work, we present a Pauli quaternion formalism of maximal component- and spin-separation for computing the Dirac-Coulomb-Gaunt Hartree-Fock ground state, with a minimal floating-point-operation count algorithm. Furthermore, this approach also allows one to explicitly separate different spin physics from the two-body interactions, such as spin-free, spin-orbit, and the spin-spin contributions. Additionally, we use this formalism to examine relativistic trends in the periodic table, and analyze the basis set dependence of atomic gold and gold dimer systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lumley's energy cascade dissipation rate model for boundary-free turbulent shear flows

True dissipation occurs mainly at the highest wavenumbers where the eddy sizes are comparatively small. These high wavenumbers receive their energy through the spectral cascade of energy starting with the largest eddies spilling energy into the smaller eddies, passing through each wavenumber until it is dissipated at the microscopic scale. However, a small percentage of the energy does not spill continuously through the cascade but is instantly passed to the higher wavenumbers. Consequently, the smallest eddies receive a certain amount of energy almost immediately. As the spectral energy cascade continues, the highest wavenumber needs a certain time to receive all the energy which has been transferred from the largest eddies. As such, there is a time delay, of the order of tau, between the generation of energy by the largest eddies and the eventual dissipation of this energy. For equilibrium turbulence at high Reynolds numbers, there is a wide range where energy is neither produced by the large eddies nor dissipated by viscosity, but is conserved and passed from wavenumber to higher wavenumbers. The rate at which energy cascades from one wavenumber to another is proportional to the energy contained within that wavenumber. This rate is constant and has been used in the past as a dissipation rate of turbulent kinetic energy. However, this is true only in steady, equilibrium turbulence. Most dissipation models contend that the production of dissipation is proportional to the production of energy and that the destruction of dissipation is proportional to the destruction of energy. In essence, these models state that the change in the dissipation rate is proportional to the change in the kinetic energy. This assumption is obviously incorrect for the case where there is no production of turbulent energy, yet energy continues to cascade from large to small eddies. If the time lag between the onset on the energy cascade to the destruction of energy at the microscale can be modeled, then there will be a better representation of the dissipation process. Development of an energy cascade time scale equation is discussed.

Duncan, B. S.↗

NGEE Arctic CO2, CH4 and Energy Eddy-Covariance (EC) Flux Tower Auxiliary Measurements, Utqiagvik (Barrow), Alaska, 2012 - 2019

Measurements began in September 2012 and are generally made throughout the snow-free growing season starting from April (pre-thawing) until November (freeze-up) when the instrumentation is taken down for maintenance and instrument calibration. CO2, CH4 and energy fluxes using the Eddy Covariance (EC) technique (Baldocchi 2003) and meteorological measurements are reported as 30-minute averages. The tower site is registered with AmeriFlux as US-NGB. See http://ameriflux.lbl.gov/sites/siteinfo/US-NGB for more information. Reported data include: (1) additional/extra meteorological and eddy covariance CO2, CH4 and energy flux data, (2) binned spectra for the three wind components, the sonic temperature and gas densities together with the binned cospectra for covariances of w (vertical wind component) and gas densities, (3) binned ogives (cumulative (co)spectra) for the three wind components, the sonic temperature and gas densities, and (4) tower site footprint matrices. Meteorological and eddy covariance CO2, CH4 and energy flux data (AmeriFlux BASE data product) are available from AmeriFlux https://doi.org/10.17190/AMF/1436326. The tower site is adjacent to the NGEE Arctic Tram system with several complementary time-series measurements https://doi.org/10.5440/1558782.This dataset supersedes the previously released dataset "Eddy-Covariance and auxiliary measurements, NGEE-Barrow, 2012-2013" DOI:10.5440/1124200 as this new product includes data for the longer time period processed using consistent methods for the entire record. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) near Utqiagvik (Barrow), Alaska and 2) multiple areas in the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets

Natural orbitals are often used to achieve a more compact representation of correlated wave-functions. Using natural orbitals computed as eigenstates of the virtual–virtual block of the state density matrix instead of the canonical Hartree–Fock orbitals results in smaller errors when the same fraction of virtual space is frozen. This strategy, termed frozen natural orbital (FNO) approach, is effective in reducing the cost of regular coupled-cluster (CC) calculations and some multistate methods, such as EOM-IP-CC (equation-of-motion CC for ionization potentials). This contribution extends the FNO approach to the EOM-SF-CC ansatz (EOM-CC with spin-flip). In contrast to EOM-IP-CCSD, EOM-SF-CCSD relies on high-spin open-shell references. Using FNOs computed for an open-shell reference leads to an erratic behavior of the EOM-SF-CC energies and properties due to an inconsistent truncation of the α and β orbital spaces. A general solution to problems arising in the EOM-CC calculations utilizing open-shell references, termed OSFNO (open-shell FNO), is proposed. Furthermore, by means of singular value decomposition (SVD) of the overlap matrix of the α and β orbitals, the OSFNO algorithm identifies the corresponding orbitals and determines virtual orbitals corresponding to the singly occupied space. This is followed by SVD of the singlet part of the state density matrix in the remaining virtual orbital subspace. The so-computed FNOs preserve the spin purity of the open-shell orbital subspace to the extent allowed by the original reference, thus facilitating a safe truncation of the virtual space. The performance of OSFNO is benchmarked for selected diradicals and triradicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automated identification and calculation of prompt effects in kinetic mechanisms using statistical models

The kinetics of prompt dissociation involves rovibrationally excited species (generally formed by exothermic reactions) which may dissociate or isomerize prior to thermalization via collisions with the bath gas. Treating such rovibrationally excited species (so-called "hot" species) with standard kinetic phenomenology may result in incorrect macroscopic representation of their reactivity. Here this work presents the first fully automated methodology for the calculation of prompt effects of a chosen species in a kinetic mechanism, including (i) reaction selection; (ii) theoretical calculation of rate constants and prompt branching fractions; and (iii) final rate constant fitting. The energy partition between hot fragments is estimated using a variety of statistical models, including a new physically sound microcanonical statistical model based on the rovibrational density of states of the fragments. The methodology is validated against literature data for the prompt dissociations of HCO and C 3 H 7 radicals. The microcanonical statistical model is in better agreement with trajectory simulations for larger species and is thus applicable for practical systems that typically involve large molecules, for which direct dynamics calculations are impractical. The automated workflow is applied to the evaluation of the effects of prompt dissociation for two isomeric radicals C 4 H 7 1-3 (1-methylallyl) and C 4 H 7 1-4 (3-buten-1-yl). Twelve H-atom abstraction reactions are selected and the corresponding rate constants are computed with first principles theory. The microcanonical statistical model predicts that prompt dissociations of C 4 H 7 1-3 and C 4 H 7 1-4 are already significant at 1000K, resulting in differences of up to an order of magnitude at 2000K with respect to the phenomenological thermal rate constants. To illustrate the effects of prompt dissociation on simulations of experimental data, the calculated prompt rate constants are implemented in both CRECK and C3MechV3.3 kinetic mechanisms. Simulations of experimental flame data illustrate the noticeable impact of prompt dissociation kinetics on the high-temperature combustion reactivity of C 4 H 8 -1 and C 4 H 8 -2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗