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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Probing the nature of electron transfer in metalloproteins on graphene-family materials as nanobiocatalytic scaffold using electrochemistry

Graphene-based nanomaterials have shown great promise not only in nanoelectronics due to ultrahigh electron mobility but also as biocatalytic scaffolds owing to irreversible protein surface adsorption and facilitating direct electron transfer. In this work, we synthesized stable dispersions of graphene using liquid-phase exfoliation approach based on non-covalent interactions between graphene and 1-pyrenesulfonic acid sodium salt (Py–1SO{sub 3}), 1-pyrenemethylamine salt (Py − Me-NH{sub 2}) and Pluronic{sup ®} P-123 surfactant using only water as solvent compatible with biomolecules. The resulting graphene nanoplatelets (Gr-LPE) are characterized by a combination of analytical (microscopy and spectroscopy) techniques revealing mono- to few-layer graphene displaying that the exfoliation efficiency strongly depends upon the type of pyrene-based salts and organic surfactants. Moreover being completely water-based approach, we build robust nanoscaffolds of graphene-family nanomaterials (GFNs) namely, monolayer graphene, Gr-LPE (the one prepared with Pluronic{sup ®} P-123), graphene oxide (GO) and its reduced form (rGO) on glassy carbon electrode surface with three important metalloproteins include cytochrome c (Cyt c) [for electron transfer], myoglobin (Mb) [for oxygen storage] and horseradish peroxidase (HRP) [for catalyzing the biochemical reaction]. In order to demonstrate the nanobiocatalytical activity of these proteins, we used electrochemical interfacial direct electron transfer (DET) kinetics and attempt to determine the rate constant (k{sub ET}) using two different analytical approaches namely, linear sweep voltammetry and Laviron’s theory. We elucidated that all of the metalloproteins retain their structural integrity (secondary structure) upon forming mixtures with GFNs confirmed through optical and vibrational spectroscopy and biological activity using electrochemistry. Among the GFNs studied, Gr-LPE, GO and rGO support the efficient electrical wiring of the redox centers (with an increase in catalytic efficiency of Cyt c and Mb in the presence of GFNs attributed partially to the surface functional (carboxyl, epoxide and hydroxyl) groups on GO and rGO facilitating rapid charge transfer.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Free energy change of a dislocation due to a Cottrell atmosphere

The free energy reduction of a dislocation due to a Cottrell atmosphere of solutes is computed using a continuum model. In this work, we show that the free energy change is composed of near-core and far-field components. The far-field component can be computed analytically using the linearized theory of solid solutions. Near the core the linearized theory is inaccurate, and the near-core component must be computed numerically. The influence of interactions between solutes in neighbouring lattice sites is also examined using the continuum model. We show that this model is able to reproduce atomistic calculations of the nickel–hydrogen system, predicting hydride formation on dislocations. The formation of these hydrides leads to dramatic reductions in the free energy. Lastly, the influence of the free energy change on a dislocation’s line tension is examined by computing the equilibrium shape of a dislocation shear loop and the activation stress for a Frank–Read source using discrete dislocation dynamics.

42 ENGINEERING↗

DOE OSTI2024

Quantum error mitigation by layerwise Richardson extrapolation

A widely used method for mitigating errors in noisy quantum computers is Richardson extrapolation, a technique in which the overall effect of noise on the estimation of quantum expectation values is captured by a single parameter that, after being scaled to larger values, is eventually extrapolated to the zero-noise limit. We generalize this approach by introducing layerwise Richardson extrapolation (LRE), an error mitigation protocol in which the noise of different individual layers (or larger chunks of the circuit) is amplified and the associated expectation values are linearly combined to estimate the zero-noise limit. The coefficients of the linear combination are analytically obtained from the theory of multivariate Lagrange interpolation. LRE leverages the flexible configurational space of layerwise unitary folding, allowing for a more nuanced mitigation of errors by treating the noise level of each layer of the quantum circuit as an independent variable. Furthermore, we provide numerical simulations demonstrating scenarios where LRE achieves superior performance compared to traditional (single-variable) Richardson extrapolation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2010

Stable structures with high topological charge in nonlinear photonic quasicrystals

Stable vortices with topological charges of 3 and 4 are examined numerically and analytically in photonic quasicrystals created by interference of five as well as eight beams, for cubic as well as saturable nonlinearities. Direct numerical simulations corroborate the analytical and numerical linear stability analysis predictions for such experimentally realizable structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Thermal control across a chain of electronic nanocavities

Here, we study a chain of alternating hot and cold electronic nanocavities—connected to one another via resonant-tunneling quantum dots—with the intent of achieving precise thermal control across the chain. This is accomplished by positioning the dots' energy levels such that a predetermined distribution of heat currents is realized across the chain in the steady state. The number of electrons in each cavity is conserved in the steady state which constrains the cavities' chemical potentials. We determine these chemical potentials analytically in the linear response regime where the energy differences between the dots' resonant levels and the neighboring chemical potentials are much smaller than the thermal energy. In this regime, the thermal control problem can be solved exactly, while, in the general case, thermal control can only be achieved in a relative sense, that is, when one only preassigns the ratios between different heat currents. We apply our results to two different cases: We first demonstrate that a “heat switch” can be easily realized with three coupled cavities, and then we show that our linear response results can provide accurate results in situations with a large number of cavities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2011

Appearance of bound states in random potentials with applications to soliton theory

We analyze the stochastic creation of a single bound state (BS) in a random potential with a compact support. We study both the Hermitian Schroedinger equation and non-Hermitian Zakharov-Shabat systems. These problems are of special interest in the inverse scattering method for Korteveg-de-Vries and the nonlinear Schroedinger equations since soliton solutions of these two equations correspond to the BSs of the two aforementioned linear eigenvalue problems. Analytical expressions for the average width of the potential required for the creation of the first BS are given in the approximation of delta-correlated Gaussian potential and additionally different scenarios of eigenvalue creation are discussed for the non-Hermitian case.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2010

Statistical analysis of dose heterogeneity in circulating blood: Implications for sequential methods of total body irradiation

Purpose: Improvements in delivery techniques for total body irradiation (TBI) using Tomotherapy and intensity modulated radiation therapy have been proven feasible. Despite the promise of improved dose conformality, the application of these ''sequential'' techniques has been hampered by concerns over dose heterogeneity to circulating blood. The present study was conducted to provide quantitative evidence regarding the potential clinical impact of this heterogeneity. Methods: Blood perfusion was modeled analytically as possessing linear, sinusoidal motion in the craniocaudal dimension. The average perfusion period for human circulation was estimated to be approximately 78 s. Sequential treatment delivery was modeled as a Gaussian-shaped dose cloud with a 10 cm length that traversed a 183 cm patient length at a uniform speed. Total dose to circulating blood voxels was calculated via numerical integration and normalized to 2 Gy per fraction. Dose statistics and equivalent uniform dose (EUD) were calculated for relevant treatment times, radiobiological parameters, blood perfusion rates, and fractionation schemes. The model was then refined to account for random dispersion superimposed onto the underlying periodic blood flow. Finally, a fully stochastic model was developed using binomial and trinomial probability distributions. These models allowed for the analysis of nonlinear sequential treatment modalities and treatment designs that incorporate deliberate organ sparing. Results: The dose received by individual blood voxels exhibited asymmetric behavior that depended on the coherence among the blood velocity, circulation phase, and the spatiotemporal characteristics of the irradiation beam. Heterogeneity increased with the perfusion period and decreased with the treatment time. Notwithstanding, heterogeneity was less than {+-}10% for perfusion periods less than 150 s. The EUD was compromised for radiosensitive cells, long perfusion periods, and short treatment times. However, the EUD was unaffected (within 10%) for perfusion periods of less than 150 s or treatment times of 20 min or greater. Treatment over six fractions improved the EUD per fraction such that all parametric combinations resulted in unaffected EUD. The stochastic models confirmed these results. Conclusions: Dose heterogeneity in circulating blood cells is clinically acceptable for typical treatment times, perfusion rates, and cell types. Development of conformal, sequential TBI treatment techniques should not be withheld based on concerns over circulating blood dose heterogeneity.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Spatiotemporal optical pulse transformation by a resonant diffraction grating

The diffraction of a spatiotemporal optical pulse by a resonant diffraction grating is considered. The pulse diffraction is described in terms of the signal (the spatiotemporal incident pulse envelope) passage through a linear system. An analytic approximation in the form of a rational function of two variables corresponding to the angular and spatial frequencies has been obtained for the transfer function of the system. A hyperbolic partial differential equation describing the general form of the incident pulse envelope transformation upon diffraction by a resonant diffraction grating has been derived from the transfer function. A solution of this equation has been obtained for the case of normal incidence of a pulse with a central frequency lying near the guided-mode resonance of a diffraction structure. The presented results of numerical simulations of pulse diffraction by a resonant grating show profound changes in the pulse envelope shape that closely correspond to the proposed theoretical description. The results of the paper can be applied in creating new devices for optical pulse shape transformation, in optical information processing problems, and analog optical computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗

DOE OSTI2017

A Hybrid Multi-Scale Model of Crystal Plasticity for Handling Stress Concentrations

Microstructural effects become important at regions of stress concentrators such as notches, cracks and contact surfaces. A multiscale model is presented that efficiently captures microstructural details at such critical regions. The approach is based on a multiresolution mesh that includes an explicit microstructure representation at critical regions where stresses are localized. At regions farther away from the stress concentration, a reduced order model that statistically captures the effect of the microstructure is employed. The statistical model is based on a finite element representation of the orientation distribution function (ODF). As an illustrative example, we have applied the multiscaling method to compute the stress intensity factor K I around the crack tip in a wedge-opening load specimen. The approach is verified with an analytical solution within linear elasticity approximation and is then extended to allow modeling of microstructural effects on crack tip plasticity.

36 MATERIALS SCIENCE↗

DOE OSTI2016

CALCULATING ROTATING HYDRODYNAMIC AND MAGNETOHYDRODYNAMIC WAVES TO UNDERSTAND MAGNETIC EFFECTS ON DYNAMICAL TIDES

To understand magnetic effects on dynamical tides, we study the rotating magnetohydrodynamic (MHD) flow driven by harmonic forcing. The linear responses are analytically derived in a periodic box under the local WKB approximation. Both the kinetic and Ohmic dissipations at the resonant frequencies are calculated, and the various parameters are investigated. Although magnetic pressure may be negligible compared to thermal pressure, the magnetic field can be important for the first-order perturbation, e.g., dynamical tides. It is found that the magnetic field splits the resonant frequency, namely the rotating hydrodynamic flow has only one resonant frequency, but the rotating MHD flow has two, one positive and the other negative. In the weak field regime the dissipations are asymmetric around the two resonant frequencies and this asymmetry is more striking with a weaker magnetic field. It is also found that both the kinetic and Ohmic dissipations at the resonant frequencies are inversely proportional to the Ekman number and the square of the wavenumber. The dissipation at the resonant frequency on small scales is almost equal to the dissipation at the non-resonant frequencies, namely the resonance takes its effect on the dissipation at intermediate length scales. Moreover, the waves with phase propagation that is perpendicular to the magnetic field are much more damped. It is also interesting to find that the frequency-averaged dissipation is constant. This result suggests that in compact objects, magnetic effects on tidal dissipation should be considered.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Self-consistent nonlinear kinetic simulations of the anomalous Doppler instability of suprathermal electrons in plasmas

Suprathermal tails in the distributions of electron velocities parallel to the magnetic field are found in many areas of plasma physics, from magnetic confinement fusion to solar system plasmas. Parallel electron kinetic energy can be transferred into plasma waves and perpendicular gyration energy of particles through the anomalous Doppler instability (ADI), provided that energetic electrons with parallel velocities v{sub ||}≥(ω+Ω{sub ce})/k{sub ||} are present; here Ω{sub ce} denotes electron cyclotron frequency, ω the wave angular frequency, and k{sub ||} the component of wavenumber parallel to the magnetic field. This phenomenon is widely observed in tokamak plasmas. Here, we present the first fully self-consistent relativistic particle-in-cell simulations of the ADI, spanning the linear and nonlinear regimes of the ADI. We test the robustness of the analytical theory in the linear regime and follow the ADI through to the steady state. By directly evaluating the parallel and perpendicular dynamical contributions to j·E in the simulations, we follow the energy transfer between the excited waves and the bulk and tail electron populations for the first time. We find that the ratio Ω{sub ce}/(ω{sub pe}+Ω{sub ce}) of energy transfer between parallel and perpendicular, obtained from linear analysis, does not apply when damping is fully included, when we find it to be ω{sub pe}/(ω{sub pe}+Ω{sub ce}); here ω{sub pe} denotes the electron plasma frequency. We also find that the ADI can arise beyond the previously expected range of plasma parameters, in particular when Ω{sub ce}>ω{sub pe}. The simulations also exhibit a spectral feature which may correspond to the observations of suprathermal narrowband emission at ω{sub pe} detected from low density tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

An improved understanding of the natural resonances of moonpools contained within floating rigid-bodies: Theory and application to oscillating water column devices

The fundamental interactions between waves, a floating rigid-body, and a moonpool that is selectively open to atmosphere or enclosed to purposefully induce pressure fluctuations are investigated. The moonpool hydrodynamic characteristics and the hydrodynamic coupling to the rigid-body are derived implicitly through reciprocity relations on an array of field points. By modeling the free surface of the moonpool in this manner, an explicit hydrodynamic coupling term is included in the equations of motion. This coupling results in the migration of the moonpool's natural resonance frequency from the piston frequency to a new frequency when enclosed in a floating rigid-body. Two geometries that highlight distinct aspects of marine vessels and oscillating water column (OWC) renewable energy devices are analyzed to reveal the coupled natural resonance migration. The power performance of these two OWCs in regular waves is also investigated. The air chamber is enclosed and a three-dimensional, linear, frequency domain performance model that links the rigid-body to the moonpool through a linear resistive control strategy is detailed. Furthermore, an analytic expression for the optimal linear resistive control values in regular waves is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

On the road to per cent accuracy – III. Non-linear reaction of the matter power spectrum to massive neutrinos

ABSTRACT We analytically model the non-linear effects induced by massive neutrinos on the total matter power spectrum using the halo model reaction framework of Cataneo et al. In this approach, the halo model is used to determine the relative change to the matter power spectrum caused by new physics beyond the concordance cosmology. Using standard fitting functions for the halo abundance and the halo mass–concentration relation, the total matter power spectrum in the presence of massive neutrinos is predicted to per cent-level accuracy, out to $k=10 \,{ h}\,{\rm Mpc}^{-1}$. We find that refining the prescriptions for the halo properties using N-body simulations improves the recovered accuracy to better than 1 per cent. This paper serves as another demonstration for how the halo model reaction framework, in combination with a single suite of standard Λ cold dark matter (ΛCDM) simulations, can recover per cent-level accurate predictions for beyond ΛCDM matter power spectra, well into the non-linear regime.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

An efficient closed-form solution for acoustic emission source location in three-dimensional structures

This paper presents an efficient closed-form solution (ECS) for acoustic emission(AE) source location in three-dimensional structures using time difference of arrival (TDOA) measurements from N receivers, N ≥ 6. The nonlinear location equations of TDOA are simplified to linear equations. The unique analytical solution of AE sources for unknown velocity system is obtained by solving the linear equations. The proposed ECS method successfully solved the problems of location errors resulting from measured deviations of velocity as well as the existence and multiplicity of solutions induced by calculations of square roots in existed close-form methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Bowtie filters for dedicated breast CT: Theory and computational implementation

Purpose: To design bowtie filters with improved properties for dedicated breast CT to improve image quality and reduce dose to the patient. Methods: The authors present three different bowtie filters designed for a cylindrical 14-cm diameter phantom with a uniform composition of 40/60 breast tissue, which vary in their design objectives and performance improvements. Bowtie design #1 is based on single material spectral matching and produces nearly uniform spectral shape for radiation incident upon the detector. Bowtie design #2 uses the idea of basis material decomposition to produce the same spectral shape and intensity at the detector, using two different materials. Bowtie design #3 eliminates the beam hardening effect in the reconstructed image by adjusting the bowtie filter thickness so that the effective attenuation coefficient for every ray is the same. All three designs are obtained using analytical computational methods and linear attenuation coefficients. Thus, the designs do not take into account the effects of scatter. The authors considered this to be a reasonable approach to the filter design problem since the use of Monte Carlo methods would have been computationally intensive. The filter profiles for a cone-angle of 0° were used for the entire length of each filter because the differences between those profiles and the correct cone-beam profiles for the cone angles in our system are very small, and the constant profiles allowed construction of the filters with the facilities available to us. For evaluation of the filters, we used Monte Carlo simulation techniques and the full cone-beam geometry. Images were generated with and without each bowtie filter to analyze the effect on dose distribution, noise uniformity, and contrast-to-noise ratio (CNR) homogeneity. Line profiles through the reconstructed images generated from the simulated projection images were also used as validation for the filter designs. Results: Examples of the three designs are presented. Initial verification of performance of the designs was done using analytical computations of HVL, intensity, and effective attenuation coefficient behind the phantom as a function of fan-angle with a cone-angle of 0°. The performance of the designs depends only weakly on incident spectrum and tissue composition. For all designs, the dynamic range requirement on the detector was reduced compared to the no-bowtie-filter case. Further verification of the filter designs was achieved through analysis of reconstructed images from simulations. Simulation data also showed that the use of our bowtie filters can reduce peripheral dose to the breast by 61% and provide uniform noise and CNR distributions. The bowtie filter design concepts validated in this work were then used to create a computational realization of a 3D anthropomorphic bowtie filter capable of achieving a constant effective attenuation coefficient behind the entire field-of-view of an anthropomorphic breast phantom. Conclusions: Three different bowtie filter designs that vary in performance improvements were described and evaluated using computational and simulation techniques. Results indicate that the designs are robust against variations in breast diameter, breast composition, and tube voltage, and that the use of these filters can reduce patient dose and improve image quality compared to the no-bowtie-filter case.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Sci—Fri PM: Topics — 05: Experience with linac simulation software in a teaching environment

Medical linear accelerator education is usually restricted to use of academic textbooks and supervised access to accelerators. To facilitate the learning process, simulation software was developed to reproduce the effect of medical linear accelerator beam adjustments on resulting clinical photon beams. The purpose of this report is to briefly describe the method of operation of the software as well as the initial experience with it in a teaching environment. To first and higher orders, all components of medical linear accelerators can be described by analytical solutions. When appropriate calibrations are applied, these analytical solutions can accurately simulate the performance of all linear accelerator sub-components. Grouped together, an overall medical linear accelerator model can be constructed. Fifteen expressions in total were coded using MATLAB v 7.14. The program was called SIMAC. The SIMAC program was used in an accelerator technology course offered at our institution; 14 delegates attended the course. The professional breakdown of the participants was: 5 physics residents, 3 accelerator technologists, 4 regulators and 1 physics associate. The course consisted of didactic lectures supported by labs using SIMAC. At the conclusion of the course, eight of thirteen delegates were able to successfully perform advanced beam adjustments after two days of theory and use of the linac simulator program. We suggest that this demonstrates good proficiency in understanding of the accelerator physics, which we hope will translate to a better ability to understand real world beam adjustments on a functioning medical linear accelerator.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2023

Analytical Large-Signal Modeling of Inverter-Based Microgrids With Koopman Operator Theory for Autonomous Control

The microgrid (MG) plays a crucial role in the energy transition, but its nonlinearity presents a significant challenge for large-signal power systems studies in the electromagnetic transient (EMT) time scale. In this paper, we develop a large-signal linear MG model that considers the detailed dynamics of the primary and zero-control levels based on the Koopman operator (KO) theory. Firstly, a set of observable functions is carefully designed to capture the nonlinear dynamics of the MG. The corresponding linear KO is then analytically derived based on these observables, resulting in the linear representation of the original nonlinear MG with observables as the new coordinate. The influence of external input on the system dynamics is also considered during the derivation, enabling control of the MG. We solve the voltage control problem using the traditional linear quadratic integrator (LQI) method to demonstrate that textbook linear control techniques can accurately control the original nonlinear MG via the developed KO-linearized MG model. Our proposed KO linearization method is generic and can be easily extended for different control objectives and MG structures using our analytical derivation procedure. Here, we validate the effectiveness of our methodology through various case studies.

25 ENERGY STORAGE↗