Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear correlation”

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 919 records · Page 51

Isotope effects on energy transport in the core of ASDEX-Upgrade tokamak plasmas: Turbulence measurements and model validation

Design and operation of future tokamak fusion reactors using a deuterium–tritium 50:50 mix requires a solid understanding of how energy confinement properties change with ion mass. This study looks at how turbulence and energy transport change in L-mode plasmas in the ASDEX Upgrade tokamak when changing ion species between hydrogen and deuterium. For this purpose, both experimental turbulence measurements and modeling are employed. Local measurements of ion-scale (with wavevector of fluctuations perpendicular to the B-field k⊥< 2 cm−1, k⊥ρs< 0.2, where ρs is the ion sound Larmor radius using the deuterium ion mass) electron temperature fluctuations have been performed in the outer core (normalized toroidal flux ρTor=0.65−0.8) using a multi-channel correlation electron cyclotron emission diagnostic. Lower root mean square perpendicular fluctuation amplitudes and radial correlation lengths have been measured in hydrogen vs deuterium. Measurements of the cross-phase angle between a normal-incidence reflectometer and an ECE signal were made to infer the cross-phase angle between density and temperature fluctuations. The magnitude of the cross-phase angle was found larger (more out-of-phase) in hydrogen than in deuterium. TRANSP power balance simulations show a larger ion heat flux in hydrogen where the electron-ion heat exchange term is found to play an important role. These experimental observations were used as the basis of a validation study of both quasilinear gyrofluid trapped gyro-Landau fluid-SAT2 and nonlinear gyrokinetic GENE codes. Linear solvers indicate that, at long wavelengths (k⊥ρs<1), energy transport in the deuterium discharge is dominated by a mixed ion-temperature-gradient (ITG) and trapped-electron mode turbulence while in hydrogen transport is exclusively and more strongly driven by ITG turbulence. The Ricci validation metric has been used to quantify the agreement between experiments and simulations taking into account both experimental and simulation uncertainties as well as four different observables across different levels of the primacy hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assimilation of satellite microwave observations over the rainbands of tropical cyclones

A novel Bayesian Monte Carlo integration (BMCI) technique was developed to retrieve geophysical variables from satellite microwave radiometer data in the presence of tropical cyclones. The BMCI technique includes three steps: generating a stochastic database, simulating satellite brightness temperatures using a radiative transfer model, and retrieving geophysical variables such as profiles of temperature, relative humidity, and cloud liquid and ice water content from real observations. The technique also provides uncertainty estimates for each retrieval and can output the error covariance matrix of selected parameters. The measurements from the Advanced Technology Microwave Sounder (ATMS) on board Suomi National Polar-Orbiting Partnership (Suomi NPP) and the Global Precipitation Measurement (GPM) Microwave Imager (GMI) were used as input. A new technique was developed to correct the ATMS and GMI observations for the beam-filling effect, which is due to small-scale variability of precipitation and clouds when compared with the instrument footprint and also the nonlinear relation between the brightness temperature and precipitation. In addition, the assimilation of the BMCI retrievals into the NASA GEOS model is discussed for Hurricane Maria. The results show that assimilating the BMCI retrievals can influence the dynamical features of the cyclone, including a stronger warm core, a symmetric eye, and vertically aligned wind columns. Two possible factors that may limit the impact of the BMCI retrievals include 1) the resolution of the model (about 25 km), which was too coarse to show the potential of the BMCI data in improving the representation of tropical storms in the model forecast, and 2) the data assimilation system not being able to consider vertically correlated observation errors.

Isaac Moradi↗

Mesoscopic-scale complexity in macroscopically-uniform plastic flow of an Al 0.3 CoCrFeNi high-entropy alloy

The concept of a smooth and homogeneous plastic flow of solids is nowadays constantly challenged by various observations of the self-organization of crystal defects on mesoscopic scales pertaining, e.g., to acoustic emission or the evolution of the local strain field. Such investigations would be of particular interest for High-Entropy Alloys (HEAs) characterized by extremely complex microstructures. However, the complexity of their deformation has been only studied in relation to the jerky flow. The present work explores a mesoscale complexity of the macroscopically-smooth plastic flow of an Al 0.3 CoCrFeNi alloy. Here, an inclusive approach involving various experimental methods and complementary nonlinear analyses allowed to reveal transitions between distinct dynamical regimes, which may be assimilated with noises of different colors, in particular the blue noise that is rarely observed in complex systems, and various kinds of reddened behavior. Such findings provide significant new insights into the micro/macro transition in the deformation behavior of HEAs.

36 MATERIALS SCIENCE↗

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗

Two-dimensional electronic–vibrational spectroscopy: Exploring the interplay of electrons and nuclei in excited state molecular dynamics

Two-dimensional electronic-vibrational spectroscopy (2DEVS) is an emerging spectroscopic technique which exploits two different frequency ranges for the excitation (visible) and detection (infrared) axes of a 2D spectrum. In contrast to degenerate 2D techniques, such as 2D electronic or 2D infrared spectroscopy, the spectral features of a 2DEV spectrum report cross correlations between fluctuating electronic and vibrational energy gaps rather than autocorrelations as in the degenerate spectroscopies. The center line slope of the spectral features reports on this cross correlation function directly and can reveal specific electronic-vibrational couplings and rapid changes in the electronic structure, for example. The involvement of the two types of transition moments, visible and infrared, makes 2DEVS very sensitive to electronic and vibronic mixing. 2DEV spectra also feature improved spectral resolution, making the method valuable for unraveling the highly congested spectra of molecular complexes. The unique features of 2DEVS are illustrated in this paper with specific examples and their origin described at an intuitive level with references to formal derivations provided. Although early in its development and far from fully explored, 2DEVS has already proven to be a valuable addition to the tool box of ultrafast nonlinear optical spectroscopy and is of promising potential in future efforts to explore the intricate connection between electronic and vibrational nuclear degrees of freedom in energy and charge transport applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast laser-driven many-body dynamics and Kondo coherence collapse

Ultrafast laser pulse has provided a systematic way to inspect the dynamics of electrons in condensed matter systems. In this paper, by means of time-dependent density matrix renormalization group, we study an ultrafast laser-driven Kondo lattice model, in which conduction electrons are strongly coupled with magnetically local moments. The single-particle spectral function due to strong correlation effects and photon emission in the nonequilibrium states under laser driving are calculated. We find laser field excited collective doublon-hole pairs and an associated transient melting of Kondo coherence phase, signifying the collapse of Kondo energy gap. Moreover, we show that the photon emission, induced by a strong laser field, exhibits a different intensity characteristics than in the equilibrium Kondo insulator, which could be explained by the Kondo collapse and related suppression of both intraband and interband contribution in Kondo melting liquid. This theoretical insight is accessible with time- and angle-resolved photoemission spectroscopy and high-harmonic generation spectroscopy and will stimulate the investigation of nonequilibrium dynamics and nonlinear phenomenon in heavy fermion systems.

36 MATERIALS SCIENCE↗

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↗

Mechanical characterization and modeling of non-linear deformation and fracture of a fiber reinforced metal matrix composite

The nonlinear anisotropic mechanical behavior of an aluminum alloy metal matrix composite reinforced with continuous alumina fibers was determined experimentally. The mechanical behavior of the composite were modeled by assuming that the composite has a periodical microstructure. The resulting unit cell problem was solved with the finite element method. Excellent agreement was found between theoretically predicted and measured stress-strain responses for various tensile and shear loadings. The stress-strain responses for transverse and inplane shear were found to be identical and this will provide a simplification of the constitutive equations for the composite. The composite has a very low ductility in transverse tension and a limited ductility in transverse shear that was correlated to high hydrostatic stresses that develop in the matrix. The shape of the initial yield surface was calculated and good agreement was found between the calculated shape and the experimentally determined shape.

Jansson, S.↗

Parametric Analysis of a Hover Test Vehicle using Advanced Test Generation and Data Analysis

Large complex aerospace systems are generally validated in regions local to anticipated operating points rather than through characterization of the entire feasible operational envelope of the system. This is due to the large parameter space, and complex, highly coupled nonlinear nature of the different systems that contribute to the performance of the aerospace system. We have addressed the factors deterring such an analysis by applying a combination of technologies to the area of flight envelop assessment. We utilize n-factor (2,3) combinatorial parameter variations to limit the number of cases, but still explore important interactions in the parameter space in a systematic fashion. The data generated is automatically analyzed through a combination of unsupervised learning using a Bayesian multivariate clustering technique (AutoBayes) and supervised learning of critical parameter ranges using the machine-learning tool TAR3, a treatment learner. Covariance analysis with scatter plots and likelihood contours are used to visualize correlations between simulation parameters and simulation results, a task that requires tool support, especially for large and complex models. We present results of simulation experiments for a cold-gas-powered hover test vehicle.

Gundy-Burlet, Karen↗

Universal Decoder for PPM of any Order

A recently developed algorithm for demodulation and decoding of a pulse-position- modulation (PPM) signal is suitable as a basis for designing a single hardware decoding apparatus to be capable of handling any PPM order. Hence, this algorithm offers advantages of greater flexibility and lower cost, in comparison with prior such algorithms, which necessitate the use of a distinct hardware implementation for each PPM order. In addition, in comparison with the prior algorithms, the present algorithm entails less complexity in decoding at large orders. An unavoidably lengthy presentation of background information, including definitions of terms, is prerequisite to a meaningful summary of this development. As an aid to understanding, the figure illustrates the relevant processes of coding, modulation, propagation, demodulation, and decoding. An M-ary PPM signal has M time slots per symbol period. A pulse (signifying 1) is transmitted during one of the time slots; no pulse (signifying 0) is transmitted during the other time slots. The information intended to be conveyed from the transmitting end to the receiving end of a radio or optical communication channel is a K-bit vector u. This vector is encoded by an (N,K) binary error-correcting code, producing an N-bit vector a. In turn, the vector a is subdivided into blocks of m = log2(M) bits and each such block is mapped to an M-ary PPM symbol. The resultant coding/modulation scheme can be regarded as equivalent to a nonlinear binary code. The binary vector of PPM symbols, x is transmitted over a Poisson channel, such that there is obtained, at the receiver, a Poisson-distributed photon count characterized by a mean background count nb during no-pulse time slots and a mean signal-plus-background count of ns+nb during a pulse time slot. In the receiver, demodulation of the signal is effected in an iterative soft decoding process that involves consideration of relationships among photon counts and conditional likelihoods of m-bit vectors of coded bits. Inasmuch as the likelihoods of all the m-bit vectors of coded bits mapping to the same PPM symbol are correlated, the best performance is obtained when the joint mbit conditional likelihoods are utilized. Unfortunately, the complexity of decoding, measured in the number of operations per bit, grows exponentially with m, and can thus become prohibitively expensive for large PPM orders. For a system required to handle multiple PPM orders, the cost is even higher because it is necessary to have separate decoding hardware for each order. This concludes the prerequisite background information. In the present algorithm, the decoding process as described above is modified by, among other things, introduction of an lbit marginalizer sub-algorithm. The term "l-bit marginalizer" signifies that instead of m-bit conditional likelihoods, the decoder computes l-bit conditional likelihoods, where l is fixed. Fixing l, regardless of the value of m, makes it possible to use a single hardware implementation for any PPM order. One could minimize the decoding complexity and obtain an especially simple design by fixing l at 1, but this would entail some loss of performance. An intermediate solution is to fix l at some value, greater than 1, that may be less than or greater than m. This solution makes it possible to obtain the desired flexibility to handle any PPM order while compromising between complexity and loss of performance.

Moision, Bruce E.↗

Wall-interference assessment and corrections for transonic NACA 0012 airfoil data from various wind tunnels

A nonlinear, four wall, post-test wall interference assessment/correction (WIAC) code was developed for transonic airfoil data from solid wall wind tunnels with flexibly adaptable top and bottom walls. The WIAC code was applied over a broad range of test conditions to four sets of NACA 0012 airfoil data, from two different adaptive wall wind tunnels. The data include many test points for fully adapted walls, as well as numerous partially adapted and unadapted test points, which together represent many different model/tunnel configurations and possible wall interference effects. Small corrections to the measured Mach numbers and angles of attack were obtained from the WIAC code even for fully adapted data; these corrections generally improve the correlation among the various sets of airfoil data and simultaneously improve the correlation of the data with calculations for a 2-D, free air, Navier-Stokes code. The WIAC corrections for airfoil data taken in fully adapted wall test sections are shown to be significantly smaller than those for comparable airfoil data from straight, slotted wall test sections. This indicates, as expected, a lesser degree of wall interference in the adapted wall tunnels relative to the slotted wall tunnels. Application of the WIAC code to this data was, however, somewhat more difficult and time consuming than initially expected from similar previous experience with WIAC applications to slotted wall data.

Green, Lawrence L.↗

Gamma-Ray Observations of the Orion Molecular Clouds with the Fermi Large Area Telescope

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between approx 100 MeV and approx 100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to approx 10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity-integrated CO intensity (W(sub CO)) at a 1 deg 1 deg pixel level. The correlation is found to be linear over a W(sub CO) range of approx 10-fold when divided in three regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The W(sub CO)-to-mass conversion factor, X(sub CO), is found to be approx 2.3 10(exp 20) / sq cm (K km/s)(exp -1) for the high-longitude part of Orion A (l > 212 deg), approx 1.7 times higher than approx 1.3 10(exp 20) found for the rest of Orion A and B. We interpret the apparent high X(sub CO) in the high-longitude region of Orion A in the light of recent works proposing a nonlinear relation between H2 and CO densities in the diffuse molecular gas.W(sub CO) decreases faster than the H2 column density in the region making the gas "darker" to W(sub CO).

Ackermann, M.↗

Structural Shape Sensing from Strain Data Using the Basis Function Method

An aircraft shape must be measured during flight to implement an active trim shape control, such as an active dihedral control for a highly flexible aircraft with a high-aspect-ratio wing, and an active sonic boom control for a supersonic transport aircraft. A basis function method for linear as well as geometrically nonlinear structural shape sensing is proposed in this study. Basis functions can be mode shapes, linear static deformation shapes, or geometrically nonlinear static deformation shapes. The proposed basis function method is validated using a high-aspect-ratio wing, a swept test plate, and the National Aeronautics and Space Administration Low Boom Flight Demonstration mission X-59 Quiet Supersonic Technology (QueSST) aircraft (Lockheed Martin, Bethesda, Maryland). A large structural deformation of a high-aspect-ratio wing is successfully captured using the proposed basis function method with less than 0.3-percentprediction error at the wing-tip section. The basis function method gives excellent deformation prediction, even with stress concentration. Performance of the basis functions is compared using the X-59 stabilator sample case. In most of the sample load cases, static deformations give a better correlation with the target deformation than do mode shapes. Wing shape sensing with sparse strain data is also demonstrated in this study using the X-59 QueSST aircraft. Predicted structural deformations match reasonably with the target deformations even without strain gauges on some structural components. The predicted deformations have a good match with the target deformations.

Chan-gi Pak↗

Field Induced Density Wave in a Kagome Superconductor

On the kagome lattice, electrons benefit from the simultaneous presence of band topology, flat electronic bands, and Van Hove singularities, forming competing or cooperating orders. Understanding the interrelation between these distinct order parameters remains a significant challenge, leaving much of the associated physics unexplored. In the kagome superconductor KV 3 Sb 5 , which exhibits a charge density wave (CDW) state below T ≃ 78 K, we uncover an unpredicted field-induced phase transition below 6 K. The observed transition is marked by a hysteretic anomaly in the resistivity, nonlinear electrical transport, and a change in the symmetry of the electronic response as probed via the angular dependence of the magnetoresistivity. These observations surprisingly suggest the emergence of an unanticipated broken symmetry state coexisting with the original CDW. To understand this experimental observation, we developed a theoretical minimal model for the normal state inside the high-temperature parent CDW phase, where an incommensurate CDW order emerges as an instability subleading to superconductivity. Here, the incommensurate CDW emerges when superconducting fluctuations become fully suppressed by large magnetic fields. Our results suggest that, in kagome superconductors, quantum states can either coexist or are nearly degenerate in energy, indicating that these are rich platforms to expose new correlated phenomena.

Charge density waves↗

Optimal Initial Perturbations for Ensemble Prediction of the Madden-Julian Oscillation during Boreal Winter

An initialization strategy, tailored to the prediction of the Madden-Julian oscillation (MJO), is evaluated using the Goddard Earth Observing System Model, version 5 (GEOS-5), coupled general circulation model (CGCM). The approach is based on the empirical singular vectors (ESVs) of a reduced-space statistically determined linear approximation of the full nonlinear CGCM. The initial ESV, extracted using 10 years (1990-99) of boreal winter hindcast data, has zonal wind anomalies over the western Indian Ocean, while the final ESV (at a forecast lead time of 10 days) reflects a propagation of the zonal wind anomalies to the east over the Maritime Continent an evolution that is characteristic of the MJO. A new set of ensemble hindcasts are produced for the boreal winter season from 1990 to 1999 in which the leading ESV provides the initial perturbations. The results are compared with those from a set of control hindcasts generated using random perturbations. It is shown that the ESV-based predictions have a systematically higher bivariate correlation skill in predicting the MJO compared to those using the random perturbations. Furthermore, the improvement in the skill depends on the phase of the MJO. The ESV is particularly effective in increasing the forecast skill during those phases of the MJO in which the control has low skill (with correlations increasing by as much as 0.2 at 20 25-day lead times), as well as during those times in which the MJO is weak.

Ham, Yoo-Geun↗

Parameter identification of material constants in a composite shell structure

One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Martinez, David R.↗

Progressive Failure Analysis Methodology for Laminated Composite Structures

A progressive failure analysis method has been developed for predicting the failure of laminated composite structures under geometrically nonlinear deformations. The progressive failure analysis uses C(exp 1) shell elements based on classical lamination theory to calculate the in-plane stresses. Several failure criteria, including the maximum strain criterion, Hashin's criterion, and Christensen's criterion, are used to predict the failure mechanisms and several options are available to degrade the material properties after failures. The progressive failure analysis method is implemented in the COMET finite element analysis code and can predict the damage and response of laminated composite structures from initial loading to final failure. The different failure criteria and material degradation methods are compared and assessed by performing analyses of several laminated composite structures. Results from the progressive failure method indicate good correlation with the existing test data except in structural applications where interlaminar stresses are important which may cause failure mechanisms such as debonding or delaminations.

Sleight, David W.↗

Nonlinear interactions of electromagnetic waves with a thermally generated cesium plasma.

The interaction of a large‐amplitude disturbing wave with a thermally generated, highly ionized cesium plasma has been investigated using a small‐amplitude sensing wave as a probe. The subsequent interaction between the disturbing and sensing waves is related to the plasma electron temperature. The plasma was contained in a magnetic field of 500–1000 G and had an electron plasma frequency near 10 GHz and an electron gyromagnetic frequency around 2 GHz. The sensing wave was a 10 GHz plane wave of approximately 1.5 mW input power and the disturbing wave was 7.5 GHz plane wave of about 1.2 W peak power. Measurements of the sensing‐wave transmitted power were made using standard microwave techniques. A notch filter in the sensing‐wave receiving circuit reduced the scattered disturbing‐wave interference to an immeasurably low level. Langmuir probes were used to detect disturbing‐wave modulation and to measure plasma temperatures and densities. Measurements of the interaction were correlated with a solution of the Boltzmann transport equation for a perturbed Maxwellian distribution function with binary Coulomb collisions.

Jack G Hewitt, Jr↗