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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

An analytical method as a preconditioning modeling for systems of linear equations

In this paper, we propose a simple and practical preconditioner based on the homotopy analysis method. Furthermore, using this modeling, we design an efficient algorithm to solve linear systems. Within the computational framework, convection–diffusion equation candidates are investigated and evaluated. Finally, computational comparisons of the standard methods against preconditioned methods based on experiments are presented which illustrate the merits of the simplicity, power, and effectiveness of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Combined mechanical loading of composite tubes

An analytical/experimental investigation was performed to study the effect of material nonlinearities on the response of composite tubes subjected to combined axial and torsional loading. The effect of residual stresses on subsequent mechanical response was included in the investigation. Experiments were performed on P75/934 graphite-epoxy tubes with a stacking sequence of (15/0/ + or - 10/0/ -15), using pure torsion and combined axial/torsional loading. In the presence of residual stresses, the analytical model predicted a reduction in the initial shear modulus. Experimentally, coupling between axial loading and shear strain was observed in laminated tubes under combined loading. The phenomenon was predicted by the nonlinear analytical model. The experimentally observed linear limit of the global shear response was found to correspond to the analytically predicted first ply failure. Further, the failure of the tubes was found to be path dependent above a critical load level.

Derstine, Mark S.↗
DOE OSTI2012

Electron cyclotron resonance near the axis of the gas-dynamic trap

Propagation of an extraordinary electromagnetic wave in the vicinity of electron cyclotron resonance surface in an open linear trap is studied analytically, taking into account inhomogeneity of the magnetic field in paraxial approximation. Ray trajectories are derived from a reduced dispersion equation that makes it possible to avoid the difficulty associated with a transition from large propagation angles to the case of strictly longitudinal propagation. Our approach is based on the theory, originally developed by Zvonkov and Timofeev [Sov. J. Plasma Phys. 14, 743 (1988)], who used the paraxial approximation for the magnetic field strength, but did not consider the slope of the magnetic field lines, which led to considerable error, as has been recently noted by Gospodchikov and Smolyakova [Plasma Phys. Rep. 37, 768-774 (2011)]. We have found ray trajectories in analytic form and demonstrated that the inhomogeneity of both the magnetic field strength and the field direction can qualitatively change the picture of wave propagation and significantly affect the efficiency of electron cyclotron heating of a plasma in a linear magnetic trap. Analysis of the ray trajectories has revealed a criterion for the resonance point on the axis of the trap to be an attractor for the ray trajectories. It is also shown that a family of ray trajectories can still reach the resonance point on the axis if the latter generally repels the ray trajectories. As an example, results of general theory are applied to the electron cyclotron resonance heating experiment which is under preparation on the gas dynamic trap in the Budker Institute of Nuclear Physics [Shalashov et al., Phys. Plasmas 19, 052503 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Magnetic islands and singular currents at rational surfaces in three-dimensional magnetohydrodynamic equilibria

Using the recently developed multiregion, relaxed MHD (MRxMHD) theory, which bridges the gap between Taylor's relaxation theory and ideal MHD, we provide a thorough analytical and numerical proof of the formation of singular currents at rational surfaces in non-axisymmetric ideal MHD equilibria. These include the force-free singular current density represented by a Dirac δ-function, which presumably prevents the formation of islands, and the Pfirsch-Schlüter 1/x singular current, which arises as a result of finite pressure gradient. An analytical model based on linearized MRxMHD is derived that can accurately (1) describe the formation of magnetic islands at resonant rational surfaces, (2) retrieve the ideal MHD limit where magnetic islands are shielded, and (3) compute the subsequent formation of singular currents. The analytical results are benchmarked against numerical simulations carried out with a fully nonlinear implementation of MRxMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Closed bipolar electrode-enabled dual-cell electrochromic detectors for chemical sensing

Bipolar electrodes (BPE) are electrically floating metallic elements placed in electrified fluids that enable the coupling of anodic and cathodic redox reactions at the opposite ends by electron transfer through the electrode. One particularly compelling application allows electron transfer reactions at one end of a closed BPE to be read out optically by inducing a redox-initiated change in the optical response function of a reporter system at the other end. Here, a BPE-enabled method for electrochemical sensing based on the electrochromic response of a methyl viologen (MV) reporter is developed, characterized, and rendered in a field-deployable format. BPE-enabled devices based on two thin-layer-cells of ITO and Pt were fabricated to couple an analytical reaction in one cell with an MV reporter reaction, producing a color change in the complementary cell. Using Fe(CN) 6 3/4- as a model analyte, the electrochemically induced color change of MV was determined initially by measuring its absorbance via a CCD camera coupled to a microscope. Then, smartphone-based detection and RGB analysis were employed to further simplify the sensing scheme. In both methods a linear relationship between the analyte concentration was produced, the quantity of MV generated, and the colorimetric response, yielding a limit of detection of 1.0 μM. Similar responses were also observed in the detection of dopamine and acetaminophen. Further evolution of the device replaced the potentiostat with batteries to control potential, demonstrating the simplicity and portability of the device. Finally, the physical separation of the reporter and analytical cells renders the device competent to detect analytes in different (e.g. non-aqueous) phases, as demonstrated by using the electrochromic behavior of aqueous MV to detect ferrocene in acetonitrile in the analytical cell.

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

RELAXATION OF WARPED DISKS: THE CASE OF PURE HYDRODYNAMICS

Orbiting disks may exhibit bends due to a misalignment between the angular momentum of the inner and outer regions of the disk. We begin a systematic simulational inquiry into the physics of warped disks with the simplest case: the relaxation of an unforced warp under pure fluid dynamics, i.e., with no internal stresses other than Reynolds stress. We focus on the nonlinear regime in which the bend rate is large compared to the disk aspect ratio. When warps are nonlinear, strong radial pressure gradients drive transonic radial motions along the disk's top and bottom surfaces that efficiently mix angular momentum. The resulting nonlinear decay rate of the warp increases with the warp rate and the warp width, but, at least in the parameter regime studied here, is independent of the sound speed. The characteristic magnitude of the associated angular momentum fluxes likewise increases with both the local warp rate and the radial range over which the warp extends; it also increases with increasing sound speed, but more slowly than linearly. The angular momentum fluxes respond to the warp rate after a delay that scales with the square root of the time for sound waves to cross the radial extent of the warp. These behaviors are at variance with a number of the assumptions commonly used in analytic models to describe linear warp dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The spin polarized linear response from density functional theory: Theory and application to atoms

Within the context of spin polarized conceptual density functional theory, the spin polarized linear response functions are introduced both in the [N, N{sub s}] and [N{sub α}, N{sub β}] representations. The mathematical relations between the spin polarized linear response functions in both representations are examined and an analytical expression for the spin polarized linear response functions in the [N{sub α}, N{sub β}] representation is derived. The spin polarized linear response functions were calculated for all atoms up to and including argon. To simplify the plotting of our results, we integrated χ(r, r′) to a quantity χ(r, r{sup ′}), circumventing the θ and ϕ dependence. This allows us to plot and to investigate the periodicity throughout the first three rows in the periodic table within the two different representations. For the first time, χ{sub αβ}(r, r{sup ′}), χ{sub βα}(r, r{sup ′}), and χ{sub SS}(r, r{sup ′}) plots have been calculated and discussed. By integration of the spin polarized linear response functions, different components to the polarisability, α{sub αα}, α{sub αβ}, α{sub βα}, and α{sub ββ} have been calculated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Direct Linearization and Adjoint Approaches to Evaluation of Atmospheric Weighting Functions and Surface Partial Derivatives: General Principles, Synergy and Areas of Application

This slide presentation reviews the observable radiances as functions of atmospheric parameters and of surface parameters; the mathematics of atmospheric weighting functions (WFs) and surface partial derivatives (PDs) are presented; and the equation of the forward radiative transfer (RT) problem is presented. For non-scattering atmospheres this can be done analytically, and all WFs and PDs can be computed analytically using the direct linearization approach. For scattering atmospheres, in general case, the solution of the forward RT problem can be obtained only numerically, but we need only two numerical solutions: one of the forward RT problem and one of the adjoint RT problem to compute all WFs and PDs we can think of. In this presentation we discuss applications of both the linearization and adjoint approaches

remote sensing↗
DOE OSTI2016

Unsteady fluid flow in a slightly curved pipe: A comparative study of a matched asymptotic expansions solution with a single analytical solution

The present work is motivated by the fact that blood flow in the aorta and the main arteries is governed by large finite values of the Womersley number α and for such values of α there is not any analytical solution in the literature. The existing numerical solutions, although accurate, give limited information about the factors that affect the flow, whereas an analytical approach has an advantage in that it can provide physical insight to the flow mechanism. Having this in mind, we seek analytical solution to the equations of the fluid flow driven by a sinusoidal pressure gradient in a slightly curved pipe of circular cross section when the Womersley number varies from small finite to infinite values. Initially the equations of motion are expanded in terms of the curvature ratio δ and the resulting linearized equations are solved analytically in two ways. In the first, we match the solution for the main core to that for the Stokes boundary layer. This solution is valid for very large values of α. In the second, we derive a straightforward single solution valid to the entire flow region and for 8 ≤ α < ∞, a range which includes the values of α that refer to the physiological flows. Each solution contains expressions for the axial velocity, the stream function, and the wall stresses and is compared to the analogous forms presented in other studies. The two solutions give identical results to each other regarding the axial flow but differ in the secondary flow and the circumferential wall stress, due to the approximations employed in the matched asymptotic expansion process. The results on the stream function from the second solution are in agreement with analogous results from other numerical solutions. The second solution predicts that the atherosclerotic plaques may develop in any location around the cross section of the aortic wall unlike to the prescribed locations predicted by the first solution. In addition, it gives circumferential wall stresses augmented by approximately 100% with respect to the matched asymptotic expansions, a factor that may contribute jointly with other pathological factors to the faster aging of the arterial system and the possible malfunction of the aorta.

42 ENGINEERING↗
DOE OSTI2018

Exploiting graphical processing units to enable quantum chemistry calculation of large solvated molecules with conductor-like polarizable continuum models

The conductor–like polarizable continuum model (C–PCM) with switching/Gaussian smooth discretization is a widely used implicit solvation model in quantum chemistry. We have previously implemented C–PCM solvation for Hartree–Fock (HF) and density functional theory (DFT) on graphical processing units (GPUs), enabling the quantum mechanical treatment of large solvated biomolecules. Here, we first propose a GPU–based algorithm for the PCM conjugate gradient linear solver that greatly improves the performance for very large molecules. The overhead for PCM–related evaluations now consumes less than 15% of the total runtime for DFT calculations on large molecules. Second, we demonstrate that our algorithms tailored for ground state C–PCM are transferable to excited state properties. Using a single GPU, our method evaluates the analytic gradient of the linear response PCM time–dependent density functional theory energy up to 80× faster than a conventional central processing unit (CPU)–based implementation. In addition, our C–PCM algorithms are transferable to other methods that require electrostatic potential (ESP) evaluations. For example, we achieve speed–ups of up to 130× for restricted ESP–based atomic charge evaluations, when compared to CPU–based codes. As a result, we also summarize and compare the different PCM cavity discretization schemes used in some popular quantum chemistry packages as a reference for both users and developers.

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

Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators

Abstract The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using numerical optimization methods. It is found that while the near-axis expansion is able to replicate the growth rates, real frequencies and perpendicular wave-vector at the inner core (both using simplified dispersion relations and first-principle gyrokinetic simulations), it leads to an overestimation of the growth rate at larger radii. An approximate analytic solution of the ITG dispersion relation for the non-resonant limit suggests growth rates could be systematically higher in quasi-axisymmetric (QA) configurations compared to quasi-helically (QH) symmetric ones. However except for very close to the axis, linear gyrokinetic simulations do not show systematic differences between QA and QH configurations.

Physics↗
NASA NTRS1985

Durability predictions of adhesively bonded composite structures using accelerated characterization methods

The utilization of adhesive bonding for composite structures is briefly assessed. The need for a method to determine damage initiation and propagation for such joints is outlined. Methods currently in use to analyze both adhesive joints and fiber reinforced plastics is mentioned and it is indicated that all methods require the input of the mechanical properties of the polymeric adhesive and composite matrix material. The mechanical properties of polymers are indicated to be viscoelastic and sensitive to environmental effects. A method to analytically characterize environmentally dependent linear and nonlinear viscoelastic properties is given. It is indicated that the methodology can be used to extrapolate short term data to long term design lifetimes. That is, the method can be used for long term durability predictions. Experimental results for near adhesive resins, polymers used as composite matrices and unidirectional composite laminates is given. The data is fitted well with the analytical durability methodology. Finally, suggestions are outlined for the development of an analytical methodology for the durability predictions of adhesively bonded composite structures.

Brinson, H. F.↗
NASA NTRS1982

Coupled flap-lag-torsional dynamics of hingeless rotor blades in forward flight

The aeroelastic stability and response problem of the coupled flap-lag-torsional dynamics of an isolated hingeless rotor blade in forward flight is considered. Linear, quasi-steady aerodynamics below stall is included. The spatial dependence of the partial differential nonliner equations of motion is discretized using a multimodal Galerkin method. The nonlinear time dependent equilibrium position (steady state response) about which the equations are linearized is obtained by solving a sequence of linear periodic response problems (quasilinearization). Results illustrating blade behavior in forward flight for both soft-in-plane and stiff-in-plane designs are presented. Quasilinearization provides a clear indication of the cases when nonlinear terms due to moderate deflections are important. The results indicate that nonlinearities affect system stability much more than system response. This implies that for obtaining blade vibrations and loads, approximate analytical models based on linearized formulations could be employed, whereas for determining blade stability a more accurate nonlinear analysis would be required.

Friedmann, P. P.↗
NASA NTRS1983

Linearized transfer between inclined circular orbits using low-thrust blow down propulsion system

Noncoplanar transfers between neighboring circular orbits are presented for spacecraft using their own low-thrust blow down propulsion system. It is assumed that the out-of-plane angle between the decaying thrust vector and the current orbit plane remains constant for each extended burn. Switching conditions are derived for the cutoff and relight of the propulsion system in order to carry out a given transfer with inclination change. Furthermore the location where the thrust acceleration is initially applied with respect to the line of nodes of the two orbits is uniquely determined. Finally an analytic derivation of the linearized coplanar motion for stationkeeping and terminal rendezvous studies is also presented and a scheme for deriving the second order correction shown.

Kechichian, J. A.↗
NASA NTRS1987

A simple parameterization for quality factor as a function of linear energy transfer

The paper presents a simple analytic approximation of the radiation quality factor (Q) as a function of linear energy transfer for use in radiation protection calculations. The paper also presents estimated quality factors in water for protons over a broad range of incident energies. It is shown that the quality factors are less than unity for all proton energies greater than 13 MeV.

Townsend, Lawrence W.↗
DOE OSTI2015

Stable static structures in models with higher-order derivatives

We investigate the presence of static solutions in generalized models described by a real scalar field in four-dimensional space–time. We study models in which the scalar field engenders higher-order derivatives and spontaneous symmetry breaking, inducing the presence of domain walls. Despite the presence of higher-order derivatives, the models keep to equations of motion second-order differential equations, so we focus on the presence of first-order equations that help us to obtain analytical solutions and investigate linear stability on general grounds. We then illustrate the general results with some specific examples, showing that the domain wall may become compact and that the zero mode may split. Moreover, if the model is further generalized to include k-field behavior, it may contribute to split the static structure itself.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A boundary condition to the Khokhlov-Zabolotskaya equation for modeling strongly focused nonlinear ultrasound fields

An equivalent source model was proposed as a boundary condition to the nonlinear parabolic Khokhlov-Zabolotskaya (KZ) equation to simulate high intensity focused ultrasound (HIFU) fields generated by medical ultrasound transducers with the shape of a spherical shell. The boundary condition was set in the initial plane; the aperture, the focal distance, and the initial pressure of the source were chosen based on the best match of the axial pressure amplitude and phase distributions in the Rayleigh integral analytic solution for a spherical transducer and the linear parabolic approximation solution for the equivalent source. Analytic expressions for the equivalent source parameters were derived. It was shown that the proposed approach allowed us to transfer the boundary condition from the spherical surface to the plane and to achieve a very good match between the linear field solutions of the parabolic and full diffraction models even for highly focused sources with F-number less than unity. The proposed method can be further used to expand the capabilities of the KZ nonlinear parabolic equation for efficient modeling of HIFU fields generated by strongly focused sources.

42 ENGINEERING↗