Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled inversion”

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 379 records · Page 21

Evaluating Terrain as a Turbulence Generation Method

When driving microscale large-eddy simulations with mesoscale model solutions, turbulence will take space to develop, known as fetch, on the microscale domain. To reduce fetch, it is common to add perturbations near the boundaries to speed up turbulence development. However, when simulating domains over complex terrain, it is possible that the terrain itself can quickly generate turbulence within the boundary layer. It is shown here that rugged terrain is able to generate turbulence without the assistance of a perturbation strategy; however, the levels of turbulence generated are improved when adding perturbations at the inlet. Flow over smoothed, but not flat, terrain fails to generate adequate turbulence throughout the boundary layer in all tests conducted herein. Sensitivities to the strength of the mean wind speed and boundary layer height are investigated and show that higher wind speeds produce turbulence over terrain features that slower wind speeds do not. Further, by increasing the height of the capping inversion, the effectiveness of topography alone to generate turbulence throughout the depth of the boundary is diminished. In all cases, the inclusion of a perturbation strategy improved simulation performance with respect to turbulence development.

17 WIND ENERGY↗

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interfacial Dzyaloshinskii-Moriya interaction of antiferromagnetic materials

The interface between a ferromagnet (FM) or antiferromagnet (AFM) and a heavy metal (HM) results in an antisymmetric exchange interaction known as the interfacial Dzyaloshinskii-Moriya interaction (iDMI) which favors non-collinear spin configurations. The iDMI is responsible for stabilizing noncollinear spin textures such as skyrmions in materials with bulk inversion symmetry. Interfacial DMI values have been previously determined theoretically and experimentally for FM/HM interfaces, and, in this work, they are calculated for the metallic AFM MnPt and the insulating AFM NiO. The heavy metals considered are W, Re, and Au. Values for the iDMI, exchange, and anisotropy constants are determined for different AFM and HM thicknesses. The iDMI values of the MnPt heterolayers are comparable to those of the common FM materials, and those of NiO are lower. In few-layer films of (001) MnPt, the high spin orbit coupling of the Pt layers can give rise to a small DMI in the absence of a HM layer.

36 MATERIALS SCIENCE↗

Multichannel Deconvolution of Vibrational Shock Signals: An Inverse Filtering Approach

When transporting critical systems of national security interest, out-of-the ordinary, impulsive events that can potentially be undetected and affect overall system performance are of great concern. Impulsive events that can occur are essentially pulse-like, transient signals of short duration that evolve from various phenomena. Here the event can be created by either the dropping of a test object, the system, subjecting it to a compact high-energy blow or being struck unintentionally during transit resulting in potential damage. The intensity and location of the strike can cause an inoperability condition that is unacceptable in a national security environment. Therefore, it is essential to detect, classify and localize damage of any test object subjected to an impulsive-event. This effort was targeted to evaluate the vibrational response of test objects that are subjected to “transport” shocks and roadway vibrations during shipping and handling. Any potential damage that could be inflicted during transportation must not only be detected, but also be evaluated to determine the operational readiness of a test object before and after transport. This event is a critical task that must be addressed as part of the Lawrence Livermore National Laboratory (LLNL) national security mission. The estimation of excitation signals from noisy data is termed the deconvolution problem in the signal processing literature. The deconvolution problem is based on recovering the input excitation signal from a system characterized by its impulse response sequence. Using this model of the system, an “inverse” representation or filter is developed to remove the system from the measured data and recover the input. Deconvolution techniques have existed for a long-time; however, transient deconvolution presents a few uncommon problems, since the signal has a finite-length time duration resulting in a limited amount of data containing information about the excitation process. The transient is wideband in the frequency domain relative to any measurement sensor implying that the smaller bandwidth sensor system “filters” the excitation eliminating some of its essential information for recovery. This fact, coupled with the filtering effect of the test object itself makes this excitation recovery (deconvolution) problem a challenge for signal processing.

42 ENGINEERING↗

Bayesian inference of fiber orientation and polymer properties in short fiber-reinforced polymer composites

Herein we present a Bayesian methodology to infer the elastic modulus of the constituent polymer and the fiber orientation state in a short-fiber reinforced polymer composite (SFRP). The properties are inversely determined using only a few experimental tests. Developing composite manufacturing digital twins for SFRP composite processes, including injection molding and extrusion deposition additive manufacturing (EDAM) requires extensive experimental material characterization. In particular, characterizing the composite mechanical properties is time consuming and therefore, micromechanics models are used to fully identify the elasticity tensor. Hence, the objective of this paper is to infer the fiber orientation and the effective polymer modulus and therefore, identify the elasticity tensor of the composite with minimal experimental tests. To that end, we develop a hierarchical Bayesian model coupled with a micromechanics model to infer the fiber orientation and the polymer elastic modulus simultaneously which we then use to estimate the composite elasticity tensor. We motivate and demonstrate the methodology for the EDAM process but the development is such that it is applicable to other SFRP composites processed via other methods. Our results demonstrate that the approach provides a reliable framework for the inference, with as few as three tensile tests, while accounting for epistemic and aleatory uncertainty. Posterior predictive checks show that the model is able to recreate the experimental data well. The ability of the Bayesian approach to calibrate the material properties and its associated uncertainties, make it a promising tool for enabling a probabilistic predictive framework for composites manufacturing digital twins.

36 MATERIALS SCIENCE↗

Image-domain DAS 3D VSP elastic transmission tomography

SUMMARY Full-wavefield elastic imaging of active-source seismic data acquired by downhole receivers commonly offers higher-resolution subsurface images in the vicinity of a borehole compared to conventional surface seismic data sets, which can lack higher-frequency wavefield components due to longer travel paths and increased attenuation. An increasingly used approach for downhole acquisition is vertical seismic profiling (VSP), which has become more attractive when coupled with distributed acoustic sensing (DAS) using optical fibres installed in wells. The main difficulty for generating high-quality images with full-wavefield imaging tools for DAS VSP data, though, is the need for an accurate velocity model. To build plausible velocity models using active-source DAS VSP data, we adopt a 3-D image-domain elastic transmission tomography technique, originally developed for surface-recorded passive (microseismic) data, by exchanging the source and receiver positions (i.e. reciprocity) to mimic a passive-seismic surface monitoring scenario. The inversion approach exploits various images for each source constructed through time-reverse imaging (TRI) of downgoing P- and S-wave first-arrival waveforms. The TRI process uses the kinetic term of the (extended) PS energy imaging condition that exhibits sufficient sensitivity to velocity model errors. The method automatically updates the P- and S-wave velocity models to optimize image focusing via adjoint-state inversion. We illustrate the efficacy of the adopted elastic inversion technique using an active-source DAS 3-D VSP field data set acquired in the North Slope of Alaska. The numerical experiments demonstrate that the inverted elastic velocity models can be further used in full-wavefield acoustic/elastic imaging algorithms to obtain accurate subsurface images.

Oren, Can (ORCID:000000019534265X)↗

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗

Anomalous electrodynamics and quantum geometry in the Dirac-Harper model for a graphene bilayer

In this work, graphene bilayers with layer antisymmetric strains are studied using the Dirac-Harper model for a pair of single layer Dirac Hamiltonians coupled by a one-dimensional moiré-periodic interlayer tunneling amplitude. This model hosts low energy, nearly dispersionless bands near charge neutrality that support anomalous polarizations of its charge multipole distributions. These are analyzed introducing a generalized Berry curvature that encodes the field-induced dynamics of multipole fields allowed in a chiral medium with time reversal symmetry. The formulation identifies a reciprocity relation between responses to layer-symmetric and layer-antisymmetric in-plane electric fields and reveals momentum-space quantum oscillations produced by a spatial pattern of band inversions on the moiré scale.

36 MATERIALS SCIENCE↗

The light-ray OPE and conformal colliders

We derive a nonperturbative, convergent operator product expansion (OPE) for null-integrated operators on the same null plane in a CFT. The objects appearing in the expansion are light-ray operators, whose matrix elements can be computed by the generalized Lorentzian inversion formula. For example, a product of average null energy (ANEC) operators has an expansion in the light-ray operators that appear in the stress-tensor OPE. An important application is to collider event shapes. The light-ray OPE gives a nonperturbative expansion for event shapes in special functions that we call celestial blocks. As an example, we apply the celestial block expansion to energy-energy correlators in N = 4 Super Yang-Mills theory. Using known OPE data, we find perfect agreement with previous results both at weak and strong coupling, and make new predictions at weak coupling through 4 loops (NNNLO).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kramers nodal lines and Weyl fermions in SmAlSi

Kramers nodal lines (KNLs) have recently been proposed theoretically as a special type of Weyl line degeneracy connecting time-reversal invariant momenta. KNLs are robust to spin orbit coupling and are inherent to all non-centrosymmetric achiral crystal structures, leading to unusual spin, magneto-electric, and optical properties. However, their existence in in real quantum materials has not been experimentally established. Here we gather the experimental evidence pointing at the presence of KNLs in SmAlSi, a non-centrosymmetric metal that develops incommensurate spin density wave order at low temperature. Using angle-resolved photoemission spectroscopy, density functional theory calculations, and magneto-transport methods, we provide evidence suggesting the presence of KNLs, together with observing Weyl fermions under the broken inversion symmetry in the paramagnetic phase of SmAlSi. We discuss the nesting possibilities regarding the emergent magnetic orders in SmAlSi. Our results provide a solid basis of experimental observations for exploring correlated topology in SmAlSi.

36 MATERIALS SCIENCE↗

Search for Non-resonant Higgs Boson Pair Production in the Four Bottom Quark Decay Channel With the CMS Experiment

This dissertation presents a search for non-resonant Higgs boson pair production, focusing on the four bottom quark decay channel. It explores the gluon fusion and vector boson fusion production mechanisms. The analysis is performed with a dataset of proton-proton collisions at a center-of-mass energy of 13 TeV, collected by the CMS detector at the LHC, corresponding to an integrated luminosity of around 138 inverse femtobarns. Innovative techniques in the areas of object identification, event categorization, signal identification, and background modeling are used to maximize the analysis sensitivity. No excess of signal events are observed relative to the background-only expectation, and 95% CL upper limits on the production cross section are set. The observed upper limit on the standard model production cross section is set at 3.6 times the theoretical expectation. At the time of writing this dissertation, it is the most stringent constraint at the LHC from an individual channel. Furthermore, the observed constraint on the coupling modifier of the Higgs boson self-interaction is set between -2.3 and 9.4. The observed constraint on the coupling modifier of the di-vector-boson-di-Higgs-boson interaction is set between -0.1 and 2.2.

Guerrero Ibarra, Guerrero Fernando↗

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Examining Cloud Feedback Components in the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM)

Cloud feedback remains the main source of uncertainty in climate sensitivity estimated by global climate models (GCMs), largely because subgrid cloud responses are parameterized in GCMs due to their coarse resolution. Here, this study examines cloud feedback in the global 3.25-km Simple Cloud-Resolving Energy Exascale Earth System Model (E3SM) Atmosphere Model (SCREAM 3 km) through a pair of 1-yr atmosphere-only simulations with control and +4-K sea surface temperature perturbations. SCREAM 3 km produces a positive cloud feedback that falls within but at the upper end of the range of Coupled Model Intercomparison Project phase 5 (CMIP5) and CMIP phase 6 (CMIP6) models and expert judgment. The positive cloud feedback arises from positive contributions from both high- and low-level clouds, with increases in high-cloud altitude and decreases in low-cloud amount and optical depth playing key roles. The stronger-than-CMIP-average feedback is mainly attributable to the high-cloud altitude feedback, owing to cloud tops rising nearly isothermally in SCREAM 3 km. The positive low-cloud amount feedback is weaker in SCREAM than in GCMs because estimated inversion strength (EIS) increases more dramatically with warming. A coarser 12-km resolution version of SCREAM exhibits a weaker positive cloud feedback than SCREAM 3 km, mainly because its low-cloud-radiative flux is more sensitive to EIS, leading to a stronger negative low-cloud amount feedback. With this process-level assessment of cloud feedback, this study reveals where SCREAM aligns with and diverges from conventional GCMs and expert assessment, providing insights to inform further model improvement and future expert assessment.

Cloud radiative effects↗

In-Operando XRD Study of the Effects of Water Vapor on Phase Formation in LSM/YSZ SOFC Cathodes

Here, in-operando XRD was conducted on anode supported SOFC button cells with LSM-YSZ cathodes operated at varying cathode air compositions under constant current conditions for over 1000 h. 1-hour XRD scans were continuously collected throughout the entire operation duration. By taking the sum of measured intensities from all of the collected patterns, the resulting cumulative XRD count times allowed the identification of minor phases present at concentrations as low as <0.1 wt%. In dry air with no contaminants, the cathode exhibited improving power output during the first couple of hundred hours, followed by stable operation. The effect of 3% H 2 O + 12% CO 2 on the LSM-YSZ cathode was very similar to the effect of 3% H 2 O alone, exhibiting performance degradation. Increasing contaminants in the cathode air were found to decrease the performance of the cells. In-operando XRD discovered an increase in MnO concentration and decrease in La 2 O 3 . A gradual expansion of the LSM lattice resulting from loss of Sr or O was discovered in the LSM/YSZ cathodes tested in humid cathode air. The inverse relationship between the unit cell volume and operating voltage suggests a possible correlation between Sr segregation and performance degradation in the LSM/YSZ cathodes.

25 ENERGY STORAGE↗

Observable proton decay in flipped SU(5)

We explore proton decay in a class of realistic supersymmetric flipped SU(5) models supplemented by a U(1) R symmetry which plays an essential role in implementing hybrid inflation. Two distinct neutrino mass models, based on inverse seesaw and type I seesaw, are identified, with the latter arising from the breaking of U(1) R by nonrenormalizable superpotential terms. Depending on the neutrino mass model an appropriate set of intermediate scale color triplets from the Higgs superfields play a key role in proton decay channels that include p → (e + , μ + ) π 0 , p → (e + , μ + ) K 0 , p → $\bar v$π + , and p → $\bar v$K + . We identify regions of the parameter space that yield proton lifetime estimates which are testable at Hyper-Kamiokande and other next generation experiments. We discuss how gauge coupling unification in the presence of intermediate scale particles is realized, and a Z 4 symmetry is utilized to show how such intermediate scales can arise in flipped SU(5). Finally, we compare our predictions for proton decay with previous work based on SU(5) and flipped SU(5).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Strong electron-phonon coupling driven pseudogap modulation and density-wave fluctuations in a correlated polar metal

Abstract There is tremendous interest in employing collective excitations of the lattice, spin, charge, and orbitals to tune strongly correlated electronic phenomena. We report such an effect in a ruthenate, Ca 3 Ru 2 O 7 , where two phonons with strong electron-phonon coupling modulate the electronic pseudogap as well as mediate charge and spin density wave fluctuations. Combining temperature-dependent Raman spectroscopy with density functional theory reveals two phonons, B 2 P and B 2 M , that are strongly coupled to electrons and whose scattering intensities respectively dominate in the pseudogap versus the metallic phases. The B 2 P squeezes the octahedra along the out of plane c -axis, while the B 2 M elongates it, thus modulating the Ru 4d orbital splitting and the bandwidth of the in-plane electron hopping; Thus, B 2 P opens the pseudogap, while B 2 M closes it. Moreover, the B 2 phonons mediate incoherent charge and spin density wave fluctuations, as evidenced by changes in the background electronic Raman scattering that exhibit unique symmetry signatures. The polar order breaks inversion symmetry, enabling infrared activity of these phonons, paving the way for coherent light-driven control of electronic transport.

36 MATERIALS SCIENCE↗

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗

Singlet fission

Singlet fission (SF) is a process whereby a photo-excited singlet exciton spontaneously splits into lower energy triplet excitons. This process has been known to occur in select organic materials for over 50 years but in the past was viewed as a curiosity and not as a problem worthy of focused attention. Over the past decade, this viewpoint has changed as the promising use of SF to improve photovoltaic efficiency has spurred vast activity, which has revealed the richness of the phenomenon. Currently, our understanding of the role played by various bright and dark states and their couplings to vibrations and to each other and by chemical composition and packing, among other factors, has enabled a large increase in the number of viable SF systems as well as a far deeper understanding of the process itself. These advances have allowed for the demonstration of proof-of-concept increases in the solar cell efficiency via the use of SF materials in simple tandem cells and as sensitizers in conjunction with conventional photovoltaics. Indeed, concepts surrounding the field have now invaded other fields, with a prominent example being the use of photon upconversion associated with triplet–triplet fusion (the inverse process to SF) in organic chemistry and photocatalysis. The 23 research articles in this special topic showcase recent cutting-edge advances in the field and highlight the diverse new directions being explored in the field of SF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗