SEARCH · Search NASA

Results for “force matrix”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A generalized 3D elastic model for nanoscale, self-assembled oxide-metal thin films with pillar-in-matrix configurations

In recent years, functional oxide-metal based vertically aligned nanocomposite (VAN) thin films have gained interest due to their intriguing physical properties and multifunctionalities stemming from the complex interactions between the two phases in the film and the substrate. In this work, we develop a model for studying the energetics of these thin film systems, including the effects of both lattice mismatch and capillary forces due to interface curvature. Each phase is incorporated into the model using a phase indicator function, and we introduce the capillary forces as body forces using a vector density representation of the interface. The model is implemented using the finite element method to study the deformation of the thin film which is composed of Au nanopillars embedded in a La 0.7 Sr 0.3 MnO 3 (LSMO) matrix on an SrTiO 3 (STO) substrate. Further, the results suggest that the total energy is lowest for random configurations of pillars compared to ordered square and hexagonal lattice configurations, consistent with the random distribution of pillars found in experiments. Furthermore, we find that the interfacial energy dominates the total energy of each configuration, suggesting that interfacial energy in the system is an important design parameter for nanocomposite growth, along with the lattice mismatch.

36 MATERIALS SCIENCE↗

NASA NTRS1994

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗

NASA NTRS2022

The Turning Point of the Aerosol Era

Over the CMIP6 historical period (1850–2014), aerosols provided the largest negative forcing compared to all other climate forcings via their ability to absorb or scatter solar radiation and alter clouds. Aerosols played an important role in counterbalancing warming by greenhouse gases (GHGs). Here we study aerosol forcing trends in the CMIP6 simulations of the NASA Goddard Institute for Space Studies (GISS) ocean-atmosphere ModelE version 2.1 (GISS-E2.1-G) using a fully coupled atmospheric composition configuration, including interactive gas-phase chemistry, and either an aerosol microphysical (MATRIX) or a mass-based aerosol (OMA) module. Simulations of the CMIP6 historical period are analyzed as well as four Shared Socioeconomic Pathway (SSP) future scenarios for 2015–2100: SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5. The main conclusion of this study is that aerosol forcing in the GISS model has reached its turning point, switching from globally increasing to a decreasing trend in the first decade of the 21st century. This result is robust, independent of which aerosol module or SSP scenario is used. Non-linear aerosol-cloud interactions dominate as a forcing agent over aerosol-radiation interactions. Aerosols' ability to counterbalance GHG forcing on the global scale is today at a level comparable to that at the beginning of the last century. In the 1980s, the decade of largest global aerosol loads, aerosols balanced up to 80% of GHG forcing. As a consequence, global warming of the last decades, which is primarily driven by greenhouse gases, has been augmented by the effect of decreasing aerosol cooling in our model. By the end of this century, following the SSP scenarios, aerosols will only counterbalance 0%–20% of GHG forcing, depending on model and on scenario.

Aerosol forcing↗

DOE OSTI2023

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)—a selected configuration interaction technique for vibrational structure theory—has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm−1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Chemistry↗

NASA NTRS1972

Dynamics of a spinning space station with a counterweight connected by multiple cables.

A unique and relatively simple approach is presented for obtaining the linearized equations of motion. A conservative system consisting of two rigid bodies connected by any number of massless cables with linear axial stiffness is assumed. The cable forces and torques are expanded in a Taylor series about the equilibrium values of the system coordinates which results in a cable stiffness matrix. The method of obtaining the equilibrium values of the coordinates is discussed and results are presented. The range of validity of the linear model is determined by comparing results with a digital simulation of the nonlinear system.

Nixon, D. D.↗

NASA NTRS1975

Flapping response characteristics of hingeless rotor blades by a gereralized harmonic balance method

Linearized equations of motion for the flapping response of flexible rotor blades in forward flight are derived in terms of generalized coordinates. The equations are solved using a matrix form of the method of linear harmonic balance, yielding response derivatives for each harmonic of the blade deformations and of the hub forces and moments. Numerical results and approximate closed-form expressions for rotor derivatives are used to illustrate the relationships between rotor parameters, modeling assumptions, and rotor response characteristics. Finally, basic hingeless rotor response derivatives are presented in tabular and graphical form for a wide range of configuration parameters and operating conditions.

Peters, D. A.↗

NASA NTRS1978

The role of rotor impedance in the vibration analysis of rotorcraft, part 4

A method for a strongly idealized case of vertical excitation and for rolling and pitching moment excitation of a four bladed hingeless rotor on an up-focussing flexible mount is developed. The aeroelastic rotor impedances are computed directly with a finite blade element method that includes aerodynamics. The rotor impedance matrix for three or more blades is determined from the root moment impedance for a single blade by a simple multiblade transformation rule. Force and moment amplitudes transferred from the rotor to support are found to be critically dependent on the support dynamics.

Hohenemser, K. H.↗

NASA NTRS1978

The role of rotor impedance in the vibration analysis of rotorcraft

In an improved method which retains the advantage of separate treatment of rotor and airframe, the rotor impedance is used to correct the input to the airframe. This improved method is illustrated for a strongly idealized case of vertical excitation and then for rolling and pitching moment excitation of a four bladed hingeless rotor on an up-focussing flexible mount. Contrary to the usual approach that represents aeroelastic blade motions by a series of normal blade modes in vacuum, the aeroelastic rotor impedances are computed directly with a finite blade element method that includes aerodynamics. The rotor impedance matrix for three or more blades is determined from the root moment impedance for a single blade by a simple multiblade transformation rule. Force and moment amplitudes transferred from the rotor to the support are found to be critically dependent on the support dynamics.

Hohenemser, K. H.↗

NASA NTRS1987

Use of the generalized transmission error in the equations of motion of gear systems

The vibratory excitation arising from a gear pair is widely recognized to be a consequence of the nonuniform transmission of motion by the gear pair. In the preceding paper in this issue, it is shown that a three-component transmission error is required to describe the nonuniform transmission of motion by bevel gears for vibration excitation characterization purposes. The expression for the three-component transmission error derived in that paper is combined in the present paper with an analysis of the mesh forces and mesh elasticity to yield an equation of constraint involving the six degree-of-freedom unknown vibratory displacements of the gear shaft centerliners, the three unknown components of the generalized force transmitted by the mesh, and the geometric deviations of the tooth running surfaces from perfect involute surfaces which are assumed known. This matrix equation can be combined with the equations of motion of a gear system to predict the vibratory response of the system to the generalized transmission error excitation arising from meshing gear pairs within the system.

Mark, W. D.↗

NASA NTRS1994

Thermoviscoplastic analysis of fibrous periodic composites by the use of triangular subvolumes

The non-linear viscoplastic behavior of fibrous periodic composites is analyzed by discretizing the unit cell into triangular subvolumes. A set of these subvolumes can be configured by the analyst to construct a representation for the unit cell of a periodic composite. In each step of the loading history the total strain increment at any point is governed by an integral equation which applies to the entire composite. A Fourier series approximation allows the incremental stresses and strains to be determined within a unit cell of the periodic lattice. The non-linearity arising from the viscoplastic behavior of the constituent materials comprising the composite is treated as a fictitious body force in the governing integral equation. Specific numerical examples showing the stress distributions in the unit cell of a fibrous tungsten/copper metal-matrix composite under viscoplastic loading conditions are given. The stress distribution resulting in the unit cell when the composite material is subjected to an overall transverse stress loading history perpendicular to the fibers is found to be highly heterogeneous, and typical homogenization techniques based on treating the stress and strain distributions within the constituent phases as homogeneous result in large errors under inelastic loading conditions.

Walker, Kevin P.↗

NASA NTRS1989

An integrating matrix formulation for buckling of rotating beams including the effects of concentrated masses

The integrating matrix technique of computational mechanics is extended to include the effects of concentrated masses. The stability of a flexible rotating beam with discrete masses is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. In this problem, the beam is subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam. To determine the eigenvalues from which stability is assessed, the differential equations of motion are solved numerically by combining the extended integrating matrix method with an eigenanalysis. Stability boundaries for a discrete mass representation of a uniform beam are shown to asymptotically approach the stability boundaries for the corresponding continuous mass beam as the number of concentrated masses is increased. An error in the literature is also noted for the discrete mass problem concerning the behavior of the critical rotational speed for inplane buckling as the radius of rotation of the clamped end of the beam is reduced.

Lakin, William D.↗

NASA NTRS1990

The influence of buoyant forces and volume fraction of particles on the particle pushing/entrapment transition during directional solidification of Al/SiC and Al/graphite composites

Directional solidification experiments in a Bridgman-type furnace were used to study particle behavior at the liquid/solid interface in aluminum metal matrix composites. Graphite or silicon-carbide particles were first dispersed in aluminum-base alloys via a mechanically stirred vortex. Then, 100-mm-diameter and 120-mm-long samples were cast in steel dies and used for directional solidification. The processing variables controlled were the direction and velocity of solidification and the temperature gradient at the interface. The material variables monitored were the interface energy, the liquid/particle density difference, the particle/liquid thermal conductivity ratio, and the volume fraction of particles. These properties were changed by selecting combinations of particles (graphite or silicon carbide) and alloys (Al-Cu, Al-Mg, Al-Ni). A model which consideres process thermodynamics, process kinetics (including the role of buoyant forces), and thermophysical properties was developed. Based on solidification direction and velocity, and on materials properties, four types of behavior were predicted. Sessile drop experiments were also used to determine some of the interface energies required in calculation with the proposed model. Experimental results compared favorably with model predictions.

Stefanescu, Doru M.↗

DOE OSTI2021

Enhanced resistance to decay of imprinted nanopatterns in thin films by bare nanoparticles compared to polymer-grafted nanoparticles

We extend a previous study on the influence of nanoparticles on the decay of nanoimprinted polymer film patterns to compare the effects of “bare” silica (SiO 2 ) nanoparticles and SiO 2 nanoparticles with grafted polymer layers having the same chemical composition as the polymer matrix. This method involves nanoimprinting substrate-supported polymer films using a pattern replicated from a digital versatile disc (DVD), and then annealing the patterned polymer nanocomposite films at elevated temperatures to follow the decay of the topographic surface pattern with time by atomic force microscopy imaging after quenching. We quantified the relaxation of the pattern height (“slumping”) and determined the relaxation time τ for this pattern decay process as a function of nanoparticle filler type and concentration to determine how nanoparticle additives influence relative film stability. Attractive interactions between the bare nanoparticles and the polymer matrix significantly enhance the thermal resilience of the nanopatterns to decay, compared to those of the particle brushes, wherein the particle core interactions are screened from the matrix via the brush layer. A novel aspect of this method is that it readily lends itself to in situ film relaxation measurements in a manufacturing context. We observe that the relaxation time of the pattern relaxation exhibits entropy–enthalpy compensation in the free energy parameters governing the pattern relaxation process as a function of temperature, irrespective of the NP system used, consistent with our previous experimental and computational studies.

36 MATERIALS SCIENCE↗

NASA NTRS2004

Combining Hard with Soft Materials in Nanoscale Under High-Pressure High-Temperature Conditions

Nano-composites with a primary nanocrystalline ceramic matrix and a secondary nanocrystalline material (metal or semiconductor) were synthesized by infiltration of an appropriate liquid into ceramic compacts under pressures of up to 8 GPa and temperatures of up to 2000 K. The purpose of our work is to obtain nanocomposites which constitute homoger?ous mixtures of two phases, both forming nano- grains of about 10 nm in size. The high pressure is used to bring the porosity of the compacted powders down to the nano-scale and force a given liquid into the nano-sized pores. The advantage of the infiltration technique is that, in a single, continuous process, we start with a nanocrystalline powder, compress it to form the matrix of the composite, and crystallize and/or synthesize a second nanomaterial in the matrix pores. The key limitation of this technology is, that the pores in the matrix need to stay open during the entire process of infiltration. Thus the initial powder should form a rigid skeleton, otherwise the so-called self-stop process can limit cr block a further flow of the liquid phase and hinder the process of the composite formation. Therefore powders of only very hard ceramic materials like diamond, Sic, or Alz03, which can withstand a substantial external load without undesired deformation, can be used as the primary phase. With this technique, using diamond and S i c ceramic powders infiltrated by liquid metals (AI, Zn, Sn, Ag, Au) and semiconductors (Si, Ge, GaAs, CdTe), we obtained nano-composites with the grain size in the range of 10 - 30 nm. Our work addresses the key problem in manufacturing bulk nanocrystalline materials, i.e. preservation of nano-scale during the fabrication process. In this paper we discuss basic technical and methodological problems associated with nano-infiltration based on the results obtained for Zn-Sic composites.

Palosz, B.↗

DOE OSTI2025

Theoretical and experimental quantification of Suzuki segregation enthalpy and strengthening mechanisms in a binary alloy

Solute segregation to planar defects in metallic alloys has been shown to drastically alter mechanical properties. While various works using first-principles and thermodynamic calculations have studied the fundamental driving forces for solute segregation via the Suzuki criterion, planar defect energy, or a comparison of energies of the HCP-like phase and FCC matrix, a quantitative experimental and computational comparison of equilibrium composition and segregation enthalpies has not yet been reported. In this work, we predict the equilibrium composition and segregation enthalpy to intrinsic stacking faults in a Ni-60Co (at.%) alloy and compare the results to two independent experimental methods. We observed that Co segregates to the innermost two planes of the intrinsic stacking fault, and we found that the experimental segregation enrichment, measured from transmission electron microscopy energy dispersive X-ray spectroscopy, of the faults is 6.8 at.% Co, which is 2.2 at.% less than the predicted value at the same temperature. We also find that the segregation enthalpy measured from the composition profile is −21.1 ± 6.4 meV/atom and separately from differential scanning calorimetry segregation enthalpy is −33.2 meV/atom, whereas the predicted enthalpy is −31 ± 1 meV/atom. Based on these results, we determine that segregation occurs very rapidly, within 8 min at temperatures as low as 36% of the homologous solidus temperature. Furthermore, this analysis provides an overview of the possible dislocation mechanisms responsible for strengthening effects due to solute segregation, and concludes that changes in room temperature hardness from local phase transformation is likely tied to post-segregation room temperature equilibrium partial separation distance.

Ab initio calculation↗

DOE OSTI2021

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Multistep Coulomb excitation of Ni 64 : Shape coexistence and nature of low-spin excitations

Here, the structure of 64 Ni, the heaviest stable Ni isotope, has been investigated via high-statistics, multistep safe Coulomb excitation to search for shape coexistence, a phenomenon recently observed in neutron-rich 66 Ni and 70 Ni as well as in doubly magic, N = 40, 68 Ni. The study was motivated by recent, state-of-the-art Monte Carlo shell-model calculations (MCSM), where a Hamiltonian with effective interactions incorporating the monopole tensor force predicts the existence of shape coexistence, also in the lower-mass 62,64 Ni isotopes. A set of transition and static E2 matrix elements for both yrast and near-yrast structures was extracted from the differential Coulomb excitation cross sections. From comparisons between the new results and MCSM as well as other shell-model calculations, a clearer picture of the structure of 64 Ni emerges. Specifically, the low-spin states are shown to be dominated by proton and neutron excitations mainly within the fp shell, with minimal contribution from the g 9/2 shape-driving neutron orbital. The agreement between experimental data and MCSM results indicates a small oblate deformation for the $0$$^{+}_{2}$ level and a spherical shape for the $0$$^{+}_{3}$ state. In addition, the small upper limit determined for the B(E2) probability of a transition associated with the decay of the recently observed 3463-keV, $0$$^{+}_{4}$ state agrees with its proposed assignment to a prolate shape, herewith providing first evidence for triple shape coexistence in a stable Ni isotope.

59 ≤ A ≤ 89↗

NASA NTRS1987

Method for the determination of the three-dimensional aerodynamic field of a rotor-stator combination to compressible flow

Using the lifting surface theory and the acceleration potential method for the flow field of an axial turbocompressor stage, a recursive and a direct method are presented that make use of the eigenfunction solutions of the isolated rotor and stator to solve for the rotor-stator interaction problem. The net pressure distribution on the rotor and stator blades is represented by modified Birnbaum series, whose coefficients are determined using a matrix procedure and satisfying the boundary conditions on the surface of the blades. The relation between the matrix operators of the recursive and the direct methods is also shown. Expressions have been given for the blade circulation, the axial and tangential forces on the blade, the rotor power required, and the induced upwash velocity of the stage.

Ramachandra, Sridhar M.↗