Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Here, this study overviews and extends a recently developed stochastic finite-temperature Kohn-Sham density functional theory to study warm dense matter using Langevin dynamics, specifically under periodic boundary conditions. The method's algorithmic complexity exhibits nearly linear scaling with system size and is inversely proportional to the temperature. Additionally, a linear-scaling stochastic approach is introduced to assess the Kubo-Greenwood conductivity, demonstrating exceptional stability for dc conductivity. Utilizing the developed tools, we investigate the equation of state, radial distribution, and electronic conductivity of hydrogen at a temperature of 30 000 K. As for the radial distribution functions, we reveal a transition of hydrogen from gaslike to liquidlike behavior as its density exceeds 4 g/cm 3 . As for the electronic conductivity as a function of the density, we identified a remarkable isosbestic point at frequencies around 7 eV, which may be an additional signature of a gas-liquid transition in hydrogen at 30 000 K.

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

Unavoidable multilevel biaxial symmetry breaking in chiral hybrid liquid crystals

Chiral nematic or cholesteric liquid crystals (LCs) are fluid mesophases with long-ranged orientational order featuring a quasilayered periodicity imparted by a helical director configuration but lacking long-range positional order. Doping molecular cholesteric LCs with strongly anisotropic uniaxial colloidal particles adds another level of complexity because of the interplay between weak surface-anchoring boundary conditions and bulk-based elastic distortions near the particle-LC interface. Using cylindrical colloidal disks and rods with different geometric shapes and surface conditions, we demonstrate that these colloidal inclusions generically exhibit biaxial orientational probability distributions which may impart anomalously strong local biaxiality onto the hybrid cholesteric LC structure. Unlike nonchiral hybrid molecular-colloidal LCs, where biaxial order emerges only at critical colloid volume fractions exceeding some uniaxial-biaxial transition value, the orientational probability of the colloidal inclusions immersed in chiral nematic hosts is unambiguously biaxial even at infinite dilution. We demonstrate that the colloids induce local biaxial perturbations within the molecular orientational order of the LC host medium which strongly enhances the weak but native biaxial order of chiral nematic LC induced by the chiral symmetry breaking of the director field. With the help of analytical modeling and computer simulations based on the Landau–de Gennes free energy of the host LC around the colloids, we rationalize the observed multilevel biaxial order and conclude that it is not only unavoidable but also strongly enhanced compared to both achiral hybrid LCs and purely molecular cholesteric LCs.

biaxial liquid crystals↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2021

Learning unknown physics of non-Newtonian fluids

Here, we use physics-informed neural networks (PINNs) to learn viscosity models of two non-Newtonian systems (polymer melts and suspensions of particles) using only velocity measurements. For synthetic velocity data generated with the power-law viscosity model, the PINN-inferred viscosity model agrees with the analytical model for shear rates with large absolute values but deviates for shear rates near zero where the analytical model has an unphysical singularity. Once the viscosity model is learned the PINN method can solve the momentum conservation equation using only the boundary conditions.

42 ENGINEERING↗
DOE OSTI2022

Dimensionless parameters for cloudy Rayleigh-Bénard convection: Supersaturation, Damköhler, and Nusselt numbers

In steady-state Rayleigh-Bénard convection, heat is transported by turbulent thermal convection from the bottom, hot surface to the top, cold surface, leading to a height-independent sensible heat flux. When water vapor is present and cloud formation occurs, there is also an additional latent heat flux. Heat transport in cloudy Rayleigh-Bénard convection depends on turbulent flow as well as the microphysical state of the clouds: specifically, whether substantial supersaturations exist and whether cloud liquid water is removed through sedimentation/precipitation. In this article we bridge between the Rayleigh-Bénard convection literature and the atmospheric literature. We express the governing equations for cloudy convection in dimensionless form, thereby explicitly identifying the governing parameters relevant to the cloudy case, including Schmidt, Damköhler, supersaturation, and sedimentation numbers. We further connect to the atmospheric literature by obtaining a Nusselt number (dimensionless heat flux) for a cloud-convection system, directly from the conservation equations for temperature and water vapor. This flux has the same form as that identified by Zhang et al. [L. Zhang, K. L. Chong, and K.-Q. Xia, J. Fluid Mech. 874, 1041 (2019)] for convection with water vapor, but is extended to the cloudy case. For equal thermal and water vapor diffusivities, the flux corresponds to the widely used atmospheric quantities equivalent temperature and moist static energy. Using large eddy simulation (LES) of an idealized cloudy Rayleigh-Bénard convection system with fixed boundary conditions, we find that the equivalent heat flux (Nusselt number) is only weakly dependent on the microphysical details of the system, such as liquid water mixing ratio and cloud droplet number concentration. Finally, from the results, we show the vertical profiles of sensible and latent heat fluxes depend on the liquid water content, whereas the equivalent heat flux remains a constant throughout the height of the chamber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Topological Defect Engineering and $\mathcal{PT}$ Symmetry in Non-Hermitian Electrical Circuits

In this work, we employ electric circuit networks to study topological states of matter in non-Hermitian systems enriched by parity-time symmetry $\mathcal{PT}$ and chiral symmetry anti-$\mathcal{PT (APT)}$. The topological structure manifests itself in the complex admittance bands which yields excellent measurability and signal to noise ratio. We analyze the impact of $\mathcal{PT}$-symmetric gain and loss on localized edge and defect states in a non-Hermitian Su-Schrieffer-Heeger (SSH) circuit. We realize all three symmetry phases of the system, including the $\mathcal{APT}$-symmetric regime that occurs at large gain and loss. We measure the admittance spectrum and eigenstates for arbitrary boundary conditions, which allows us to resolve not only topological edge states, but also a novel $\mathcal{PT}$-symmetric $\mathbb{Z}_2$ invariant of the bulk. We discover the distinct properties of topological edge states and defect states in the phase diagram. In the regime that is not $\mathcal{PT}$ symmetric, the topological defect state disappears and only reemerges when $\mathcal{APT}$ symmetry is reached, while the topological edge states always prevail and only experience a shift in eigenvalue. Our findings unveil a future route for topological defect engineering and tuning in non-Hermitian systems of arbitrary dimension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Convert Widespread Paraelectric Perovskite to Ferroelectrics

While nature provides a plethora of perovskite materials, only a few exhibit large ferroelectricity and possibly multiferroicity. The majority of perovskite materials have the nonpolar CaTiO 3 (CTO) structure, limiting the scope of their applications. Based on the effective Hamiltonian model as well as first-principles calculations, we propose a general thin-film design method to stabilize the functional BiFeO 3 (BFO)- type structure, which is a common metastable structure in widespread CTO-type perovskite oxides. It is found that the improper antiferroelectricity in CTO-type perovskite and ferroelectricity in BFO-type perovskite have distinct dependences on mechanical and electric boundary conditions, both of which involve oxygen octahedral rotation and tilt. Thus the above difference can be used to stabilize the highly polar BFO-type structure in many CTO-type perovskite materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Global Data-Driven Determination of Baryon Transition Form Factors

Hadronic resonances emerge from strong interactions encoding the dynamics of quarks and gluons. The structure of these resonances can be probed by virtual photons parametrized in transition form factors. Here, in this study, twelve N* and Δ transition form factors at the pole are extracted from data with the center-of-mass energy from πN threshold to 1.8 GeV, and the photon virtuality 0 ≤ Q 2 /GeV 2 ≤ 8. For the first time, these results are determined from a simultaneous analysis of more than one state, i.e., ~10 5 π⁢N, η⁢N, and K⁢Λ electroproduction data. In addition, about 5×10 4 data in the hadronic sector as well as photoproduction serve as boundary conditions. For the Δ⁡(1232) and N⁡(1440) states our results are in qualitative agreement with previous studies, while the transition form factors at the poles of some higher excited states are estimated for the first time. Realistic uncertainties are determined by further exploring the parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Domain fluctuations in a ferroelectric low-strain BaTiO 3 thin film

In this work, a ferroelectric BaTiO 3 thin film grown on a NdScO 3 substrate was studied using x-ray photon correlation spectroscopy (XPCS) to characterize thermal fluctuations near the a/b to a/c domain structure transformation present in this low-strain material, which is absent in the bulk. XPCS studies provide a direct comparison of the role of domain fluctuations in first- and second-order phase transformations. The a/b to a/c domain transformation is accompanied by a decrease in fluctuation timescales, and an increase in intensity and correlation length. Surprisingly, domain fluctuations are observed up to 25 °C above the transformation, concomitant with the growth of a/c domains and coexistence of both domain types. After a small window of stability, as the Curie temperature is approached, a/c domain fluctuations are observed, albeit slower, potentially due to the structural transformation associated with the ferroelectric to paraelectric transformation. The observed time evolution and reconfiguration of domain patterns highlight the role played by phase coexistence and elastic boundary conditions in altering fluctuation timescales in ferroelectric thin films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗
DOE OSTI2024

Predicting Thermochemical Equilibria with Interacting Defects: Sr 1 − x Ce x Mn O 3 − δ Alloys for Water Splitting

Solar thermochemical hydrogen is one of the few potential routes towards direct fuel production from renewable energy sources, but the thermodynamic boundary conditions for efficient and economic energy conversion are challenging. Success or failure of a given oxide working material depends on the subtle balance between enthalpy and entropy contributions in the redox processes. Developing a mechanistic understanding of the behavior of materials on the basis of atomistic models and first-principles calculations is an important part of advancing the technology. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox-active defects start to interact with each other, thereby impeding the formation of additional defects. This problem is of more general importance to applications that rely on high levels of off-stoichiometry or doping, including, for example, batteries, thermoelectrics, and ceramic fuel cells. To account for such repulsive defect interactions, we introduce a statistical mechanics approach, defining an expression for the free energy of defect interaction based on limited sampling of defect configurations in density functional theory supercell calculations. The parameterization of this energy contribution as a function of defect concentration and temperature allows on-the-fly simulation of thermochemical equilibria. The approach consistently incorporates finite temperature effects by including the leading contributions to the temperature-dependent free energy for the case at hand, i.e., the ideal gas and configurational enthalpies and entropies. We demonstrate the capability and utility of the approach by simulating the water splitting redox processes for Sr 1 − x Ce x Mn O 3 − δ alloys. Published by the American Physical Society 2024

08 HYDROGEN↗
DOE OSTI2026

Quantum simulation of QED in Coulomb gauge

A recent work considered quantum simulation of Quantum Electrodynamics on a lattice in the Coulomb gauge with gauge degrees of freedom represented in the occupation basis in momentum space. In this work, we consider the more efficient representation of the gauge degrees of freedom in field basis in position space and develop a quantum algorithm for real-time simulation. We show that the continuum Coulomb gauge Hamiltonian is equivalent to the temporal gauge Hamiltonian when acting on physical states consisting of fermion and transverse gauge fields. The Coulomb gauge Hamiltonian is discretized by using the Green’s function of the discrete Laplacian operator under the Dirichlet boundary conditions. Both the continuum Coulomb gauge Hamiltonian and the discretized one proposed here guarantee that the unphysical longitudinal gauge fields are decoupled and commute with the corresponding Hamiltonian. Thus there is no need to impose any constraint. The local gauge field basis and the canonically conjugate variable basis are swapped efficiently using the quantum Fourier transform. We prove that the qubit cost to represent physical states and the gate count for real-time simulation scale polynomially with the lattice size, energy, time, accuracy, and Hamiltonian parameters in lattice units. The gate cost here for implementing the time evolution of the gauge field is reduced at least by a factor on the order of 10 8 for modest lattice size and accuracy level compared with the previous work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Asymptotic gauge symmetry and UV extension of the nonperturbative coupling in holographic QCD

We extend our recent analytic study of the strong coupling 𝛼 eff in the nonperturbative and near-perturbative regimes [Phys. Rev. Lett. 133, 181901 (2024)] by imposing rigorous renormalization-group constraints from asymptotically free gauge theories at 𝑄 2 → ∞. The asymptotic boundary conditions modify the scaling properties of 𝛼eff at large values of the momentum transfer 𝑄 2 and lead to a scale-dependent confinement strength 𝜅⁡(𝑄 2 ). This requires that both 𝜅⁡(𝑄 2 ) and 𝛼 eff⁡ (𝑄 2 , 𝜅⁡(𝑄 2 )) remain holomorphic in the complex 𝑄 2 plane, except at the physical cuts associated with the heavy-quark thresholds and the singularity flow trajectory studied in our previous Letter. For color SU(3), a precise connection is found between the scaling exponent of 𝜅⁡(𝑄 2 ) in the ultraviolet, the value of the infrared fixed point of the strong coupling, and the number of flavors in agreement with observations. The nonperturbative analytic model gives an accurate description of the strong coupling across all scales, up to the highest available data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Irreversible thermodynamics of curved lipid membranes. II. Permeability and osmosis

We present a theory that combines the framework of irreversible thermodynamics with modified integral theorems to model arbitrarily curved and deforming membranes immersed in bulk fluid solutions. We study the coupling between the mechanics and permeability of a viscous and elastically bendable membrane, and a multicomponent bulk fluid solution. An equation for the internal entropy production for irreversibilities at the membrane is derived, determining the generalized thermodynamic forces and fluxes, from which we identify the deviatoric stress as a novel driving force for permeability. A complete set of equations of motion, constitutive laws, and boundary conditions to model the lipid membrane and bulk fluid system are provided.

Alkadri, Ahmad M↗
DOE OSTI2025

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗