Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field 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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Preliminary calibration plan for the Advanced Particles and Field Observatory (APAFO) magnetometer experiment

Prelaunch and postlaunch calibration plans for the APAFO magnetometer experiment are presented. A study of tradeoffs between boom length and spacecraft field is described; the results are summarized. The prelaunch plan includes: calibration of the Vector Fluxgate Magnetometer (VFM), Star Sensors, and Scalar Helium Magnetometer (SHM); optical bench integration; and acquisition of basic spacecraft field data. Postlaunch calibration has two phases. In phase one, SHM data are used to calibrate the VFM, total vector magnetic field data are used to calibrate a physical model of the spacecraft field, and both calibrations are refined by iteration. In phase two, corrected vector data are transformed into geocentric coordinates, previously undetected spacecraft fields are isolated, and initial geomagnetic field models are computed. Provided the SHM is accurate to the required 1.0 nT and can be used to calibrate the VFM to the required 3.0- nT accuracy, the tradeoff study indicates that a 12 m boom and a spacecraft field model uncertainty of 5 percent together allow the 1.0 nT spacecraft field error requirement to be met.

Voorhies, C. V.↗
DOE OSTI2019

Field-induced quantum spin liquid in the Kitaev-Heisenberg model and its relation to α – RuCl 3

Recently considerable excitement has arisen due to the experimental observation of a field-induced spin liquid phase in the compound α– RuCl 3 . Yet, the nature of this putative spin liquid phase and the relevant microscopic model Hamiltonian remain still unclear. In this work, we address these questions by performing large-scale numerical simulations of a generalized Kitaev-Heisenberg model proposed to describe the physics of α– RuCl 3 . Although an intermediate phase does not appear for in-plane magnetic fields, our results strongly suggest that a stable intermediate spin liquid phase, sandwiched between a magnetically ordered phase at low fields and a high-field polarized phase, can be induced by out-of-plane magnetic fields. Moreover, we show that this field-induced spin liquid phase can be smoothly connected to a spin liquid possessing a spinon Fermi surface as proposed recently for the Kitaev model. The relevance of our findings to α– RuCl 3 is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2012

Broadband Spectroscopy Using Two Suzaku Observations of the HMXB GX 301-2

We present the analysis of two Suzaku observations of GX 301-2 at two orbital phases after the periastron passage. Variations in the column density of the line-of-sight absorber are observed, consistent with accretion from a clumpy wind. In addition to a CRSF, multiple fluorescence emission lines were detected in both observations. The variations in the pulse profiles and the CRSF throughout the pulse phase have a signature of a magnetic dipole field. Using a simple dipole model we calculated the expected magnetic field values for different pulse phases and were able to extract a set of geometrical angles, loosely constraining the dipole geometry in the neutron star. From the variation of the CRSF width and energy, we found a geometrical solution for the dipole, making the inclination consistent with previously published values.

Suchy, Slawomir↗
NASA NTRS1991

Ordering and phase transitions in random-field Ising systems

An exact analysis of the Ising model with infinite-range interactions in a random field and a local mean-field theory in three dimensions is carried out leading to a phase diagram with several coexistence surfaces and lines of critical points. The results show that the phase diagram depends crucially on whether the distribution of random fields is symmetric or not. Thus, Ising-like phase transitions in a porous medium (the asymmetric case) are in a different universality class from the conventional random-field model (symmetric case).

Maritan, Amos↗
NASA NTRS1989

Novel local symmetries and chiral-symmetry-broken phases in S = 1/2 triangular-lattice Heisenberg model

Using a nonmean-field approach the triangular-lattice S = 1/2 Heisenberg antiferromagnet with nearest- and next-nearest-neighbor couplings is shown undergo an Ising-type phase transition into a chiral-symmetry-broken phase (Kalmeyer-Laughlin-like state) at small T. Removal of next-nearest-neighbor coupling introduces a local Z2 symmetry, thereby suppressing any finite-T chiral order.

Baskaran, G.↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Theory of the Kitaev model in a [111] magnetic field

Recent numerical studies indicate that the antiferromagnetic Kitaev honeycomb lattice model undergoes a magnetic-field-induced quantum phase transition into a new spin-liquid phase. This intermediate-field phase has been previously characterized as a gapless spin liquid. By implementing a recently developed variational approach based on the exact fractionalized excitations of the zero-field model, we demonstrate that the field-induced spin liquid is gapped and belongs to Kitaev’s 16-fold way. Specifically, the low-field non-Abelian liquid with Chern number C = ±1 transitions into an Abelian liquid with C = ±4. The critical field and the field-dependent behaviors of key physical quantities are in good quantitative agreement with published numerical results. Furthermore, we derive an effective field theory for the field-induced critical point which readily explains the ostensibly gapless nature of the intermediate-field spin liquid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2023

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗
DOE OSTI2020

Grand potential sintering simulations of doped UO 2 accident-tolerant fuel concepts

Chromium-doped UO is a widely-studied near-term deployable accident-tolerant fuel concept because it results in a dense, large-grain structure that increases the fuel resistance to densification, swelling, and fission gas release. A new charged-interstitial mechanism was recently proposed to describe the behavior of dopants like chromium in sintered UO 2 . Based on that mechanism, manganese was suggested as an even stronger dopant than chromium. In the current work we use mesoscale sintering simulations in an effort to validate the new mechanism. We compare the relative behavior of Cr-doped and undoped UO 2 against experimental data in the literature. We also make predictions of the relative behavior of Mn-doped UO 2 . Simulations are done using the phase field-based grand potential sintering model. Dopants have two effects on sintered UO 2 . They marginally increase the densification rate and greatly increase the average grain size. Both of these effects are individually tested in small-scale simulations. Then large 3D sintering simulations are used to test the combined behavior of both effects. The results for the densification rate simulations are consistent with experiments for Cr-doped fuels. However, we note the grain growth rates are lower than what is found in the literature for Cr-doped fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling impacts of ice-nucleating particles from marine aerosols on mixed-phase orographic clouds during 2015 ACAPEX field campaign

Abstract. A large fraction of annual precipitation over the western United States comes from wintertime orographic clouds associated with atmospheric rivers (ARs). Transported African and Asian dust and marine aerosols from the Pacific Ocean may act as ice-nucleating particles (INPs) to affect cloud and precipitation properties over the region. Here we explored the effects of INPs from marine aerosols on orographic mixed-phase clouds and precipitation at different AR stages for an AR event observed during the 2015 ACAPEX field campaign under low dust (<0.02 cm−3) conditions. Simulations were conducted using the chemistry version of the Weather Research and Forecasting Model coupled with the spectral-bin microphysics at 1 km grid spacing, with ice nucleation connected with dust and marine aerosols. By comparing against airborne and ground-based observations, accounting for marine INP effects improves the simulation of AR-precipitation. The marine INPs enhance the formation of ice and snow, leading to less shallow warm clouds but more mixed-phase and deep clouds, as well as to a large spillover effect of precipitation after AR landfall. The responses of cloud and precipitation to marine INPs vary with the AR stages, with more significant effects before AR landfall and post-AR than after AR landfall, mainly because the moisture and temperature conditions change with the AR evolution. This work suggests weather and climate models need to consider the impacts of marine INPs since their contribution is notable under low dust conditions despite the much lower relative ice nucleation efficiency of marine INPs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2019

Phase space quantum–classical hybrid model

Highlights: •Semiclassical model including back-action reaction beyond mean-field approximation. •Quantum–classical hybridization holds for all times as well as all intensities. •Wigner–Weyl phase-space correspondence mixing quantum with classical systems. -- Abstract: In this work we provide a complete model of semiclassical theories by including back-reaction and correlation into the picture. We specially aim at the interaction between light and a two-level atom, and we also illustrate it via the coupling of two harmonic oscillators. Quantum and classical systems are treated on the same grounds via the Wigner–Weyl phase-space correspondence of the quantum theory. We show that this model provides a suitable mixture of the quantum and classical degrees of freedom, including the fact that the evolution transfers nonclassical features to the classical subsystem and nonquantum behaviour to the quantum subsystem. In that sense, we can no longer distinguish between classical and quantum variables and we need to talk about a hybrid model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Charge Order of Strongly Bounded Electron Pairs on the Triangular Lattice: the Zero-Bandwidth Limit of the Extended Hubbard Model with Strong Onsite Attraction

The inhomogeneous distribution of charged particles is known as charge-order phenomena. Recent experimental studies of (quasi-) two-dimensional materials found different ordered patterns on the triangular lattice. In this work, the preliminary studies of the extended Hubbard model (with the intersite Coulomb density-density repulsion) in the atomic limit on the triangular lattice are presented. The model is investigated within the mean-field approach that treats the onsite term exactly. The finite temperature phase diagrams including metastable phases and some characteristics of the model are determined for the limit of strong onsite attraction. It is shown that in the limit considered all transitions between stable phases are discontinuous (for fixed chemical potential), which is totally different from the model on the alternate lattices. Continuous transitions between metastable phases are also found.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗
DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Simulating Heisenberg interactions in the Ising model with strong drive fields

The time evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. Finally, these approximate equivalences, valid above a critical driving field strength set by dynamical phase transitions in the Ising model, are expected to enable quantum devices that natively evolve qubits according to the Ising model to simulate more complex systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗