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 181 records · Page 10

Linear response of molecular polaritons

In this article, we show that the collective light–matter strong coupling regime where N molecular emitters couple to the photon mode of an optical cavity can be mapped to a quantum impurity model where the photon is the impurity that is coupled to a bath of anharmonic transitions. In the thermodynamic limit where N ≫ 1, we argue that the bath can be replaced with an effective harmonic bath, leading to a dramatic simplification of the problem into one of the coupled harmonic oscillators. We derive simple analytical expressions for linear optical spectra (transmission, reflection, and absorption) where the only molecular input required is the molecular linear susceptibility. This formalism is applied to a series of illustrative examples, showing the role of temperature, disorder, vibronic coupling, and optical saturation of the molecular ensemble, explaining that it is useful even when describing an important class of nonlinear optical experiments. For completeness, we provide Appendixes A–C that include a self-contained derivation of the relevant spectroscopic observables for arbitrary anharmonic systems (for both large and small N) within the rotating-wave approximation. While some of the presented results herein have already been reported in the literature, we provide a unified presentation of the results as well as new interpretations that connect powerful concepts in open quantum systems and linear response theory with molecular polaritonics.

Chemistry↗

Effects of linear trends on estimation of noise in GNSS position time-series

A thorough understanding of time-dependent noise in Global Navigation Satellite System (GNSS) position time-series is necessary for computing uncertainties in any signals found in the data. However, estimation of time-correlated noise is a challenging task and is complicated by the difficulty in separating noise from signal, the features of greatest interest in the time-series. In this study, we investigate how linear trends affect the estimation of noise in daily GNSS position time-series. We use synthetic time-series to study the relationship between linear trends and estimates of time-correlated noise for the six most commonly cited noise models. We find that the effects of added linear trends, or conversely de-trending, vary depending on the noise model. The commonly adopted model of random walk (RW), flicker noise (FN) and white noise (WN) is the most severely affected by de-trending, with estimates of low-amplitude RW most severely biased. FN plus WN is least affected by adding or removing trends. Non-integer power-law noise estimates are also less affected by de-trending, but are very sensitive to the addition of trend when the spectral index is less than one. We derive an analytical relationship between linear trends and the estimated RW variance for the special case of pure RW noise. Finally, overall, we find that to ascertain the correct noise model for GNSS position time-series and to estimate the correct noise parameters, it is important to have independent constraints on the actual trends in the data.

58 GEOSCIENCES↗

Turbofan forced mixer lobe flow modeling. 1: Experimental and analytical assessment

A joint analytical and experimental investigation of three-dimensional flowfield development within the lobe region of turbofan forced mixer nozzles is described. The objective was to develop a method for predicting the lobe exit flowfield. In the analytical approach, a linearized inviscid aerodynamical theory was used for representing the axial and secondary flows within the three-dimensional convoluted mixer lobes and three-dimensional boundary layer analysis was applied thereafter to account for viscous effects. The experimental phase of the program employed three planar mixer lobe models having different waveform shapes and lobe heights for which detailed measurements were made of the three-dimensional velocity field and total pressure field at the lobe exit plane. Velocity data was obtained using Laser Doppler Velocimetry (LDV) and total pressure probing and hot wire anemometry were employed to define exit plane total pressure and boundary layer development. Comparison of data and analysis was performed to assess analytical model prediction accuracy. As a result of this study a planar mixed geometry analysis was developed. A principal conclusion is that the global mixer lobe flowfield is inviscid and can be predicted from an inviscid analysis and Kutta condition.

Barber, T.↗

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗

Assessing the structure of classical molecular optimal control landscapes

Here we consider laser-driven optimal control landscape of a molecule from a classical mechanical perspective. The goal of optimal control in the present work is to steer the molecule from an initial state to a target state, denoted by two distinct points in phase space. Thus, a particular control objective is given as the difference between the final achieved phase space point and the target. The corresponding control landscape is defined as the latter control objective as a functional of the control field. While previous examination of the landscape critical points (i.e., a suboptimal point on the landscape where there is a zero gradient) has shown that the landscape topology is generally trap-free, the structure of the landscape away from these critical points is not well understood. We explore the landscape structure by examining an underlying metric defined as the ratio R of the gradient-based optimization path length of the control field evolution to the Euclidean distance between a given initial control field and the resultant optimal control field, where the latter field corresponds to a point at the top of the landscape. We analyze the path length-to-distance ratio R analytically for a linear forced harmonic oscillator and numerically for a nonlinear forced Morse oscillator. For the linear forced harmonic oscillator, we find that R ≤ $\sqrt{2}$ and reaches its minimum value of 1 (i.e., corresponding to “a straight shot” through control space) in the large target time limit, as well as at special finite target times. The ratio R is similarly small for Morse oscillator simulations when following a steepest-ascent path to the top of the landscape, implying that the landscape is quite smooth and devoid of gnarled features. This conclusion is exemplified for a path discovered with R≃ 1.0 where simply following the initial gradient direction takes the climb very close to the top of the landscape. These findings are consistent with a variety of previous like simulations examining R in quantum control scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

System Identification of X-33 Neural Network

Modern flight control research has improved spacecraft survivability as its goal. To this end we need to have a failure detection system on board. In case the spacecraft is performing imperfectly, reconfiguration of control is needed. For that purpose we need to have parameter identification of spacecraft dynamics. Parameter identification of a system is called system identification. We treat the system as a black box which receives some inputs that lead to some outputs. The question is: what kind of parameters for a particular black box can correlate the observed inputs and outputs? Can these parameters help us to predict the outputs for a new given set of inputs? This is the basic problem of system identification. The X33 was supposed to have the onboard capability of evaluating the current performance and if needed to take the corrective measures to adapt to desired performance. The X33 is comprised of both rocket and aircraft vehicle design characteristics and requires, in general, analytical methods for evaluating its flight performance. Its flight consists of four phases: ascent, transition, entry and TAEM (Terminal Area Energy Management). It spends about 200 seconds in ascent phase, reaching an altitude of about 180,000 feet and a speed of about 10 to 15 Mach. During the transition phase which lasts only about 30 seconds, its altitude may increase to about 190,000 feet but its speed is reduced to about 9 Mach. At the beginning of this phase, the Main Engine is Cut Off (MECO) and the control is reconfigured with the help of aerosurfaces (four elevons, two flaps and two rudders) and reaction control system (RCS). The entry phase brings down the altitude of X33 to about 90,000 feet and its speed to about Mach 3. It spends about 250 seconds in this phase. Main engine is still cut off and the vehicle is controlled by complex maneuvers of aerosurfaces. The last phase TAEM lasts for about 450 seconds and the altitude and speed, both are reduced to zero. The present attempt, as a start, focuses only on the entry phase. Since the main engine remains cut off in this phase, there is no thrust acting on the system. This considerably simplifies the equations of motion. We introduce another simplification by assuming the system to be linear after some non-linearities are removed analytically from our consideration. Under these assumptions, the problem could be solved by Classical Statistics by employing the least sum of squares approach. Instead we chose to use the Neural Network method. This method has many advantages. It is modern, more efficient, can be adapted to work even when the assumptions are diluted. In fact, Neural Networks try to model the human brain and are capable of pattern recognition.

Aggarwal, Shiv↗

Electron cloud considerations for HIF drivers

Here, we review some previous results on electron-cloud dynamics and accumulation and the impact of electron clouds on ion beam dynamics, and assess these in the context of a magnetic-quadrupole-based heavy-ion fusion driver. We present a new analysis which exploits analytic solution of linearized envelope equations between accelerating gaps to derive a set of mapping equations which we use to extrapolate from previously obtained particle-simulation results for a beam transport system to an HIF driver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diffraction-controlled backscattering threshold and application to Raman gap

In most classic analytical models of linear stimulated scatter, light diffraction is omitted, a priori. However, modern laser optic typically includes a variant of the random phase plate [Y. Kato et al., Phys. Rev. Lett. 53, 1057 (1984)], resulting in diffraction limited laser intensity fluctuations - or localized speckles - which may result in explosive reflectivity growth as the average laser intensity approaches a critical value [H. A. Rose and D. F. DuBois, Phys. Rev. Lett. 72, 2883 (1994)]. Among the differences between stimulated Raman scatter (SRS) and stimulated Brillouin scatter is that the SRS scattered light diffracts more strongly than the laser light with increase of electron density. This weakens the tendency of the SRS light to closely follow the most amplified paths, diminishing gain. Let G{sub 0} be the one-dimensional power gain exponent of the stimulated scatter. In this paper we show that differential diffraction gives rise to an increase of G{sub 0} at the SRS physical threshold with increase of electron density up to a drastic disruption of SRS as electron density approaches one fourth of its critical value from below. For three wave interaction lengths not small compared to a speckle length, this is a physically robust Raman gap mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tsunami and acoustic-gravity waves in water of constant depth

A study of wave radiation by a rather general bottom displacement, in a compressible ocean of otherwise constant depth, is carried out within the framework of a three-dimensional linear theory. Simple analytic expressions for the flow field, at large distance from the disturbance, are derived. Realistic numerical examples indicate that the Acoustic-Gravity waves, which significantly precede the Tsunami, are expected to leave a measurable signature on bottom-pressure records that should be considered for early detection of Tsunami.

58 GEOSCIENCES↗

Dissipative solitons in pair-ion plasmas

The effects of ion-neutral collisions on the dynamics of the nonlinear ion acoustic wave in pair-ion plasma are investigated. The standard perturbative approach leads to a Korteweg-de Vries equation with a linear damping term for the dynamics of the finite amplitude wave. The ion-neutral collision induced dissipation is responsible for the linear damping. The analytical solution and numerical simulation reveal that the nonlinear wave propagates in the form of a weakly dissipative compressive solitons. Furthermore, the width of the soliton is proportional to the amplitude of the wave for fixed soliton velocity. Results are discussed in the context of the fullerene pair-ion plasma experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗