A multipoint stress mixed finite element method for elasticity on quadrilateral grids
Not provided.
SEARCH · Search NASA
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.
Not provided.
In the present work we propose a nonlinear anti- P T -symmetric dimer, that at the linear level has been experimentally created in the realm of electric circuit resonators. We find four families of solutions, the so-called upper and lower branches, both in a symmetric and in an asymmetric (symmetry-broken) form. We unveil analytically and confirm numerically the critical thresholds for the existence of such branches and explore the bifurcations (such as saddle-node ones) that delimit their existence, as well as transcritical ones that lead to their potential exchange of stability. We find that out of the four relevant branches, only one, the upper symmetric branch, corresponds to a spectrally and dynamically robust solution. We subsequently leverage detailed direct numerical computations in order to explore the dynamics of the different states, corroborating our spectral analysis results.
Synchrotron accelerator facilities such as the NSLS II require extreme stability, both transient (short-term) and quasi-static (long-term) to achieve the desired resolution performance. Consequently, even μm-level movements, particularly differential movements between locations in the concrete structure supporting the accelerator electron beam lattice (storage ring) or high sensitivity experiments (experimental floor) will lead to serious degradation of its performance. Differential settlement in the overall structure or structural movement exceeding anticipated levels will inevitably degrade the performance and will require intervention. Presented in this paper are the design philosophy of the NSLS II ring structure favoring a monolithic ring, the observed cracking behavior of the young NSLS II concrete following casting and in combination with extreme ambient temperature fluctuation, the results of a non-linear, high-fidelity numerical analysis used to emulate the observed cracking and establish the driving mechanism, the numerical analysis-based identification of the crack-arresting solution and finally the implementation of the remediation solution and the long-term performance of the adopted engineering solution. The multi-stage process revealed that computational methods such as non-linear finite element methods have the potential of providing engineering guidance even when complex structures and in combination with non-linear materials, such as steel reinforcement and concrete are involved.
This study reports the design and development of a novel chemiresistor (CR) sensor using ion imprinted polymer (IIP)-functionalized reduced graphene oxide (rGO) [IIP/rGO-CR] for cadmium ions (Cd(II)) determination in water. The sensor consisted of a CR transducer made of rGO channel bridging source and drain electrodes prepared by self-assembly and thermal reduction of graphene oxide (GO) on Au interdigitated electrodes chip fabricated on Si/SiO 2 substrate. The IIP was then grafted on rGO using surface-initiated reversible addition-fragmentation chain transfer (RAFT) polymerization with polyethylenimine (PEI) and methylacrylic acid (MAA) as dual functional monomers and Cd(II) ions as template through UV light-initiated copolymerization. The IIP functionalized on rGO acted as an effective recognition element that modulated the resistance of rGO-CR upon binding of Cd(II), enabling Cd(II) detection at ppb level in aqueous solutions. The prepared IIP/rGO-CR sensor worked effectively in the linear range of 2~200ppb and achieved a limit of detection (LOD) of 0.83 ppb, which is lower than the World Health Organization guidelines of 3ppb for drinking water quality. The developed sensor of IIP/rGO-CR showed a high selectivity against a variety of trace and heavy metal ions found in water and good stability for up to 60 days when stored at room temperature for Cd(II) determination in water. Further, the sensor was successfully applied to analyzing Cd(II) spiked in tap, lake and river waters with a 94.5%–113.5% recovery, demonstrating a high degree of accuracy even in complex water samples. Finally, our results illustrated that the CR sensor of IIP functionalized rGO provides a potential platform for sensitive, robust and low-cost environmental analysis of Cd(II) in water.
This study employs an explainable machine learning (ML) framework (XGBoost-SHapley Additive exPlanations analysis) to investigate controlling factors on cloud liquid water path (LWP) using EPCAPE observations near the California coast. Aerosols are found to be the dominant factor explaining LWP variability, surpassing meteorological factors (MFs). By isolating aerosol effects from meteorological influences, the ML reveals a negative linear relationship between LWP and cloud droplet number concentration (Nd) in log space, likely driven by entrainment drying via evaporation-entrainment feedback. This aligns with the negative regime of the inverted-V relationship reported in previous studies, while no positive LWP responses are found due to a limited number of precipitating cases in EPCAPE. Furthermore, the sensitivity of LWP to Nd shows a non-linear dependence on MFs like moisture contrast between surface and free troposphere and lower-tropospheric stability. This occurs due to the interplay between the MFs' direct effects on entrainment drying and indirect effects through LWP adjustments.
We develop an analytic theory to describe the interaction between electrons and K phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K phonon (160meV) of graphene is responsible for the intervalley electron-phonon interaction. Here, this phonon has recently been found in angular-resolved photoemission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of the phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest T c pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy-fermion mapping of TBG which has allowed us to obtain an analytic form of a phonon-mediated attractive interaction and to analytically solve the linearized and T = 0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat-band limit. While the flat-band Coulomb interaction can be screened sufficiently enough—around the Van Hove singularities—to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states.
In earlier papers (Leblond et al., 2011; 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a “double” propagation criterion combining (Griffith, 1920)’s energetic condition and (Goldstein and Salganik, 1974)’s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the “threshold” value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.
In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.
Poly(4–4′-oxydiphenylene-pyromellitimide)-based polyimides─trade name Kapton─have wide-ranging engineering applications owing to their thermal and mechanical stability, but little is known about underlying chain-level characteristics. While theoretical models have conceptualized Kapton as inflexible polycyclic rods separated by freely rotating diphenyl ether hinge groups, the model’s core predictions remain untested and subtleties of the relaxation behavior are missed, which atomistic modeling can resolve. To these ends, we examine all-atom Kapton structures in crystalline and glassy amorphous configurations using a DFT-validated class II force field. Constructing amorphous configurations is challenging, as the fused-ring-containing backbone has slow relaxation dynamics and scaling suggestive of entanglements even in oligomers. In conclusion, we find larger backbone rearrangements of the linear polycyclic segments about ether groups that are consistent with the rod-hinge picture on the monomer scale, whereas a ring rotation analysis suggests partially flexible rod-like segments and involves multiple facile rotational relaxation modes.
Metal silicon phosphides composed of earth-abundant Si and P tend to exhibit semiconducting properties and adopt diverse crystal structures with relatively small additions of structure-directing elements. The potential of silicon phosphide materials in nonlinear optical applications has been hindered by the inability to systematically produce noncentrosymmetric structures with such a flexible framework. Here, in this work, two isostructural compounds with a novel noncentrosymmetric structure were made possible by the inclusion of elements with stereochemically active lone pairs (Sn 2+ and Pb 2+ ). The structures were determined through single-crystal and synchrotron powder X-ray diffraction. Analysis of chemical bonding in real space through the electron localization function revealed stereochemically active Pb 2+ and Sn 2+ species in a trigonal pyramidal coordination with {Pb/Sn}–P bonds. Such covalent bonding between Pb and P is quite uncommon in extended solids and has been reported in a few rare instances. Band structure calculations and linear optical measurements confirm the semiconducting nature of Cs X Si 15 P 21 ( X = Sn or Pb). The synthesis was optimized to yield high-purity polycrystalline samples. The nonlinear optical properties show promising second-harmonic generation (SHG) coefficients from the Kurtz–Perry method. First-principles calculations of the nonlinear optical properties support the experimentally determined SHG values and provide moderate values of birefringence, suggesting Cs X Si 15 P 21 could be phase-matchable and practical nonlinear optical materials in the mid-IR region.
Understanding the structural diversity of honeybee-infecting viruses is critical to maintain pollinator health and manage the spread of diseases in ecology and agriculture. We determine cryo-EM structures of $T$ = 4 and $T$ = 3 capsids of virus-like particles (VLPs) of Lake Sinai virus (LSV) 2 and delta-N48 LSV1, belonging to tetraviruses, at resolutions of 2.3–2.6 Å in various pH environments. Structural analysis shows that the LSV2 capsid protein (CP) structural features, particularly the protruding domain and C-arm, differ from those of other tetraviruses. The anchor loop on the central β-barrel domain interacts with the neighboring subunit to stabilize homo-trimeric capsomeres during assembly. Delta-N48 LSV1 CP interacts with ssRNA via the rigid helix α1’, α1’–α1 loop, β-barrel domain, and C-arm. Cryo-EM reconstructions, combined with X-ray crystallographic and small-angle scattering analyses, indicate that pH affects capsid conformations by regulating reversible dynamic particle motions and sizes of LSV2 VLPs. C-arms exist in all LSV2 and delta-N48 LSV1 VLPs across varied pH conditions, indicating that autoproteolysis cleavage is not required for LSV maturation. The observed linear domino-scaffold structures of various lengths, made up of trapezoid-shape capsomeres, provide a basis for icosahedral $T$ = 4 and $T$ = 3 architecture assemblies. These findings advance understanding of honeybee-infecting viruses that can cause Colony Collapse Disorder.
We report that biofuels are particularly attractive and play an increasingly important role in sustainable energy. However, biofuels originating from lignocellulosic biomass (LCB) are extremely challenging because of their low carbon content, low stability, and high oxygen content. This work evaluates the supercritical methanol (scMeOH) co-liquefaction of rice straw and linear low-density polyethylene (LLDPE) at temperature range of 240-340 °C for 0-2 h, to obtain hydrocarbons (HCs)-rich oil and carbon-rich solid product. Results show that reaction temperature dominated the yield and properties of products, but not the holding time. Among parameters tested, 30.07 wt% oil yield with 75.79% HCs content and 33.05 wt% oil yield with 70.91% HCs content were obtained at 300 °C for 1 h and 1.5 h, respectively. Simultaneously, the remaining solid products were still as high as 53.85 wt% with a carbon content of 79.59% and 48.28 wt% with carbon content of 81.34% under 300 °C for 1 h and 1.5 h, respectively. Ultimate analysis, FT-IR, TGA, and SEM show that solid products could be used as sustainable carbon resources, and solid fuel rather than soil amendment because of risk of micro plastic or adsorbent due to smooth surface without pores.
The emergent contaminant family, per- and poly-fluorinated alkyl substances (PFAS) has gained research attention due to their widespread detection and stability within the environment. Despite the growing amount of research on perfluorooctanesulfonic acid (PFOS) and perfluoro-n-octanoic acid (PFOA) in aquatic organisms, investigations detailing behavioral and physiological effects of aquatic organisms exposed to a mixture of PFAS analytes in the wild have been limited. The objective of this study was to evaluate the potential behavioral and histological effects of environmental exposure to PFAS compounds within multiple trophic levels of aquatic ecosystems. The current study investigates effects of environmentally relevant PFAS concentration exposures in crayfish (Faxonius immunis, F. rusticus, F. virilis) and bluegill (Lepomis macrochirus) sourced from four water bodies in Northern Michigan. Antipredator response and foraging behavioral assays were used to investigate potential effects on crayfish; a swimming speed behavioral assay and liver and gill histology analysis were used to investigate potential effects on fish. Linear mixed model and multiple regression analyses resulted in significant relationships between tissue accumulation levels of long chain PFAS compounds and crayfish foraging and fish critical swimming speed responses. Crayfish foraging decreased and fish critical swim speeds increased with PFAS exposure which may lead to energetic and population concerns. Antipredator response in crayfish and liver and gill histology in fish were not significantly related to PFAS tissue or water concentrations. The sensitivity of crayfish and bluegill behavior contributes to the growing body of research regarding the differential toxicity of short-chain and long-chain PFAS compounds. The sensitivity of some aquatic organism behaviors to PFAS accumulated in tissue may have implications for PFAS transfer and alterations to ecosystem functioning; based on the results of this field study, further laboratory research is recommended to further evaluate these relationships.
Here, a straightforward in situ detection method for dengue infection was demonstrated through the molecular imprinting of a dengue nonstructural protein 1 (NS1) epitope into an electropolymerized molecularly imprinted polyterthiophene (E-MIP) film sensor. The key enabling step in the sensor fabrication is based on an epitope imprinting strategy, in which short peptide sequences derived from the original target molecules were employed as the main template for detection and analysis. The formation of the E-MIP sensor films was facilitated using cyclic voltammetry (CV) and monitored in situ by electrochemical quartz crystal microbalance (EC-QCM). Surface properties were analyzed using different techniques including atomic force microscopy (AFM), X-ray photoelectron spectroscopy (XPS), and polarization modulation-infrared reflection-adsorption (PM-IRRAS). The standard calibration curve (R = 0.9830) was generated for the detection of the epitope, Ac-VHTWTEQYKFQ-NH 2 , with a linear range of 0.2 to 30 μg/mL and detection limit of 0.073 μg/mL. A separate calibration curve (R = 0.9786) was obtained using spiked buffered solutions of dengue NS1 protein, which resulted in a linear range of 0.2 to 10 μg/mL and a detection limit of 0.056 μg/mL. The fabricated E-MIP sensor exhibited long-term stability, high sensitivity, and good selectivity towards the targeted molecules. These results indicated that the formation of the exact and stable cavity imprints in terms of size, shape, and functionalities was successful. In our future work, we aim to use our E-MIP sensors for NS1 detection in real-life samples such as serum and blood.
The influence of elastic strains on the adsorption of H, O, and OH on the (111) surfaces of 8 fcc (Ni, Cu, Pd, Ag, Pt, Au, Rh, Ir) and on the (0001) surfaces of 3 hcp (Co, Zn, Cd) transition metals was analyzed by means of density functional theory calculations. To this end, surface slabs were subjected to different strain states (uniaxial, biaxial, shear, and a combination of them) up to strains dictated by the mechanical stability limits indicated by phonon calculations. It was found that the adsorption energy followed the predictions of the d-band theory but – surprisingly – the variations in the adsorption energy only depended on the area of the adsorption hole and not on the particular elastic strain tensor applied to achieve this area. The analysis of the electronic structure showed that the applied strains did not modify the shape of the Projected Density of States (PDOS) of the d-orbitals of the transition metals but only led to a shift in the energy levels. Moreover, the presence of the adsorbates on the surfaces led to negligible changes in the PDOS. Thus, the adsorption energies were a function of the Fermi energy which in turn was associated with the change of the area of the adsorption through a general linear law that was valid for all metals. Furthermore, the information in this paper allows the immediate and accurate estimation of the effect of any elastic strain on the adsorption energies of H, O, and OH in 11 transition metals with more than half-filled d-orbitals.
The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.
Bastnäsites (LnCO 3 (F,OH)) are a group of common rare earth elements (REE)-bearing minerals and are one of the primary global sources of REE. Due to the chemical similarities among REE, bastnӓsites tend to occur as solid solutions instead of end members in REE containing ores. To better understand the processes and the mechanisms of formation of such deposits, it is essential to determine the thermodynamic properties of bastnӓsites, including hydroxylbastnӓsite (LnCO 3 OH) solid solutions. In this work, we performed detailed structural and calorimetric investigations on synthetic hexagonal La–Nd hydroxylbastnӓsite (La 1-$\chi$ NdCO 3 OH, $\chi$ = 0, 0.25, 0.5, 0.75, 1) solid solutions. X-ray diffraction confirms the crystal structure of the solid solution series in the $\ P$6 space group, and pair distribution function (PDF) analysis reveals local bonding environments characterized by three different types of 9-coordinated metal-oxygen polyhedra. Unit cell parameters of La1–xNdxCO3OH exhibit a nearly linear relation with the Nd content x, suggesting a random distribution of La and Nd in the structure. Their standard enthalpies of formation (Δ$\ H$° f ) were determined by high temperature oxide melt drop solution calorimetry, from which the enthalpies of mixing (Δ$\ H$ mix ) were derived. The Δ$\ H$ mix can be fitted by a regular solution model with an interaction parameter of 12.58 ± 0.16 kJ/mol, suggesting enthalpic metastability of La1–xNdxCO3OH relative to the two endmembers. Combining entropy and enthalpy, we further estimated the Gibbs free energies of mixing (Δ$\ G$ mix ) at relevant temperatures, revealing favorable temperatures under which the intermediate La 1-$\chi$ NdCO 3 OH phases can be stabilized. Such entropy-driven stabilization, as is consistent with our geochemical modeling results, may explain the enhancement of thermal stability of the solid solutions in nature. Additionally, the temperature range constrained from this study may be used to estimate the thermal history of REE bastnӓsite deposit.
Here, we introduce a new framework to study the dynamics of open quantum systems with linearly coupled Gaussian baths. Our approach replaces the continuous bath with an auxiliary discrete set of pseudomodes with dissipative dynamics, but we further relax the complete positivity requirement in the Lindblad master equation and formulate a quasi-Lindblad pseudomode theory. We show that this quasi-Lindblad pseudomode formulation directly leads to a representation of the bath correlation function in terms of a complex weighted sum of complex exponentials, an expansion that is known to be rapidly convergent in practice and thus leads to a compact set of pseudomodes. The pseudomode representation is not unique and can differ by a gauge choice. When the global dynamics can be simulated exactly, the system dynamics is unique and independent of the specific pseudomode representation. However, the gauge choice may affect the stability of the global dynamics, and we provide an analysis of why and when the global dynamics can retain stability despite losing positivity. We showcase the performance of this formulation across various spectral densities in both bosonic and fermionic problems, finding significant improvements over conventional pseudomode formulations.