Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Molecular Architecture Directs Linear–Bottlebrush–Linear Triblock Copolymers to Self-Assemble to Soft Reprocessable Elastomers

Linear–bottlebrush–linear (LBBL) triblock copolymers represent an emerging system for creating multifunctional nanostructures. Their self-assembly depends on molecular architecture but remains poorly explored. We synthesize polystyrene- block -bottlebrush polydimethylsiloxane- block -polystyrene triblock copolymers with controlled molecular architecture and use them as a model system to study the self-assembly of LBBL polymers. Unlike classical stiff rod-flexible linear block copolymers that are prone to form highly ordered nanostructures such as lamellae, at small weight fractions of the linear blocks, LBBL polymers self-assemble to a disordered sphere phase, regardless of the bottlebrush stiffness. Microscopically, characteristic lengths increase with the bottlebrush stiffness by a power of 2/3, which is captured by a scaling analysis. Macroscopically, the formed nanostructures are ultrasoft, reprocessable elastomers with shear moduli of about 1 kPa, two orders of magnitude lower than that of conventional polydimethylsiloxane elastomers. Our results provide insights on exploiting the self-assembly of LBBL polymers to create soft functional nanostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Exciton Absorption Spectra by Linear Response Methods: Application to Conjugated Polymers

The theoretical description of the time-evolution of excitons requires, as an initial step, the calculation of their spectra, which has been inaccessible to most users due to the high computational scaling of conventional algorithms and accuracy issues caused by common density functionals. Previously (J. Chem. Phys.2016, 144, 204105), we developed a simple method that resolves these issues. Our scheme is based on a two-step calculation in which a linear-response TDDFT calculation is used to generate orbitals perturbed by the excitonic state, and then a second linear-response TDDFT calculation is used to determine the spectrum of excitations relative to the excitonic state. Herein, we apply this theory to study near-infrared absorption spectra of excitons in oligomers of the ubiquitous conjugated polymers poly(3-hexylthiophene) (P3HT), poly(2-methoxy-5-(2-ethylhexyloxy)-1,4-phenylenevinylene) (MEH-PPV), and poly(benzodithiophene-thieno[3,4-b]thiophene) (PTB7). For P3HT and MEH-PPV oligomers, the calculated intense absorption bands converge at the longest wavelengths for 10 monomer units, and show strong consistency with experimental measurements. The calculations confirm that the exciton spectral features in MEH-PPV overlap with those of the bipolaron formation. In addition, our calculations identify the exciton absorption bands in transient absorption spectra measured by our group for oligomers (1, 2, and 3 units) of PTB7. For all of the cases studied, we report the dominant orbital excitations contributing to the optically active excited state–excited state transitions, and suggest a simple rule to identify absorption peaks at the longest wavelengths. We suggest our methodology could be considered for further developments in theoretical transient spectroscopy to include nonadiabatic effects, coherences, and to describe the formation of species such as charge-transfer states and polaron pairs.

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

Linear Response Function Reveals the Most Effective Remote Forcing in Causing September Arctic Sea Ice Melting in CESM

We apply the linear response function method to investigate the most excitable mode of the September Arctic sea ice and its corresponding remote oceanic forcing in climate models. The method is useful in identifying the fundamental, dynamical and causal relationship between the remote forcing and sea ice response. We find that the most excitable mode of the September Arctic sea ice preferentially takes place over the Pacic side of the Arctic and its remote forcing corresponds to a dipole pattern of precipitation anomaly in the tropics with an increase of precipitation over the western and central tropical and subtropical Pacic ocean while a decrease over the Maritime Continent. The tropical precipitation anomaly likely drives a Rossby wave train propagating toward higher latitudes and leads to a ridge anomaly over the Pacic side of the Arctic, resulting in poleward atmospheric heat transport, enhanced downward longwave radiation and thus melting of the sea ice. In addition, a good agreement is found with the leading tropical Arctic teleconnection mode detected in a pre-industrial control simulation, which supports the usefulness and robustness of the linear response function method and the importance of tropical forcing in shaping the sea ice variability.

Wu, Yutian↗
DOE OSTI2014

Enhanced localized energetic ion losses resulting from first-orbit linear and non-linear interactions with Alfvén eigenmodes in DIII-D

Two key insights into interactions between Alfvén eigenmodes (AEs) and energetic particles in the plasma core are gained from measurements and modeling of first-orbit beam-ion loss in DIII-D. First, the neutral beam-ion first-orbit losses are enhanced by AEs and a single AE can cause large fast-ion displacement. The coherent losses are from born trapped full energy beam-ions being non-resonantly scattered by AEs onto loss orbits within their first poloidal transit. The loss amplitudes scale linearly with the mode amplitude but the slope is different for different modes. The radial displacement of fast-ions by individual AEs can be directly inferred from the measurements. Second, oscillations in the beam-ion first-orbit losses are observed at the sum, difference, and harmonic frequencies of two independent AEs. These oscillations are not plasma modes and are absent in magnetic, density, and temperature fluctuations. The origin of the non-linearity as a wave-particle coupling is confirmed through bi-coherence analysis, which is clearly observed because the coherences are preserved by the first-orbit loss mechanism. An analytic model and full orbit simulations show that the non-linear features seen in the loss signal can be explained by a non-linear interaction between the fast ions and the two independent AEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Linear calculations of edge current driven kink modes with BOUT++ code

This work extends previous BOUT++ work to systematically study the impact of edge current density on edge localized modes, and to benchmark with the GATO and ELITE codes. Using the CORSICA code, a set of equilibria was generated with different edge current densities by keeping total current and pressure profile fixed. Based on these equilibria, the effects of the edge current density on the MHD instabilities were studied with the 3-field BOUT++ code. For the linear calculations, with increasing edge current density, the dominant modes are changed from intermediate-n and high-n ballooning modes to low-n kink modes, and the linear growth rate becomes smaller. The edge current provides stabilizing effects on ballooning modes due to the increase of local shear at the outer mid-plane with the edge current. For edge kink modes, however, the edge current does not always provide a destabilizing effect; with increasing edge current, the linear growth rate first increases, and then decreases. In benchmark calculations for BOUT++ against the linear results with the GATO and ELITE codes, the vacuum model has important effects on the edge kink mode calculations. By setting a realistic density profile and Spitzer resistivity profile in the vacuum region, the resistivity was found to have a destabilizing effect on both the kink mode and on the ballooning mode. With diamagnetic effects included, the intermediate-n and high-n ballooning modes can be totally stabilized for finite edge current density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The linear and non-linear characterization of dust ion acoustic mode in complex plasma in presence of dynamical charging of dust

A systematic theoretical investigation has been carried out on the role of dust charging dynamics on the nature and stability of DIA (Dust Ion Acoustic) mode in complex plasma. The study has been made for both linear and non-linear scale regime of DIA mode. The observed results have been characterized in terms of background plasma responses towards dust surface responsible for dust charge fluctuation, invoking important dusty plasma parameters, especially the ion flow speed and dust size. The linear analyses confirm the nature of instability in DIA mode in presence of dust charge fluctuation. The instability shows a damping of DIA mode in subsonic flow regime followed by a gradual growth in instability in supersonic limit of ion flow. The strength of non-linearity and their existence domain is found to be driven by different dusty plasma parameters. As dust is ubiquitous in interstellar medium with plasma background, the study also addresses the possible effect of dust charging dynamics in gravito-electrostatic characterization and the stability of dust molecular clouds especially in proto-planetary disc. The observations are influential and interesting towards the understanding of dust settling mechanism and formation of dust environments in different regions in space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Novel linear analysis for a gyrotron oscillator based on a spectral approach

With the aim of gaining a better physical insight into linear regimes in gyrotrons, a new linear model was developed. This model is based on a spectral approach for solving the self-consistent system of equations describing the wave-particle interaction in the cavity of a gyrotron oscillator. Taking into account the wall-losses self-consistently and including the main system inhomogeneities in the cavity geometry and in the magnetic field, the model is appropriate to consider real system parameters. The main advantage of the spectral approach, compared with a time-dependent approach, is the possibility to describe all of the stable and unstable modes, respectively, with negative and positive growth rates. This permits to reveal the existence of a new set of eigenmodes, in addition to the usual eigenmodes issued from cold-cavity modes. The proposed model can be used for studying other instabilities such as, for instance, backward waves potentially excited in gyrotron beam tunnels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Gyrokinetic neoclassical study of the bootstrap current in the tokamak edge pedestal with fully non-linear Coulomb collisions

As a follow-up on the drift-kinetic study of the non-local bootstrap current in the steep edge pedestal of tokamak plasma by Koh et al. [Phys. Plasmas 19, 072505 (2012)], a gyrokinetic neoclassical study is performed with gyrokinetic ions and drift-kinetic electrons. Besides the gyrokinetic improvement of ion physics from the drift-kinetic treatment, a fully non-linear Fokker-Planck collision operator—that conserves mass, momentum, and energy—is used instead of Koh et al.'s linearized collision operator in consideration of the possibility that the ion distribution function is non-Maxwellian in the steep pedestal. An inaccuracy in Koh et al.'s result is found in the steep edge pedestal that originated from a small error in the collisional momentum conservation. The present study concludes that (1) the bootstrap current in the steep edge pedestal is generally smaller than what has been predicted from the small banana-width (local) approximation [e.g., Sauter et al., Phys. Plasmas 6, 2834 (1999) and Belli et al., Plasma Phys. Controlled Fusion 50, 095010 (2008)], (2) the plasma flow evaluated from the local approximation can significantly deviate from the non-local results, and (3) the bootstrap current in the edge pedestal, where the passing particle region is small, can be dominantly carried by the trapped particles in a broad trapped boundary layer. A new analytic formula based on numerous gyrokinetic simulations using various magnetic equilibria and plasma profiles with self-consistent Grad-Shafranov solutions is constructed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic reconnection in 3D fusion devices: non-linear reduced equations and linear current-driven instabilities

Abstract Magnetic reconnection in 3D fusion devices is investigated. With the use of Boozer co-ordinates, we reduce the non-linear resistive magnetohydrodynamic equations in the limit of large aspect ratio and finite pressure fluctuations, to obtain a set of non-linear equations suitable for magnetic reconnection studies in stellarators. Magnetic flux unfreezing due to a finite electron mass is also considered. Equations that govern the linear regime and some of their general properties are given. We emphasise the role of magnetic geometry and identify how some aspects of stellarator optimisation could have an impact on reconnecting instabilities, in particular by exacerbating those enabled by electron inertia. The effect of 3D coupling on the linear reconnection rates and the mode structure is quantitatively addressed in the case in which the equilibrium rotational transform has one specific resonant location for which one mode can reconnect while coupled to an arbitrary number of non-resonant harmonics. The full problem is rigorously reduced to an equivalent cylindrical one, by introducing some geometrically modified plasma inertial and dissipative scales. The 3D scalings for the growth rates of reconnection instabilities and their destabilisation criteria are given.

Physics↗
DOE OSTI2016

The linear -- non-linear frontier for the Goldstone Higgs

The minimal $SO(5)/SO(4)$ sigma model is used as a template for the ultraviolet completion of scenarios in which the Higgs particle is a low-energy remnant of some high-energy dynamics, enjoying a (pseudo) Nambu-Goldstone boson ancestry. Varying the $\sigma$ mass allows to sweep from the perturbative regime to the customary non-linear implementations. The low-energy benchmark effective non-linear Lagrangian for bosons and fermions is obtained, determining as well the operator coefficients including linear corrections. At first order in the latter, three effective bosonic operators emerge which are independent of the explicit soft breaking assumed. The Higgs couplings to vector bosons and fermions turn out to be quite universal: the linear corrections are proportional to the explicit symmetry breaking parameters. Furthermore, we define an effective Yukawa operator which allows a simple parametrization and comparison of different heavy fermion ultraviolet completions. In addition, one particular fermionic completion is explored in detail, obtaining the corresponding leading low-energy fermionic operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Designing linear lattices for round beam in electron storage rings using the solution by linear matrices analysis

For some synchrotron light source beamline applications, a round beam is preferable to a flat one. A conventional method of obtaining round beam in an electron storage ring is to shift its tune close to a linear difference resonance. The linearly coupled beam dynamics is analyzed with perturbation theories, which have certain limitations. In this paper, we adopt the Solution by LInear Matrices (SLIM) analysis to calculate exact beam sizes to design round beam lattices. The SLIM analysis can deal with a generally linearly coupled accelerator lattice. The effects of various coupling sources on beam emittances and sizes can be studied within a self-consistent frame. Both the on- and off-resonance schemes to obtain round beams are explained with examples. The SLIM formalism for two widely used magnet models: combined-function bending magnets, and planar wigglers and undulators, is also derived.

43 PARTICLE ACCELERATORS↗
NASA NTRS1987

User's manual for LINEAR, a FORTRAN program to derive linear aircraft models

This report documents a FORTRAN program that provides a powerful and flexible tool for the linearization of aircraft models. The program LINEAR numerically determines a linear system model using nonlinear equations of motion and a user-supplied nonlinear aerodynamic model. The system model determined by LINEAR consists of matrices for both state and observation equations. The program has been designed to allow easy selection and definition of the state, control, and observation variables to be used in a particular model.

Duke, Eugene L.↗
DOE OSTI2015

A simple, stable, and accurate linear tetrahedral finite element for transient, nearly, and fully incompressible solid dynamics: A dynamic variational multiscale approach [A simple, stable, and accurate tetrahedral finite element for transient, nearly incompressible, linear and nonlinear elasticity: A dynamic variational multiscale approach]

Here, we propose a new approach for the stabilization of linear tetrahedral finite elements in the case of nearly incompressible transient solid dynamics computations. Our method is based on a mixed formulation, in which the momentum equation is complemented by a rate equation for the evolution of the pressure field, approximated with piece-wise linear, continuous finite element functions. The pressure equation is stabilized to prevent spurious pressure oscillations in computations. Incidentally, it is also shown that many stabilized methods previously developed for the static case do not generalize easily to transient dynamics. Extensive tests in the context of linear and nonlinear elasticity are used to corroborate the claim that the proposed method is robust, stable, and accurate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Linear solvers for power grid optimization problems: A review of GPU-accelerated linear solvers

The linear equations that arise in interior methods for constrained optimization are sparse symmetric indefinite, and they become extremely ill-conditioned as the interior method converges. These linear systems present a challenge for existing solver frameworks based on sparse LU or LDL T decompositions. Here, we benchmark five well known direct linear solver packages on CPU- and GPU-based hardware, using matrices extracted from power grid optimization problems. The achieved solution accuracy varies greatly among the packages. None of the tested packages delivers significant GPU acceleration for our test cases. For completeness of the comparison we include results for MA57, which is one of the most efficient and reliable CPU solvers for this class of problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Linear optical response of finite systems using multishift linear system solvers

We discuss the application of multishift linear system solvers to linear-response time-dependent density functional theory. Using this technique the complete frequency-dependent electronic density response of finite systems to an external perturbation can be calculated at the cost of a single solution of a linear system via conjugate gradients. We show that multishift time-dependent density functional theory yields excitation energies and oscillator strengths in perfect agreement with the standard diagonalization of the response matrix (Casida's method), while being computationally advantageous. We present test calculations for benzene, porphin, and chlorophyll molecules. We argue that multishift solvers may find broad applicability in the context of excited-state calculations within density-functional theory and beyond.

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

Investigations of the linear and non-linear flow harmonics using the a multi-phase transport model

The multi-phase transport model (AMPT) is used to study the effects of the parton-scattering cross-sections (σ pp ) and hadronic re-scattering on the linear contributions to the flow harmonic v 4 , the non-linear response coefficients, and the correlations between different order flow symmetry planes in Au + Au collisions at 200 GeV. The model results, which agree with current experimental measurements, indicate that the higher-order flow harmonics are sensitive to the σ pp variations. However, the non-linear response coefficients and the correlations between different order flow symmetry planes are σ pp independent. Furthermore, these results suggest that further detailed experimental measurements which span a broad range of collision systems and beam energies could serve as an additional constraint for the theoretical models' calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗