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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Conjugacy Analysis of the Finite-Difference Equations of the Method of Characteristics on the Basis of a Linear-Algebraic Formulation

A linear-algebraic form of the equations of the method of characteristics, which is used to approximate the neutron transport equation, is obtained. It is shown on the basis of the obtained linear-algebraic form that the discrete form of the conjugate equation differs from the algebraically discrete problem constructed by linear-algebraic transformations of the discrete form of the normal problem. The reason for the discrepancy lies in the approximation of the volumes of the spatial cells in covering the working region by a network of characteristics. It is shown by means of test calculations that when the network of characteristics is refined the solution of the conjugate transport equation converges to the solution of the algebraically conjugate problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Excited State and Non-linear Optical Properties of NIR Absorbing β-Thiophene-Fused BF{sub 2}-Azadipyrromethene Dyes—Computational Investigation

Density functional theory and time-dependent density functional theory computations were used to understand the electronic and photophysical parameters of NIR β-thiophene-fused BF{sub 2}-azadipyrromethene dyes. The computed data are in good agreement with those obtained experimentally and they provide insights into the origin of red shifted optical spectra compared to the parent aza-BODIPY, low Stokes shift, non-linear optical responses and quantitative description of the singlet–triplet energy gap. The resultant decrease in the HOMO − LUMO energy gap is responsible for the red shift. The possible use as non-linear optical materials is supported by large enhancement in the non-linear optical properties. On the basis of vertical triplet energies, their possible potential therapeutic use as a photosensitizer in photodynamic therapy is proposed. The singlet–triplet energetic gaps suggest that the β-thiophene-fused BF{sub 2}-azadipyrromethene dyes can act as a sensitizer to produce an efficient generation of singlet oxygen. Their optimal use as an efficient singlet fission materials has been proposed on the basis of excitation energies in the ground, lowest singlet and triplet excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Piecewise Linear Dynamical System Modeling Gene Network with Variable Feedback

We construct a discretization of the phase portrait of a 3-dimensional dynamical system of biochemical kinetics with piecewise linear right-hand sides. We describe geometry of the phase portrait and construct an invariant piecewise linear surface bounded by a stable cycle of this system composed of eight linear segments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Linear Reciprocating Tribometer for In Situ Neutron Reflectometry of Soft Matter

Neutron reflectometry is a technique for measuring structure near planar interfaces that has been previously used to non-destructively characterize the polymer density of hydrated, dilute, and soft materials. Previous investigations have conducted neutron reflectometry measurements of liquids, gels, emulsion, and polymer solutions at rest, in compression, and subject to shear stress. However, correlating structure with tribological properties of soft materials presents significant experimental challenges for prior instruments due to wall slip, sample thickness, and structural heterogeneity (e.g., depth-wise gradients). A linear reciprocating tribometer offers several advantages for in situ neutron reflectometry studies, including uniform velocity profiles, constant shear stress over large regions of interest, and independent control of normal force and sliding velocity during measurements. This work outlines basic considerations for the design of a custom linear reciprocating tribometer that operates in a neutron beamline and includes commissioning measurements. The tribometer is designed to compress soft and hydrated materials against linearly reciprocating silicon disks. The three key design considerations for this tribometer are (1) safety, (2) neutron transmission, and (3) sample positioning. This instrument design will enable in situ studies of soft matter and illuminate the role of interfacial structure on tribological phenomena.

42 ENGINEERING↗
DOE OSTI2018

Determination of Collision Cross Sections Using a Fourier Transform Electrostatic Linear Ion Trap Mass Spectrometer

Collision cross sections (CCSs) were determined from the frequency-domain linewidths in a Fourier transform electrostatic linear ion trap. With use of an ultrahigh-vacuum precision leak valve and nitrogen gas, transients were recorded as the background pressure in the mass analyzer chamber was varied between 4× 10{sup -8} and 7 × 10{sup -7} Torr. The energetic hard-sphere ion–neutral collision model, described by Xu and coworkers, was used to relate the recorded image charge to the CCS of the molecule. In lieu of our monoisotopically isolating the mass of interest, the known relative isotopic abundances were programmed into the Lorentzian fitting algorithm such that the linewidth was extracted from a sum of Lorentzians. Although this works only if the isotopic distribution is known a priori, it prevents ion loss, preserves the high signal-to-noise ratio, and minimizes the experimental error on our homebuilt instrument. Six tetraalkylammonium cations were used to correlate the CCS measured in the electrostatic linear ion trap with that measured by drift-tube ion mobility spectrometry, for which there was an excellent correlation (R{sup 2} ≈ 0.9999). Although the absolute CCSs derived with our method differ from those reported, the extracted linear correlation can be used to correct the raw CCS. With use of [angiotensin II]{sup 2+} and reserpine, the corrected CCSs (334.9 ± 2.1 and 250.1 ± 0.5, respectively) were in good agreement with the reported ion mobility spectrometry CCSs (335 and 254.3, respectively). With sufficient signal-to-noise ratio, the CCSs determined are reproducible to within a fraction of a percent, comparable to the uncertainties reported on dedicated ion mobility instruments. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

A fast eigenvalue approach for solving the trust region subproblem with an additional linear inequality

In this paper, we study the extended trust region subproblem (eTRS) in which the trust region intersects the Euclidean ball with a single linear inequality constraint. By reformulating the Lagrangian dual of eTRS as a two-parameter linear eigenvalue problem, we state a necessary and sufficient condition for its strong duality in terms of an optimal solution of a linearly constrained bivariate concave maximization problem. This results in an efficient algorithm for solving eTRS of large size whenever the strong duality is detected. Finally, some numerical experiments are given to show the effectiveness of the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Automated Model Updating Effect on the Linear and Nonlinear Dynamic Responses of Historical Masonry Structures

Finite element models supply the simulation of the problems to obtain the structural behaviors under different loading conditions. The updating processes based on changing the uncertain parameters such as material properties, boundary conditions, meshing size etc. are made to create actual finite element models in order to represent the real conditions of the problems. This procedure can be made using two methods, namely manual model updating and global/local automatic model updating. The global/local automatic model updating procedures are very popular owing to the applicability to all kind of engineering structures, minimizing the differences nearly to zero, availability of damage localization and structural health monitoring. In this paper a historical masonry armory building selected as an application and its finite element analyses, experimental measurements and automated model updating processes are carried out. The finite element model of the structures was developed with ANSYS software and first three natural frequencies are attained. The experimental dynamic characteristics such as natural frequencies, damping ratios and mode shapes are obtained using ambient vibration method. The differences between experimental and numerical frequencies are minimized up to the range of 0–1%, by using local updating procedure. The linear and nonlinear seismic analyses of the building are conducted on initial and updated models. At the end of the analyses, maximum displacements, principal stresses and strains are given with detail. In addition, crack distributions diagrams are attained to validate the stresses accumulation points. It can be seen that there is an increasing trend in the displacement, stress and strain values after finite element model updating both linear and nonlinear analysis. In addition, linear and nonlinear analysis results are compared with each other to evaluate the nonlinear analysis effect. It is concluded that the displacement and strain values are increased after nonlinear analysis. But, there is not any agreement in stress values.

42 ENGINEERING↗
DOE OSTI2022

Non-intrusive reduced order modeling of natural convection in porous media using convolutional autoencoders: Comparison with linear subspace techniques

Natural convection in porous media is a highly nonlinear multiphysical problem relevant to many engineering applications (e.g., the process of CO 2 sequestration). Here, we extend and present a non-intrusive reduced order model of natural convection in porous media employing deep convolutional autoencoders for the compression and reconstruction and either radial basis function (RBF) interpolation or artificial neural networks (ANNs) for mapping parameters of partial differential equations (PDEs) on the corresponding nonlinear manifolds. To benchmark our approach, we also describe linear compression and reconstruction processes relying on proper orthogonal decomposition (POD) and ANNs. Further, we present comprehensive comparisons among different models through three benchmark problems. The reduced order models, linear and nonlinear approaches, are much faster than the finite element model, obtaining a maximum speed-up of 7 × 10 6 because our framework is not bound by the Courant–Friedrichs–Lewy condition; hence, it could deliver quantities of interest at any given time contrary to the finite element model. Our model’s accuracy still lies within a relative error of 7% in the worst-case scenario. We illustrate that, in specific settings, the nonlinear approach outperforms its linear counterpart and vice versa. We hypothesize that a visual comparison between principal component analysis (PCA) and t-Distributed Stochastic Neighbor Embedding (t-SNE) could indicate which method will perform better prior to employing any specific compression strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Linear response of tripartite entanglement to infinitesimal noise

Recent experimental progress in prolonging the coherence time of a quantum system prompts us to explore the behavior of quantum entanglement at the beginning of the decoherence process. The response of the entanglement under an infinitesimal noise can serve as a signature of the robustness of entangled states. A crucial problem of this topic in multipartite systems is to compute the degree of entanglement in a mixed state. We find a family of global noise in three-qubit systems, which is composed of four W states. Under its influence, the linear response of the tripartite entanglement of a symmetrical three-qubit pure state is studied. A lower bound of the linear response is found to depend completely on the initial tripartite and bipartite entanglement. This result shows that the decay of tripartite entanglement is hastened by the bipartite one. - Highlights: • We study a set of W-type noise and its linear effect on symmetric pure states. • Its effect on two-qubit entanglement depends only on the initial concurrence. • A lower bound of the effect on 3-tangle is found in terms of initial entanglements. • We obtain the time of three-tangle sudden death for two families of typical states. • These reveal that the bipartite entanglement speeds up the decay of the tripartite one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Deterministic implementations of single-photon multi-qubit Deutsch–Jozsa algorithms with linear optics

It is very important to seek an efficient and robust quantum algorithm demanding less quantum resources. We propose one-photon three-qubit original and refined Deutsch–Jozsa algorithms with polarization and two linear momentums degrees of freedom (DOFs). Our schemes are constructed by solely using linear optics. Compared to the traditional ones with one DOF, our schemes are more economic and robust because the necessary photons are reduced from three to one. Our linear-optic schemes are working in a determinate way, and they are feasible with current experimental technology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Anthraquinones quinizarin and danthron unwind negatively supercoiled DNA and lengthen linear DNA

Highlights: • Anthraquinones quinizarin and danthron unwind negatively supercoiled DNA. • Anthraquinones quinizarin and danthron lengthen linear DNA. • Anthraquinones quinizarin and danthron possess middle binding affinity to DNA. • Anthraquinones quinizarin and danthron interact with DNA by intercalating mode. - Abstract: The intercalating drugs possess a planar aromatic chromophore unit by which they insert between DNA bases causing the distortion of classical B-DNA form. The planar tricyclic structure of anthraquinones belongs to the group of chromophore units and enables anthraquinones to bind to DNA by intercalating mode. The interactions of simple derivatives of anthraquinone, quinizarin (1,4-dihydroxyanthraquinone) and danthron (1,8-dihydroxyanthraquinone), with negatively supercoiled and linear DNA were investigated using a combination of the electrophoretic methods, fluorescence spectrophotometry and single molecule technique an atomic force microscopy. The detection of the topological change of negatively supercoiled plasmid DNA, unwinding of negatively supercoiled DNA, corresponding to appearance of DNA topoisomers with the low superhelicity and an increase of the contour length of linear DNA in the presence of quinizarin and danthron indicate the binding of both anthraquinones to DNA by intercalating mode.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Predicting non-linear stress–strain response of mesostructured cellular materials using supervised autoencoder

Recent breakthroughs in advanced manufacturing capabilities have made it possible to design and print sophisticated topologies of cellular structures using diverse engineering materials such as metals, polymers, and ceramics. In these architectured materials, it is often desirable to tailor the mechanical properties by altering the unit cell topology. This necessitates an in-depth understanding of how the topology of the unit cell structure affects the macroscopic behavior of the material in both the linear and the non-linear regimes encountered under large compression. Here, we have developed a machine learning (ML) approach capable of accelerating the prediction of the stress–strain response of a polymer-based cellular structure under uniaxial confined compression. As part of generating the training data for ML, 60,000 mesostructures were generated using a relatively novel approach based on cellular automata, and their corresponding stress–strain responses were obtained from the finite element simulations. Principal component analysis (PCA) was used to reduce the dimensionality of the stress–strain curves. With only 20 principal components, PCA captured 99.89% of the variance in the stress–strain curves while reducing the dimensionality by 5X. ML using supervised autoencoder was able to successfully speed up the prediction of the non-linear stress–strain response of a unit cell by up to 4600X. The proposed method can serve as an efficient data generation tool and a rapid means for predicting the structure–property relationship through accelerated forward modeling of cellular materials under compaction, in cases where the macroscopic stress–strain response is governed by the unit-cell topology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluation of alkali–silica reaction damage in concrete using linear and nonlinear resonance techniques

Here in this study, the initiation and development of alkali–silica reaction (ASR) in concrete prisms containing reactive coarse (RCA) and reactive fine (RFA) aggregates were monitored over one year using linear and nonlinear resonance techniques. Although the linear resonance test could detect initiation and development of ASR damage in RCA specimens, it failed in the RFA specimens. Petrographic analysis confirmed the presence of early-stage ASR damage and the presence of different microcracking patterns in the RFA and RCA specimens, which may explain differences in the linear resonance data. For both RCA and RFA specimens, the nonlinear resonance test results were able to detect the initiation of ASR damage and showed a strong correlation with the ASR expansion for all specimens from different mix designs. The results demonstrate the high sensitivity of the nonlinear resonance test and its potential for quantitative evaluation of the ASR damage state of concrete specimens without a baseline measurement.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Linearized self-consistent quasiparticle GW method: Application to semiconductors and simple metals

We present a code implementing the linearized self-consistent quasiparticle GW method (QSGW) in the LAPW basis. Our approach is based on the linearization of the self-energy around zero frequency which differs it from the existing implementations of the QSGW method. The linearization allows us to use Matsubara frequencies instead of working on the real axis. This results in efficiency gains by switching to the imaginary time representation in the same way as in the space time method. The all electron LAPW basis set eliminates the need for pseudopotentials. We discuss the advantages of our approach, such as its N 3 scaling with the system size N, as well as its shortcomings. We apply our approach to study the electronic properties of selected semiconductors, insulators, and simple metals and show that our code produces the results very close to the previously published QSGW data. Our implementation is a good platform for further many body diagrammatic resummations such as the vertex-corrected GW approach and the GW+DMFT method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Subsurface failure in spherical bodies. A formation scenario for linear troughs on Vesta’s surface

Many asteroids in the Solar System exhibit unusual, linear features on their surface. The Dawn mission recently observed two sets of linear features on the surface of the asteroid 4 Vesta. Geologic observations indicate that these features are related to the two large impact basins at the south pole of Vesta, though no specific mechanism of origin has been determined. Furthermore, the orientation of the features is offset from the center of the basins. Experimental and numerical results reveal that the offset angle is a natural consequence of oblique impacts into a spherical target. We demonstrate that a set of shear planes develops in the subsurface of the body opposite to the point of first contact. Moreover, these subsurface failure zones then propagate to the surface under combined tensile-shear stress fields after the impact to create sets of approximately linear faults on the surface. Comparison between the orientation of damage structures in the laboratory and failure regions within Vesta can be used to constrain impact parameters (e.g., the approximate impact point and likely impact trajectory).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal model

In this paper, we develop a series of linear, unconditionally energy stable numerical schemes for solving the classical phase field crystal model. The temporal discretizations are based on the first order Euler method, the second order backward differentiation formulas (BDF2) and the second order Crank–Nicolson method, respectively. The schemes lead to linear elliptic equations to be solved at each time step, and the induced linear systems are symmetric positive definite. We prove that all three schemes are unconditionally energy stable rigorously. Various classical numerical experiments in 2D and 3D are performed to validate the accuracy and efficiency of the proposed schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

The PAC-MAN model: Benchmark case for linear acoustics in computational physics

Benchmark cases in the field of computational physics, on the one hand, have to contain a certain complexity to test numerical edge cases and, on the other hand, require the existence of an analytical solution, because an analytical solution allows the exact quantification of the accuracy of a numerical simulation method. This dilemma causes a need for analytical sound field formulations of complex acoustic problems. A well known example for such a benchmark case for harmonic linear acoustics is the “Cat's Eye model”, which describes the three-dimensional sound field radiated from a sphere with a missing octant analytically. In this paper, a benchmark case for two-dimensional (2D) harmonic linear acoustic problems, viz., the “PAC-MAN model”, is proposed. The PAC-MAN model describes the radiated and scattered sound field around an infinitely long cylinder with a cut out sector of variable angular width. While the analytical calculation of the 2D sound field allows different angular cut-out widths and arbitrarily positioned line sources, the computational cost associated with the solution of this problem is similar to a 1D problem because of a modal formulation of the sound field in the PAC-MAN model. - Highlights: • A benchmark case for 2D harmonic linear acoustic problems is proposed. • The PAC-MAN geometry is used as emitter/scatterer in this benchmark. • An analytic formulation of radiated and scattered sound is derived. • Line and disk sources, plane waves, and surface vibrations are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

An efficient explicit implementation of a near-optimal quantum algorithm for simulating linear dissipative differential equations

We propose an efficient block-encoding technique for the implementation of the Linear Combination of Hamiltonian Simulations (LCHS) for simulating dissipative initial-value problems. This algorithm approximates a target nonunitary operator as a weighted sum of Hamiltonian evolutions, thereby emulating a dissipative problem by mixing various time scales. We introduce an efficient encoding of the LCHS into a quantum circuit based on a simple coordinate transformation that turns the dependence on the summation index into a trigonometric function. Classically, this method is equivalent to the use of a highly accurate Fejér-Clenshaw-Curtis quadrature formula. Quantumly, this significantly simplifies block-encoding of a dissipative problem and allows one to perform an exponential number of Hamiltonian simulations by a single Quantum Signal Processing (QSP) circuit. The resulting LCHS circuit has high success probability and the selector scales logarithmically with the number of terms in the LCHS sum and linearly with time. Careful analysis of error convergence proves that this method is more efficient than other LCHS circuits that have recently appeared in the literature. We verify the quantum circuit and its scaling by simulating it on a digital emulator of fault-tolerant quantum computers and, as a test problem, solve the advection-diffusion equation. The proposed algorithm can be used for simulating a wide class of nonunitary initial-value problems including the Liouville equation with added dissipation and linear embeddings of nonlinear systems, such as the Koopman-von Neumann and Carleman embeddings.

Novikau, I [Lawrence Livermore National Laboratory↗