Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution theory”

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 361 records · Page 20

Precise 3D reactor core calculation using spherical harmonics and discontinuous Galerkin finite element methods

We study the use of P{sub N} method in angle and discontinuous Galerkin is space to solve 3D neutron transport problem. P{sub N} method consists in developing the angular flux on truncated spherical harmonics basic. In this paper, we couple this method with the discontinuous finite elements in space to obtain a complete discretization of the multigroup neutron transport equation. To investigate its precision, the method was applied to Takeda and C5G7 benchmark problems. These calculations point out that the proposed P{sub N}-DG method is capable of producing accurate solutions in small computational time, and that it is able to handle complex 3D geometries. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

TCC in the interior of moduli space and its implications for the string landscape and cosmology

We consider the classical Friedmann-Robertson-Walker solutions that describe a universe undergoing a transition from an accelerating expansion phase in the past to an eternal decelerating expansion phase in the future, driven by a scalar field evolving in a potential energy landscape. We show that any solution for which the accelerating phase violates the Trans-Planckian Censorship Conjecture (TCC), even in the interior of moduli space, never approaches the asymptotic vacuum with zero particles. Based on the assumption that the effective field theory must be valid for the vacuum on the asymptotic boundary, as motivated by holography and string theory, we argue that (multi-field) scalar potentials with such solutions are disallowed, thus strengthening the case for TCC. In particular, assuming the regularity of the future vacuum state in the string landscape, we derive results that imply a new set of highly-nonlinear constraints across the string landscape which in the absence of certain meta-stable vacua make realizing inflation impossible.

Cosmological models↗

Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

42 ENGINEERING↗

Hydrogen Bond Network Disruption by Hydration Layers in Water Solutions with Salt and Hydrogen-Bonding Polymers (PEO)

A mean field theory model describing the interaction of ion hydration layers with the network of hydrogen bonds of both water and the nonionic polymer poly(ethylene oxide) (PEO) is presented. The predictions of the model for types and statistics of hydrogen bonds, the number of water molecules bound to PEO, or their dependence on temperature are successfully verified from all-atom simulations at different NaCl and PEO concentrations. Furthermore, our simulations show that the binding of cations to PEO increases monotonically with salt concentration, in agreement with recent experimental results, through a mechanism in which the sum of the number of bound water and cations is independent of salt concentration. Here, the model introduced is general and can describe any salt or hydrogen-bond-forming polymer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Einstein-Maxwell theory and the Weyl double copy

The Weyl double copy relates vacuum solutions in general relativity to Abelian gauge fields in Minkowski spacetime. In a previous work, we showed how the Weyl double copy can be extended to provide a treatment of external gravitational sources consistent with the classical Kerr-Schild double copy. Using this generalization, here we provide a complete double copy analysis of electrovacuum Petrov type D spacetimes. This includes the first analysis of the charged C-metric, whose single copy interpretation invokes the two-potential formalism of electrodynamics. Here, we also present the first double copy prescription for the Ricci spinor, which for non-accelerating spacetimes, takes a form similar to the original double copy relation for the Weyl spinor.

79 ASTRONOMY AND ASTROPHYSICS↗

The efficacy of Lewis affinity scale metrics to represent solvent interactions with reagent salts in all-inorganic metal halide perovskite solutions

Solvents employed in the solution processing of metal halide perovskites are known to play a key role in defining the morphology and properties of the resulting thin film, and thus the performance of perovskite solar cell devices. Accurate metrics are needed that are capable of differentiating among candidates, finding solvents that adequately solubilize the various precursor species in solution and facilitate the nucleation and growth of these materials. Existing metrics such as the unsaturated Mayer bond order (UMBO) and the Gutmann donor number (DN) have been tested for lead iodide perovskite systems; but there has yet to be a comprehensive study on their transferability to lead-free perovskite solutions. Here, we use ab initio methods (density functional theory) and regression analysis tools to study the usefulness of DN and BF 3 affinity scales in this regard. We compared the relative effectiveness of these scales to describe interactions between solvents and BXn perovskite salts of lead (Pb 2+ ), tin (Sn 2+ and Sn 4+ ), germanium (Ge 2+ ), bismuth (Bi 3+ ), and antimony (Sb 3+ and Sb 5+ ). The DN proved to be a better representation than the BF 3 of such interactions, reflecting the closer similarity of these species to the “parent” SbCl 5 Lewis acid than to BF 3 . In addition, we have uncovered the usefulness of the lithium cation affinity metric (LCA) to describe the strength of interactions between solvents and A-site cations (e.g. Na + , K + , Rb + and Cs + ) in all-inorganic metal halide perovskite solutions. We find that the coordination strengths of solvents towards species in all-inorganic metal halide perovskite solutions are best described by two different metrics with distinct modes of action: DN differentiates among BX n salt complexes, and LCA among A-site cation species. This revelation can help guide the choice of solvent to optimize processing conditions. It also emphasizes the importance of selecting solvents whose DN and LCA optimize coordination to key Lewis acid species in all-inorganic perovskite solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photon phase-space dynamics in a plasma wakefield accelerator

Frequency up-shifting of laser light in a beam-driven plasma wakefield has the potential to provide high-intensity sources of short wavelength radiation. Simulations have demonstrated that a laser pulse can undergo large frequency shifts, limited only by the drive beam energy, when the plasma density is tailored to match the accelerating phase of the wake to the group velocity of the pulse. Here, we study the dynamical evolution of photons in the 1D1P phase-space vicinity of the plasma wake-phase matching condition. Numerical calculations using a photon kinetic model are validated by direct comparison with 1D and quasi-3D particle-in-cell simulations. These calculations form the basis of a linear theory of the photon dynamics which reveals several important results, including scalings for the properties of the witness pulse and a self-similar solution for the photon phase-space dynamics. One prediction of the analytic theory is that the pulse can be compressed indefinitely, though the amount of compression would ultimately be limited by practical constraints. These results suggest that photon acceleration can provide a novel source of sub-femtosecond, short wavelength radiation.

XUV generation↗

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

74 ATOMIC AND MOLECULAR PHYSICS↗

Designing linear lattices for round beam in electron storage rings using the solution by linear matrices analysis

For some synchrotron light source beamline applications, a round beam is preferable to a flat one. A conventional method of obtaining round beam in an electron storage ring is to shift its tune close to a linear difference resonance. The linearly coupled beam dynamics is analyzed with perturbation theories, which have certain limitations. In this paper, we adopt the Solution by LInear Matrices (SLIM) analysis to calculate exact beam sizes to design round beam lattices. The SLIM analysis can deal with a generally linearly coupled accelerator lattice. The effects of various coupling sources on beam emittances and sizes can be studied within a self-consistent frame. Both the on- and off-resonance schemes to obtain round beams are explained with examples. The SLIM formalism for two widely used magnet models: combined-function bending magnets, and planar wigglers and undulators, is also derived.

43 PARTICLE ACCELERATORS↗

Sensitivity measurements for a 250 MHz quartz shear-horizontal surface acoustic wave biosensor under liquid viscous loading

Surface acoustic wave (SAW) devices have been used in biochemical assays due to their high sensitivity. The device sensitivity is a function of changes in the density and viscosity of the liquid. Here, we studied the effect of fluid viscosity using a 250 MHz quartz shear-horizontal (SH)-SAW biosensor by monitoring different concentrations of binary aqueous/glycerol solutions. In this study, the sensitivity of the biosensor was determined by fitting the data to models derived from perturbation theory. Measurements in water were used as the reference. For a 0% to 50% glycerol solution, an 87°–204° separation in the phase shift was observed. The slope of the plot of the phase shift vs (ηρ)0.5 was used to indicate the sensor’s sensitivity. The sensitivity for our 250 MHz quartz SH-SAW sensors was calculated to be 3.7×10−3m2sKg. The corresponding mass sensitivity was determined to be 9.25 × 105m2Kg. The limit of detection was calculated to be 36 picograms (pg), while the limit of quantification or LOQ was calculated to be 109 pg. Traditionally, liquid phase measurements have been challenging for SAW devices because liquids dampen the vibrating sensors severely. This problem has been largely solved using a transverse (shear) wave instead of the more popular longitudinal or Rayleigh waves. Liquid measurements are now possible using transverse waves, also known as shear waves, because transverse waves are only minimally attenuated by liquids. Shear-horizontal SAW sensors (SH-SAW) show great promise as label-free biosensors because of their ability to handle liquid samples. However, the viscosity of the liquid still induces loading effects and can be measured when the liquid is loaded onto the SH-SAW propagating surface (delay line). When the liquid above the delay line is perturbed by physical or chemical changes, such as binding to a receptor, it alters the propagating acoustic wave. The SH-SAW device can measure these changes in liquid properties as a change in the wave’s phase compared to the original wave. The device’s phase shift was recorded as a function of the changes in the density and viscosity of the binary glycerol solution and used to determine the sensitivity in the linear dynamic range of responses.

Materials Science↗

Solving Stochastic Inverse Problems for Property–Structure Linkages Using Data-Consistent Inversion and Machine Learning

Determining process–structure–property linkages is one of the key objectives in material science, and uncertainty quantification plays a critical role in understanding both process–structure and structure–property linkages. In this work, we seek to learn a distribution of microstructure parameters that are consistent in the sense that the forward propagation of this distribution through a crystal plasticity finite element model matches a target distribution on materials properties. This stochastic inversion formulation infers a distribution of acceptable/consistent microstructures, as opposed to a deterministic solution, which expands the range of feasible designs in a probabilistic manner. Furthermore, to solve this stochastic inverse problem, we employ a recently developed uncertainty quantification framework based on push-forward probability measures, which combines techniques from measure theory and Bayes’ rule to define a unique and numerically stable solution. This approach requires making an initial prediction using an initial guess for the distribution on model inputs and solving a stochastic forward problem. To reduce the computational burden in solving both stochastic forward and stochastic inverse problems, we combine this approach with a machine learning Bayesian regression model based on Gaussian processes and demonstrate the proposed methodology on two representative case studies in structure–property linkages.

36 MATERIALS SCIENCE↗

Predictive Dirac neutrino spectrum with strong CP solution in SU(5)L × SU(5)R unification

Abstract We develop a grand unified theory of matter and forces based on the gauge symmetry SU(5) L × SU(5) R with parity interchanging the two factor groups. Our main motivation for such a construction is to realize a minimal GUT embedding of left-right symmetric models that provide a parity solution to the strong CP problem without the axion. We show how the gauge couplings unify with an intermediate gauge symmetry SU(3) cL × SU(2) 2L × U(1) L × SU(5) R , and establish its consistency with proton decay constraints. The model correctly reproduces the observed fermion masses and mixings and leads tonaturally light Dirac neutrinoswith their Yukawa couplings suppressed by a factorM I /M G , the ratio of the intermediate scale to the GUT scale. We call this mechanism type II-Dirac seesaw. Furthermore, the model predictsδ CP = ±(130.4±1.2)° and$${m}_{{\nu }_{1}}$$= (4.8 – 8.4) meV for the Dirac CP phase and the lightest neutrino mass. We demonstrate how the model solves the strong CP problem via parity symmetry.

Physics↗

Design of super-strong and thermally stable nanotwinned Al alloys via solute synergy

Al alloys have widespread industrial applications. However, their mechanical strength is often much lower than steels. Here, we investigate the influence of solutes on achieving ultrahigh strength and thermal stability of nanotwinned Al alloys. In situ micropillar compression tests show the addition of a small amount of Ti can significantly increase the mechanical strength of Al–Ni alloys to 2 GPa. Deformation induced detwinning, Ni segregation and grain coarsening as discovered in binary Al–Ni alloys are mostly absent in the ternary Al–Ni–Ti alloys. Moreover, the ternary Al–Ni–Ti alloys have outstanding thermal stability. Density function theory calculations reveal the synergetic pinning effect of Ni–Ti solute pairs on incoherent twin boundaries. This study demonstrates that the proper selection of synergistic solute pairs is critical to improve the thermal stability and mechanical properties of nanotwinned Al alloys.

36 MATERIALS SCIENCE↗

The role of excess vacancies in stabilizing solute clusters in low-alloy steels

The formation of solute clusters in irradiated low-alloy steels, such as reactor pressure vessel steels, is a critical cause of radiation hardening and embrittlement. However, the chemical interactions among solute elements and excess vacancies remain to be fully understood. This study employs density functional theory, cluster expansion, and lattice-based Monte Carlo simulations to examine the stability and morphology of nano-size coherent solute clusters with excess vacancies. The findings reveal that excess vacancies are crucial in stabilizing and promoting the growth of Mn-Ni-Si-vacancy clusters. A minimum of seven vacancies is required for stable nucleation and growth of these clusters. These Mn-Ni-Si clusters act as defect sinks, effectively trapping and absorbing mobile vacancies generated under irradiation. Furthermore, phosphorus (P) preferentially dissolves in Mn-Ni-Si clusters due to chemical coupling with vacancies. In conclusion, this study offers new insights into solute-vacancy interactions, enhancing our understanding of solute-defect cluster formation and embrittlement in low-alloy steels.

36 - MATERIALS SCIENCE↗

Direct Extraction of Sodium Hydroxide by Calix[4]pyrrole-Based Ion-Pair Receptors

Here, described in this work, are calix[4]pyrrole-based ion-pair receptors, cis/trans-1 and cis/trans-2, designed for the extraction of sodium hydroxide. An X-ray diffraction analysis of a single crystal of the cis-1·NaOH isomer isolated from a mixture of cis/trans-1 revealed a unique dimeric supramolecular structure. An average dimer in toluene-d 8 solution was inferred on the basis of diffusion-ordered spectroscopy (DOSY). Support for the proposed stoichiometry came from density functional theory (DFT) calculations. The structural stability of the dimeric cis-1·NaOH complex in toluene solution was further confirmed by ab initio molecular dynamics (AIMD) simulation with explicit representation of solvent. Under conditions of liquid–liquid extraction (LLE), purified receptors cis- and trans-2 were both found to remove NaOH from a pH 11.01 aqueous source phase into toluene with extraction efficiencies (E%) of 50–60% when used equimolar to NaOH. However, in all cases, precipitation was observed. Complexities associated with precipitation could be avoided by immobilization of the receptors onto a chemically inert poly(styrene) resin by means of solvent impregnation. The use of solvent-impregnated resins (SIRs) eliminated precipitation in solution while retaining the extraction efficiency toward NaOH. This allowed both the pH and salinity of the alkaline source phase to be lowered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗

Nuclear quantum effects on the quasiparticle properties of the chloride anion aqueous solution within the GW approximation

Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral abinitio molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared with the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the Cl – anion and water. Furthermore, our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗