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,099 records · Page 61

Stellar pulsations in beyond Horndeski gravity theories

Theories of gravity in the beyond Horndeski class recover the predictions of general relativity in the solar system whilst admitting novel cosmologies, including late-time de Sitter solutions in the absence of a cosmological constant. Deviations from Newton's law are predicted inside astrophysical bodies, which allow for falsifiable, smoking-gun tests of the theory. In this work we study the pulsations of stars by deriving and solving the wave equation governing linear adiabatic oscillations to find the modified period of pulsation. Using both semi-analytic and numerical models, we perform a preliminary survey of the stellar zoo in an attempt to identify the best candidate objects for testing the theory. Brown dwarfs and Cepheid stars are found to be particularly sensitive objects and we discuss the possibility of using both to test the theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulation of MHD instabilities with fluid runaway electron model in M3D- C 1

Runaway electrons are generated in a tokamak during the start up, during normal operation and during a plasma disruption. During a disruption, runaway electrons can be accelerated to high energies, potentially damaging the first wall. To predict the consequences of runaway generation during a disruption, it is necessary to consider resonant interactions of runaways with the bulk plasma. Here we consider the interactions of runaways on low mode number tearing modes. We have developed a fluid runaway electron model for the 3D MHD code M3D-C1[Jardin,et al. J Comput. Sci Discovery 6 014002 (2012)]. To benchmark, we have reproduced the MHD linear tearing mode results (with runaway electrons) in a circular cylinder presented in previous analytic studies[[Helander, P., et al, Phys. Plasmas 14 144102 (2007)] and extended here with a numerical eigenvalue calculation. Furthermore, we find that the low mode number tearing mode has a rotation caused by the MHD - runaways interaction, and the toroidal current scale length is much smaller with runaways than that for without and decreases as the runaway speed increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator

Here, a conservative discontinuous Galerkin scheme for a non-linear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it has a non-decreasing entropy functional, and conserves particles, momentum and energy. Discretely these conservative properties are maintained exactly as well, independent of numerical resolution. In this work the phase space discretization is performed using a novel version of the discontinuous Galerkin scheme, carefully constructed using concepts of weak equality and recovery. Discrete time advancement is carried out with an explicit time-stepping algorithm, whose stability limits we explore. The formulation and implementation within the long-wavelength gyrokinetic solver of Gkeyll are validated with relaxation tests, collisional Landau-damping benchmarks and the study of 5D gyrokinetic turbulence on helical, open field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stability of new exact solutions of the nonlinear Schrödinger equation in a Pöschl–Teller external potential

Here, we discuss the stability properties of the solutions of the general nonlinear Schrödinger equation in 1 + 1 dimensions in an external potential derivable from a parity-time ($ \newcommand{\PT}{\mathcal{PT}} \PT$ ) symmetric superpotential $W(x)$ that we considered earlier in Kevrekidis et al (2015 Phys. Rev. E 92 042901). In particular we consider the nonlinear partial differential equation $ \{i \, \partial_t + \partial_x^2 - V(x) + g \vert \psi(x, t) \vert ^{2\kappa} \} \, \psi(x, t) = 0 \>, $ for arbitrary nonlinearity parameter κ, where $g= \pm1$ and V is the well known Pöschl–Teller potential which we allow to be repulsive as well as attractive. Using energy landscape methods, linear stability analysis as well as a time dependent variational approximation, we derive consistent analytic results for the domains of instability of these new exact solutions as a function of the strength of the external potential and κ. For the repulsive potential we show that there is a translational instability which can be understood in terms of the energy landscape as a function of a stretching parameter and a translation parameter being a saddle near the exact solution. In this case, numerical simulations show that if we start with the exact solution, the initial wave function breaks into two pieces traveling in opposite directions. If we explore the slightly perturbed solution situations, a 1% change in initial conditions can change significantly the details of how the wave function breaks into two separate pieces. For the attractive potential, changing the initial conditions by 1% modifies the domain of stability only slightly. For the case of the attractive potential and negative g perturbed solutions merely oscillate with the oscillation frequencies predicted by the variational approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An idealized model for collapse of a molecular core within a filament

Here, we present an analytic treatment of the self-similar collapse of a uniform density ellipsoid to linear order in the deviations from sphericity. First we obtain a self-consistent closed-form solution for the evolution of an isolated ellipsoid and then impose the effects of an external forcing. This model describes the evolution under gravity of a pre-stellar core of molecular gas embedded in a much larger and lower-density filament. We compare with numerical solutions for the collapse without the limitation of small deviations. These show how the external perturbing force producing the initial deviation from sphericity is eventually surpassed by the anisotropic forces generated by the collapsing ellipsoid itself. This model should be useful in interpreting the way in which environment shapes the evolution of pre-stellar cores.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum critical collapse abhors a naked singularity

Classical critical collapse yields naked singularities from smooth initial data, challenging cosmic censorship, and shaping the spectrum of primordial black holes. We show that one-loop vacuum polarization near the threshold qualitatively changes this outcome by dressing the singularity with a horizon within a controlled semiclassical regime. In analytically tractable Einstein-scalar critical spacetimes, a one-loop $s$-wave treatment linearized around self-similar backgrounds shows that regularity uniquely selects an asymptotically Minkowskian, vacuum-polarization state. Its renormalized stress tensor carries a universal quantum growing mode that competes with the classical unstable mode, shifts the critical point, and generates a trapped surface along with a finite mass gap at the new threshold, thereby enforcing horizon formation even under arbitrary fine-tuning. In primordial collapse, the threshold shift enters exponentially into the formation fraction, while the mass gap truncates the low-mass tail, suggesting potentially important consequences for the predicted mass spectrum. Furthermore, these results provide a self-consistent semiclassical treatment of critical collapse and yield sharp predictions within the one-loop, near-critical, linearized regime.

Anomalies↗

SU-E-T-555: A Protontherapy Inverse Treatment Planning System Prototype with Linear Energy Transfer (LET) Optimization

Purpose: Develop and benchmark an inverse treatment planning system (TPS) for proton radiotherapy integrating fast analytical dose and LET calculations in patient geometries and a dual objective function with both dose and LET components, enabling us to apply optimization techniques to improve the predicted outcome of treatments based on radiobiological models. Methods: The software package was developed in MATLAB and implements a fluence-dose calculation technique based on a pencil beam model for dose calculations and a 3D LET model based on the extension of the LET in the radial direction as a function of the predicted radiological pathway. Both models were benchmarked against commissioning data from our institution, dose calculations performed with a commercial treatment planning system and Monte Carlo simulations. The optimization is based on the adaptive simulated annealing approach . Results: The dose and LET calculations were tested in a water phantom and several real patient treatments. The pass rate for the gamma index analysis (3%/3mm) test was above 90% for all test cases analyzed, and the calculation time was of the order of seconds. The inverse planning module produced plans with a significantly higher mean LET in the target compared to traditional plans, without any loss of target coverage. The clinical relevance of this improvement is under consideration . Conclusion: The developed treatment planning system is a valuable clinical and research tool that enables us to incorporate LET effects into proton radiotherapy planning in a streamlined fashion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A fluorescence-based method for rapid and direct determination of polybrominated diphenyl ethers in water

A new method was developed for rapid and direct measurement of polybrominated diphenyl ethers (PBDEs) in aqueous samples using fluorescence spectroscopy. The fluorescence spectra of tri- to deca-BDE (BDE 28, 47, 99, 153, 190, and 209) commonly found in environment were measured at variable emission and excitation wavelengths. The results revealed that the PBDEs have distinct fluorescence spectral profiles and peak positions that can be exploited to identify these species and determine their concentrations in aqueous solutions. The detection limits as determined in deionized water spiked with PBDEs are 1.71-5.82 ng/L for BDE 28, BDE 47, BDE 190, and BDE 209 and 45.55–69.95 ng/L for BDE 99 and BDE 153. The effects of environmental variables including pH, humic substance, and groundwater chemical composition on PBDEs measurements were also investigated. These environmental variables affected fluorescence intensity, but their effect can be corrected through linear additivity and separation of spectral signal contribution. Compared with conventional GC-based analytical methods, the fluorescence spectroscopy method is more efficient as it only uses a small amount of samples (2-4 mL), avoids lengthy complicated concentration and extraction steps, and has a low detection limit of a few ng/L.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear data covariances are critical input to determine upper sub-critical limits and to design experiments to increase it [Slides]

This presentation discusses how Upper Subcritical Limits (USL) are key parameters to determine operational limits in nuclear criticality safety evaluations. It also discusses an example of plutonium casting operation using tantalum at LANL PF-4. The Whisper tool at Los Alamos relies on many inputs, including covariance data, leading the presentation to ask if an existing benchmark data be used in Whisper to adjust nuclear data and covariances to justify a higher USL. If not, Whisper can be used to help design an optimal new benchmark experiment. The presentation also seeks to determine what the possible impacts are on USL and operational limits for plutonium casting. In conclusion, nuclear data covariances are used for by Whisper for: GSSL adjustment of nuclear data and covariances, identification of most similar existing benchmark experiments to application, simulation of Upper Subcritical Limit of application, and input to optimization techniques for designing most appropriate new benchmark experiment(s) to meet requirements. This requires a complete set of nuclear data covariances, benchmarks and k-effective sensitivity profiles (for both benchmarks and applications). The presentation concludes by asking if end users should trust results that depend on current covariance data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytical prediction of digital signal crosstalk of FCC

The results are presented of study effort whose aim was the development of accurate means of analyzing and predicting signal cross-talk in multi-wire digital data cables. A complete analytical model is developed n + 1 wire systems of uniform transmission lines with arbitrary linear boundary conditions. In addition, a minimum set of parameter measurements required for the application of the model are presented. Comparisons between cross-talk predicted by this model and actual measured cross-talk are shown for a six conductor ribbon cable.

Belleisle, A. P.↗

Analytical notes - Electrochemical method for early detection and monitoring of coliforms

An electrochemical method for detecting bacteria, based on a linear relationship between inoculum size and the time of hydrogen evolution, was tested for the early detection and monitoring of coliforms in naturally contaminated estuarine and fresh water samples. Standard methods for coliform analysis were performed on each sample, and membrane filtration counts were used to construct dose-response curves; relationships and results are discussed herein.

Wilkins, J. R.↗

Flat-sided rectilinear trough as a solar concentrator - An analytical study

Formulas are derived for the concentration factor and irradiance distribution at the base of the flat-sided linear trough. Performance is affected by the number of reflections the solar rays undergo before reaching the base, the cone apex angle, and the coefficient of reflection. Results are presented graphically in such a way that one can choose the optimum configuration, which is the minimum material required, to achieve a given concentration factor. Practical concentration factors range from 1.5 to 4 depending on the geometry and coefficient of reflection.

Burkhard, D. G.↗

Nonlinear spiral density waves - An inviscid theory

It is pointed out that the theory of spiral density waves, invented to explain the spiral structure of disk galaxies, has also been found useful for the study of planetary rings. The linear theory is by now well developed, while the nonlinear theory is less complete. Analytical calculations which include self-gravitation have, so far, obtained results only in the slightly nonlinear regime, or have concentrated on partial effects which are not of primary importance to the physical problem at hand. In the present paper, it is attempted to remedy these shortcomings. The simplest asymptotic ordering which can still yield useful results is adopted. Attention is given to the reduction to a nonlinear integral equation in a single variable, the use of the Wentzel-Kramers-Brillouin-Jeffreys theory, and the replacement of an equation by another which is easier to handle numerically.

Shu, F. H.↗

Toward comparing experiment and theory for corroborative research on hingeless rotor stability in forward flight

For flap lag stability of isolated rotors, experimental and analytical investigations were conducted in hover and forward flight on the adequacy of a linear quasisteady aerodynamics theory with dynamic flow. Forward flight effects on lag regressing mode were emphasized. A soft inplane hingeless rotor with three blades was tested at advance ratios as high as 0.55 and at shaft angles as high as 20 deg. The 1.62 m model rotor was untrimmed with an essentially unrestricted tilt of the tip path plane. In combination with lag natural frequencies, collective pitch settings and flap lag coupling parameters, the data base comprises nearly 1200 test points (damping and frequency) in forward flight and 200 test points in hover. By computerized symbolic manipulation, a linear model was developed in substall to predict stability margins with mode identification. To help explain the correlation between theory and data it also predicted substall and stall regions of the rotor disk from equilibrium values. The correlation showed both the strengths and weaknesses of the theory in substall ((angle of attack) equal to or less than 12 deg).

Gaonkar, G.↗

International Reference Ionosphere: Plasma densities - Status 1988

An account is given of the changes proposed in 1988 for the International Reference Ionosphere electron density profile, as well as the status of their implementation. The fully analytical profile function under development for the entire ionosphere can be achieved with a linear combination of several LAY functions. Although four LAY functions are required to describe the density features of the middle ionosphere, three LAY functions suffice to reproduce electron densities in both the topside ionosphere and lower ionosphere. Attention is given to the way in which the LAY parameters are computationally derivable from characteristic profile points.

Rawer, K.↗

Numerical Asymptotic Solutions Of Differential Equations

Numerical algorithms derived and compared with classical analytical methods. In method, expansions replaced with integrals evaluated numerically. Resulting numerical solutions retain linear independence, main advantage of asymptotic solutions.

Thurston, Gaylen A.↗

Integrated Aeroservoelastic Optimization: Status and Direction

The interactions of lightweight flexible airframe structures, steady and unsteady aerodynamics, and wide-bandwidth active controls on modern airplanes lead to considerable multidisciplinary design challenges. More than 25 years of mathematical and numerical methods' development, numerous basic research studies, simulations and wind-tunnel tests of simple models, wind-tunnel tests of complex models of real airplanes, as well as flight tests of actively controlled airplanes, have all contributed to the accumulation of a substantial body of knowledge in the area of aeroservoelasticity. A number of analysis codes, with the capabilities to model real airplane systems under the assumptions of linearity, have been developed. Many tests have been conducted, and results were correlated with analytical predictions. A selective sample of references covering aeroservoelastic testing programs from the 1960s to the early 1980s, as well as more recent wind-tunnel test programs of real or realistic configurations, are included in the References section of this paper. An examination of references 20-29 will reveal that in the course of development (or later modification), of almost every modern airplane with a high authority active control system, there arose a need to face aeroservoelastic problems and aeroservoelastic design challenges.

Livne, Eli↗