Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational modes”

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 955 records · Page 53

Anharmonic Rovibrational Partition Functions at High Temperatures: Tests of Reduced-Dimensional Models for Systems with up to Three Fluxional Modes

In this work, curvilinear coordinate Monte Carlo phase space integration and a series of full-dimensional fitted potential energy surfaces are used to study the effectiveness of reduced-dimensional models for predicting rovibrational anharmonicity at high temperatures. Fully coupled and fully anharmonic, but classical, rovibrational partition functions Q are computed for 14 species with two or three fluxional modes (inversions or torsions) and as many as 30 degrees of freedom. These results are converted to semiclassical anharmonicity correction factors f and are analyzed alongside results obtained previously for 22 species with up to two fluxional modes. As expected, fluxional species show considerable variation in f at high temperatures; f is as small as 0.2 for acetone and is as large as 9 for methylene glycol at 2500 K. This set of full-dimensional results is used to test the accuracy of reduced-dimensional models where fluxional modes are treated as coupled to one another but as separable from the remaining nonfluxional modes. For most systems, we find that such an approximation is accurate at high temperatures, with average errors in Q of just similar to 25%. For some systems, however, larger errors are found, and these are attributed to strong coupling of the fluxional modes to one or more nonfluxional modes. In particular, we identify strong coupling to low-frequency bends for some systems, and we show that by comparing curvilinear and rectilinear harmonic frequencies for the fluxional modes, we can estimate the effect of this coupling on rovibrational anharmonicity. We also quantify the accuracy of the more severe but common assumption of treating fluxional modes as separable from one another, that is, as sets of uncoupled one-dimensional inversions and torsions. This approach can work well for methyl and alkyl rotors, but it is shown to have errors as large as a factor of 7 at high temperatures for more complex systems. Lastly, we note that while the present analysis focuses on the treatment of fluxional modes, the collective anharmonicity correction associated with the more numerous nonfluxional modes, although simpler to describe, comprises a significant fraction of the overall anharmonicity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The transmission of sound in nonuniform ducts

The method of weighted residuals in the form of a modified Galerkin method with boundary residuals was developed for the study of the transmission of sound in nonuniform ducts carrying a steady, compressible flow. In this development, the steady flow was modeled as essentially one dimensional but with a kinematic modification to force tangency of the flow at the duct walls. Three forms of the computational scheme were developed using for basis functions (1) the no-flow uniform duct modes, (2) positive running uniform duct modes, with flow, and (3) positive and negative running uniform duct modes, with flow. The formulation using the no-flow modes was the most highly developed, and has advantages primarily due to relative computational simplicity. Results using the three methods are shown to converge to known solutions for several special cases, and the most significant check case is against low frequency, one dimensional results over the complete subsonic Mach number range. Development of the method is continuing, with emphasis on assessing the relative accuracy and efficiency of the three implementations.

Eversman, W.↗

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗

Simultaneous expansion and orthogonalization of measured modes for structure identification

Tests of large structures on-orbit will be performed with measurements at a relatively few structure points. Values for the unmeasured degrees of freedom (dofs) can be estimated based on measured dofs and analytical model dynamic information. These 'expanded' mode shapes are useful for optimal-update identification and damage location as well as test/analysis correlation. A new method of expansion for test mode shape vectors is developed from the orthogonal Procrustes problem from computational linear algebra. A subspace defined by the set of measured dofs is compared to a subspace defined by mode shapes from an analytical model of the structure. The method simultaneously expands and orthogonalizes the mode shape vectors. Two demonstration problems are used to compare the new method to current expansion techniques. One demonstration uses test data from a laboratory scale-model truss structure. Performance of the new method is comparable or superior to that of the previous expansion methods which require separate orthogonalization.

Smith, Suzanne Weaver↗

Suppression of Composition g -modes in Chemically Equilibrating Warm Neutron Stars

We investigate the impact of chemical equilibration and the resulting bulk viscosity on nonradial oscillation modes of warm neutron stars at temperatures up to T ≈ 5 MeV, relevant for protoneutron stars and neutron star postmerger remnants. In this regime, the relaxation rate of weak interactions becomes comparable to the characteristic frequencies of composition g-modes in the core, resulting in resonant damping. To capture this effect, we introduce the dynamical sound speed, a complex, frequency-dependent generalization of the adiabatic sound speed that encodes both the restoring force and the dissipative effects of bulk compression. Using realistic weak reaction rates and three representative equations of state, we compute the complex frequencies of composition g-modes with finite-temperature profiles. We find that bulk viscous damping becomes increasingly significant with temperature and can completely suppress composition g-modes. In contrast, the f-mode remains largely unaffected by bulk viscosity due to its nearly divergence-free character. Our results highlight the sensitivity of g-mode behavior to thermal structure, weak reaction rates, and the equation of state, and establish the dynamical sound speed as a valuable descriptor characterizing oscillation properties in dissipative neutron star matter.

Neutron stars↗

Running SINDA '85/FLUINT interactive on the VAX

Computer software as engineering tools are typically run in three modes: Batch, Demand, and Interactive. The first two are the most popular in the SINDA world. The third one is not so popular, due probably to the users inaccessibility to the command procedure files for running SINDA '85, or lack of familiarity with the SINDA '85 execution processes (pre-processor, processor, compilation, linking, execution and all of the file assignment, creation, deletions and de-assignments). Interactive is the mode that makes thermal analysis with SINDA '85 a real-time design tool. This paper explains a command procedure sufficient (the minimum modifications required in an existing demand command procedure) to run SINDA '85 on the VAX in an interactive mode. To exercise the procedure a sample problem is presented exemplifying the mode, plus additional programming capabilities available in SINDA '85. Following the same guidelines the process can be extended to other SINDA '85 residence computer platforms.

Simmonds, Boris↗

Structural and functional features of lysine acetylation of plant and animal tubulins

Abstract The study of the genome and the proteome of different species and representatives of distinct kingdoms, especially detection of proteome via wide‐scaled analyses has various challenges and pitfalls. Attempts to combine all available information together and isolate some common features for determination of the pathway and their mechanism of action generally have a highly complicated nature. However, microtubule (MT) monomers are highly conserved protein structures, and microtubules are structurally conserved from Homo sapiens to Arabidopsis thaliana . The interaction of MT elements with microtubule‐associated proteins and post‐translational modifiers is fully dependent on protein interfaces, and almost all MT modifications are well described except acetylation. Crystallography and interactome data using different approaches were combined to identify conserved proteins important in acetylation of microtubules. Application of computational methods and comparative analysis of binding modes generated a robust predictive model of acetylation of the ϵ‐amino group of Lys40 in α‐tubulins. In turn, the model discarded some probable mechanisms of interaction between elements of interest. Reconstruction of unresolved protein structures was carried out with modeling by homology to the existing crystal structure (PDBID: 1Z2B) from B. taurus using Swiss‐model server, followed by a molecular dynamics simulation. Docking of the human tubulin fragment with Lys40 into the active site of α‐tubulin acetyltransferase, reproduces the binding mode of peptidomimetic from X‐ray structure (PDBID: 4PK3).

Rayevsky, Alexey V.↗

A Study of the Terahertz Spectra of Crystalline Materials (Polyethylene, Poly(Vinylidene Fluoride) form II, and α-D-Glucose) Using NDDO Semiempirical Methods

Semiempirical quantum chemistry methods offer a very interesting compromise between accuracy and computational complexity. The low-frequency vibration modes of three crystalline materials, namely, polyethylene, poly(vinylidene fluoride) form II, and α-D-glucose, have been studied using the PM6 and PM7 Hamiltonians in order to assess the performance of NDDO methods in the interpretation of terahertz spectra. The results have been compared with the available experimental data and former calculations performed using alternative methods. NDDO calculations are in good qualitative or semiquantitative agreement with the experimentally observed terahertz spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Plastic shrinkage of mortars cured with a paraffin-based compound – Bimodal neutron/X-ray tomography study

Early-age drying (immediately after casting) of mortars and the corresponding plastic shrinkage were studied using bimodal neutron/X-ray computed tomography. This novel, correlative 3D imaging mode enabled studying simultaneously and without any source of spurious perturbation the water migration and loss processes together with the corresponding deformations due to plastic shrinkage. Bimodal imaging opens up new possibilities for studying dynamic processes of coupled water transport and deformations in porous solids. The measurements were carried out on model systems (cylindrical mortar specimens with height of 19 mm). The study focused on the effect of a paraffin-based curing compound. Our results confirm that when the curing compound was applied directly onto the drying surface in a sufficient amount, both the evaporation rate and the rate of vertical displacement (settlement) were substantially reduced. The results shed a new light on the mechanisms of plastic shrinkage and the action of curing compounds.

36 MATERIALS SCIENCE↗

Parallel exponential time differencing methods for geophysical flow simulations

Two ocean models are considered for geophysical flow simulations: the multilayer shallow water equations and the multilayer primitive equations. For the former, we investigate the parallel performance of exponential time differencing (ETD) methods, including exponential Rosenbrock–Euler, ETD2wave, and B-ETD2wave. For the latter, we take advantage of the splitting of barotropic and baroclinic modes and propose a new two-level method in which an ETD method is applied to solve the fast barotropic mode. Furthermore, these methods could improve the computational efficiency of numerical simulations because ETD methods allow for much larger time step sizes than traditional explicit time-stepping techniques that are commonly used in existing computational ocean models. Several standard benchmark tests for ocean modeling are performed and comparison of the numerical results demonstrates a great potential of applying the parallel ETD methods for simulating real-world geophysical flows.

54 ENVIRONMENTAL SCIENCES↗

Predicting Tensile Stretchability of Trimmed AA6111-T4 Sheets

An integrated manufacturing process simulation framework has been developed to predict the trimmed edge tensile stretchability of AA6111-T4 sheets by incorporating the burr geometry, damage, and plastic strain from trimming simulations into subsequent tensile stretchability simulations. The influence of the trimming die clearances on the predicted tensile stretching ductility (stretchability) is studied and quantitatively compared with experimental measurements. Stretchability is found to decrease with increasing cutting clearances, and simulation results have successfully captured experimentally observed edge crack initiation and failure mode variations for different trimming clearances. Subsequent computational sensitivity studies reveal that while deburring of previously trimmed edges has little influence on tensile stretchability, removal of trimmed edge initial plastic strain may significantly enhance the subsequent trimmed edge stretchability.

Aluminum alloys↗

Combined experimental and theoretical assessments of the lattice dynamics and optoelectronics of TaON and Ta{sub 3}N{sub 5}

Presented herein is a detailed discussion of the properties of the lattice dynamic and optoelectronic properties of tantalum(V) oxynitride (TaON) and tantalum(V) nitride (Ta{sub 3}N{sub 5}), from experimental and theoretical standpoint. The active Raman and infra red (IR) frequencies of TaON and Ta{sub 3}N{sub 5} were measured using confocal Raman and Fourier Transform Infrared spectroscopies (FTIR) and calculated using the linear response method within the density functional perturbation theory (DFPT). The detailed study leads to an exhaustive description of the spectra, including the symmetry of the vibrational modes. Electronic structures of these materials were computed using DFT within the range-separated hybrid HSE06 exchange–correlation formalism. Electronic and ionic contributions to the dielectric constant tensors of these materials were obtained from DFPT within the linear response method using the PBE functional. Furthermore, effective mass of photogenerated holes and electrons at the band edges of these compounds were computed from the electronic band structure obtained at the DFT/HSE06 level of theory. The results suggest that anisotropic nature in TaON and Ta{sub 3}N{sub 5} is present in terms of dielectric constant and effective masses. - Graphical abstract: Detailed investigation has been conducted from combined experimental and theoretical approaches on Raman and IR spectroscopies, electronic structures, dielectric constants and effective masses of TaON and Ta{sub 3}N{sub 5}. - Highlights: • Crystal structures of TaON and Ta{sub 3}N{sub 5} are discussed based on XRD and DFT calculation. • Raman and IR spectra of TaON and Ta{sub 3}N{sub 5} materials are measured and computed by DFPT method. • Optoelectronic properties of TaON and Ta{sub 3}N{sub 5} are discussed. • Dielectric constant and effective masses of TaON and Ta{sub 3}N{sub 5} are calculated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Internal Conversion and Vibrational Energy Redistribution in Chlorophyll A

For this research, we have computationally investigated the role of intramolecular vibrational modes in determining nonradiative relaxation pathways of photoexcited electronic states in isolated chlorophyll A (ChlA) molecules. To simulate the excited state relaxation from the initially excited Soret state to the lowest excited state Q y , the approach of nonadiabatic excited state molecular dynamics has been adopted. The intramolecular vibrational energy relaxation and redistribution that accompany the electronic internal conversion process is followed by analyzing the excited state trajectories in terms of the ground state equilibrium normal modes. The time dependence of the normal mode velocities is determined by projecting instantaneous Cartesian velocities onto the normal mode vectors. Our analysis of the time evolution of the average mode energies uncovers that only a small subset of the medium-to-high frequency normal modes actively participate in the electronic relaxation processes. These active modes are characterized by the highest overlap with the nonadiabatic coupling vectors (NACRs) during the electronic transitions. Further statistical analysis of the nonadiabatic transitions reveals that the electronic and vibrational energy relaxation occurs via two distinct pathways with significantly different time scales on which the hopping from Soret to Q x occurs thereby dictating the overall dynamics. Furthermore, the NACRs corresponding to each of the transitions have been consistently found to be predominantly similar to a set of normal modes that vary depending upon the transition and the identified categories. Each pathway exhibits a differential time scale of energy transfer and also a differential set of predominant active modes. Lastly, our present analysis can be considered as a general approach allowing identification of a reduced subset of specific vibrational coordinates associated with nonradiative relaxation pathways. Therefore, it represents an adequate prior strategy that can particularly facilitates mixed quantum-classical approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct simulation Monte Carlo modeling of relaxation processes in polyatomic gases

Relaxation processes of polyatomic molecules are modeled and implemented in an in-house Direct Simulation Monte Carlo code in order to enable the simulation of atmospheric entry maneuvers at Mars and Saturn’s Titan. The description of rotational and vibrational relaxation processes is derived from basic quantum-mechanics using a rigid rotator and a simple harmonic oscillator, respectively. Strategies regarding the vibrational relaxation process are investigated, where good agreement for the relaxation time according to the Landau-Teller expression is found for both methods, the established prohibiting double relaxation method and the new proposed multi-mode relaxation. Differences and applications areas of these two methods are discussed. Consequently, two numerical methods used for sampling of energy values from multi-dimensional distribution functions are compared. The proposed random-walk Metropolis algorithm enables the efficient treatment of multiple vibrational modes within a time step with reasonable computational effort. The implemented model is verified and validated by means of simple reservoir simulations and the comparison to experimental measurements of a hypersonic, carbon-dioxide flow around a flat-faced cylinder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time-dependent, multimode interaction analysis of the gyroklystron amplifier

In this paper, a time-dependent multimode nonlinear analysis for the gyroklystron amplifier has been developed by extending the analysis of gyrotron oscillators by employing the self-consistent approach. The nonlinear analysis developed here has been validated by taking into account the reported experimental results for a 32.3 GHz, three cavity, second harmonic gyroklystron operating in the TE{sub 02} mode. The analysis has been used to estimate the temporal RF growth in the operating mode as well as the nearby competing modes. Device gain and bandwidth have been computed for different drive powers and frequencies. The effect of various beam parameters, such as beam voltage, beam current, and pitch factor, has also been studied. The computational results have estimated the gyroklystron saturated RF power ∼319 kW at 32.3 GHz with efficiency ∼23% and gain ∼26.3 dB with device bandwidth ∼0.027% (8 MHz) for a 70 kV, 20 A electron beam. The computed results are found to be in agreement with the experimental values within 10%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear response of a Hall magnetic drift wave for verification of Hall MHD algorithms

Numerical implementations of Hall magnetohydrodynamics (Hall MHD) can be challenging due to the nonlinear multidimensional nature of the Hall term. Here, a model problem is presented that couples the hydrodynamic motion of the plasma to Hall MHD evolution of the magnetic field. The Hall MHD equations are linearized about unperturbed solutions in both cylindrical and Cartesian coordinates in two dimensions. The magnetic field is assumed to lie in the ignorable direction, and the linear response about the unperturbed solution is considered. The resulting ordinary differential equation is used to numerically compute the eigenfunctions and eigenfrequencies of the mode. The resulting eigenfunctions do not make the local wave approximation but are instead global solutions that depend on the spatial dependence of the unperturbed Alfvén speed. Hall MHD simulations are then performed in the Ares multiphysics code and shown to agree with the predicted phase velocities of the wave, and the simulated solutions are shown to numerically converge to the semianalytic modes. By varying the background density of the plasma (and correspondingly, the ion inertial length), the importance of Hall physics can be varied. This allows the test problem to transition from the classical MHD limit to the extreme Hall MHD limit. Furthermore, this problem is a useful tool for the verification of Hall MHD routines implemented in various codes, and the robustness of a routine can be tested in regimes in which Hall physics is dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PRIMORDIAL GRAVITATIONAL WAVE DETECTABILITY WITH DEEP SMALL-SKY COSMIC MICROWAVE BACKGROUND EXPERIMENTS

We use the Bayesian estimation on direct T - Q - U cosmic microwave background (CMB) polarization maps to forecast errors on the tensor-to-scalar power ratio r, and hence on primordial gravitational waves, as a function of sky coverage f{sub sky}. This map-based likelihood filters the information in the pixel-pixel space into the optimal combinations needed for r detection for cut skies, providing enhanced information over a first-step linear separation into a combination of E, B, and mixed modes, and ignoring the latter. With current computational power and for typical resolutions appropriate for r detection, the large matrix inversions required are accurate and fast. Our simulations explore two classes of experiments, with differing bolometric detector numbers, sensitivities, and observational strategies. One is motivated by a long duration balloon experiment like Spider, with pixel noise {proportional_to}{radical}(f{sub sky}) for a specified observing period. This analysis also applies to ground-based array experiments. We find that, in the absence of systematic effects and foregrounds, an experiment with Spider-like noise concentrating on f{sub sky} {approx} 0.02-0.2 could place a 2{sigma}{sub r} Almost-Equal-To 0.014 boundary ({approx}95% confidence level), which rises to 0.02 with an l-dependent foreground residual left over from an assumed efficient component separation. We contrast this with a Planck-like fixed instrumental noise as f{sub sky} varies, which gives a Galaxy-masked (f{sub sky} = 0.75) 2{sigma}{sub r} Almost-Equal-To 0.015, rising to Almost-Equal-To 0.05 with the foreground residuals. Using as the figure of merit the (marginalized) one-dimensional Shannon entropy of r, taken relative to the first 2003 WMAP CMB-only constraint, gives -2.7 bits from the 2012 WMAP9+ACT+SPT+LSS data, and forecasts of -6 bits from Spider (+ Planck); this compares with up to -11 bits for CMBPol, COrE, and PIXIE post-Planck satellites and -13 bits for a perfectly noiseless cosmic variance limited experiment. We thus confirm the wisdom of the current strategy for r detection of deeply probed patches covering the f{sub sky} minimum-error trough with balloon and ground experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-linear realizations of conformal symmetry and effective field theory for the pseudo-conformal universe

The pseudo-conformal scenario is an alternative to inflation in which the early universe is described by an approximate conformal field theory on flat, Minkowski space. Some fields acquire a time-dependent expectation value, which breaks the flat space so(4,2) conformal algebra to its so(4,1) de Sitter subalgebra. As a result, weight-0 fields acquire a scale invariant spectrum of perturbations. The scenario is very general, and its essential features are determined by the symmetry breaking pattern, irrespective of the details of the underlying microphysics. In this paper, we apply the well-known coset technique to derive the most general effective lagrangian describing the Goldstone field and matter fields, consistent with the assumed symmetries. The resulting action captures the low energy dynamics of any pseudo-conformal realization, including the U(1)-invariant quartic model and the Galilean Genesis scenario. We also derive this lagrangian using an alternative method of curvature invariants, consisting of writing down geometric scalars in terms of the conformal mode. Using this general effective action, we compute the two-point function for the Goldstone and a fiducial weight-0 field, as well as some sample three-point functions involving these fields.

79 ASTRONOMY AND ASTROPHYSICS↗