Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Investigation of non-linear contact for a clearance-fit bolt in a graphite/epoxy laminate

Numerous analytical studies have been published for the nonlinear load-contact variations in clearance-fit bolted joints. In these studies, stress distributions have been obtained and failure predictions have been made. However, very little experimental work has been reported regarding the contact or the stresses. This paper describes a fiber-optic technique for measuring the angle of contact in a clearance-fit bolt-loaded hole. Measurements of the contact angle have been made in a quasi-isotropic graphite-epoxy laminate by the optical as well as an electrical technique, and the results have been compared with those obtained from a finite-element analysis. The results from the two experimental techniques show excellent agreement; the finite-element results show some discrepancy, probably due to the interfacial frictions.

Prabhakaran, R.↗
DOE OSTI2016

Strengthening of reinforced concrete beams with basalt-based FRP sheets: An analytical assessment

In this paper the effectiveness of the flexural strengthening of RC beams through basalt fiber-reinforced sheets is investigated. The non-linear flexural response of RC beams strengthened with FRP composites applied at the traction side is described via an analytical formulation. Validation results and some comparative analyses confirm soundness and consistency of the proposed approach, and highlight the good mechanical performances (in terms of strength and ductility enhancement of the beam) produced by basalt-based reinforcements in comparison with traditional glass or carbon FRPs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Flux-driven algebraic damping of m = 1 diocotron mode

Recent experiments with pure electron plasmas in a Malmberg–Penning trap have observed the algebraic damping of m = 1 diocotron modes. Transport due to small field asymmetries produces a low density halo of electrons moving radially outward from the plasma core, and the mode damping begins when the halo reaches the resonant radius r = R{sub w} at the wall of the trap. The damping rate is proportional to the flux of halo particles through the resonant layer. The damping is related to, but distinct from, spatial Landau damping, in which a linear wave-particle resonance produces exponential damping. This paper explains with analytic theory the new algebraic damping due to particle transport by both mobility and diffusion. As electrons are swept around the “cat's eye” orbits of the resonant wave-particle interaction, they form a dipole (m = 1) density distribution. From this distribution, the electric field component perpendicular to the core displacement produces E × B-drift of the core back to the axis, that is, damps the m = 1 mode. The parallel component produces drift in the azimuthal direction, that is, causes a shift in the mode frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On non-autonomous dynamical systems

In usual realistic classical dynamical systems, the Hamiltonian depends explicitly on time. In this work, a class of classical systems with time dependent nonlinear Hamiltonians is analyzed. This type of problems allows to find invariants by a family of Veronese maps. The motivation to develop this method results from the observation that the Poisson-Lie algebra of monomials in the coordinates and momenta is clearly defined in terms of its brackets and leads naturally to an infinite linear set of differential equations, under certain circumstances. To perform explicit analytic and numerical calculations, two examples are presented to estimate the trajectories, the first given by a nonlinear problem and the second by a quadratic Hamiltonian with three time dependent parameters. In the nonlinear problem, the Veronese approach using jets is shown to be equivalent to a direct procedure using elliptic functions identities, and linear invariants are constructed. For the second example, linear and quadratic invariants as well as stability conditions are given. Explicit solutions are also obtained for stepwise constant forces. For the quadratic Hamiltonian, an appropriated set of coordinates relates the geometric setting to that of the three dimensional manifold of central conic sections. It is shown further that the quantum mechanical problem of scattering in a superlattice leads to mathematically equivalent equations for the wave function, if the classical time is replaced by the space coordinate along a superlattice. The mathematical method used to compute the trajectories for stepwise constant parameters can be applied to both problems. It is the standard method in quantum scattering calculations, as known for locally periodic systems including a space dependent effective mass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Nonlinear sensitivity coefficients and corrections in system identification

The evaluation of sensitivity coefficients which are the changes in the eigencharacteristics with respect to changes in the structural parameters is important in the identification of structural systems. In such studies the commonly used linearized sensitivity coefficients for system identification have occasionally resulted in poor convergence or even divergence during the iteration process to update analytical models for better correlation with the test data. Nonlinear sensitivity coefficients and nonlinear correction terms, usually eliminated during the linearization process, have been developed to evaluate sensitivity coefficients of linear systems. Application to computer simulated example problems indicates a more rapid and stable convergence of the iteration processes to update the analytical model for the improvement of the correlation with test data. The successful application of combining the nonlinear sensitivity coefficients and nonlinear correction terms with the Multiple Boundary Condition Tests enhances the ability to validate large space structures by ground tests.

Kuo, C. P.↗
DOE OSTI2018

Analytical and Computational Studies of Correlations of Hydrodynamic Fluctuations in Shear Flow

Correlations of hydrodynamic fluctuations in shear flow analytically and also by dissipative particle dynamics (DPD) simulations are studied. The hydrodynamic equations are linearized around the macroscopic velocity field and then solved by a perturbation method in Fourier-transformed space. The autocorrelation functions (ACFs) from the analytical method are compared with results obtained from DPD simulations under the same shear-flow conditions. Up to a moderate shear rate, various ACFs from the two approaches agree with each other well. At large shear rates, discrepancies between the two methods are observed, hence revealing strong additional coupling between different fluctuating variables, which is not considered in the analytical approach. In addition, the results at low and moderate shear rates can serve as benchmarks for developing multiscale algorithms for coupling of heterogeneous solvers, such as a hybrid simulation of molecular dynamics and fluctuating hydrodynamics solver, where thermal fluctuations are indispensable.

Physics↗
NASA NTRS1972

Relative motion of orbiting satellites

The relative motion problem is analyzed, as a linearized case, and as a numerically determined solution to provide a time history of the geometries representing the motion state. The displacement history and the hodographs for families of solutions are provided, analytically and graphically, to serve as an aid to understanding this problem area. Linearized solutions to relative motion problems of orbiting particles are presented for the impulsive and fixed thrust cases. Second order solutions are described to enhance the accuracy of prediction. A method was developed to obtain accurate, numerical solutions to the intercept and rendezvous problem; and, special situations are examined. A particular problem related to relative motions, where the motion traces develop a cusp, is examined in detail. This phenomenon is found to be dependent on a particular relationship between orbital eccentricity and the inclination between orbital planes. These conditions are determined, and, example situations are presented and discussed.

Eades, J. B., Jr.↗
DOE OSTI2020

Dynamical friction in a Fuzzy Dark Matter universe

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass m ~ O(10 -22 ) eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that (1) for FDM masses m & 10 -21 eV c -2 , the timing problem of the Fornax dwarf spheroidal’s globular clusters is no longer solved and (2) the effects of FDM on the process of dynamical friction for satellites of total mass M and relative velocity vrel should require detailed numerical simulations for M/109 M m/10 -22 eV 100 km s -1 /vrel ~ 1, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regi

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

The relations between the linear and bilinear atmospheric tidal fields in geometric height and in log-pressure coordinates

Theoretical models of planetary-atmosphere tidal fields are examined analytically, comparing models based on geometric-height coordinates with those employing log-pressure coordinates. The relationship between the linearized meteorological variables in the two systems is explored for classical tidal theory and its reduction to the Laplace tidal equation, and it is shown that identical horizontal and vertical equations are obtained. Also considered are the tidal zonal-mean bilinear flux convergences and their Eliassen-Palm formulations. Numerical results for problems involving the earth and Mars atmospheres are presented in graphs and discussed.

Zurek, Richard W.↗
DOE OSTI2014

Effects of initial radius of the interface and Atwood number on nonlinear saturation amplitudes in cylindrical Rayleigh-Taylor instability

Nonlinear saturation amplitudes (NSAs) of the first two harmonics in classical Rayleigh-Taylor instability (RTI) in cylindrical geometry for arbitrary Atwood numbers have been analytically investigated considering nonlinear corrections up to the fourth-order. The NSA of the fundamental mode is defined as the linear (purely exponential) growth amplitude of the fundamental mode at the saturation time when the growth of the fundamental mode (first harmonic) is reduced by 10% in comparison to its corresponding linear growth, and the NSA of the second harmonic can be obtained in the same way. The analytic results indicate that the effects of the initial radius of the interface (r{sub 0}) and the Atwood number (A) play an important role in the NSAs of the first two harmonics in cylindrical RTI. On the one hand, the NSA of the fundamental mode first increases slightly and then decreases quickly with increasing A. For given A, the smaller the r{sub 0}/λ (with λ perturbation wavelength) is, the larger the NSA of the fundamental mode is. When r{sub 0}/λ is large enough (r{sub 0}≫λ), the NSA of the fundamental mode is reduced to the prediction of previous literatures within the framework of third-order perturbation theory [J. W. Jacobs and I. Catton, J. Fluid Mech. 187, 329 (1988); S. W. Haan, Phys. Fluids B 3, 2349 (1991)]. On the other hand, the NSA of the second harmonic first decreases quickly with increasing A, reaching a minimum, and then increases slowly. Furthermore, the r{sub 0} can reduce the NSA of the second harmonic for arbitrary A at r{sub 0}≲2λ while increase it for A ≲ 0.6 at r{sub 0}≳2λ. Thus, it should be included in applications where the NSA has a role, such as inertial confinement fusion ignition target design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Transverse low frequency wave in a two fluid solar wind

Investigation is made of the properties of low frequency transverse waves in a two-fluid solar wind having a radial magnetic field and radial streaming velocity. In order to examine what effects this streaming medium has on the waves, linearly polarized waves are decomposed into left and right circularly polarized waves. Computation is made of analytic expressions valid to first order for the radial amplitude and phase dependence of these constituent waves. It is shown that after travelling a given distance r, these waves have different amplitudes and phases. The former result causes their superposition to become elliptical rather than linear. The latter causes the axis of the ellipse of polarization to rotate through a well-defined angle. Analytic expressions are obtained for the eccentricity of the ellipse and for the angle of rotation. In analogy with regular Faraday rotation, in which the plane of polarization of a linear polarized wave rotates, the effect is denoted as generalized Faraday rotation.

Solodyna, G. V.↗
DOE OSTI2018

On the geometric dependence of the molecular dipole polarizability in water: A benchmark study of higher-order electron correlation, basis set incompleteness error, core electron effects, and zero-point vibrational contributions

In this work, we investigate how geometric changes influence the static dipole polarizability ($α$) of a water molecule by explicitly computing the corresponding dipole polarizability surface (DPS) across 3125 total (1625 symmetry-unique) geometries using linear response coupled cluster theory including single, double, and triple excitations (LR-CCSDT) and the doubly augmented triple-$ζ$ basis set (d-aug-cc-pVTZ). Analytical formulae based on power series expansions of this $ab initio$ surface are generated using linear least-squares analysis and provide highly accurate estimates of this quantity as a function of molecular geometry (i.e., bond and angle variations) in a computationally tractable manner. An additional database, which consists of 25 representative molecular geometries and incorporates a more thorough treatment of both basis sets and core electron effects, is provided as a current benchmark for this quantity and the corresponding leading-order $C_6$ dispersion coefficient. This database has been utilized to assess the importance of these effects as well as the relative accuracy that can be obtained using several quantum chemical methods and a library of density functional approximations. In addition to high-level electron correlation methods (like CCSD) and our analytical least-squares formulae, we find that the SCAN0, PBE0, MN15, and B97-2 hybrid functionals yield the most accurate descriptions of the molecular polarizability tensor in H 2 O. Using first-order perturbation theory, we compute the zero-point vibrational correction to $α$ at the CCSDT/d-aug-cc-pVTZ level and find that this correction contributes approximately 3% to the isotropic ($α$ iso ) and nearly 50% to the anisotropic ($α$ aniso ) polarizability values. In doing so, we find that $α$ iso = 9.8307 bohr 3 , which is in excellent agreement with the experimental value of 9.83 ± 0.02 bohr 3 provided by Russell and Spackman. The DPS reported herein provides a benchmark-quality quantum mechanical estimate of this fundamental quantity of interest and should find extensive use in the development (and assessment) of next-generation force fields and machine-learning based approaches for modeling water in complex condensed-phase environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

MO-F-16A-02: Simulation of a Medical Linear Accelerator for Teaching Purposes

Purpose: Detailed functioning of linear accelerator physics is well known. Less well developed is the basic understanding of how the adjustment of the linear accelerator's electrical components affects the resulting radiation beam. Other than the text by Karzmark, there is very little literature devoted to the practical understanding of linear accelerator functionality targeted at the radiotherapy clinic level. The purpose of this work is to describe a simulation environment for medical linear accelerators with the purpose of teaching linear accelerator physics. Methods: Varian type lineacs were simulated. Klystron saturation and peak output were modelled analytically. The energy gain of an electron beam was modelled using load line expressions. The bending magnet was assumed to be a perfect solenoid whose pass through energy varied linearly with solenoid current. The dose rate calculated at depth in water was assumed to be a simple function of the target's beam current. The flattening filter was modelled as an attenuator with conical shape, and the time-averaged dose rate at a depth in water was determined by calculating kerma. Results: Fifteen analytical models were combined into a single model called SIMAC. Performance was verified systematically by adjusting typical linac control parameters. Increasing klystron pulse voltage increased dose rate to a peak, which then decreased as the beam energy was further increased due to the fixed pass through energy of the bending magnet. Increasing accelerator beam current leads to a higher dose per pulse. However, the energy of the electron beam decreases due to beam loading and so the dose rate eventually maximizes and the decreases as beam current was further increased. Conclusion: SIMAC can realistically simulate the functionality of a linear accelerator. It is expected to have value as a teaching tool for both medical physicists and linear accelerator service personnel.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

Stability of Brillouin flow in the presence of slow-wave structure

Including a slow-wave structure (SWS) on the anode in the conventional, planar, and inverted magnetron, we systematically study the linear stability of Brillouin flow, which is the prevalent flow in crossed-field devices. The analytic treatment is fully relativistic and fully electromagnetic, and it incorporates the equilibrium density profile, flow profile, and electric field and magnetic field profiles in the linear stability analysis. Using parameters similar to the University of Michigan's recirculating planar magnetron, the numerical data show that the resonant interaction of the vacuum circuit mode and the corresponding smooth-bore diocotron-like mode is the dominant cause for instability. This resonant interaction is far more important than the intrinsic negative (positive) mass property of electrons in the inverted (conventional) magnetron geometry. It is absent in either the smooth-bore magnetron or under the electrostatic assumption, one or both of which was almost always adopted in prior analytical formulation. This resonant interaction severely restricts the wavenumber for instability to the narrow range in which the cold tube frequency of the SWS is within a few percent of the corresponding smooth bore diocotron-like mode in the Brillouin flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Morphology-Driven Control of Metabolite Selectivity Using Nanostructure-Initiator Mass Spectrometry

Nanostructure-initiator mass spectrometry (NIMS) is a laser desorption/ionization analysis technique based on the vaporization of a nanostructure-trapped liquid "initiator" phase. Here in this work, we report an intriguing relationship between NIMS surface morphology and analyte selectivity. Scanning electron microscopy and spectroscopic ellipsometry were used to characterize the surface morphologies of a series of NIMS substrates generated by anodic electrochemical etching. Mass spectrometry imaging was applied to compare NIMS sensitivity of these various surfaces toward the analysis of diverse analytes. The porosity of NIMS surfaces was found to increase linearly with etching time where the pore size ranged from 4 to 12 nm with corresponding porosities estimated to be 7-70%. Surface morphology was found to significantly and selectively alter NIMS sensitivity. The small molecule (<2k Da) sensitivity was found to increase with increased porosity, whereas low porosity had the highest sensitivity for the largest molecules examined. Estimation of molecular sizes showed that this transition occurs when the pore size is <3× the maximum of molecular dimensions. While the origins of selectivity are unclear, increased signal from small molecules with increased surface area is consistent with a surface area restructuring-driven desorption/ionization process where signal intensity increases with porosity. In contrast, large molecules show highest signal for the low-porosity and small-pore-size surfaces. We attribute this to strong interactions between the initiator-coated pore structures and large molecules that hinder desorption/ionization by trapping large molecules. Lastly, this finding may enable us to design NIMS surfaces with increased specificity to molecules of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗