Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORTHOGONAL TRANSFORMATIONS”

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 199 records · Page 11

Wideband Spectroscopy: The Design and Implementation of a 3 GHz Bandwidth, 8192 Channel, Polyphase Digital Spectrometer

A family of state-of-the-art digital Fourier transform spectrometers has been developed, with a combination of high bandwidth and fine resolution unavailable elsewhere. Analog signals consisting of radiation emitted by constituents in planetary atmospheres or galactic sources are downconverted and subsequently digitized by a pair of interleaved Analog-to-Digital Converters, (ADC). This 6 Gsps (giga-sample per second) digital representation of the analog signal is then processed through an FPGA-based streaming Fast Fourier Transform (FFT), the key development described below. Digital spectrometers have many advantages over previously used analog spectrometers, especially in terms of accuracy and resolution, both of which are particularly important for the type of scientific questions to be addressed with next-generation radiometers. the implementation, results and underlying math for this spectrometer, as well as, potential for future extension to even higher bandwidth, resolution and channel orthogonality, needed to support proposed future advanced atmospheric science and radioastronomy, are discussed.

field programmable gate array↗

Wideband Spectroscopy: The Design and Implementation of a 3 GHz, 2048 Channel Digital Spectrometer

A state-of-the-art digital Fourier Transform spectrometer has been developed, with a combination of high bandwidth and fine resolution unavailable elsewhere. Analog signals consisting of radiation emitted by constituents in planetary atmospheres or galactic sources are downconverted and subsequently digitized by a pair of interleaved Analog-to-Digital Converters (ADC). This 6 Gsps (giga sample per second) digital representation of the analog signal is then processed through an FPGA-based streaming Fast Fourier Transform (FFT), the key development described below. Digital spectrometers have many advantages over previously used analog spectrometers, especially in terms of accuracy and resolution, both of which are particularly important for the type of scientific questions to be addressed with next-generation radiometers. The implementation, results and underlying math for this spectrometer, as well as potential for future extension to even higher bandwidth, resolution and channel orthogonality, needed to support proposed future advanced atmospheric science and radioastronomy, are discussed.

field programmable gate array↗

A Novel Partitioned Approach for Reduced Order Model—Finite Element Model (ROM-FEM) and ROM-ROM Coupling

Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.

de Castro, Amy↗

Implementation of a 622 Mbps Digital Modem

This paper presents the implementation and initial test results of an Orthogonal Frequency Division Multiplexing (OFDM) digital modem (modulator and demodulator) with an aggregate information throughput of 622 megabits per second (Mbps). The OFDM waveform is constructed by dividing an incoming data stream into four channels, each channel using either a 16-ary Quadrature Amplitude Modulation (16QAM) scheme or an 8-Phase Shift Keying (8PSK) scheme. The generation and detection of the composite waveform are performed using Discrete Fourier Transform (DFT) and polyphase filtering, to digitally stack and band-limit the individual carriers respectively. The four-channel OFDM approach enables the implementation of a modem that can be both power and bandwidth efficient, with sufficient parallelism to meet higher data rate goals. As a result, the OFDM modem requires only a 240 MHz bandwidth to transmit 622 Mbps. Hardware and simulation results in the form of spectrum diagrams and bit-error-rate (BER) curves are also presented in this paper.

Kifle, Muli↗

Inference of response functions with the help of machine-learning algorithms

Response functions are a key quantity to describe the near-equilibrium dynamics of strongly interacting many-body systems. Recent techniques that attempt to overcome the challenges of calculating these ab initio have employed expansions in terms of orthogonal polynomials. We employ a neural network prediction algorithm to reconstruct a response function 𝑆⁡(𝜔) defined over a range in frequencies 𝜔. Here, we represent the calculated response function as a truncated Chebyshev series whose coefficients can be optimized to reduce the representation error. We compare the quality of response functions obtained using coefficients calculated using a neural network (NN) algorithm with those computed using the Gaussian integral transform (GIT) method. In the regime where only a small number of terms in the Chebyshev series are retained, we find that the NN scheme outperforms the GIT method.

Kurkcuoglu, Doga Murat [Fermi National Accelerator↗

Imine Reductase-Catalyzed, Radical-Mediated Asymmetric Cyano Group Migration

Functional group migration (FGM) reactions represent a fundamental class of transformations in organic chemistry, enabling the repositioning of functional moieties in nonobvious ways. However, catalytic asymmetric radical-mediated FGMs remain rare due to the inherent challenges of achieving catalyst-controlled enantioselectivity over free radical intermediates. Herein, we repurpose imine reductases (IREDs), a class of biotechnologically important enzymes known for their substrate promiscuity, to enable the first examples of catalytic asymmetric cyano group migration via a radical mechanism. An orthogonal set of radical enzymes, including PbaIREDCym and SmiIREDCym, was engineered, allowing both 1,4- and 1,5-cyano group migration reactions to occur in an enantiodivergent fashion. The use of the nonionic surfactant TPGS-1000 was found to improve both the yield and enantioselectivity of these cyano migration reactions. Furthermore, this biocatalytic process exhibited a broad substrate scope and is readily scalable, affording a rare example of chiral nonamine product assembly with imine reductases. More broadly, stereoselective radical biocatalysis with engineered IREDs and other versatile enzymes provides a potentially general solution to challenging asymmetric FGM reactions.

Biocatalysis↗

Spatiotemporally resolved measurements of electric field around a piezoelectric transformer using electric-field induced second harmonic (E-FISH) generation

When a piezoelectric transformer (PT) is actuated at its second harmonic frequency by a low input voltage, the generated electric field at the distal end can be sufficient to breakdown the surrounding gas, making them attractive power sources for non-equilibrium plasma generation. Understanding the potential and electric fields produced in the surrounding medium by the PT is important for effectively designing and using PT plasma devices. In this work, the spatiotemporally resolved characteristics of the electric field generated by a PT operating in open air have been investigated using the femtosecond electric field-induced second harmonic generation (E-FISH) method. Here, electric field components were determined by simultaneously conducting E-FISH measurements with the incident laser polarized in two orthogonal directions relative to the PT crystal. Results of this work demonstrate the spatial distribution of electric field around the PT’s output distal end and how it evolves as a function of time. Notably, the strongest electric field appears on the face of the PT’s distal surface, near the top and bottom edges and decreases by approximately 70% over 3 mm. The time delay between the PT’s input voltage and measured electric field indicates that there is an about 0.45 phase difference between the PT’s input voltage and output signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Shifted Block Lanczos Algorithm 1: The Block Recurrence

In this paper we describe a block Lanczos algorithm that is used as the key building block of a software package for the extraction of eigenvalues and eigenvectors of large sparse symmetric generalized eigenproblems. The software package comprises: a version of the block Lanczos algorithm specialized for spectrally transformed eigenproblems; an adaptive strategy for choosing shifts, and efficient codes for factoring large sparse symmetric indefinite matrices. This paper describes the algorithmic details of our block Lanczos recurrence. This uses a novel combination of block generalizations of several features that have only been investigated independently in the past. In particular new forms of partial reorthogonalization, selective reorthogonalization and local reorthogonalization are used, as is a new algorithm for obtaining the M-orthogonal factorization of a matrix. The heuristic shifting strategy, the integration with sparse linear equation solvers and numerical experience with the code are described in a companion paper.

Grimes, Roger G.↗

Formation of Vertically Oriented Channels during Calcination of Surfactant-Templated Titania-Doped Mesoporous Silica Thin Films

Here, the mesoporous structure development during calcination of titania-doped mesoporous silica thin films is characterized by in situ grazing-incidence small angle x-ray scattering (GISAXS). Varying amounts of titania (≤ 6 mol%) are incorporated in cetyltrimethylammonium bromide (CTAB) templated film by two methods in this study: dispersion of titania in the silica matrix, and fixing titania at the pore surface by complexation of the titania precursor with a sugar-based co-surfactant, dodecyl maltoside (C 12 G 2 ). Before calcination, films templated with only CTAB display a mesoporous structure with a mixture of in-plane 2D hexagonal and randomly oriented porous structure, while C 12 G 2 -complexed film possessed randomly oriented mesostructure only. Both methods yield a final structure of hexagonal close-packed pores orthogonally oriented to the substrate after calcination at 500 °C. This is attributed to weak Si-O-Ti bonds which allow sintering to occur at 500 °C, combined with unidirectional thermal contraction of the film in the vertical direction. Anisotropic stress and annealing of the films allows the randomly oriented pores to merge vertically and form channels with pore diameter of 2.3 ± 0.3 nm. Titania doping greater than 1% is required for this transformation in films with no C 12 G 2 surfactant used while the transformation was observed for films with C 12 G 2 complexation even at low (0.005%) titania doping, perhaps accentuated by orientational disorder introduced by the sugar surfactant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A new approach to the solution of large, full matrix equations: A two-dimensional potential flow feasibility study

An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.

James, R. M.↗

9+ Years of CALIOP PSC Data: An Evolving Climatology

Polar stratospheric clouds (PSCs) play key roles in the springtime chemical depletion of ozone at high latitudes. PSC particles provide sites for heterogeneous chemical reactions that transform stable chlorine and bromine reservoir species into highly reactive ozone-destructive forms. Furthermore, large nitric acid trihydrate (NAT) PSC particles can irreversibly redistribute odd nitrogen through gravitational sedimentation, which prolongs the ozone depletion process by slowing the reformation of the stable chlorine reservoirs. However, there are still significant gaps in our understanding of PSC processes, particularly concerning the details of NAT particle formation. Spaceborne observations from the CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) lidar on the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) satellite are providing a rich new dataset for studying PSCs on unprecedented vortex-wide scales. In this paper, we examine the vertical and spatial distribution of PSCs in the Antarctic and Arctic on vortex-wide scales for entire PSC seasons over the more than nine-year data record.

Pitts, Michael C.↗

A Class of Sparse Johnson–Lindenstrauss Transforms and Analysis of their Extreme Singular Values

The Johnson–Lindenstrauss (JL) lemma is a powerful tool for dimensionality reduction in modern algorithm design. The lemma states that any set of high-dimensional points in a Euclidean space can be projected into lower dimensions while approximately preserving pairwise Euclidean distances. Random matrices satisfying this lemma are called JL transforms (JLTs). Inspired by existing $s$-hashing JLTs with exactly $s$ nonzero elements on each column, the present work introduces an ensemble of sparse matrices encompassing so-called $s$-hashing-like matrices whose expected number of nonzero elements on each column is $s$. The independence of the sub-Gaussian entries of these matrices and the knowledge of their exact distribution play an important role in their analyses. Using properties of independent sub-Gaussian random variables, these matrices are demonstrated to be JLTs, and their smallest nontrivial singular values and largest singular values are estimated nonasymptotically using a technique from geometric functional analysis. As the dimensions of the matrix grow to infinity, these singular values are proved to converge almost surely to fixed quantities (by using the universal Bai–Yin law) and in distribution to the Gaussian orthogonal ensemble Tracy–Widom law after proper rescalings. Understanding the behaviors of extreme singular values is important in general because they are often used to define a measure of stability of matrix algorithms. For example, JLTs were recently used in derivative-free optimization algorithmic frameworks to select random subspaces in which are constructed random models or poll directions to achieve scalability, and hence estimating their smallest singular value in particular helps determine the dimension of these subspaces.

97 MATHEMATICS AND COMPUTING↗

Real space iterative reconstruction for vector tomography (RESIRE-V)

Tomography has had an important impact on the physical, biological, and medical sciences. To date, most tomographic applications have been focused on 3D scalar reconstructions. However, in some crucial applications, vector tomography is required to reconstruct 3D vector fields such as the electric and magnetic fields. Over the years, several vector tomography methods have been developed. Here, we present the mathematical foundation and algorithmic implementation of REal Space Iterative REconstruction for Vector tomography, termed RESIRE-V. RESIRE-V uses multiple tilt series of projections and iterates between the projections and a 3D reconstruction. Each iteration consists of a forward step using the Radon transform and a backward step using its transpose, then updates the object via gradient descent. Incorporating with a 3D support constraint, the algorithm iteratively minimizes an error metric, defined as the difference between the measured and calculated projections. The algorithm can also be used to refine the tilt angles and further improve the 3D reconstruction. To validate RESIRE-V, we first apply it to a simulated data set of the 3D magnetization vector field, consisting of two orthogonal tilt series, each with a missing wedge. Our quantitative analysis shows that the three components of the reconstructed magnetization vector field agree well with the ground-truth counterparts. We then use RESIRE-V to reconstruct the 3D magnetization vector field of a ferromagnetic meta-lattice consisting of three tilt series. Our 3D vector reconstruction reveals the existence of topological magnetic defects with positive and negative charges. We expect that RESIRE-V can be incorporated into different imaging modalities as a general vector tomography method. To make the algorithm accessible to a broad user community, we have made our RESIRE-V MATLAB source codes and the data freely available at https://github.com/minhpham0309/RESIRE-V.

47 OTHER INSTRUMENTATION↗

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗

An efficient virus–induced gene silencing (VIGS) system for gene functional studies in Miscanthus

Virus-induced gene silencing (VIGS) is a powerful tool for transient gene functional analysis in plants, especially for monocot species (e.g., grasses) that are recalcitrant to transformation. Despite various VIGS systems that have been developed in different plant species, none was previously available for the bioenergy crop Miscanthus. Here, we report the establishment of an efficient and robust VIGS system mediated by Tobacco Rattle Virus (TRV) in Miscanthus. We first investigated the impact of various factors that may affect gene silencing efficiency using the Miscanthus sinensis Phytoene Desaturase (MsPDS) gene as a visual indicator of photobleaching. Then, we optimized the TRV-elicited VIGS procedure using an orthogonal experimental design with four factors (sprout size, Agrobacterium concentration, vacuum infiltration time, and co-incubation time) each at three levels. The following led to the highest silencing efficiency (~76%): inoculation of germinating seedlings (1.0–2.0 mm), Agrobacterium tumefaciens culture grown to optical density at 600 nm (OD 600 ) of 0.4, vacuum infiltration for 90 min, and co-incubation for 5 h. The VIGS system established was applicable for both M. sinensis and M. lutarioriparius, with comparable gene silencing efficiency. We verified the efficacy of the VIGS system via the functional characterization of the role of a MYB transcription factor, MsMYB112, in salt stress tolerance. Expression of MsMYB112 was successfully knocked down using the VIGS system, and this led to compromised salt tolerance in the silenced Miscanthus plants. The TRV-based VIGS system established may, therefore, substantially facilitate functional genomic studies in Miscanthus.

60 APPLIED LIFE SCIENCES↗

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗

Modular flow in JT gravity and entanglement wedge reconstruction

It has been shown in recent works that JT gravity with matter with two boundaries has a type II ∞ algebra on each side. As the bulk spacetime between the two boundaries fluctuates in quantum nature, we can only define the entanglement wedge for each side in a pure algebraic sense. As we take the semiclassical limit, we will have a fixed long wormhole spacetime for a generic partially entangled thermal state (PETS), which is prepared by inserting heavy operators on the Euclidean path integral. Under this limit, with appropriate assumptions of the matter theory, geometric notions of the causal wedge and entanglement wedge emerge in this background. In particular, the causal wedge is manifestly nested in the entanglement wedge. Different PETS are orthogonal to each other, and thus the Hilbert space has a direct sum structure over sub-Hilbert spaces labeled by different Euclidean geometries. The full algebra for both sides is decomposed accordingly. From the algebra viewpoint, the causal wedge is dual to an emergent type III 1 subalgebra, which is generated by boundary light operators. To reconstruct the entanglement wedge, we consider the modular flow in a generic PETS for each boundary. We show that the modular flow acts locally and is the boost transformation around the global RT surface in the semiclassical limit. It follows that we can extend the causal wedge algebra to a larger type III 1 algebra corresponding to the entanglement wedge. Within each sub-Hilbert space, the original type II ∞ reduces to type III 1 .

2D Gravity↗