Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear embedding”

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 109 records · Page 6

Application of the locally self-consistent embedding approach to the Anderson model with non-uniform random distributions

Highlights: • Typical Medium Theory (TMT) for the Anderson Localization. • Locally Self-Consistent Multiple Scattering Method (LSMS) for Random Disordered Systems. • Linear Scaling Computational Method for Random Disordered Systems. We apply the recently developed embedding scheme for the locally self-consistent method to random disorder electrons systems. The method is based on the locally self-consistent multiple scattering theory and the typical medium theory. The locally self-consistent multiple scattering theory divides a system into many small designated local interaction zones. The subsystem within each local interaction zone is embedded in a self-consistent field from the typical medium theory. This approximation allows the study of random systems with large numbers of sites. We present results for the three dimensional Anderson model with different random disorder potential distributions. Using the typical density of states as an indicator of Anderson localization, we find that the method can capture the localization for commonly studied disorder potentials. These include the uniform distribution, the Gaussian distribution, and even the unbounded Cauchy distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large-scale harmonic balance simulations with Krylov subspace and preconditioner recycling

The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. Further, the method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.

97 MATHEMATICS AND COMPUTING↗

Design of Cylindrical Implosion Experiments to Demonstrate Scale-Invariant Rayleigh-Taylor Instability Growth

Radiation-hydrodynamics simulations are used to design laser-driven cylindrical implosion experiments to directly measure hydrodynamic instability growth in convergent geometry. Designs for two different size targets, varying in radial dimension by a factor of three, are presented. A set of beam pointings and powers are identified for each scale design that result in a nearly axially uniform implosion of an embedded marker layer. The implosion trajectories are shown to be scale-invariant between designs, with nearly identical scaled acceleration profiles. Finally, linear theory and radiation-hydrodynamics simulations predict that Rayleigh-Taylor instability growth of an azimuthal perturbation, machined on the inner surface of the embedded marker, is also scale-invariant between designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Finite-strain large-deflection elastic-viscoplastic finite-element transient response analysis of structures

A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.

Rodal, J. J. A.↗

Detection of Delaminations in Composite Beams Using Piezoelectric Sensors

This paper investigates the feasibility of a proposed technique for detecting delamination using piezoelectric layers or patches embedded or bonded to a composite structure. Variations in the voltage generated in the piezoelectric layers indicates the presence and location of delamination, when the structure is excited either externally or via actuators. The theoretical foundations of a method for predicting the dynamic response of delaminated composite beams with piezoelectric layers are described. The governing equations are presented for the case of external vibroacoustic excitation, as well as, for the case of locally induced vibrations by some of the embedded piezoelectric elements. An exact solution is developed within the limits of linear laminate theory. Applications illustrate the feasibility of delamination detection in cantilever beams. The results illustrate that the proposed technique may provide accurate detection of the presence, size, and location of a delamination.

Saravanos, Dimitris A.↗

Growth and Morphology of Phase Separating Supercritical Fluids

The scientific objective is to study the relation between the morphology and the growth kinetics of domains during phase separation. We know from previous experiments performed near the critical point of pure fluids and binary liquids that there are two simple growth laws at late times. The 'fast' growth appears when the volumes of the phases are nearly equal and the droplet pattern is interconnected. In this case the size of the droplets grows linearly in time. The 'slow' growth appears when the pattern of droplets embedded in the majority phase is disconnected. In this case the size of the droplets increases in proportion to time to the power 1/3. The volume fraction of the minority phase is a good candidate to determine this change of behavior. All previous attempts to vary the volume fraction in a single experimental cell have failed because of the extreme experimental difficulties.

Hegseth, John↗

Ultrafast vibrational dynamics of a trigonal planar anionic probe in ionic liquids (ILs): A two-dimensional infrared (2DIR) spectroscopic investigation

A major impediment limiting the widespread application of ionic liquids (ILs) is their high shear viscosity. Incorporation of a tricyanomethanide (TCM - ) anion in ILs leads to low shear viscosity and improvement of several characteristics suitable for large scale applications. However, properties including interactions of TCM - with the local environment and dynamics of TCM - have not been thoroughly investigated. Herein, we have studied the ultrafast dynamics of TCM - in several imidazolium ILs using linear IR and two-dimensional infrared spectroscopy techniques. The spectral diffusion dynamics of the CN stretching modes of TCM - in all ILs exhibit a nonexponential behavior with a short time component of ~2 ps and a long time component spanning ~9 ps to 14 ps. The TCM - vibrational probe reports a significantly faster relaxation of ILs compared to those observed previously using linear vibrational probes, such as thiocyanate and selenocyanate. Our results indicate a rapid relaxation of the local ion-cage structure embedding the vibrational probe in the ILs. The faster relaxation suggests that the lifetime of the local ion-cage structure decreases in the presence of TCM - in the ILs. Linear IR spectroscopic results show that the hydrogen-bonding interaction between TCM - and imidazolium cations in ILs is much weaker. Shorter ion-cage lifetimes together with weaker hydrogen-bonding interactions account for the low shear viscosity of TCM - based ILs compared to commonly used ILs. In addition, this work demonstrates that TCM - can be used as a potential vibrational reporter to study the structure and dynamics of ILs and other molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the conformal symmetry of exceptional scalar theories

The DBI and special galileon theories exhibit a conformal symmetry at unphysical values of the spacetime dimension. We find the Lagrangian form of this symmetry. The special conformal transformations are non-linearly realized on the fields, even though conformal symmetry is unbroken. Commuting the conformal transformations with the extended shift symmetries, we find new symmetries, which when taken together with the conformal and shift symmetries close into a larger algebra. For DBI this larger algebra is the conformal algebra of the higher dimensional bulk in the brane embedding view of DBI. For the special galileon it is a real form of the special linear algebra. We also find the Weyl transformations corresponding to the conformal symmetries, as well as the necessary improvement terms to make the theories Weyl invariant, to second order in the coupling in the DBI case and to lowest order in the coupling in the special galileon case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

Comparison of Damage Models for Predicting the Non-Linear Response of Laminates Under Matrix Dominated Loading Conditions

Five models for matrix damage in fiber reinforced laminates are evaluated for matrix-dominated loading conditions under plane stress and are compared both qualitatively and quantitatively. The emphasis of this study is on a comparison of the response of embedded plies subjected to a homogeneous stress state. Three of the models are specifically designed for modeling the non-linear response due to distributed matrix cracking under homogeneous loading, and also account for non-linear (shear) behavior prior to the onset of cracking. The remaining two models are localized damage models intended for predicting local failure at stress concentrations. The modeling approaches of distributed vs. localized cracking as well as the different formulations of damage initiation and damage progression are compared and discussed.

Schuecker, Clara↗

Control of a Thrust Alignment Table for Modeling the Coning Dynamics of a Spinning Spacecraft with a Follower Force

This document presents a system controlling the motion of a spherical air bearing used in the modeling of spacecraft dynamics and controls in a laboratory environment. The system is part of the Spinning Rocket Simulator (SRS), used to simulate the coning of spacecraft during a thrusting stage. The reaction force at the spherical air bearing supporting the spacecraft model must coincide with the thrust axis of the model for proper simulation. Therefore, the bearing is translated in a circular path to introduce a centrifugal force. This horizontal force along with the gravitational reaction force at the bearing combines to simulate the direction of the spacecraft's thrust force. The control system receives attitude information from the spacecraft model via a laser beam embedded in the model that impinges on a photosensitive array. The non-linear system is controlled using high-speed lookup tables and digital techniques. A vector-controlled motor and a stepper motor are given the necessary signals to accurately control the turntable and platform supporting the air bearing. Preliminary performance data is presented. Mechanical elements of the table and platform are described in detail. A wireless (RF) data path for all devices on the spacecraft model to an off-table command computer is also described.

Halsmer, Dominic↗

NIR Imaging Spectroscopy of the Inner Few Arcseconds of NGC 4151 with OSIRIS at Keck

We present H- and K-band data from the inner arcsecond of the Seyfert 1.5 galaxy NGC 4151 obtained with the adaptive optics assisted near-infrared imaging field spectrograph OSIRIS at the Keck Observatory. The angular resolution is about a few parsecs on-site and thus competes easily with optical images taken previously with the Hubble Space Telescope. We present the morphology and dynamics of most species detected but focus on the morphology and dynamics of the narrow line region (as traced by emission of [FeII]λ1.644 μm), the interplay between plasma ejected from the nucleus (as traced by 21 cm continuum radio data) and hot H2 gas and characterize the detected nuclear HeIλ2.058 μm absorption feature as a narrow absorption line (NAL) phenomenon. Emission from the narrow line region (NLR) as traced by [FeII] reveals a biconical morphology and we compare the measured dynamics in the [FeII] emission line with models proposing acceleration of gas in the NLR and simple ejection of gas into the NLR. In the inner 2.5 arcseconds the acceleration model reveals a better fit to our data than the ejection model.We also see evidence that the jet very locally enhances emission in [FeII] at certain positions in our field-of-view such that we were able to distinct the kinematics of these clouds from clouds generally accelerated in the NLR. Further, the radio jet is aligned with the bicone surface rather than the bicone axis such that we assume that the jet is not the dominant mechanism responsible for driving the kinematics of clouds in the NLR. The hot H2 gas is thermal with a temperature of about 1700 K. We observe a remarkable correlation between individual H2 clouds at systemic velocity with the 21 cm continuum radio jet. We propose that the radio jet is at least partially embedded in the galactic disk of NGC 4151 such that deviations from a linear radio structure are invoked by interactions of jet plasma with H2 clouds that are moving into the path of the jet because of rotation of the galactic disk of NGC 4151. Additionally, we observe a correlation of the jet as traced by the radio data, with gas as traced in Brγ and H2, at velocities between systemic and +/- 200 km/s at several locations along the path of the jet. The HeIλ2.058 μm line in NGC 4151 appears in emission with a blueshifted absorption component from an outflow. The emission (absorption) component has a velocity offset of 10 km/s (-280 km/s) with a Gaussian (Lorentzian) full-width (half-width) at half maximum of 160 km/s (440 km/s). The absorption component remains spatially unresolved and its kinematic measures differ from that of UV resonance absorption lines. From the amount of absorption we derive a lower limit of the HeI 2S column density of 1 × 10(exp 14) cm−2 with a covering factor along the line-of-sight of C(sub los) approximately equal to 0.1.

K-band data↗

Efficient solutions to the Euler equations for supersonic flow with embedded subsonic regions

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two dimensions is described. Convergence of the basic algorithm to the steady state is quadratic for fully supersonic flows and is linear for other flows. This is in contrast to the block alternating direction implicit methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented herein is easily coupled with methods to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, and yields a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, Robert W.↗

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗

Double photoionization of atomic oxygen: Feshbach resonances in the two-electron continuum

In this work, we describe a joint experimental and theoretical investigation on oxygen double photoionization—the emission of two electrons from atomic oxygen following single photon absorption. High-resolution experimental measurements were performed at the Advanced Light Source, revealing sharp resonance structure superimposed on the more familiar Wannier-like, nearly-linear background. These resonance features are attributed to ionization-plus-excitation Feshbach resonances embedded in the double ionization continuum, doubly-excited states that lie above the double-ionization threshold. Such features are absent in the double photoionization cross section of He, or other quasi-two-electron systems, for which the doubly-ionized atomic core remains inert. For a corresponding theoretical analysis, the R-matrix with pseudostates (RMPS) method was invoked by calculating final-state, two-electron resonances-plus-continua wavefunctions and corresponding single-photon absorption cross sections. Overall agreement is found in the direct, background double photoionization cross section. However, the RMPS method, using a small basis due to practical computational limitations, was unable to reproduce quantitatively the smooth background or the sharper resonance features observed in the measurements, showing instead large-scale oscillations about the experimental background, and characteristic pseudoresonance jitter, associated with an insufficient convergence of the pseudostate representation to the true two-electron infinite series of Feshbach resonances embedded in the two-electron continuum. The prominent resonance structure observed highlights the need to consider multiple excitation processes in atoms more complex than He or quasi-two-electron systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Electric dipole radiation at VLF in a uniform warm magneto-plasma.

Use of a linear full electromagnetic wave theory to calculate the input impedance of an electric antenna embedded in a uniform, lossless, unbounded warm magnetoplasma, which is assumed to consist of warm electrons and cold ions. In calculating the dipole radiation resistance for the thermal modes and the thermally modified whistler mode the analysis includes the finite temperature only for the electrons. In deriving the formal solution of the warm plasma dipole input impedance a full-wave analysis is used and two antenna orientations are considered, parallel and perpendicular to the static magnetic field. A general dispersion equation governing the modes of propagation is derived and a detailed analysis is made of the propagation characteristics of these modes.

Wang, T. N. C.↗

Forcing function effects on unsteady aerodynamic gust response. II - Low solidity airfoil row response

The paper investigates the fundamental gust modeling assumption on the basis of a series of experiments performed in the Purdue Annular Cascade Research Facility. The unsteady period flow field is generated by rotating flows of perforated plates and airfoil cascades, with the resulting unsteady periodic chordwise pressure response of a downstream low solidity stator row determined by miniature pressure transducers embedded within selected airfoils. When the forcing function exhibited the characteristics of a linear-theory gust, the resulting response on the downstream stator airfoils was in excellent agreement with the linear-theory models. When the forcing function did not exhibit linear-theory gust characteristics, the resulting unsteady aerodynamic response of the downstream stators was much more complex and correlated poorly with the linear-theory gust predictions. It is shown that the forcing function generator significantly affects the resulting gust response, with the complexity of the response characteristics increasing from the perforated-plate to the airfoil-cascade forcing functions.

Henderson, Gregory H.↗

Depth dependence of absorbed dose, dose equivalent and linear energy transfer spectra of galactic and trapped particles in polyethylene and comparison with calculations of models

A matched set of five tissue-equivalent proportional counters (TEPCs), embedded at the centers of 0 (bare), 3, 5, 8 and 12-inch-diameter polyethylene spheres, were flown on the Shuttle flight STS-81 (inclination 51.65 degrees, altitude approximately 400 km). The data obtained were separated into contributions from trapped protons and galactic cosmic radiation (GCR). From the measured linear energy transfer (LET) spectra, the absorbed dose and dose-equivalent rates were calculated. The results were compared to calculations made with the radiation transport model HZETRN/NUCFRG2, using the GCR free-space spectra, orbit-averaged geomagnetic transmission function and Shuttle shielding distributions. The comparison shows that the model fits the dose rates to a root mean square (rms) error of 5%, and dose-equivalent rates to an rms error of 10%. Fairly good agreement between the LET spectra was found; however, differences are seen at both low and high LET. These differences can be understood as due to the combined effects of chord-length variation and detector response function. These results rule out a number of radiation transport/nuclear fragmentation models. Similar comparisons of trapped-proton dose rates were made between calculations made with the proton transport model BRYNTRN using the AP-8 MIN trapped-proton model and Shuttle shielding distributions. The predictions of absorbed dose and dose-equivalent rates are fairly good. However, the prediction of the LET spectra below approximately 30 keV/microm shows the need to improve the AP-8 model. These results have strong implications for shielding requirements for an interplanetary manned mission.

Flight Experiment↗