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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

New non-linear photovoltaic effect in uniform bipolar semiconductor

A linear theory of the new non-linear photovoltaic effect in the closed circuit consisting of a non-uniformly illuminated uniform bipolar semiconductor with neutral impurities is developed. The non-uniform photo-excitation of impurities results in the position-dependant current carrier mobility that breaks the semiconductor homogeneity and induces the photo-electromotive force (emf). As both the electron (or hole) mobility gradient and the current carrier generation rate depend on the light intensity, the photo-emf and the short-circuit current prove to be non-linear functions of the incident light intensity at an arbitrarily low illumination. The influence of the sample size on the photovoltaic effect magnitude is studied. Physical relations and distinctions between the considered effect and the Dember and bulk photovoltaic effects are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Beam envelope calculations in general linear coupled lattices

The envelope equations and Twiss parameters (β and α) provide important bases for uncoupled linear beam dynamics. For sophisticated beam manipulations, however, coupling elements between two transverse planes are intentionally introduced. The recently developed generalized Courant-Snyder theory offers an effective way of describing the linear beam dynamics in such coupled systems with a remarkably similar mathematical structure to the original Courant-Snyder theory. In this work, we present numerical solutions to the symmetrized matrix envelope equation for β which removes the gauge freedom in the matrix envelope equation for w. Furthermore, we construct the transfer and beam matrices in terms of the generalized Twiss parameters, which enables calculation of the beam envelopes in arbitrary linear coupled systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Stochastic linear multistep methods for the simulation of chemical kinetics

In this paper, we introduce the Stochastic Adams-Bashforth (SAB) and Stochastic Adams-Moulton (SAM) methods as an extension of the τ-leaping framework to past information. Using the Θ-trapezoidal τ-leap method of weak order two as a starting procedure, we show that the k-step SAB method with k ≥ 3 is order three in the mean and correlation, while a predictor-corrector implementation of the SAM method is weak order three in the mean but only order one in the correlation. These convergence results have been derived analytically for linear problems and successfully tested numerically for both linear and non-linear systems. A series of additional examples have been implemented in order to demonstrate the efficacy of this approach.

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

Tunnel magnetoresistance and linear conductance of double quantum dots strongly coupled to ferromagnetic leads

We analyze the spin-dependent linear-response transport properties of double quantum dots strongly coupled to external ferromagnetic leads. By using the numerical renormalization group method, we determine the dependence of the linear conductance and tunnel magnetoresistance on the degree of spin polarization of the leads and the position of the double dot levels. We focus on the transport regime where the system exhibits the SU(4) Kondo effect. It is shown that the presence of ferromagnets generally leads the suppression of the linear conductance due to the presence of an exchange field. Moreover, the exchange field gives rise to a transition from the SU(4) to the orbital SU(2) Kondo effect. We also analyze the dependence of the tunnel magnetoresistance on the double dot levels' positions and show that it exhibits a very nontrivial behavior.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Internal transport barrier triggered by non-linear lower hybrid wave deposition under condition of beam-driven toroidal rotation

By using gyro-Landau fluid transport model (GLF23), time-dependent integrated modeling is carried out using TRANSP to explore the dynamic process of internal transport barrier (ITB) formation in the neutral beam heating discharges. When the current profile is controlled by LHCD (lower hybrid current drive), with appropriate neutral beam injection, the nonlinear interplay between the transport determined gradients in the plasma temperature (T{sub i,e}) and toroidal velocity (V{sub ϕ}) and the E×B flow shear (including q-profile) produces transport bifurcations, generating spontaneously a stepwise growing ITB. In the discharge, the constraints imposed by the wave propagation condition causes interplay of the LH driven current distribution with the plasma configuration modification, which constitutes non-linearity in the LH wave deposition. The non-linear effects cause bifurcation in LHCD, generating two distinct quasi-stationary reversed magnetic shear configurations. The change of current profile during the transition period between the two quasi-stationary states results in increase of the E×B shearing flow arising from toroidal rotation. The turbulence transport suppression by sheared E×B flow during the ITB development is analysed, and the temporal evolution of some parameters characterized the plasma confinement is examined. Ample evidence shows that onset of the ITB development is correlated with the enhancement of E×B shearing rate caused by the bifurcation in LHCD. It is suggested that the ITB triggering is associated with the non-linear effects of the LH power deposition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Synchronization, non-linear dynamics and low-frequency fluctuations: Analogy between spontaneous brain activity and networked single-transistor chaotic oscillators

In this paper, the topographical relationship between functional connectivity (intended as inter-regional synchronization), spectral and non-linear dynamical properties across cortical areas of the healthy human brain is considered. Based upon functional MRI acquisitions of spontaneous activity during wakeful idleness, node degree maps are determined by thresholding the temporal correlation coefficient among all voxel pairs. In addition, for individual voxel time-series, the relative amplitude of low-frequency fluctuations and the correlation dimension (D{sub 2}), determined with respect to Fourier amplitude and value distribution matched surrogate data, are measured. Across cortical areas, high node degree is associated with a shift towards lower frequency activity and, compared to surrogate data, clearer saturation to a lower correlation dimension, suggesting presence of non-linear structure. An attempt to recapitulate this relationship in a network of single-transistor oscillators is made, based on a diffusive ring (n = 90) with added long-distance links defining four extended hub regions. Similarly to the brain data, it is found that oscillators in the hub regions generate signals with larger low-frequency cycle amplitude fluctuations and clearer saturation to a lower correlation dimension compared to surrogates. The effect emerges more markedly close to criticality. The homology observed between the two systems despite profound differences in scale, coupling mechanism and dynamics appears noteworthy. These experimental results motivate further investigation into the heterogeneity of cortical non-linear dynamics in relation to connectivity and underline the ability for small networks of single-transistor oscillators to recreate collective phenomena arising in much more complex biological systems, potentially representing a future platform for modelling disease-related changes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Non-perturbative aspects of particle acceleration in non-linear electrodynamics

We undertake an investigation of particle acceleration in the context of non-linear electrodynamics. We deduce the maximum energy that an electron can gain in a non-linear density wave in a magnetised plasma, and we show that an electron can “surf” a sufficiently intense Born-Infeld electromagnetic plane wave and be strongly accelerated by the wave. The first result is valid for a large class of physically reasonable modifications of the linear Maxwell equations, whilst the second result exploits the special mathematical structure of Born-Infeld theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Response analysis of a laminar premixed M-flame to flow perturbations using a linearized compressible Navier-Stokes solver

The response of a laminar premixed methane-air flame subjected to flow perturbations around a steady state is examined experimentally and using a linearized compressible Navier-Stokes solver with a one-step chemistry mechanism to describe combustion. The unperturbed flame takes an M-shape stabilized both by a central bluff body and by the external rim of a cylindrical nozzle. This base flow is computed by a nonlinear direct simulation of the steady reacting flow, and the flame topology is shown to qualitatively correspond to experiments conducted under comparable conditions. The flame is then subjected to acoustic disturbances produced at different locations in the numerical domain, and its response is examined using the linearized solver. This linear numerical model then allows the componentwise investigation of the effects of flow disturbances on unsteady combustion and the feedback from the flame on the unsteady flow field. It is shown that a wrinkled reaction layer produces hydrodynamic disturbances in the fresh reactant flow field that superimpose on the acoustic field. This phenomenon, observed in several experiments, is fully interpreted here. The additional perturbations convected by the mean flow stem from the feedback of the perturbed flame sheet dynamics onto the flow field by a mechanism similar to that of a perturbed vortex sheet. The different regimes where this mechanism prevails are investigated by examining the phase and group velocities of flow disturbances along an axis oriented along the main direction of the flow in the fresh reactant flow field. It is shown that this mechanism dominates the low-frequency response of the wrinkled shape taken by the flame and, in particular, that it fully determines the dynamics of the flame tip from where the bulk of noise is radiated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Improved linearity in AlGaN/GaN metal-insulator-semiconductor high electron mobility transistors with nonlinear polarization dielectric

We demonstrate highly improved linearity in a nonlinear ferroelectric of Pb(Zr{sub 0.52}Ti{sub 0.48})-gated AlGaN/GaN metal-insulator-semiconductor high electron mobility transistor (MIS-HEMT). Distinct double-hump feature in the transconductance-gate voltage (g{sub m}-V{sub g}) curve is observed, yielding remarkable enhancement in gate voltage swing as compared to MIS-HEMT with conventional linear gate dielectric. By incorporating the ferroelectric polarization into a self-consistent calculation, it is disclosed that in addition to the common hump corresponding to the onset of electron accumulation, the second hump at high current level is originated from the nonlinear polar nature of ferroelectric, which enhances the gate capacitance by increasing equivalent dielectric constant nonlinearly. This work paves a way for design of high linearity GaN MIS-HEMT by exploiting the nonlinear properties of dielectric.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations

Fundamental understanding of complex dynamics in many-particle systems on the atomistic level is of utmost importance. Often the systems of interest are of macroscopic size but can be partitioned into a few important degrees of freedom which are treated most accurately and others which constitute a thermal bath. Particular attention in this respect attracts the linear generalized Langevin equation, which can be rigorously derived by means of a linear projection technique. Within this framework, a complicated interaction with the bath can be reduced to a single memory kernel. This memory kernel in turn is parametrized for a particular system studied, usually by means of time-domain methods based on explicit molecular dynamics data. Here, we discuss that this task is more naturally achieved in frequency domain and develop a Fourier-based parametrization method that outperforms its time-domain analogues. Very surprisingly, the widely used rigid bond method turns out to be inappropriate in general. Importantly, we show that the rigid bond approach leads to a systematic overestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.

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

Lattice model of linear telechelic polymer melts. I. Inclusion of chain semiflexibility in the lattice cluster theory

The lattice cluster theory (LCT) for the thermodynamics of polymer systems has recently been reformulated to treat strongly interacting self-assembling polymers composed of fully flexible linear telechelic chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. 136, 064902 (2012)]. Here, we further extend the LCT for linear telechelic polymer melts to include a description of chain semiflexibility, which is treated by introducing a bending energy penalty whenever a pair of consecutive bonds from a single chain lies along orthogonal directions. An analytical expression for the Helmholtz free energy is derived for the model of semiflexible linear telechelic polymer melts. The extension provides a theoretical tool for investigating the influence of chain stiffness on the thermodynamics of self-assembling telechelic polymers, and for further exploring the influence of self-assembly on glass formation in such systems.

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

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in quantum chemistry and beyond.

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

Linear and nonlinear heavy ion-acoustic waves in a strongly coupled plasma

A theoretical study on the propagation of linear and nonlinear heavy ion-acoustic (HIA) waves in an unmagnetized, collisionless, strongly coupled plasma system has been carried out. The plasma system is assumed to contain adiabatic positively charged inertial heavy ion fluids, nonextensive distributed electrons, and Maxwellian light ions. The normal mode analysis is used to study the linear behaviour. On the other hand, the well-known reductive perturbation technique is used to derive the nonlinear dynamical equations, namely, Burgers equation and Korteweg-de Vries (K-dV) equation. They are also numerically analyzed in order to investigate the basic features of shock and solitary waves. The adiabatic effects on the HIA shock and solitary waves propagating in such a strongly coupled plasma are taken into account. It has been observed that the roles of the adiabatic positively charged heavy ions, nonextensivity of electrons, and other plasma parameters arised in this investigation have significantly modified the basic features (viz., polarity, amplitude, width, etc.) of the HIA solitary/shock waves. The findings of our results obtained from this theoretical investigation may be useful in understanding the linear as well as nonlinear phenomena associated with the HIA waves both in space and laboratory plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Various notions of positivity for bi-linear maps and applications to tri-partite entanglement

We consider bi-linear analogues of s-positivity for linear maps. The dual objects of these notions can be described in terms of Schmidt ranks for tri-tensor products and Schmidt numbers for tri-partite quantum states. These tri-partite versions of Schmidt numbers cover various kinds of bi-separability, and so we may interpret witnesses for those in terms of bi-linear maps. We give concrete examples of witnesses for various kinds of three qubit entanglement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A widefield fluorescence microscope with a linear image sensor for image cytometry of biospecimens: Considerations for image quality optimization

Linear image sensors have been widely used in numerous research and industry applications to provide continuous imaging of moving objects. Here, we present a widefield fluorescence microscope with a linear image sensor used to image translating objects for image cytometry. First, a calibration curve was characterized for a custom microfluidic chamber over a span of volumetric pump rates. Image data were also acquired using 15 μm fluorescent polystyrene spheres on a slide with a motorized translation stage in order to match linear translation speed with line exposure periods to preserve the image aspect ratio. Aspect ratios were then calculated after imaging to ensure quality control of image data. Fluorescent beads were imaged in suspension flowing through the microfluidics chamber being pumped by a mechanical syringe pump at 16 μl min{sup −1} with a line exposure period of 150 μs. The line period was selected to acquire images of fluorescent beads with a 40 dB signal-to-background ratio. A motorized translation stage was then used to transport conventional glass slides of stained cellular biospecimens. Whole blood collected from healthy volunteers was stained with 0.02% (w/v) proflavine hemisulfate was imaged to highlight leukocyte morphology with a 1.56 mm × 1.28 mm field of view (1540 ms total acquisition time). Oral squamous cells were also collected from healthy volunteers and stained with 0.01% (w/v) proflavine hemisulfate to demonstrate quantifiable subcellular features and an average nuclear to cytoplasmic ratio of 0.03 (n = 75), with a resolution of 0.31 μm pixels{sup −1}.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Linearly scaling and almost Hamiltonian dielectric continuum molecular dynamics simulations through fast multipole expansions

Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved—up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.

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

Linear-scaling time-dependent density-functional theory beyond the Tamm-Dancoff approximation: Obtaining efficiency and accuracy with in situ optimised local orbitals

We present a solution of the full time-dependent density-functional theory (TDDFT) eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspaces with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate gradient algorithm that is very memory-efficient. The algorithm is validated on a small test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of bacteriochlorophyll in solution. The largest systems considered in this work are of the same order of magnitude as a variety of widely studied pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.

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

Efficient linear-scaling second-order Møller-Plesset perturbation theory: The divide–expand–consolidate RI-MP2 model

The Resolution of the Identity second-order Møller-Plesset perturbation theory (RI-MP2) method is implemented within the linear-scaling Divide-Expand-Consolidate (DEC) framework. In a DEC calculation, the full molecular correlated calculation is replaced by a set of independent fragment calculations each using a subset of the total orbital space. The number of independent fragment calculations scales linearly with the system size, rendering the method linear-scaling and massively parallel. The DEC-RI-MP2 method can be viewed as an approximation to the DEC-MP2 method where the RI approximation is utilized in each fragment calculation. The individual fragment calculations scale with the fifth power of the fragment size for both methods. However, the DEC-RI-MP2 method has a reduced prefactor compared to DEC-MP2 and is well-suited for implementation on massively parallel supercomputers, as demonstrated by test calculations on a set of medium-sized molecules. The DEC error control ensures that the standard RI-MP2 energy can be obtained to the predefined precision. The errors associated with the RI and DEC approximations are compared, and it is shown that the DEC-RI-MP2 method can be applied to systems far beyond the ones that can be treated with a conventional RI-MP2 implementation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗