Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Polar system modeling”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Structure of an Atmospheric River over Australia and the Southern Ocean. Part I: Tropical and Midlatitude Water Vapor Fluxes

An atmospheric river (AR) impacting Australia-Tasmania and the Southern Ocean during the Austral summer on 28-29 January 2018 during the SOCRATES campaign is analyzed using both a modeling and observational approach. Gulfstream-V dropsonde measurements and GPM radar analyses were used in conjunction with Weather Research and Forecasting model simulations with water vapor tracers to investigate the relative contributions of tropical and mid-latitude moisture sources to the AR. Moisture associated with a monsoonal tropical depression became entrained into a mid-latitude frontal system that extended to 60°S over the Southern Ocean reaching the associated low pressure system 850 km off the east coast of Antarctica – effectively connecting the tropics and the polar region. Tropical moisture contributed to about 50% of the precipitable water within the AR. The tropical contribution to precipitation decreased with latitude, from > 70% over the Australian continent, to ~50% off the Australian coast, to less than 10% poleward of 60°S. Precipitation in the AR conformed to a seeder-feeder model, the seeder precipitation forming primarily from tropical sourced moisture above the altitude of the 0°C isotherm, and the feeder precipitation sourced from mid-latitude sourced moisture below the altitude of the 0°C isotherm. The integrated vapor transport (IVT) through the core of the AR reached above 500 kg m-1s-1 during 1200 UTC 28 January to 0600 UTC 29 January 29, 1.29 times the average amount of water carried by the world’s largest river, the Amazon. The high IVT strength might be attributed to the higher water vapor content associated with the warmer temperatures in Austral summer.

Rauber, Robert M.↗
DOE OSTI2021

Unraveling Strain Gradient Induced Electromechanical Coupling in Twisted Double Bilayer Graphene Moiré Superlattices

Moiré superlattices of 2D materials with a small twist angle are thought to exhibit appreciable flexoelectric effect, though unambiguous confirmation of their flexoelectricity is challenging due to artifacts associated with commonly used piezoresponse force microscopy (PFM). For example, unexpectedly small phase contrast (≈8°) between opposite flexoelectric polarizations is reported in twisted bilayer graphene (tBG), though theoretically predicted value is 180°. We report a methodology is developed to extract intrinsic moiré flexoelectricity using twisted double bilayer graphene (tDBG) as a model system, probed by lateral PFM. For small twist angle samples, it is found that a vectorial decomposition is essential to recover the small intrinsic flexoelectric response at domain walls from a large background signal. The obtained threefold symmetry of commensurate domains with significant flexoelectric response at domain walls is fully consistent with the theoretical calculations. Incommensurate domains in tDBG with relatively large twist angles can also be observed by this technique. A general strategy is provided here for unraveling intrinsic flexoelectricity in van der Waals moiré superlattices while providing insights into engineered symmetry breaking in centrosymmetric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Perturbative theoretical model of electronic transient circular dichroism spectroscopy of molecular aggregates

A chiral analog of transient absorption spectroscopy, transient circular dichroism (TCD) spectroscopy is an emerging time-resolved method. Both spectroscopic methods can probe the electronic transitions of a sample, and TCD is additionally sensitive to the dynamic aspects of chirality, such as those induced by molecular excitons. Here, we develop a theoretical description of TCD for electronic multi-level models in which the pump pulse is linearly polarized and probe pulse is alternately left- and right-circularly polarized. We derive effective response functions analogous to those often used to describe other four-wave mixing methods and then simulate and analyze TCD spectra for three representative multi-level electronic model systems. We elaborate on the presence and detection of the spectral signatures of electronic coherences.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Atomic-Scale Imaging of Polarization Switching in an (Anti-)Ferroelectric Memory Material: Zirconia (ZrO2)

Direct, atomic-scale visualization of polarization switching in a functional, polycrystalline, binary oxide via in-situ high-resolution transmission electron microscopy (HRTEM) biasing is reported for the first time. Antiferroelectric (AFE) ZrO2 was used as the model system, which is important for commercial DRAMs and as emerging NVMs (through work-function engineering). We observed (1) clear shifting and coalescing of domains within a single grain, and (2) dramatic changes of the atomic arrangements and crystalline phases—both at voltages above the critical voltage measured for AFE switching. Similar synergistic in-situ structural-electrical characterization can pave the way to understand and engineer microscopic mechanisms for retention, fatigue, variability, sub-coercive switching and analog states in ferroelectric and AFE-based memory devices.

Lombardo, Sarah↗
DOE OSTI2024

Data and scripts associated with the manuscript "Estimating soil moisture and salinity response to simulated coastal flooding using time-lapse electrical resistivity and induced polarization monitoring”

This package contains geophysics datasets generated during the simulated ecosystem-flooding experiment – TEMPEST (Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatment). These geophysics data consist of Electrical Resistivity Imaging (ERI) and Induced Polarization (IP) datasets collected to estimate the soil moisture and salinity response to the simulated coastal flooding experiment. In addition to using geophysics datasets to estimate soil moisture and salinity response, petrophysical relationship between the measured resistivity and moisture content and salinity using the multisalinity. This package consists of two folders: “MultiSalinity Experiment” folder and “TEMPEST Experiment and petrophysical conversion” folders. The “MultiSalinity Experiment” folder contains the datasets generated from the multisalinity experiment used to develop the petrophysical reslationship. The “TEMPEST Experiment and petrophysical conversion” folder contains ERI, IP, Ground Penetrating Radar (GPR) and Electromagnetic Induction (EMI) datasets collected during the TEMPEST experiment. This study shows the ability to used geophysical datasets in estimating moisture content and salinity at scale will be useful in calibrating Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2024

Higher‐Order Gravity Waves and Traveling Ionospheric Disturbances From the Polar Vortex Jet on 11–15 January 2016: Modeling With HIAMCM‐SAMI3 and Comparison With Observations in the Thermosphere and Ionosphere

Abstract In Vadas et al. (2024, https://doi.org/10.1029/2024ja032521 ), we modeled the atmospheric gravity waves (GWs) during 11–14 January 2016 using the HIAMCM, and found that the polar vortex jet generates medium to large‐scale, higher‐order GWs in the thermosphere. In this paper, we model the traveling ionospheric disturbances (TIDs) generated by these GWs using the HIAMCM‐SAMI3 and compare with ionospheric observations from ground‐based Global Navigation Satellite System (GNSS) receivers, Incoherent Scatter Radars (ISR) and the Super Dual Auroral Radar Network (SuperDARN). We find that medium to large‐scale TIDs are generated worldwide by the higher‐order GWs from this event. Many of the TIDs over Europe and Asia have concentric ring/arc‐like structure, and most of those over North/South America have planar wave structure and occur during the daytime. Those over North/South America propagate southward and are generated by higher‐order GWs from Europe/Asia which propagate over the Arctic. These latter TIDs can be misidentified as arising from geomagnetic forcing. We find that the higher‐order GWs that propagate to Africa and Brazil from Europe may aid in the formation of equatorial plasma bubbles (EPBs) there. We find that the simulated GWs, TIDs and EPBs agree with EISCAT, PFISR, GNSS, and SuperDARN measurements. We find that the higher‐order GWs are concentrated at N at 200 km, in agreement with GOCE and CHAMP data. Thus the polar vortex jet is important for generating TIDs in the northern winter ionosphere via multi‐step vertical coupling through GWs.

Vadas, Sharon L. [Northwest Research Associates Bo↗
DOE OSTI2024

Bayesian Conavigation: Dynamic Designing of the Material Digital Twins via Active Learning

Scientific advancement is universally based on the dynamic interplay between theoretical insights, modeling, and experimental discoveries. However, this feedback loop is often slow, including delayed community interactions and the gradual integration of experimental data into theoretical frameworks. This challenge is particularly exacerbated in domains dealing with high-dimensional object spaces, such as molecules and complex microstructures. Hence, the integration of theory within automated and autonomous experimental setups, or theory in the loop-automated experiment, is emerging as a crucial objective for accelerating scientific research. The critical aspect is to use not only theory but also on-the-fly theory updates during the experiment. Furthermore, we introduce a method for integrating theory into the loop through Bayesian conavigation of theoretical model space and experimentation. Our approach leverages the concurrent development of surrogate models for both simulation and experimental domains at the rates determined by latencies and costs of experiments and computation, alongside the adjustment of control parameters within theoretical models to minimize epistemic uncertainty over the experimental object spaces. This methodology facilitates the creation of digital twins of material structures, encompassing both the surrogate model of behavior that includes the correlative part and the theoretical model itself. While being demonstrated here within the context of functional responses in ferroelectric materials, our approach holds promise for broader applications, such as the exploration of optical properties in nanoclusters, microstructure-dependent properties in complex materials, and properties of molecular systems.

Microscopy↗
DOE OSTI2020

Determining Activity Coefficients of SOA from Isothermal Evaporation in a Laboratory Chamber

Nonideal molecular interactions in aerosol particles influence the partitioning of semivolatile organic compounds (SVOCs). However, few direct measurements exist that determine activity coefficients (γ), which quantify nonideality, for individual organic compounds in different atmospheric mixtures. By measuring the isothermal evaporation of SOA formed from multifunctional SVOCs in an environmental chamber, this study determines the mole-fraction-based γ of the bulk SOA in a variety of pre-existing chemically complex particles, approximated as if the bulk SOA behaved like a single SVOC. These multicomponent mixtures contain molecular structures commonly found in the atmosphere. When treated as a single liquid phase, γ of the bulk SOA tends to increase from ~1 to ~5 as the organic seeds and the SOA have more differing polarities. A high computed γ value of 74 for a wet ammonium sulfate-SOA system indicates phase separation. The γ for some individual species was also quantified based on gas-phase measurements and AIOMFAC model predictions. However, the bulk SOA γ cannot be explained by the simplified speciated SVOC–seed interactions, suggesting incomplete compositional understanding of the bulk SOA. In this work, these results demonstrate a method to quantify nonideal behavior and show that it can occur in multicomponent mixtures and therefore influence OA formation, evolution, and lifetimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Manipulating Growth and Propagation of Correlations in Dipolar Multilayers: From Pair Production to Bosonic Kitaev Models

Here, we study the nonequilibrium dynamics of dipoles confined in multiple stacked two-dimensional layers realizing a long-range interacting quantum spin 1/2 XXX model. We demonstrate that strong in-plane interactions can protect a manifold of collective layer dynamics. This then allows us to map the many-body spin dynamics to bosonic models. In a bilayer configuration we show how to engineer the paradigmatic two-mode squeezing Hamiltonian known from quantum optics, resulting in exponential production of entangled pairs and generation of metrologically useful entanglement from initially prepared product states. In multilayer configurations we engineer a bosonic variant of the Kitaev model displaying chiral propagation along the layer direction. Our study illustrates how the control over interactions, lattice geometry, and state preparation in interacting dipolar systems uniquely afforded by AMO platforms such as Rydberg and magnetic atoms, polar molecules, or trapped ions allows for the control over the temporal and spatial propagation of correlations for applications in quantum sensing and quantum simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Valley-spin polarization at zero magnetic field induced by strong hole-hole interactions in monolayer WSe 2

Monolayer transition metal dichalcogenides have emerged as prominent candidates to explore the complex interplay between spin and valley degrees of freedom. Their strong spin-orbit interaction and broken inversion symmetry lead to the spin-valley locking effect, in which carriers occupying theKandK′ valleys of the reciprocal space must have opposite spins. This effect is particularly strong for holes due to a larger spin-orbit gap in the valence band. By reducing the dimensionality of a monolayer of WSe 2 to 1D via electrostatic confinement, we demonstrate that spin-valley locking and strong hole-hole interactions lead to a ferromagnetic state where hole transport is spin-valley polarized, even without an applied magnetic field. A massive Dirac fermion model in the Hartree-Fock approximation reveals that many-body hole-exchange interactions lead to this polarized ground-state. This observation opens the possibility of implementing a robust and stable valley-polarized system, essential in valleytronic applications.

Science & Technology - Other Topics↗
DOE OSTI2024

Comparing the influence of cation order and composition in simulated Zn(Sn, Ge)N 2 on structure, elastic moduli, and polarization for solid state lighting

Alloying and site ordering play complementary roles in dictating a material’s properties. However, deconvolving the impacts of these separate phenomena can be challenging. In this work, we simulate structures of Zn(Sn,Ge)N 2 with varied Sn content and site ordering to determine the impacts of order and composition on structural and electronic properties. We assess the formation enthalpy, lattice parameters, elastic constants, spontaneous polarization, and piezoelectric coefficients. In mostly disordered structures (order parameters ranging from 0.2 to 0.4), the formation enthalpy exhibits local extrema as a function of the order parameter, deviating from the more linear trends seen in both fully disordered and fully ordered systems. This anomalous deviation from the otherwise linear trend in formation enthalpy with order manifests in each of the other properties calculated. This range of order parameters of interest may be caused by a transition in the ordering of the quaternary material similar to phase changes seen in ternary compounds but stretched over a region inclduing 20% of the order parameter range. Most parameters calculated are more sensitive to order than to composition in the limited composition range tested; however, the lattice parameter c, piezoelectric coefficient e 33 , and elastic moduli C 12 , C 13 , and C 23 are more sensitive to composition. Of the properties compared, the piezoelectric coefficients are influenced most significantly by changes in both the composition and order parameter. Lattice parameters undergo the smallest changes with order and composition, but these small differences appear to impart large trends in the other properties. Better understanding the effects of disorder and group IV alloying in Zn(Sn,Ge)N 2 allows for more accurate modeling of characteristics of this material system for solid state lighting and other applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Hydration and Restructuring of Polar Polymer Interfaces: Implications in Antifouling and Responsive Materials

Manipulating polymer interfaces is crucial for understanding how structure influences function in applications spanning biofouling prevention to energy storage. Moreover, observing how polymers adapt their microscopic structure to changes in their local environment can reveal essential properties that govern their performance in such applications, providing key insights into how to design more effective interfaces. Here, in this study, a series of “grafting-from” polymer brushes with side chains varying in charge, are probed by sum frequency generation (SFG) and modeled using all-atom molecular dynamics (MD) simulations to elucidate how chemical makeup and charge mediate interfacial restructuring in dry versus hydrated states. Results show that charge, in progressing from nonpolar to cationic to zwitterionic, results in dramatic changes in interfacial structure and overall hydration. While net neutral systems, regardless of bulk phase polarity, show minimal interfacial water structuring, the cationic species exhibits strong bulk water signals from the surface potential. Meanwhile, the polymer brushes themselves restructure in water differently independent of charge, impacting the functional groups that are presented to the aqueous phase. Nonpolar and cationic species for instance undergo a change in alkyl group orientations to accommodate hydrating water molecules, whereas the zwitterionic polymer becomes completely disordered in water. Overall, the structure-based behavior trends presented herein have implications in antifouling applications and responsive material interfaces.

biointerface↗
DOE OSTI2021

Maximizing selectivity: An analysis of isoporous membranes

Membranes are compelling candidates for many separation-based industrial processes. However, all synthetic membranes are subject to well-recognized permeability-selectivity trade-offs that hinder their use in advanced applications. Having uniform nanoscale pore size, isoporous nanofiltration and ultrafiltration membranes have attracted wide research interest. However, quantitative analysis of isoporous membrane performance is lacking - particularly at module and system-level scales. In this study, an analytical model is built to investigate how membrane pore-size distribution affects the solute separation efficiency (i.e., rejection and selectivity). The model accounts for convective hindrance as well as concentration polarization and explores effects on separation from the solute-to-pore-size ratio and percent recovery (i.e., fraction of feed flow that passes through the membrane). At the module level, our analysis suggests that a notably sharper rejection vs. solute size curve can be obtained as the pore-size distribution becomes narrower. Additionally, as the recovery increases, the rejection of a solute is maintained at a high level using isoporous membranes while the rejection drops rapidly for membranes that possess a disperse pore-size distribution. Furthermore, in the separation of mixed solutes, isoporous membranes achieve a large selectivity with a modest loss of permeability. Overall, this study elaborates the unique role that isoporous membranes play within the permeability-selectivity trade-off relations and identifies the most exciting opportunities for their future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Dynamics of Electric Polarization and Relaxation of Ions at Humid Calcite Surfaces

Mobile ions at mineral surfaces can respond to an applied electric field, adopting a new distribution that effectively represents polarization of the electrical double layer. When the field is released, the ions relax to their equilibrium distribution. In both cases, the dynamics are characteristic of the interface. However, current models of electrokinetic phenomena are not sufficiently robust to accurately predict collective ion dynamics at structurally and chemically specific mineral–water interfaces. Here, in this study, we use electrostatic force microscopy (EFM) to investigate the dynamics of ion relaxation at hydrated calcite (104) surfaces at controlled relative humidity (RH). Electrically biased probes are used to polarize the distributions of calcium and carbonate ions that are intrinsic to this interface across a range of RH values. Polarization kinetics are tracked by monitoring the tip–sample force gradient during charging, and EFM imaging is used to characterize the spatial relaxation dynamics after the applied field is released. Electrostatic finite element modeling of the sample/probe system across length-scales from nanometers to millimeters reproduces the observed stretched exponential charging response. Together, these results allow us to estimate the ion diffusivities at the interface across a wide range of RH values. These diffusivities increase by roughly 5 orders of magnitude as the RH is increased from 5 to 90%, highlighting the critical role of adsorbed water for surface ion solvation that enables ion mobility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Antarctic Water Masses and Ice Shelves: Visualizing the Physics

High-resolution simulation of global climate physics enables us to model how the climate may change under a variety of future scenarios. Such simulations produce vast amounts of information and dense datasets. If interrogated in tandem, these datasets can provide holistic, vital information on Earth’s many integrated systems by revealing the manifold interrelated properties of the atmosphere, ocean, and polar ice, framed by real-world terrain in three-dimensional space as they vary over time. To accomplish this, climate scientists have joined with computer scientists and an artist to develop techniques enabling scientists to see these relationships. Finally, the impact of ocean water properties on Antarctic ice shelves illustrates the benefit of this analysis in understanding land ice melt rates and thus sea-level rise.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Non-adiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (wherethe electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase- space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Furthermore, the F-PSSH algorithm sets the stage for modeling nonadiabatic dynamics under strong externally pumped circular polarization in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Magnetized laser–plasma interactions in high-energy-density systems: Parallel propagation

In this study, we investigate parametric processes in magnetized plasmas, driven by a large-amplitude pump light wave. Our focus is on laser–plasma interactions relevant to high-energy-density (HED) systems, such as the National Ignition Facility and the Sandia MagLIF concept. We present a self-contained derivation of a “parametric” dispersion relation for magnetized three-wave interactions, meaning the pump wave is included in the equilibrium, similar to the unmagnetized work of Drake et al., Phys. Fluids 17, 778 (1974). For this, we use a multi-species plasma fluid model and Maxwell's equations. The application of an external B field causes right- and left-polarized light waves to propagate with differing phase velocities. This leads to Faraday rotation of the polarization, which can be significant in HED conditions. Phase-matching and linear wave dispersion relations show that Raman and Brillouin scattering have modified spectra due to the background B field, though this effect is usually small in systems of current practical interest. We study a scattering process we call stimulated whistler scattering, where a light wave decays to an electromagnetic whistler wave (ω≲ω ce ) and a Langmuir wave. This only occurs in the presence of an external B field, which is required for the whistler wave to exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗