Search NASA⌕ Search

SEARCH · Search NASA

Results for “matching complex”

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 217 records · Page 12

Pressure-driven density match nucleates metastable r8 phases from amorphous Si and Ge

The pressure–temperature phase behavior of covalent disordered solids such as amorphous silicon and germanium is complex. Questions remain on possible glass transitions, on polyamorphism via amorphous–amorphous transitions, on connections with liquid–liquid transitions, on structure-behavior relationships, and on their potential as precursor for novel methods for material discovery. Here we demonstrate experimentally the nucleation of a metastable, four-fold coordinated rhombohedral r8 phase from pure amorphous silicon and germanium upon room temperature compression at pressures below 10 GPa. Accompanying theory reveals a strong pressure-driven distortion of the bond angle transforming the starting tetrahedral low-density amorphous network to a distorted four-fold coordinated medium-density state. This state is of lower density than metallic high-density networks, resembles the crystalline r8 phase and initiates its nucleation. Our finding shows that polyamorphism is not the only possible transformation mode for these amorphous solids and that instead nucleation of interesting functional phases at potentially useful pressures is possible. Such novel access modes to metastable structures are critical for future exploitability and could be useful for other tetrahedral materials including carbon, where the related (bc8) post-diamond phase remains elusive. Our observed density match between an amorphous and a metastable crystalline phase clearly allows for a new phase transition pathway, while corresponding theory demonstrates how carefully validated atomistic simulations can guide prediction, discovery and synthesis of novel material structures.

Materials discovery↗

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

V‐Groove Si Nanopatterning for the Direct Epitaxy of Orientation‐Patterned III–V Nonlinear Optical Crystals

Orientation‐patterned (OP) III–V semiconductors—used as quasi‐phase‐matched crystals for nonlinear optics applications—are typically epitaxially grown on expensive III–V substrates using a complex process requiring three separate epitaxy steps. In this work, a method is demonstrated for growing orientation‐patterned III–V crystals on Si substrates through V‐groove nanopatterning and a single epitaxial growth. V‐groove Si allows for suppression of random antiphase domain formation that is typical of III–V growth on (001)‐oriented Si through the use of (111)‐faceted trenches patterned on (001)‐oriented Si substrates. By alternating the directions of the V‐groove trenches between [110] and [10], antiphase boundaries can be selectively induced at the boundaries between the two directions of trenches due to the difference in symmetry between the III–V material and Si. This approach allows for a greatly simplified process for growing OP‐III‐Vs and more broadly is a new, robust approach for precisely patterning antiphase boundaries of arbitrary shape and length scales down to 50 nm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗

PDBspheres: a method for finding 3D similarities in local regions in proteins

Abstract We present a structure-based method for finding and evaluating structural similarities in protein regions relevant to ligand binding. PDBspheres comprises an exhaustive library of protein structure regions (‘spheres’) adjacent to complexed ligands derived from the Protein Data Bank (PDB), along with methods to find and evaluate structural matches between a protein of interest and spheres in the library. PDBspheres uses the LGA (Local–Global Alignment) structure alignment algorithm as the main engine for detecting structural similarities between the protein of interest and template spheres from the library, which currently contains >2 million spheres. To assess confidence in structural matches, an all-atom-based similarity metric takes side chain placement into account. Here, we describe the PDBspheres method, demonstrate its ability to detect and characterize binding sites in protein structures, show how PDBspheres—a strictly structure-based method—performs on a curated dataset of 2528 ligand-bound and ligand-free crystal structures, and use PDBspheres to cluster pockets and assess structural similarities among protein binding sites of 4876 structures in the ‘refined set’ of the PDBbind 2019 dataset.

59 BASIC BIOLOGICAL SCIENCES↗

General-Purpose Unsupervised Cyber Anomaly Detection via Non-Negative Tensor Factorization

Distinguishing malicious anomalous activities from unusual but benign activities is a fundamental challenge for cyber defenders. Prior studies have shown that statistical user behavior analysis yields accurate detections by learning behavior profiles from observed user activity. These unsupervised models are able to generalize to unseen types of attacks by detecting deviations from normal behavior, without knowledge of specific attack signatures. However, approaches proposed to date based on probabilistic matrix factorization are limited by the information conveyed in a two-dimensional space. Non-negative tensor factorization, on the other hand, is a powerful unsupervised machine learning method that naturally models multi-dimensional data, capturing complex and multi-faceted details of behavior profiles. Herein, our new unsupervised statistical anomaly detection methodology matches or surpasses state-of-the-art supervised learning baselines across several challenging and diverse cyber application areas, including detection of compromised user credentials, botnets, spam e-mails, and fraudulent credit card transactions.

97 MATHEMATICS AND COMPUTING↗

Galileo SSI Observations of Volcanic Activity at Tvashtar Catena, Io

Introduction: We report on the analysis of the Galileo SSI's observations of the volcanic activity at Tvashtar Catena, Io as discussed by Milazzo et al. Galileo's Solid State Imager (SSI) observed Tvashtar Catena (63 deg N, 120 deg W) four times between November 1999 and October 2001, providing a unique look at the distinctive high latitude volcanism on Io. The November 1999 observation spatially resolved, for the first time, an active extraterrestrial fissure eruption. The brightness temperature of the lavas at the November 1999 fissure eruption was 1300 K. The second observation (orbit I27, February 2000) showed a large (approx. 500 sq km) region with many, small spots of hot, active lava. The third observation was taken in conjunction with a Cassini observation in December 2000 and showed a Pele-like plume deposition ring, while the Cassini images revealed a 400 km high Pele-type plume above the Catena. The final Galileo SSI observation of Tvashtar was acquired in October 2001, and all obvious (to SSI) activity had ceased, although data from Galileo's Near Infrared Mapping Spectrometer (NIMS) indicated that there was still significant thermal emission from the Tvashtar region. We have concentrated on analyzing the style of eruption during orbit I27 (February 2000). Comparison with a lava flow cooling model indicates that the behavior of the Tvashtar eruption during I27 does not match that of "simple" advancing lava flows. Instead, it may be an active lava lake or a complex set of lava flows with episodic, overlapping (in time and space) eruptions.

Milazzo, M. P.↗

Digital Radar-Signal Processors Implemented in FPGAs

High-performance digital electronic circuits for onboard processing of return signals in an airborne precipitation- measuring radar system have been implemented in commercially available field-programmable gate arrays (FPGAs). Previously, it was standard practice to downlink the radar-return data to a ground station for postprocessing a costly practice that prevents the nearly-real-time use of the data for automated targeting. In principle, the onboard processing could be performed by a system of about 20 personal- computer-type microprocessors; relative to such a system, the present FPGA-based processor is much smaller and consumes much less power. Alternatively, the onboard processing could be performed by an application-specific integrated circuit (ASIC), but in comparison with an ASIC implementation, the present FPGA implementation offers the advantages of (1) greater flexibility for research applications like the present one and (2) lower cost in the small production volumes typical of research applications. The generation and processing of signals in the airborne precipitation measuring radar system in question involves the following especially notable steps: The system utilizes a total of four channels two carrier frequencies and two polarizations at each frequency. The system uses pulse compression: that is, the transmitted pulse is spread out in time and the received echo of the pulse is processed with a matched filter to despread it. The return signal is band-limited and digitally demodulated to a complex baseband signal that, for each pulse, comprises a large number of samples. Each complex pair of samples (denoted a range gate in radar terminology) is associated with a numerical index that corresponds to a specific time offset from the beginning of the radar pulse, so that each such pair represents the energy reflected from a specific range. This energy and the average echo power are computed. The phase of each range bin is compared to the previous echo by complex conjugate multiplication to obtain the mean Doppler shift (and hence the mean and variance of the velocity of precipitation) of the echo at that range.

Berkun, Andrew↗

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

boundary layer transition↗

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

Boundary Layer Transition↗

A Three-Dimensional Coupled Internal/External Simulation of a Film-Cooled Turbine Vane

A three-dimensional Navier-Stokes simulation has been performed for a realistic film-cooled turbine vane using the LeRC-HT code. The simulation includes the flow regions inside the coolant plena and film cooling holes in addition to the external flow. The vane is the subject of an upcoming NASA Glenn Research Center experiment and has both circular cross-section and shaped film cooling holes. This complex geometry is modeled using a multi-block grid which accurately discretizes the actual vane geometry including shaped holes. The simulation matches operating conditions for the planned experiment and assumes periodicity in the spanwise direction on the scale of one pitch of the film cooling hole pattern. Two computations were performed for different isothermal wall temperatures, allowing independent determination of heat transfer coefficients and film effectiveness values. The results indicate separate localized regions of high heat transfer coefficient values, while the shaped holes provide a reduction in heat flux through both parameters. Hole exit data indicate rather simple skewed profiles for the round holes, but complex profiles for the shaped holes with mass fluxes skewed strongly toward their leading edges.

Heidmann, James D.↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Biconjugate Impedance Matching Control of the HERO WEC: Preprint

Effective control strategies are essential for maximizing the power output of wave energy converters (WECs). In this study, we explore the methodology of port-based modelling which can be used for bi-conjugate impedance matching, a type of reactive control that involves maximizing power transfer by setting the load (PTO) impedance to be the complex conjugate of the source impedance (WEC). Despite its foundational role in electrical engineering, this method has seen limited application in wave energy and has not yet been demonstrated in wave-powered desalination systems. We represent our complex multidisciplinary system as an equivalent circuit to explore frequency-based relations. The HERO WEC, a small-scale, modular point absorber (PA) designed for reverse osmosis desalination is used as a case study to demonstrate port-based modelling for control applications. Improving the performance of wave-powered water desalination systems could play a critical role in addressing freshwater scarcity for coastal communities, especially those in remote or disaster-prone regions.

16 TIDAL AND WAVE POWER↗

Neutral–neutral synthesis of organic molecules in cometary comae

Remote and in situ observations of cometary gases have revealed the presence of a wealth of complex organic molecules, including carbon chains, alcohols, imines, and the amino acid glycine. Such chemical complexity in cometary material implies that impacts by comets could have supplied reagents for prebiotic chemistry to young planetary surfaces. However, the assumption that some of the molecules observed in cometary comae at millimetre wavelengths originate from ices stored inside the nucleus has not yet been proven. In fact, the comae of moderately-active comets reach sufficient densities within a few thousand kilometres of the nucleus for an active (solar radiation-driven) photochemistry to ensue. Here, we present results from our latest chemical-hydrodynamic models incorporating an updated reaction network, and show that the commonly-observed HC3N (cyanoacetylene) and NH2CHO (formamide) molecules can be efficiently produced in cometary comae as a result of two-body, neutral–neutral, gas-phase reactions involving well-known coma species. In the presence of a near-nucleus distributed source of CN (similar to that observed by the Rosetta spacecraft at comet 67P), we find that sufficient HC3N and NH2CHO can be synthesized to match the abundances of these molecules observed previously in Oort cloud comets. The precise coma origins of these (and other) complex organic molecules can be verified through radio interferometric mapping observations, for example using the Atacama Large Millimeter/submillimeter Array.

Comets↗

Inference in fuzzy rule bases with conflicting evidence

Inference based on fuzzy 'If ... then' rules has played a very important role since when Zadeh proposed the Compositional Rule of Inference and, especially, since the first successful application presented by Mamdani. From the mid-1980's when the 'fuzzy boom' started in Japan, numerous industrial applications appeared, all using simplified techniques because of the high levels of computational complexity. Another feature is that antecedents in the rules are distributed densely in the input space, so the conclusion can be calculated by some weighted combination of the consequents of the matching (fired) rules. The CRI works in the following way: If R is a rule and A* is an observation, the conclusion is computed by B* = R o A* (o stands for the max-min composition). Algorithms implementing this idea directly have an exponential time complexity (maybe the problem is NP-hard) as the rules are relations in X x Y, a k1 x k2 dimensional space, if X is k1, Y is k2 dimensional. The simplified techniques usually decompose the relation into k1 projections in X(sub i) and measure in some way the degree of similarity between observation and antecedent by some parameter of the overlapping. These parameters are aggregated to a single value in (0,1) which is applied as a resulting weight for the given rule. The projections of rules in dimensions Y(sub i) are weighted by these aggregated values and then they are combined in order to obtain a resulting conclusion separately in every dimension. This method is unapplicable with sparse bases as there is no guarantee that an arbitrary observation matches with any of the antecedents. Then, the degree of similarity is 0 and all consequents are weighted by 0. Some considerations for such a situation are summarized in the next sections.

Koczy, Laszlo T.↗

Mixed-precision iterative refinement using tensor cores on GPUs to accelerate solution of linear systems

Double-precision floating-point arithmetic (FP64) has been the de facto standard for engineering and scientific simulations for several decades. Problem complexity and the sheer volume of data coming from various instruments and sensors motivate researchers to mix and match various approaches to optimize compute resources, including different levels of floating-point precision. In recent years, machine learning has motivated hardware support for half-precision floating-point arithmetic. A primary challenge in high-performance computing is to leverage reduced-precision and mixed-precision hardware. We show how the FP16/FP32 Tensor Cores on NVIDIA GPUs can be exploited to accelerate the solution of linear systems of equations Ax = b without sacrificing numerical stability. The techniques we employ include multiprecision LU factorization, the preconditioned generalized minimal residual algorithm (GMRES), and scaling and auto-adaptive rounding to avoid overflow. We also show how to efficiently handle systems with multiple right-hand sides. On the NVIDIA Quadro GV100 (Volta) GPU, we achieve a 4×-5× performance increase and 5× better energy efficiency versus the standard FP64 implementation while maintaining an FP64 level of numerical stability.

GMRES↗

Duplicating Curved Tile Surfaces for Pull Testing

Strength of adhesive bonds to fragile objects with complex shapes tested easily with vacuum chucks. In bond-strength test, contact surface of chuck pressed against matching surface of bonded tile. Vacuum line extending through fine bore in chuck creates vacuum at chuck-and-tile interface. Thus, when chuck pulled away, atmospheric pressure tends to push tile away with chuck. If adhesive resists predetermined pull exerted by chuck on tile, bond acceptable. Originally developed for vacuum chucks used on silica tiles attached to outer surface of Space Shuttle orbiter, used to custom-fabricate chuck for specific tile.

Sebenick, R. H.↗

Applying transfer matrix method to the estimation of the modal characteristics of the NASA Mini-Mass Truss

It is beneficial to use a distributed parameter model for large space structures because the approach minimizes the number of model parameters. Holzer's transfer matrix method provides a useful means to simplify and standardize the procedure for solving the system of partial differential equations. Any large space structures can be broken down into sub-structures with simple elastic and dynamical properties. For each single element, such as beam, tether, or rigid body, we can derive the corresponding transfer matrix. Combining these elements' matrices enables the solution of the global system equations. The characteristics equation can then be formed by satisfying the appropriate boundary conditions. Then natural frequencies and mode shapes can be determined by searching the roots of the characteristic equation at frequencies within the range of interest. This paper applies this methodology, and the maximum likelihood estimation method, to refine the modal characteristics of the NASA Mini-Mast Truss by successively matching the theoretical response to the test data of the truss. The method is being applied to more complex configurations.

Shen, Ji-Yao↗