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 1,405 records · Page 78

Cosmological perturbation theory in 1+1 dimensions

Many recent studies have highlighted certain failures of the standard Eulerian-space cosmological perturbation theory (SPT). Its problems include (1) not capturing large-scale bulk flows [leading to an O( 1) error in the 1-loop SPT prediction for the baryon acoustic peak in the correlation function], (2) assuming that the Universe behaves as a pressureless, inviscid fluid, and (3) treating fluctuations on scales that are non-perturbative as if they were. Recent studies have highlighted the successes of perturbation theory in Lagrangian space or theories that solve equations for the effective dynamics of smoothed fields. Both approaches mitigate some or all of the aforementioned issues with SPT. We discuss these physical developments by specializing to the simplified 1D case of gravitationally interacting sheets, which allows us to substantially reduces the analytic overhead and still (as we show) maintain many of the same behaviors as in 3D. In 1D, linear-order Lagrangian perturbation theory ('the Zeldovich approximation') is exact up to shell crossing, and we prove that n{sup th}-order Eulerian perturbation theory converges to the Zeldovich approximation as narrow ∞. In no 1D cosmology that we consider (including a CDM-like case and power-law models) do these theories describe accurately the matter power spectrum on any mildly nonlinear scale. We find that theories based on effective equations are much more successful at describing the dynamics. Finally, we discuss many topics that have recently appeared in the perturbation theory literature such as beat coupling, the shift and smearing of the baryon acoustic oscillation feature, and the advantages of Fourier versus configuration space. Our simplified 1D case serves as an intuitive review of these perturbation theory results.

79 ASTRONOMY AND ASTROPHYSICS↗

Two-peaked and flat-top perfect bright solitons in nonlinear metamaterials with epsilon near zero

We investigate analytically transverse-magnetic spatial bright solitons, as exact solutions of Maxwell's equations, propagating through nonlinear metamaterials whose linear dielectric permittivity is very close to zero and whose effective nonlinear Kerr parameters can be tailored to achieve values not available in standard materials. Exploiting the fact that, in the medium considered, linear and nonlinear polarization can be comparable at feasible and realistic optical intensities, we identify two self-trapping mechanisms able to support two-peaked and flat-top solitons, respectively. Specifically, these two mechanisms are based on the occurrence of critical points at which the effective nonlinear permittivity vanishes, the two mechanisms differing in the way the compensation between linear and nonlinear polarization is achieved through the nonstandard values of the nonlinear parameters.

36 MATERIALS SCIENCE↗

Crosstalk Compensation for a Rapid, Higher Resolution Impedance Spectrum Measurement

Batteries and other energy storage devices are playing larger roles in various industries (e.g., military, automotive, electric utilities, etc.) as the U.S. seeks to reduce its dependence on foreign energy resources. As such, there exists a significant need for accurate, robust state-of-health assessment techniques. Present techniques tend to focus on simple, passive monitoring of voltage and current at a given ambient temperature. However, this approach has the disadvantage of ignoring key elements of health, that is, changes in resistance growth and power fade. Impedance spectroscopy is considered a useful laboratory tool in gauging changes in the resistance and power performance, but it has not been widely considered as an onboard diagnostic tool due to the length of time required to complete the measurement. Cross-Talk Compensation (CTC) is a novel approach that enables rapid, high resolution impedance spectra measurements using a hardware platform that could be designed as an embedded system. This input signal consists of a sum-of-sines excitation current that has a known frequency spread and a duration of one period of the lowest frequency. The voltage response is then captured at a sufficiently fast sample rate. Previously developed rapid impedance spectrum measurement techniques either required a longer excitation signal or a sum-of-sines signal that was separated by harmonic frequencies to reduce or eliminate, respectively, the cross-talk interference in the calculated results. The distinct advantage of CTC, however, is that non-harmonic frequencies can now be included within the excitation signal while still keeping the signal duration at one period of the lowest frequency. Since the frequency spread of the input signal is known, the crosstalk interference between sinusoidal signals within the sum-of-sines at a given frequency of interest can be pre-determined and assigned to an error matrix. Consequently, the real and imaginary components of the impedance at each frequency of interest can be calculated using simple linear algebra based on the error matrix and measured response from the energy storage device given the excitation signal. Analytical validation of CTC over a frequency range between 2000 and 0.1 Hz (i.e., a ten-second input signal duration) was performed using a standardized battery lumped parameter model. The results indicated that the CTC was able to successfully resolve more than 45 frequencies within a sum-of-sines excitation signal, whereas previous techniques could only resolve up to 15 frequencies. A simplified derivation of the CTC technique and its corresponding analytical validation studies using the lumped-parameter model will be presented.

25 ENERGY STORAGE↗

SU-E-J-85: Leave-One-Out Perturbation (LOOP) Fitting Algorithm for Absolute Dose Film Calibration

Purpose: To introduce an outliers-recognition fitting routine for film dosimetry. It cannot only be flexible with any linear and non-linear regression but also can provide information for the minimal number of sampling points, critical sampling distributions and evaluating analytical functions for absolute film-dose calibration. Methods: The technique, leave-one-out (LOO) cross validation, is often used for statistical analyses on model performance. We used LOO analyses with perturbed bootstrap fitting called leave-one-out perturbation (LOOP) for film-dose calibration . Given a threshold, the LOO process detects unfit points (“outliers”) compared to other cohorts, and a bootstrap fitting process follows to seek any possibilities of using perturbations for further improvement. After that outliers were reconfirmed by a traditional t-test statistics and eliminated, then another LOOP feedback resulted in the final. An over-sampled film-dose- calibration dataset was collected as a reference (dose range: 0-800cGy), and various simulated conditions for outliers and sampling distributions were derived from the reference. Comparisons over the various conditions were made, and the performance of fitting functions, polynomial and rational functions, were evaluated. Results: (1) LOOP can prove its sensitive outlier-recognition by its statistical correlation to an exceptional better goodness-of-fit as outliers being left-out. (2) With sufficient statistical information, the LOOP can correct outliers under some low-sampling conditions that other “robust fits”, e.g. Least Absolute Residuals, cannot. (3) Complete cross-validated analyses of LOOP indicate that the function of rational type demonstrates a much superior performance compared to the polynomial. Even with 5 data points including one outlier, using LOOP with rational function can restore more than a 95% value back to its reference values, while the polynomial fitting completely failed under the same conditions. Conclusion: LOOP can cooperate with any fitting routine functioning as a “robust fit”. In addition, it can be set as a benchmark for film-dose calibration fitting performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Octavius1500 2D ion chamber array and its associated phantoms: Dosimetric characterization of a new prototype

Purpose: The purpose of the study is to characterize the prototype of the new Octavius1500 (PTW, Freiburg, Germany) 2D ion chamber array, covering its use in different phantom setups, from the most basic solid water sandwich setup to the more complex cylindrical Octavius{sup ®} 4D (Oct4D) (PTW) phantom/detector combination. The new detector houses nearly twice the amount of ion chambers as its predecessors (Seven29 and Octavius729), thereby tackling one of the most important limitations of ion chamber (or diode) arrays, namely the limited detector density. The 0.06 cm{sup 3} cubic ion chambers are now arranged in a checkerboard pattern, leaving no lines (neither longitudinally nor laterally) without detectors. Methods: All measurements were performed on a dual energy (6 MV and 18 MV) iX Clinac (Varian Medical Systems, Palo Alto, CA) and all calculations were done in the Eclipse treatment planning system (Varian) with the Anisotropic Analytical Algorithm. First, the basic characteristics of the 2D array, such as measurement stability, dose rate dependence and dose linearity were investigated in the solid water sandwich setup. Second, the directional dependence was assessed to allow the evaluation of the new Octavius2D phantom (Oct2D{sup 1500}) for planar verification measurements of composite plans. Third, measurements were performed in the Oct4D phantom to evaluate the impact of the increased detector density on the accuracy of the volumetric dose reconstruction. Results: While showing equally good dose linearity and dose rate independence, the Octavius1500 outperforms the previous models because of its instantaneous measurement stability and its twofold active area coverage. Orthogonal field-by-field measurements immediately benefit from the increased detector density. The 3.9 cm wide compensation cavity in the new Oct2D{sup 1500} phantom prototype adequately corrects for directional dependence from the rear, resulting in good agreement within the target dose. Discrepancies may arise towards the sides of the array because of uncompensated lateral beam incidence. The beneficial impact of the detector density is most prominent in the Oct4D system, for which the average pass rate (PR) is now nearly 100% (99.31 ± 0.37) when using gamma criteria of 2%G,2 mm (10% dose threshold). In search of gamma analysis criteria that are not too lenient to detect possibly relevant deviations, the authors conclude that for our radiotherapy environment, the authors choose to adopt 3%L,3 mm PR97% (threshold 10%) criteria for the Oct2D{sup 1500}/Octavius1500 system and 2%L,3 mm PR97% (threshold 10%) for the Oct4D/Octavius1500 system. These are first line pass/check criteria and plans that fail are not necessarily rejected, but submitted to a more detailed investigation. Conclusions: When irradiated from the front, the Octavius1500 array has two main advantages over its 729 predecessors: its instantaneous measurement stability and-–most importantly—its twofold detector density. In the Oct2D1500 phantom, these advantages are counterbalanced by the more pronounced directional dependence. The measurement-based 3D dose reconstruction in the Oct4D system, however, benefits considerably from the higher detector density in the checkerboard panel design.

07 ISOTOPE AND RADIATION SOURCES↗

Nonlinear quantum theory of stimulated Cherenkov radiation of transverse electromagnetic waves from a low-density relativistic electron beam in a dielectric medium

A nonlinear quantum theory of stimulated Cherenkov radiation of transverse electromagnetic waves from a low-density relativistic electron beam in an isotropic dielectric medium is presented. A quantum model based on the Klein-Gordon equation is used. The growth rates of beam instabilities caused by the effect of stimulated Cherenkov radiation have been determined in the linear approximation. Mechanisms of the nonlinear saturation of relativistic quantum Cherenkov beam instabilities have been analyzed and the corresponding analytical solutions have been obtained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Microstructural evolution of Q12{sup TM} alloy irradiated in PWRs and comparison with other Zr base alloys - 2016-0061

Based on the M5{sup R}* alloy metallurgy, the Q12{sup TM} alloy (Zr-1Nb-0.5Sn-0.1Fe) was developed by AREVA NP for structural components, with ultra-low tin addition and slightly increased iron content. The behavior of this alloy was tested under irradiation in a pressurized water reactor (PWR) and has shown improvement in irradiation creep strength and similar free growth compared with M5. This paper provides results on dimensional stability and details the microstructural evolution of the Q12 alloy under neutron irradiation in PWRs. The Q12 microstructural evolution under irradiation was studied for fast neutron fluences up to 13 X 10{sup 25} n/m{sup 2} (E>1 MeV) with analytical transmission electron microscopy observations. We focused on radiation-enhanced needle-like particles, Laves phases, and the linear density of -component loops. These results are compared with other quaternary zirconium-niobium-tin-iron alloys and with M5. All these results allow a general discussion about microstructural evolution and behavior under irradiation of quaternary-type alloys compared to M5 alloy. This study, in agreement with previous works on Zr-1Nb and quaternary alloys, seems to show that increasing the iron content with the presence of niobium and tin will decrease the -component loop linear density and delay the growth breakaway. (authors)

36 MATERIALS SCIENCE↗

Analytical studies of ultrashort optical pulse propagation in a resonant absorbing medium

Resonant absorption in a laser pulse can be reduced to zero if the trailing part induces stimulated emission and returns the absorbed energy back to the pulse. This kind of effect is called Self-Induced Transparency (SIT). Analytic descriptions of short pulse propagation in a resonant medium, based on Slowly Varying Envelope Approximation (SVEA) and in linear media, have been extended in studies of SIT along two principal directions: First, to improve the SVEA by introducing Short Pulse Approximation (SPA); and second, to consider the effect of nonlinear refractive index.

Lee, C. T.↗

A higher order panel method for linearized supersonic flow

The basic integral equations of linearized supersonic theory for an advanced supersonic panel method are derived. Methods using only linear varying source strength over each panel or only quadratic doublet strength over each panel gave good agreement with analytic solutions over cones and zero thickness cambered wings. For three dimensional bodies and wings of general shape, combined source and doublet panels with interior boundary conditions to eliminate the internal perturbations lead to a stable method providing good agreement experiment. A panel system with all edges contiguous resulted from dividing the basic four point non-planar panel into eight triangular subpanels, and the doublet strength was made continuous at all edges by a quadratic distribution over each subpanel. Superinclined panels were developed and tested on s simple nacelle and on an airplane model having engine inlets, with excellent results.

Ehlers, F. E.↗

LANDSAT-D MSS/TM tuned orbital jitter analysis model LDS900

The final LANDSAT-D orbital dynamic math model (LSD900), comprised of all test validated substructures, was used to evaluate the jitter response of the MSS/TM experiments. A dynamic forced response analysis was performed at both the MSS and TM locations on all structural modes considered (thru 200 Hz). The analysis determined the roll angular response of the MSS/TM experiments to improve excitation generated by component operation. Cross axis and cross experiment responses were also calculated. The excitations were analytically represented by seven and nine term Fourier series approximations, for the MSS and TM experiment respectively, which enabled linear harmonic solution techniques to be applied to response calculations. Single worst case jitter was estimated by variations of the eigenvalue spectrum of model LSD 900. The probability of any worst case mode occurrence was investigated.

Pollak, T. E.↗

Differential rotation in a solar-driven quasi-axisymmetric circulation

The concept of a quasi-axisymmetric circulation is used to explore the global scale dynamics of planetary atmospheres. A numerical circulation model applicable to Jupiter is presented, and an analytical study is performed elucidating the conditions leading to differential rotation in an atmosphere which is convectively unstable. A linear system forced by solar differential heating is considered, with nonlinear effects arising from advection being represented in the form of eddy diffusion. An empirical, latitudinal spectrum of the observed zonal wind field on Jupiter is discussed. Numerical solutions are presented which reveal banded wind fields with alternating and equatorial zonal jets and a multicellular Ferrel-Thomson meridional circulation consistent with the observed cloud striations on Jupiter. The vertical derivatives are parameterized to construct a simplified one-layer model.

Mayr, H. G.↗

Re-examination of cumulative fatigue damage analysis: An engineering perspective

A method which has evolved in our laboratories for the past 20 yr is re-examined with the intent of improving its accuracy and simplicity of application to engineering problems. Several modifications are introduced both to the analytical formulation of the Damage Curve Approach, and to the procedure for modifying this approach to achieve a Double Linear Damage Rule formulation which immensely simplifies the calculation. Improvements are also introduced in the treatment of mean stress for determining fatigue life of the individual events that enter into a complex loading history. While the procedure is completely consistent with the results of numerous two level tests that have been conducted on many materials, it is still necessary to verify applicability to complex loading histories. Caution is expressed that certain phenomena can also influence the applicability - for example, unusual deformation and fracture modes inherent in complex loading - especially if stresses are multiaxial. Residual stresses at crack tips, and metallurgical factors are also important in creating departures from the cumulative damage theories; examples of departures are provided.

Manson, S. S.↗

Simplex finite element analysis of viscous incompressible flow with penalty function formulation

Viscous flow calculations are important for the determination of separated flows, recirculating flows, secondary flows and so on. This paper presents a penalty function approach for the finite element analysis of steady incompressible viscous flow. A simplex element is used with linear velocity and constant pressure in contrast to other works which usually employ higher order elements. Simplex elements yield analytical expressions for the element matrices which in turn lead to efficient solutions. Earlier works have partially indicated how constraint and lock-up problems might be avoided for simplex elements. This paper extends the earlier works by indicating the approach in detail and verifying that it is successful for several applications not discussed in the literature so far. Solution times and accuracy considerations are discussed for Couette flow, plane Poiseuille flow, a driven cavity problem, and laminar and turbulent flow over a step.

Allaire, P. E.↗

A perturbation approach to cosmic ray transients in interplanetary space

A perturbation approach is used to model the linear response of the cosmic ray distribution as a function of perturbations in the solar wind transport parameters. The analytical technique permits examination of the effects of the different solar wind parameters, i.e., velocity, drift velocity and the diffusion tensor, which are additive. Cosmic ray changes in the solar wind are neglected. Variations in the diffusion coefficient and the convection-diffusion equation are applied to describing Forbush decreases and the 11-yr solar cycle variations. The treatment is limited to energies above 100 MeV, yet is considered valid enough to have identified a hysteresis effect in the 11-yr variation, wherein high-energy particles lead low-energy particles to a magnitude that increases with increasing heliocentric distance within the inner heliosphere. The calculations also indicate that the Forbush decreases are a cumulative effect of precipitous precursors accompanied by slow recoveries from travelling solar wind perturbations.

Chih, P. P.↗

Elevated temperature crack growth

The objective of the Elevated Temperature Crack Growth Program is to evaluate proposed nonlinear fracture mechanics methods for application to hot section components of aircraft gas turbine engines. Progress during the past year included linear-elastic fracture mechanics data reduction on nonlinear crack growth rate data on Alloy 718. The bulk of the analytical work centered on thermal gradient problems and proposed fracture mechanics parameters. Good correlation of thermal gradient experimental displacement data and finite element prediction was obtained.

Malik, S. N.↗

Re-examination of cumulative fatigue damage analysis - An engineering perspective

A method which has evolved in the laboratories for the past 20 yr is re-examined with the intent of improving its accuracy and simplicity of application to engineering problems. Several modifications are introduced both to the analytical formulation of the Damage Curve Approach, and to the procedure for modifying this approach to achieve a Double Linear Damage Rule formulation which immensely simplifies the calculation. Improvements are also introduced in the treatment of mean stress for determining fatigue life of the individual events that enter into a complex loading history. While the procedure is completely consistent with the results of numerous two level tests that have been conducted on many materials, it is still necessary to verify applicability to complex loading histories. Caution is expressed that certain phenomenon can also influence the applicability - for example, unusual deformation and fracture modes inherent in complex loading especially if stresses are multiaxial. Residual stresses at crack tips, and metallurgical factors are also important in creating departures from the cumulative damage theories; examples of departures are provided.

Manson, S. S.↗

Stability analysis of double-peaked proton distribution functions in the solar wind

Twelve proton velocity distribution functions as measured for different solar wind speeds by the plasma analyzer at the Helios spacecraft are modeled analytically by a superposition of two bi-Maxwellians. The stability of these distributions is investigated by solving the linear dispersion equation for parallel propagation of both left- and right-hand polarized modes. Due to the temperature anisotropy of the main proton component the presence of the resonant left-hand ion cyclotron instability is found to be typical for high-speed streams and is independent of an ion beam. The results indicate that the ratio of the beam density to total particle density and the beam drift velocity are generally close to the threshold values necessary to drive a highly unstable nonresonant right-hand mode. This instability is found for one proton distribution and is a consequence of the beam particles. It is concluded that both the ion cyclotron instability and the ion beam instability act as a regulating mechanism on the shape of solar wind double humped proton distribution functions.

Leubner, Manfred P.↗

Effects of unequal particle number densities on Alfven waves

Analytic plasma theory and numerical solutions of the dispersion equation are used to show that the assumption that the linear properties of the waves are determined by a charge-neutral plasma in the absence of the nonthermal particles, while the nonthermal particles cause growth or additional damping superposed onto the background, is seriously flawed even for stable plasmas. Even when the nonthermal particles do not contribute significantly to the dispersion equation, unequal thermal electron and ion number densities (due to the presence of the nonthermal particles) may cause fundamental low wave number modifications to the Alfven modes, including the creation of a new resonance and severely modified dispersion. These results are found for both cold and warm plasmas. Previous work on Alfven waves should be reevaluated in view of these results.

Cairns, I. H.↗