Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constant”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mineral Scale Formation during Crossflow Reverse Osmosis at Constant Flux and Constant Transmembrane Pressure Conditions

Mineral scale formation on membrane surfaces is a significant challenge in reverse osmosis water purification. Laboratory fouling experiments are typically run such that the transmembrane pressure (TMP) is fixed, and the permeate flux decreases over time as scales accumulate on the membrane surface. However, this change in flux means that the hydrodynamic conditions at the membrane surface are continuously changing, which could affect crystallization and foulant deposition processes. Operating under constant permeate flux conditions, in contrast, is advantageous because it keeps the hydrodynamic conditions relatively consistent, making it possible to compare how membrane properties (e.g., surface chemistry) affect fouling propensity. Industrial reverse osmosis operations are not run strictly in either constant TMP or constant flux mode; while they may start at a constant TMP, feed pressure may be periodically adjusted to maintain permeate water production within a specified range. Here, the scarcity of constant permeate flux reverse osmosis scaling experiments reported in the literature frustrates efforts to compare membrane fouling processes under constant TMP and constant flux conditions. For the first time, the evolution of the fouling layer resistance was compared as a function of cumulative permeate volume per membrane area during constant TMP and constant flux reverse osmosis filtrations. Scaling experiments were conducted by challenging commercial reverse osmosis membranes with a model feed solution nearly saturated with calcium sulfate dihydrate (gypsum). At low fluxes, the increase in fouling layer resistance was quantitatively similar for the two operational modes. In contrast, at high fluxes, the fouling layer resistance increased more rapidly in constant flux filtration than in constant TMP filtration. The mechanism of scale formation in constant TMP and constant flux operation was self-limiting and self-reinforcing, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Second‐ and Third‐Order Elastic Constants of Inert and Energetic Molecular Crystals From Density Functional Theory

Complete tensors of the second- and third-order elastic constants of the organic molecular crystals acetaminophen, pentaerythritol tetranitrate (PETN), cyclotrimethylene trinitramine (RDX), cyclotetramethylene tetranitramine (HMX), 1,1-diamino-2,2-dinitroethylene (FOX-7), hexanitrohexaazaisowurtzitane (CL-20), and erythritol tetranitrate (ETN) have been calculated using dispersion-corrected density functional theory. The sets of second- and third-order elastic constants are expected to provide a more accurate and reliable description of the behavior of these materials under nonhydrostatic loads than pressure- and volume-dependent second-order elastic constants. The tensors of second-order constants have been compared with experimental data and/or other calculations when possible, and with the exception of results for CL-20 from Brillouin scattering experiments, we find good agreement. The calculated third-order elastic constants of PETN are in very good agreement with the subset of third-order constants derived from experimental wave speed measurements. The elastic anisotropies of the crystals have been estimated using the universal elastic anisotropy index, which shows that the crystals fall into three groups with low anisotropy (PETN, RDX, and CL-20), moderate anisotropy (acetaminophen, HMX, and ETN), and high elastic anisotropy (FOX-7).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗
DOE OSTI2026

An Infrared Database of n/k Optical Constants for Calculating Aerosol Spectra for the PICARD Program

The objective of the PICARD Program is to develop fieldable sensing platforms for the rapid chemical identification of aerosol particles in plumes. Standoff detection involves interrogating the aerosol cloud from a distance using optical methods and probing the signal returned from direct backscattering from the aerosol particles or transmitted through the plume after reflection from a retroreflector or surface of opportunity (SOO). The identification of chemical species, however, in aerosols is complicated by their complex compositions and morphologies, chemical interferants, and non-uniform particle sizes. To advance standoff detection of aerosols, modelling the infrared transmittance, reflectance and scattering spectra of aerosolized liquid and solid chemical compounds is required, and then testing that model experimentally via laboratory and field experiments. To perform accurate modeling, the infrared optical constants (n/k), i.e., the complex refractive index, of the compounds of interest are required. Thus, PNNL was tasked to provide the optical constants for a set of analytes relevant to the PICARD program. This report describes the experimental techniques used to derive the optical constants of both liquid and solid compounds using established “gold-standard” protocols, how the experimental data are processed to produce the wavenumber dependent optical constant vectors, and how the data are used in the aerosol absorption spectra modeling.

complex refractive index↗
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2025

Beyond leading twist: 𝜌 meson decay constants and distribution amplitudes in a self-consistent light-front quark model

In this study, we present a comprehensive analysis of decay constants and chiral-even and chiral-odd distribution amplitudes (DAs) up to twist 4 for the 𝜌 meson in the standard light-front quark model (LFQM) based on the Bakamjian-Thomas (BT) construction. For the 𝜌 meson, which possesses both longitudinal (ℎ=0) and transverse (ℎ =±1) polarizations, two types of decay constants, 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$, arises accordingly. We demonstrate that these decay constants can be self-consistently extracted from both local (𝑧𝜇 =0) and nonlocal (𝑧𝜇 ≠0) matrix elements ⟨0⁢|$\overline{𝑞}$(𝑧)⁢Γ⁢𝑞⁡(−𝑧)|⁢𝜌⁡(𝑃,ℎ)⟩, with Γ=(𝛾 𝜇 ,𝜎 𝜇⁢𝜈 ,𝛾 𝜇 ⁢𝛾 5 ,𝟏), in a manner independent of current components, polarizations, and reference frames. In particular, we emphasize the role of nonlocal matrix elements involving axial-vector and scalar currents, where mixing between 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$ occurs. We show that this mixing is consistently resolved through the BT construction, ensuring the proper extraction of these decay constants.

Electroweak interaction↗
DOE OSTI2024

Constant-Depth Preparation of Matrix Product States with Adaptive Quantum Circuits

Adaptive quantum circuits, which combine local unitary gates, midcircuit measurements, and feedforward operations, have recently emerged as a promising avenue for efficient state preparation, particularly on near-term quantum devices limited to shallow-depth circuits. Matrix product states (MPS) comprise a significant class of many-body entangled states, efficiently describing the ground states of one-dimensional gapped local Hamiltonians and finding applications in a number of recent quantum algorithms. Recently, it has been shown that the Affleck-Kennedy-Lieb-Tasaki state—a paradigmatic example of an MPS—can be exactly prepared with an adaptive quantum circuit of constant depth, an impossible feat with local unitary gates alone due to its nonzero correlation length [Smith , PRX Quantum 4, 020315 (2023)]. In this work, we broaden the scope of this approach and demonstrate that a diverse class of MPS can be exactly prepared using constant-depth adaptive quantum circuits, outperforming theoretically optimal preparation with unitary circuits. We show that this class includes short- and long-ranged entangled MPS, symmetry-protected topological (SPT) and symmetry-broken states, MPS with finite Abelian, non-Abelian, and continuous symmetries, resource states for MBQC, and families of states with tunable correlation length. Moreover, we illustrate the utility of our framework for designing constant-depth sampling protocols, such as for random MPS or for generating MPS in a particular SPT phase. We present sufficient conditions for particular MPS to be preparable in constant time, with global on-site symmetry playing a pivotal role. Altogether, this work demonstrates the immense promise of adaptive quantum circuits for efficiently preparing many-body entangled states and provides explicit algorithms that outperform known protocols to prepare an essential class of states. Published by the American Physical Society 2024

Smith, Kevin C. (ORCID:0000000223971518)↗
DOE OSTI2026

A Solid State Zwitterionic Plastic Crystal With High Static Dielectric Constant

Solid materials with a high dielectric constant have a wide range of applications in the energy storage field. In this research, an imidazolium-based zwitterion is designed, synthesized, and confirmed to have a plastic crystal phase based on the following experimental and computational evidence: (i) the presence of long-range order with weak intermolecular forces and competing attractive-repulsive interactions along different crystallographic directions; (ii) the observation of more than one endotherm on heating including a solid-solid phase transition at T s-s = −26 °C and melting of the plastic crystal at T m = 72°C; (iii) a low entropy of fusion at melting (2.1 JK −1 mol −1 ); (iv) a strongly anisotropic morphology; (v) relatively fast dynamics originating from short-range degrees of freedom. Furthermore, it exhibits a very high dielectric constant in the plastic crystal solid state (147 at −10°C and 103 at 70°C) due to the rotational degrees of freedom of plastic crystals that arise from weak net intermolecular interactions of zwitterions due to only having two carbons between the anion and cation. This material conveniently remains in the plastic crystal phase within 50 K of ambient temperature. This discovery opens new opportunities in the search for solid-state high dielectric constant materials.

dielectric relaxation↗
DOE OSTI2026

Probing fundamental constant oscillation in the Galactic Center with S-Star Spectroscopy

Astrophysical spectroscopy provides a powerful probe of spacetime variations of fundamental constants, as atomic and ionic emission and absorption lines depend sensitively on the fine-structure constant. In particular, coherent temporal oscillations induced by an ultralight scalar background produce characteristic, time-resolved signatures that can be robustly disentangled from intrinsic variability. In the Galactic Center, such scalar backgrounds can be substantially enhanced, either through the formation of dense scalar clouds powered by black hole rotational energy extraction or as ultralight scalar dark matter forming a soliton-like core. These scalar configurations generically induce oscillations of the fine-structure constant, with periods set by the scalar mass and spatial profiles determined by the scalar wavefunction and its coupling to the electromagnetic sector. We show that precise, time-resolved spectroscopy of S-stars orbiting the supermassive black hole Sgr A^* provides a sensitive test of these effects, enabling constraints on quadratic scalar-photon couplings in the exceptionally high boson-density environment of the Galactic Center.

Bai, Zhaoyu [Weizmann Inst.] (ORCID:00000002075815↗
DOE OSTI2025

Comparing Constant and Transient Membrane Transport Parameters for Use in Wave Desalination Models

Directly pressurizing seawater for desalination with reverse osmosis membranes via wave motion is a promising and sustainable method for producing freshwater in coastal regions. However, such a system could result in significant pressure fluctuations and a departure from conventional steady-state desalination operations. This study sought to assess if membrane transport parameters (apparent water and salt permeability) should be modeled as transient or constant in solution–diffusion-based modeling efforts of dynamically operated desalination systems, such as those coupled to wave power. Two approaches were used to model membrane transport parameters: one considered each parameter to be a function of the net driving pressure of the system, and the other assumed they were constant across all conditions. A pilot-scale system was used to conduct steady-state and controlled ramping experiments. Data from steady-state experiments were used to calculate transient and constant transport parameters. Parameter combinations were used in a simulation model to predict water flux and effective permeate salinity, and simulation outcomes were compared against experimental ramping results. The transient relationships for both water and salt permeability produced the most accurate results for water flux and comparable results for effective permeate salinity. Development of such relationships would be unique to a specific system but could be valuable in modeling wave-driven desalination systems across the wide range of operating conditions they experience.

13 HYDRO ENERGY↗
DOE OSTI2025

Redshifting the Cosmological Constant in Unimodular Gravity via Nonlinear Quantum Mechanics

The cosmological constant problem represents a profound conflict between quantum field theory and general relativity. Unimodular gravity offers a compelling starting point by de-gravitating the vacuum energy of the Standard Model, but this framework traditionally trades the problem of vacuum energy for a fine-tuning of initial conditions, which manifest as a ``shadow" cosmological constant. In this paper, we resolve this initial conditions problem by proposing a novel modification to gravity based on nonlinear quantum mechanics. We introduce specific state-dependent terms to the Hamiltonian, constructed from expectation values of the metric such as the average Ricci scalar. These terms alter the dynamical equations of gravity such that the shadow energy density associated with unconstrained initial conditions redshifts away with cosmic expansion, rendering it negligible at late times. The resulting cosmology is naturally dominated by matter and radiation without fine-tuning. We demonstrate that this significant infrared modification of gravity is consistent with local and cosmological tests of gravity. We comment on the possibility of testing this solution in cosmological measurements of Newton's constant.

Kaplan, David E. [Johns Hopkins U.; Tokyo U., IPMU↗
DOE OSTI2024

Free-Charge Carrier Generation in Homojunction Non-Polymeric Organic Semiconductor Films - The Role of the Optical Frequency Dielectric Constant

Engineering the dielectric constant (..epsilon..) to lower the exciton binding energy of the light-absorbing semiconductor can improve organic photovoltaic (OPV) device performance. Here, a series of materials are reported with 2-(3-oxo-2,3-dihydro-1H-inden-1-ylidene)malononitrile (INCN) acceptor end groups and a central glycolated bis(4H-cyclopenta[2,1-b:3,4-b']dithiophene) unit with large low-frequency (..epsilon..lf = 7.4-7.9 at 0.1-0.2 MHz) and optical-frequency (..epsilon..opt up to 6.6 at 2 x 1014 Hz) dielectric constants. The INCN end groups differed in whether they were protonated, chlorinated, or fluorinated, with the latter having the highest ..epsilon..opt. An ..epsilon..opt of 6.6 is predicted to lead to a low exciton binding energy of ~0.04 eV. Time-resolved microwave conductivity measurements showed a temperature-dependent yield-mobility product, with it increasing linearly from 340 K. The onset temperature was near that required to overcome the calculated exciton binding energy and indicates increased free charge generation in a homojunction film. Room temperature transient absorption spectroscopy revealed that photoexcitation rapidly converted to a lower energy state that was consistent with the formation of polarons or a charge transfer state. This work provides experimental evidence of the importance of ..epsilon..opt for the generation of free charges, and a strategy for development of efficient single chromophore homojunction OPV devices.

charge generation↗
DOE OSTI2025

Modeling of Inductive Constant Power Load for Electromagnetic-Transient Simulations-Part II

This paper improves the dynamic constant power (CP) load model that was published in Part I, which is appropriate for electromagnetic-transient (EMT). The improved model conserves all features of its predecessor. For instance, it maintains a fixed power consumption (both active and reactive parts) and a fixed power factor for loads that are predominantly inductive. Furthermore, as the proposed model is a time-dependent system, it is applicable to both sinusoidal and non-sinusoidal case studies. However, the previous model cannot be easily integrated with numeric solvers because it simulated load data over one cycle all together, not sequentially in a time-step manner, due to the limitation involved with the power factor. The improved version, on the contrary, allows the load to be simulated at every time step, which would facilitate its integration with numeric solvers. The model's validity is confirmed by comparing its response with data that is synthesized from constant impedance load, and the result is satisfactory.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Complete coverage fouling model for constant flux crossflow ultrafiltration and experimental validation

A complete coverage model is proposed to describe fouling in constant flux crossflow ultrafiltration. Constant flux crossflow fouling experiments were conducted using dilute latex bead suspensions and commercial poly(ether sulfone) flat sheet ultrafiltration membranes to investigate the influence of operating conditions on evolution of transmembrane pressure (ΔP, TMP) with time. Changes in permeate flux or crossflow rate had little influence on the normalized TMP profile at high latex bead concentration (i.e., above 25 ppm) because the membrane surface was covered with latex beads. At low concentration (i.e., below 25 ppm), increases in permeate flux or foulant concentration increased normalized ΔP. However, this increase in normalized ΔP with permeate flux or foulant concentration diminishes when the permeate flux/concentration is high enough to overwhelm particle removal due to crossflow. Furthermore, these results are in good agreement with the new complete coverage model which describes the influence of operating parameters on fouling better than the previous model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2024

Tackling the Giants: Applying Smart Labs Principles to Constant Air Volume Lab Buildings

Laboratories typically consume 3 to 10 times more energy than similarly sized commercial buildings, and as much as 50% of that energy is wasted by inefficient and poorly operating fume hoods and ventilation systems. One challenge faced by older laboratory buildings is the heating, ventilation, and air-conditioning systems serving many of these buildings. The older systems are usually constant air volume (CAV) systems that maintain constant ventilation rates that cause excess airflow and inefficient energy use. Variable air volume systems can be more efficient systems with sensors to detect the need for a change in volumetric flow rate; however; renovation of ventilation systems can create disruption to ongoing research and operations along with considerable up-front costs. When a Smart Labs program is implemented, an organization has a systems-based management approach that yields a high-performing laboratory building. As decarbonization continues as a priority for sites, buildings with CAV systems are difficult to address. This work centers around practical guidance for improving lab buildings with CAV. In conjunction with industry input on top technology solutions and best practices, recommendations will include performing a laboratory ventilation risk assessment in conjunction with robust retro-commissioning work, which is a crucial step in the Smart Lab process. By applying Smart Labs principles, the aging laboratory building stock of 153,343 (Lawrence Berkeley National Laboratory [LBNL] 2017), comprising roughly 500,000 lab spaces in the United States, can be brought to safe and high-performance operations.

building↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗