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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-D-17A-03: 5D Respiratory Motion Model Based Iterative Reconstruction Method for 4D Cone-Beam CT

Purpose: The purpose of this work is to develop a new iterative reconstruction method for 4D cone-beam CT (CBCT) based on a published time-independent 5D respiratory motion model. The proposed method will offer a single high-resolution image at a user-selected breathing phase and the 5D motion model parameters, which could be used to generate the breathing pattern during the CT acquisition. Methods: 5D respiratory motion model was proposed for accurately modeling the motion of lung and lung tumor tissues. 4D images are then parameterized by a reference image, measured breathing amplitude, breathing rate, two time-independent vector fields that describe the 5D model parameters, and a scalar field that describes the change in HU as a function of breathing amplitude. In contrast with the traditional method of reconstructing multiple temporal image phases to reduce respiratory artifact, 5D model based method simplify the problem into the reconstruction of a single reference image and the 5D motion model parameters. The reconstruction formulation of the reference image and scalar and vector fields is a nonlinear least-square optimization problem that consists of solving the reference image and fields alternately, in which the reference image is regularized with the total variation sparsity transform and the vector fields are solved through linearizations regularized by the H1 norm. 2D lung simulations were performed in this proof-of-concept study. Results: The breathing amplitude, its rate, and the corresponding scalar and vector fields were generated from a patient case. Compared with filtered backprojection method and sparsity regularized iterative method for the phase-by-phase reconstruction, the proposed 5D motion model based method yielded improved image quality. Conclusion: Based on 5D respiratory motion model, we have developed a new iterative reconstruction method for 4D CBCT that has the potential for improving image quality while providing needed on-table motion information.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Simulation Framework for Regional Geologic CO 2 Storage Along Arches Province of Midwestern United States

This report presents final technical results for the project Simulation Framework for Regional Geologic CO 2 Storage Infrastructure along Arches Province of the Midwest United States . The Arches Simulation project was a three year effort designed to develop a simulation framework for regional geologic carbon dioxide (CO 2 ) storage infrastructure along the Arches Province through development of a geologic model and advanced reservoir simulations of large-scale CO 2 storage. The project included five major technical tasks: (1) compilation of geologic, hydraulic and injection data on Mount Simon, (2) development of model framework and parameters, (3) preliminary variable density flow simulations, (4) multi-phase model runs of regional storage scenarios, and (5) implications for regional storage feasibility. The Arches Province is an informal region in northeastern Indiana, northern Kentucky, western Ohio, and southern Michigan where sedimentary rock formations form broad arch and platform structures. In the province, the Mount Simon sandstone is an appealing deep saline formation for CO 2 storage because of the intersection of reservoir thickness and permeability. Many CO 2 sources are located in proximity to the Arches Province, and the area is adjacent to coal fired power plants along the Ohio River Valley corridor. Geophysical well logs, rock samples, drilling logs, and geotechnical tests were evaluated for a 500,000 km2 study area centered on the Arches Province. Hydraulic parameters and historical operational information was also compiled from Mount Simon wastewater injection wells in the region. This information was integrated into a geocellular model that depicts the parameters and conditions in a numerical array. The geologic and hydraulic data were integrated into a three-dimensional grid of porosity and permeability, which are key parameters regarding fluid flow and pressure buildup due to CO 2 injection. Permeability data were corrected in locations where reservoir tests have been performed in Mount Simon injection wells. The geocellular model was used to develop a series of numerical simulations designed to support CO 2 storage applications in the Arches Province. Variable density fluid flow simulations were initially run to evaluate model sensitivity to input parameters. Two dimensional, multiple-phase simulations were completed to evaluate issues related to arranging injection fields in the study area. A basin-scale, multiple-phase model was developed to evaluate large scale injection effects across the region. Finally, local scale simulations were also completed with more detailed depiction of the Eau Claire formation to investigate to the potential for upward migration of CO 2 . Overall, the technical work on the project concluded that injection large-scale injection may be achieved with proper field design, operation, siting, and monitoring. Records from Mount Simon injection wells were compiled, documenting more than 20 billion gallons of injection into the Mount Simon formation in the Arches Province over the past 40 years, equivalent to approximately 60 million metric tons CO 2 . The multi-state team effort was useful in delineating the geographic variability in the Mount Simon reservoir properties. Simulations better defined potential well fields, well field arrangement, CO 2 pipeline distribution system, and operational parameters for large-scale injection in the Arches Province. Multiphase scoping level simulations suggest that injection fields with arrays of 9 to 50+ wells may be used to accommodate large injection volumes. Individual wells may need to be separated by 3 to 10 km. Injection fields may require spacing of 25 to 40 km to limit pressure and saturation front interference. Basin-scale multiple-phase simulations in STOMP reflect variability in the Mount Simon. While simulations suggest a total injection rate of 100 million metric tons per year (approximately to a 40% reduction of CO 2 emissions from large point sources across the Arches Province) may be feasible, some areas are more suitable due to favorable geology. Sustainable injection rates were higher in areas where there was higher thickness and reservoir permeability. Distribution of injection across well fields is an effective method to reduce pressure and CO 2 saturation front interference.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Universe acceleration and fine structure constant variation in BSBM theory

In this work we investigate the utility of using SNe Ia observations in constraining the cosmological parameters in BSBM theory where a scalar field is responsible for both fine structure constant variation and late time universe acceleration. The model is discussed in the presence of an exponential self potential for the scalar field. Stability and phase space analysis of the solutions are studied. The model is tested against observational data for Hubble parameter and quasar absorption spectra. With the best fitted model parameters, the theory predicts a good match with the experimental results and exhibits fine structure constant variation. The analysis also predicts the recent universe acceleration and possible phantom crossing in future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Phantom behavior bounce with tachyon and non-minimal derivative coupling

The bouncing cosmology provides a successful solution of the cosmological singularity problem. In this paper, we study the bouncing behavior of a single scalar field model with tachyon field non-minimally coupled to itself, its derivative and to the curvature. By utilizing the numerical calculations we will show that the bouncing solution can appear in the universe dominated by such a quintom matter with equation of state crossing the phantom divide line. We also investigate the classical stability of our model using the phase velocity of the homogeneous perturbations of the tachyon scalar field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Regularized Finite Element Formulations for Shear Band Instabilities in Metals

In my research I have addressed several fundamental issues of dynamic fracture: 1. established the role of heat conduction as a natural regularizing mechanism in shearband problems (i.e. using the correct physics for this problem), 2. developed a unified model to incorporate fracture through a modified phase field formulation, 3. developed stability analysis tools to determine the onset of localization, 4. developed discretization methods and solvers to solve this problem (on parallel computers) that lead to mesh insensitive results and fast convergence. For such difficult problems, the physics and numerics are strongly coupled. To accurately model the physics, one must develop reliable numerics but reliable numerics depends on correct description of the physics. To this end, we have already made significant progress in the field showing that a fully-coupled implicit solution with heat conduction is the key to obtain reliable results for shear bands (see the next sections for more details). This research fits well within the mission of DOE ASCR program on all aspects: from governing equations to discretization to solvers to parallel computers to applications problems. It also relies on the computational tools developed through the SciDAC institutes FASTMath and SDAV, such as PETSc (and SLEPc), ParaView, Trilinos and Masquite.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Two- and three-color two-quark states in warm, dense quark matter

We study the Gaussian fluctuations of the two-flavor, meson-diquark bosonized NJL model for two and three colors at the color superconducting phase transition. The analysis is based on analytical properties of the polarization matrix. Pionic fluctuations are shown to be stabilised in the 2SC phase compared to the two-color result where they are right on threshold.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Transformer quantum state: A multipurpose model for quantum many-body problems

Here, inspired by the advancements in large language models based on transformers, we introduce the transformer quantum state (TQS): a versatile machine learning model for quantum many-body problems. In sharp contrast to Hamiltonian/task specific models, TQS can generate the entire phase diagram, predict field strengths with experimental measurements, and transfer such a knowledge to new systems it has never been trained on before, all within a single model. With specific tasks, fine-tuning the TQS produces accurate results with small computational cost. Versatile by design, TQS can be easily adapted to new tasks, thereby pointing towards a general-purpose model for various challenging quantum problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

A solvable problem in statistical mechanics: The dipole-type Hamiltonian mean field model

The present study documents a type of mean field approximation inspired by the dipole interaction model, which is analytically solved in the canonical and microcanonical ensembles. The current calculations were derived from the Hamiltonian mean field model developments. After describing the canonical partition function, the free and internal energies, magnetization, and specific heat are derived and graphically depicted. In the microcanonical ensemble, the entropy is calculated as well as other thermodynamic properties. The system shows a second-order phase transition emphasizing that both methods coincide, which is only valid only in equilibrium. In addition, the current model represents a nonsymmetric Hamiltonian mean field model that shows a phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

SIMULATION STUDY OF HEMISPHERIC PHASE-ASYMMETRY IN THE SOLAR CYCLE

Observations of the Sun suggest that solar activities systematically create north–south hemispheric asymmetries. For instance, the hemisphere in which sunspot activity is more active tends to switch after the early half of each solar cycle. Svalgaard and Kamide recently pointed out that the time gaps of polar field reversal between the northern and southern hemispheres are simply consequences of the asymmetry of sunspot activity. However, the mechanism underlying the asymmetric feature in solar cycle activity is not yet well understood. In this paper, in order to explain the cause of the asymmetry from the theoretical point of view, we investigate the relationship between the dipole- and quadrupole-type components of the magnetic field in the solar cycle using the mean-field theory based on the flux transport dynamo model. As a result, we found that there are two different attractors of the solar cycle, in which either the north or the south polar field is first reversed, and that the flux transport dynamo model explains well the phase-asymmetry of sunspot activity and the polar field reversal without any ad hoc source of asymmetry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Dependence of neoclassical toroidal viscosity on the poloidal spectrum of applied nonaxisymmetric fields

This paper presents a single mode model that accurately predicts the coupling of applied nonaxisymmetric fields to the plasma response that induces neoclassical toroidal viscosity (NTV) torque in DIII-D H-mode plasmas. The torque is measured and modeled to have a sinusoidal dependence on the relative phase of multiple nonaxisymmetric field sources, including a minimum in which large amounts of nonaxisymmetric drive is decoupled from the NTV torque. This corresponds to the coupling and decoupling of the applied field to a NTV-driving mode spectrum. Modeling using the perturbed equilibrium nonambipolar transport (PENT) code confirms an effective single mode coupling between the applied field and the resultant torque, despite its inherent nonlinearity. Lastly, the coupling to the NTV mode is shown to have a similar dependence on the relative phasing as that of the IPEC dominant mode, providing a physical basis for the efficacy of this linear metric in predicting error field correction optima in NTV dominated regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Statistics of particle time-temperature histories

Particles in non-isothermal turbulent flow are subject to a stochastic environment that produces a distribution of particle time-temperature histories. This distribution is a function of the dispersion of the non-isothermal (continuous) gas phase and the distribution of particles relative to that gas phase. In this work we extend the one-dimensional turbulence (ODT) model to predict the joint dispersion of a dispersed particle phase and a continuous phase. The ODT model predicts the turbulent evolution of continuous scalar fields with a model for the cascade of fluctuations to smaller scales (the ‘triplet map’) at a rate that is a function of the fully resolved one-dimensional velocity field. Stochastic triplet maps also drive Lagrangian particle dispersion with finite Stokes numbers including inertial and eddy trajectory-crossing effects included. Two distinct approaches to this coupling between triplet maps and particle dispersion are developed and implemented along with a hybrid approach. An ‘instantaneous’ particle displacement model matches the tracer particle limit and provides an accurate description of particle dispersion. A ‘continuous’ particle displacement model translates triplet maps into a continuous velocity field to which particles respond. Particles can alter the turbulence, and modifications to the stochastic rate expression are developed for two-way coupling between particles and the continuous phase. Each aspect of model development is evaluated in canonical flows (homogeneous turbulence, free-shear flows and wall-bounded flows) for which quality measurements are available. ODT simulations of non-isothermal flows provide statistics for particle heating. These simulations show the significance of accurately predicting the joint statistics of particle and fluid dispersion. Inhomogeneous turbulence coupled with the influence of the mean flow fields on particles of varying properties alters particle dispersion. The joint particle-temperature dispersion leads to a distribution of temperature histories predicted by the ODT. Predictions are shown for the lower moments and the full distributions of the particle positions, particle-observed gas temperatures and particle temperatures. An analysis of the time scales affecting particle-temperature interactions covers Lagrangian integral time scales based on temperature autocorrelations, rates of temperature change associated with particle motion relative to the temperature field and rates of diffusional change of temperatures. These latter two time scales have not been investigated previously; they are shown to be strongly intermittent having peaked distributions with long tails. The logarithm of the absolute value of these time scales exhibits a distribution closer to normal.

Hewson, John C.↗
DOE OSTI2015

Luttinger liquid behavior in the alternating spin-chain system copper nitrate

We determine the phase diagram of copper nitrate Cu(NO 3 ) 2 ∙2.5D 2 O in the context of quantum phase transitions and novel states of matter. We establish this compound as an ideal candidate to study quasi-1D Luttinger liquids, 3D Bose-Einstein-Condensation of triplons, and the crossover between 1D and 3D physics. Magnetocaloric effect, magnetization, and neutron scattering data provide clear evidence for transitions into a Luttinger liquid regime and a 3D long-range ordered phase as a function of field and temperature. Theoretical simulations of this model material allow us to fully establish the phase diagram and to discuss it in the context of dimerized spin systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Unit-density SU(3) Fermi-Hubbard model with spin-flavor imbalance

The advent of ultracold alkaline-earth atoms in optical lattices has established a platform for investigating correlated quantum matter with SU⁡(𝑁) symmetry, offering highly tunable model parameters that allow experiments to access phenomena that are unavailable in conventional materials. Understanding the ground-state physics of SU⁡(𝑁) Fermi-Hubbard models away from the Heisenberg limit and from the spin-flavor balanced setting is important, as examining the flavor imbalance reveals new physics in Fermi-Hubbard models and shows how SU⁡(𝑁) phases react to practical experimental imperfections in optical lattices. Here, in this study, mean-field phase diagrams are presented for the unit-density SU(3) Fermi-Hubbard model at two sets of flavor densities, ($\frac{1}{3}$ − 𝛿, $\frac{1}{3}$ + 𝛿, $\frac{1}{3}$) and ($\frac{1}{4}$ − 𝛿, $\frac{1}{4}$ +𝛿, $\frac{1}{2}$), with the flavor imbalance introduced as 𝛿. Novel phases are identified at moderate interaction strengths for both densities, and their robustness is investigated in the presence of flavor imbalance. Furthermore, we provide microscopic explanations of the phases found and their stability. Analysis of thermal ensembles of random mean-field solutions indicate that, at temperatures accessible in state-of-the-art cold atom experiments, some spin orders are hard for conventional scattering or local observable measurements to detect, but can be more accessible with quantum gas microscopy in optical lattice experiments. This work also shows that nesting and Mottness, intertwined in the usual SU(2) Hubbard model in stark contrast to generic materials, can be tuned in the SU(3) model and play distinct roles. The resulting phase diagrams not only deepen our understanding of SU⁡(𝑁) Fermi-Hubbard models but also inform future experimental search for new phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

Thin current sheets in the magnetotail during substorms: CDAW 6 revisited

The global magnetic field configuration during the growth phase of the Coordinated Data Analysis Workshop (CDAW) 6 substorm (March 22, 1979, 1054 UT) is modeled using data from two suitably located spacecraft and temporally evolving variations of the Tsyganenko magnetic field model. These results are compared with a local calculation of the current sheet location and thickness carried out by McPherron et al. (1987) and Sanny et al. (this issue). Both models suggest that during the growth phase the current sheet rotated away from its nominal location, and simultaneously thinned strongly. The locations and thickness obtained from the two models are in good agreement. The global model suggests that the peak current density is approximately 120 nA/sq m and that the cross-tail current almost doubled its intensity during this very strong growth phase. The global model predicts a field configuration that is sufficiently stretched to scatter thermal electrons, which may be conducive to the onset of ion tearing in the tail. The electron plasma data further support this scenario, as the anisotropy present in the low-energy electrons disappears close to the substorm onset. The electron contribution to the intensifying current in this case is of the order of 10% before the isotropization of the distribution.

Pulkkinen, T. I.↗
DOE OSTI2024

Cohesive phase-field chemo-mechanical simulations of inter- and trans- granular fractures in polycrystalline NMC cathodes via image-based 3D reconstruction

The optimal design and durable utilization of lithium-ion batteries necessitates an objective modeling approach to understand fracture and failure mechanisms. This paper presents a comprehensive chemo-mechanical modeling study focused on elucidating fracture-induced damage and degradation phenomena in the polycrystalline Li $\mathcal{x}$ Ni 0.5 Mn 0.3 Co 0.2 O 2 (NMC532) cathode. An innovative approach that utilizes image-based reconstructed 3D geometry as finite element (FE) mesh input is employed to enhance the precision in capturing the convoluted architecture and morphological features. For accurately representing the intricate crack configurations within the polycrystalline system, we adopted the cohesive phase-field fracture (CPF) model. Through the integration of advanced image-based geometry reconstruction technique and the promising CPF modeling approach, lithium (de)intercalation induced crack evolution (e.g., nucleation, propagation, branching and diverse modes including inter-/trans-(intra-) granular patterns) and the resulting chemical degradation can be precisely captured, which is also compared and validated with numerical predictions using a continuum damage model. In particular, this model predicts fracture induced degradation under varying fracture properties of grain boundaries and charging rates; the conclusion that NMC particles comprised of larger grains are predicted to have less degradation than those with smaller grains can also be drawn. This comprehensive analysis provides valuable insights into the fracture and degradation within polycrystalline NMC cathodes.

25 ENERGY STORAGE↗
NASA NTRS1983

Spinning mode acoustic radiation from the flight inlet

A mathematical model was developed for spinning mode acoustic radiation from a thick wall duct without flow. This model is based on a series of experiments (with and without flow). A nearly pure azimuthal spinning mode was isolated and then reflection coefficients and far field pressure (amplitude and phase) were measured. In our model the governing boundary value problem for the Helmholtz equation is first converted into an integral equation for the unknown acoustic pressure over a disk, S1, near the mouth of the duct and over the exterior surface, S2, of the duct. Assuming a pure azimuthal mode excitation, the azimuthal dependence is integrated out which yields an integral equation over the generator C1 of S1 and the generator C2 of S2. The sound pressure on C1 was approximated by a truncated modal expansion of the interior acoustic pressure. Piecewise linear spline approximation on C2 was used.

Moss, W. F.↗
DOE OSTI2014

Model-based feasibility assessment and evaluation of prostate hyperthermia with a commercial MR-guided endorectal HIFU ablation array

Purpose: Feasibility of targeted and volumetric hyperthermia (40–45 °C) delivery to the prostate with a commercial MR-guided endorectal ultrasound phased array system, designed specifically for thermal ablation and approved for ablation trials (ExAblate 2100, Insightec Ltd.), was assessed through computer simulations and tissue-equivalent phantom experiments with the intention of fast clinical translation for targeted hyperthermia in conjunction with radiotherapy and chemotherapy. Methods: The simulations included a 3D finite element method based biothermal model, and acoustic field calculations for the ExAblate ERUS phased array (2.3 MHz, 2.3 × 4.0 cm{sup 2}, ∼1000 channels) using the rectangular radiator method. Array beamforming strategies were investigated to deliver protracted, continuous-wave hyperthermia to focal prostate cancer targets identified from representative patient cases. Constraints on power densities, sonication durations and switching speeds imposed by ExAblate hardware and software were incorporated in the models. Preliminary experiments included beamformed sonications in tissue mimicking phantoms under MR temperature monitoring at 3 T (GE Discovery MR750W). Results: Acoustic intensities considered during simulation were limited to ensure mild hyperthermia (T{sub max} < 45 °C) and fail-safe operation of the ExAblate array (spatial and time averaged acoustic intensity I{sub SATA} < 3.4 W/cm{sup 2}). Tissue volumes with therapeutic temperature levels (T > 41 °C) were estimated. Numerical simulations indicated that T > 41 °C was calculated in 13–23 cm{sup 3} volumes for sonications with planar or diverging beam patterns at 0.9–1.2 W/cm{sup 2}, in 4.5–5.8 cm{sup 3} volumes for simultaneous multipoint focus beam patterns at ∼0.7 W/cm{sup 2}, and in ∼6.0 cm{sup 3} for curvilinear (cylindrical) beam patterns at 0.75 W/cm{sup 2}. Focused heating patterns may be practical for treating focal disease in a single posterior quadrant of the prostate and diffused heating patterns may be useful for heating quadrants, hemigland volumes or even bilateral targets. Treatable volumes may be limited by pubic bone heating. Therapeutic temperatures were estimated for a range of physiological parameters, sonication duty cycles and rectal cooling. Hyperthermia specific phasing patterns were implemented on the ExAblate prostate array and continuous-wave sonications (∼0.88 W/cm{sup 2}, 15 min) were performed in tissue-mimicking material with real-time MR-based temperature imaging (PRFS imaging at 3.0 T). Shapes of heating patterns observed during experiments were consistent with simulations. Conclusions: The ExAblate 2100, designed specifically for thermal ablation, can be controlled for delivering continuous hyperthermia in prostate while working within operational constraints.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2010

BEC-BCS crossover driven by the axial anomaly in the NJL model

We study the QCD phase structure in the three-flavor Nambu-Jona-Lasinio model, incorporating the chiral-diquark interplay due to the axial anomaly. We demonstrate that for a certain range of model parameters, the low temperature critical point predicted by a Ginzburg-Landau analysis appears in the phase diagram. In addition, we show that the axial anomaly presents a new scenario for a possible BEC-BCS crossover in the color-flavor locked phase of QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗