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 991 records · Page 55

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗

Overview of the Field Phase of the NASA Tropical Cloud Systems and Processes (TCSP)Experiment

The Tropical Cloud Systems and Processes experiment is sponsored by the National Aeronautics and Space Administration (NASA) to investigate characteristics of tropical cyclone genesis, rapid intensification and rainfall using a three-pronged approach that emphasizes satellite information, suborbital observations and numerical model simulations. Research goals include demonstration and assessment of new technology, improvements to numerical model parameterizations, and advancements in data assimilation techniques. The field phase of the experiment was based in Costa Rica during July 2005. A fully instrumented NASA ER-2 high altitude airplane was deployed with Doppler radar, passive microwave instrumentation, lightning and electric field sensors and an airborne simulator of visible and infrared satellite sensors. Other assets brought to TCSP were a low flying uninhabited aerial vehicle, and a surface-based radiosonde network. In partnership with the Intensity Forecasting Experiment of the National Oceanic and Atmospheric Administration (NOAA) Hurricane Research Division, two NOAA P-3 aircraft instrumented with radar, passive microwave, microphysical, and dropsonde instrumentation were also deployed to Costa Rica. The field phase of TCSP was conducted in Costa Rica to take advantage of the geographically compact tropical cyclone genesis region of the Eastern Pacific Ocean near Central America. However, the unusual 2005 hurricane season provided numerous opportunities to sample tropical cyclone development and intensification in the Caribbean Sea and Gulf of Mexico as well. Development of Hurricane Dennis and Tropical Storm Gert were each investigated over several days in addition to Hurricane Emily as it was close to Saffir-Simpson Category 5 intensity. An overview of the characteristics of these storms along with the pregenesis environment of Tropical Storm Eugene in the Eastern Pacific will be presented.

Hood, Robbie E.↗

Low-energy phenomenology of scalarless standard-model extensions with high-energy Lorentz violation

We consider renormalizable standard model extensions that violate Lorentz symmetry at high energies, but preserve CPT, and do not contain elementary scalar fields. A Nambu-Jona-Lasinio mechanism gives masses to fermions and gauge bosons and generates composite Higgs fields at low energies. We study the effective potential at the leading order of the large-N{sub c} expansion, prove that there exists a broken phase, and study the phase space. In general, the minimum may break invariance under boosts, rotations, and CPT, but we give evidence that there exists a Lorentz invariant phase. We study the spectrum of composite bosons and the low-energy theory in the Lorentz phase. Our approach predicts relations among the parameters of the low-energy theory. We find that such relations are compatible with the experimental data within theoretical errors. We also study the mixing among generations, the emergence of the CKM matrix, and neutrino oscillations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Numerical and experimental characterization of twin transmission across grain boundaries along the forward and lateral directions

Pervasive deformation twinning and transmission events across grain boundaries (GBs) affect the strength and failure of hexagonal close-packed (HCP) magnesium. A three-dimensional twin can transmit along the twinning shear direction, (forward), and along the direction perpendicular to both the twinning plane normal and the shear direction, (lateral). For the first time, phase-field calculations and electron backscatter diffraction (EBSD)-based statistical analysis are combined to investigate the effect of the twinned grain boundary characteristics on twin transmission (TT) along the forward and lateral directions. This combined analysis reveals that TT propensity decreases with increasing misorientation angle across the GB for both forward and lateral directions. Also, the TT is more favorable along the lateral than along the forward direction. Twin transmission seems harder across GBs with a misorientation axis closer to the twin -direction than the other directions ( and ). Further, the PF calculations reveal that, at the onset of a transmission process, the crystallography tends to be preserved in the case of lateral transmission, whereas, in the forward transmission case, the transmitted twin punches straight through the GBs and its morphology prevails. Additionally, the EBSD analysis finds that pure forward and lateral transmissions do not occur often, yet reveals a preference for lateral propagation consistent with PF simulations. Further, the local twin transmission configurations observed in the actual material do not correspond to pure tilt or twist GBs, which are most commonly considered as model cases.

36 MATERIALS SCIENCE↗

Simulations of phase space distributions of storm time proton ring current

We use results of guiding-center simulations of ion transport to map phase space densities of the stormtime proton ring current. We model a storm as a sequence of substorm-associated enhancements in the convection electric field. Our pre-storm phase space distribution is an analytical solution to a steady-state transport model in which quiet-time radial diffusion balances charge exchange. This pre-storm phase space spectra at L approximately 2 to 4 reproduce many of the features found in observed quiet-time spectra. Using results from simulations of ion transport during model storms having main phases of 3, 6, and 12 hr, we map phase space distributions from the pre-storm distribution in accordance with Liouville's theorem. We find stormtime enhancements in the phase space densities at energies E approximately 30-160 keV for L approximately 2.5 to 4. These enhancements agree well with the observed stormtime ring current. For storms with shorter main phases (approximately 3 hr), the enhancements are caused mainly by the trapping of ions injected from open night side trajectories, and diffusive transport of higher-energy (greater than or approximately 160 keV) ions contributes little to the stormtime ring current. However, the stormtime ring current is augmented also by the diffusive transport of higher-energy ions (E greater than or approximately 160 keV) durinng stroms having longer main phases (greater than or approximately 6 hr). In order to account for the increase in Dst associated with the formation of the stormtime ring current, we estimate the enhancement in particle-energy content that results from stormtime ion transport in the equatorial magnetosphere. We find that transport alone cannot account for the entire increase in absolute value of Dst typical of a major storm. However, we can account for the entire increase in absolute value of Dst by realistically increasing the stormtime outer boundary value of the phase space density relative to the quiet-time value. We compute the magnetic field produced by the ring current itself and find that radial profiles of the magnetic field depression resemble those obtained from observational data.

Chen, Margaret W.↗

Persistent ocean anomalies as a response to Northern Hemisphere heating induced by biomass burning variability

Biomass burning aerosol (BBA) emissions in the Coupled Model Intercomparison Project Phase 6 (CMIP6) historical forcing fields have enhanced temporal variability during the years 1997–2014 compared to earlier periods. Recent studies document that the corresponding inhomogeneous shortwave forcing over this period can cause changes in clouds, permafrost, and soil moisture, which contribute to a net terrestrial Northern Hemisphere warming relative to earlier periods. Here, we investigate the ocean response to the hemispherically asymmetric warming, using a 100-member ensemble of the Community Earth System Model version 2 Large Ensemble forced by two different BBA emissions (CMIP6 default and temporally smoothed over 1990–2020). Differences between the two subensemble means show that ocean temperature anomalies occur during periods of high BBA variability and subsequently persist over multiple decades. In the North Atlantic, surface warming is efficiently compensated for by decreased northward oceanic heat transport due to a slowdown of the Atlantic Meridional Overturning Circulation. In the North Pacific, surface warming is compensated for by an anomalous cross-equatorial cell (CEC) that reduces northward oceanic heat transport. The heat that converges in the South Pacific through the anomalous CEC is shunted into the subsurface and contributes to formation of long-lasting ocean temperature anomalies. The anomalous CEC is maintained through latitude-dependent contributions from narrow western boundary currents and basin-wide near-surface Ekman transport. Furthermore, these results indicate that interannual variability in forcing fields may significantly change the background climate state over long timescales, presenting a potential uncertainty in CMIP6-class climate projections forced without interannual variability.

54 ENVIRONMENTAL SCIENCES↗

Wavefront Sensing for WFIRST with a Linear Optical Model

In this paper we develop methods to use a linear optical model to capture the field dependence of wavefront aberrations in a nonlinear optimization-based phase retrieval algorithm for image-based wavefront sensing. The linear optical model is generated from a ray trace model of the system and allows the system state to be described in terms of mechanical alignment parameters rather than wavefront coefficients. This approach allows joint optimization over images taken at different field points and does not require separate convergence of phase retrieval at individual field points. Because the algorithm exploits field diversity, multiple defocused images per field point are not required for robustness. Furthermore, because it is possible to simultaneously fit images of many stars over the field, it is not necessary to use a fixed defocus to achieve adequate signal-to-noise ratio despite having images with high dynamic range. This allows high performance wavefront sensing using in-focus science data. We applied this technique in a simulation model based on the Wide Field Infrared Survey Telescope (WFIRST) Intermediate Design Reference Mission (IDRM) imager using a linear optical model with 25 field points. We demonstrate sub-thousandth-wave wavefront sensing accuracy in the presence of noise and moderate undersampling for both monochromatic and polychromatic images using 25 high-SNR target stars. Using these high-quality wavefront sensing results, we are able to generate upsampled point-spread functions (PSFs) and use them to determine PSF ellipticity to high accuracy in order to reduce the systematic impact of aberrations on the accuracy of galactic ellipticity determination for weak-lensing science.

Jurling, Alden S.↗

The QCD Phase Transition in Strong Magnetic Fields

We investigate the properties of the deconfining/chiral restoring transition for two flavor QCD in presence of a uniform background magnetic field. We adopt a standard staggered discretization of the fermion action, different values of the bare quark mass corresponding to m{sub {pi}} ranging from 200 to 500 MeV, and magnetic fields up to eB 1 GeV{sup 2}. We present first results regarding the dependence of the deconfinement and chiral transition temperature and strength on the magnetic field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fluctuations and the QCD phase diagram

In this contribution the role of quantum fluctuations for the QCD phase diagram is discussed. This concerns in particular the importance of the matter back-reaction to the gluonic sector. The impact of these fluctuations on the location of the confinement/deconfinement and the chiral transition lines as well as their interrelation are investigated. Consequences of our findings for the size of a possible quarkyonic phase and location of a critical endpoint in the phase diagram are drawn.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Progress on HL-LHC Nb 3 Sn Magnets

The high-luminosity Large Hadron Collider (HL-LHC) project aims at allowing to increase the collisions in the LHC by a factor of ten in the decade 2025-2035. One essential element is the superconducting magnet around the interaction region points, where the large aperture magnets will be installed to allow to further reduce the beam size in the interaction point. The core of this upgrade is the Nb 3 Sn triplet, made up of 150-mm aperture quadrupoles in the range of 7-8 m. The project is being shared between the European Organization for Nuclear Research and the US Accelerator Upgrade Program, based on the same design, and on the two strand technologies. The project is ending the short model phase, and entering the prototype construction. We will report on the main results of the short model program, including the quench performance and field quality. A second important element is the 11 T dipole that replaces a standard dipole making space for additional collimators. The magnet is also ending the model development and entering the prototype phase. A critical point in the design of this magnet is the large current density, allowing increase of the field from 8 to 11 T with the same coil cross section as in the LHC dipoles. This is also the first two-in-one Nb 3 Sn magnet developed so far. We will report the main results on the test and the critical aspects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Connecting Pressure-Saturation and Relative Permeability Models to Fracture Properties: The Case of Capillary-Dominated Flow of Supercritical CO 2 and Brine

Fractures are potential pathways for subsurface fluids. In the context of geologic CO 2 sequestration, fractures could threaten the security of reservoirs since they are pathways for leakage. However, the constitutive equations describing the pressure-saturation and relative permeability-saturation, key inputs for flow modeling and prediction, are poorly understood for two-phase fracture flow at the field scale where each fracture might be an element of a complex fracture network or dual-porosity domain. This knowledge is required for the safe storage of CO 2 under potentially fractured caprocks. To address this, we numerically simulated two-phase viscous flow through horizontal heterogeneous fractures across a broad range of roughness and aperture correlation length. The numerical experiments, which solved the modified local cubic law, showed a weak to strong phase interference during the drainage process; that is, supercritical CO 2 displaces brine if the imposed pressure + local viscous force > local capillary force. The drainage process ended when no new brine cells can be further invaded by CO 2 . Afterward, we froze the saturation field for each individual phase and calculated the relative permeability by single-phase flow modeling. Thousands of simulated pressure-saturation and relative permeability curves enabled us to estimate and connect the parameters of the Van-Genuchten model and of the generalized ν-type model to fracture roughness and aperture correlation length. In conclusion, this empirical connection will be useful for assessing and predicting immiscible two-phase flow in horizontal rough fractures at the large scale.

54 ENVIRONMENTAL SCIENCES↗

Interactions Enhance Ramp Reversal Memory in Locally Phase Separated Materials

The ramp-reversal memory (RRM) effect in metal–insulator transition metal oxides (TMOs), a non-volatile resistance change induced by repeated temperature cycling, has attracted considerable interest in neuromorphic computing and non-volatile memory devices. Our previous defect motion model successfully explained RRM in vanadium dioxide (VO 2 ), capturing observed critical temperature shifts and memory accumulation throughout the sample. However, this approach lacked interactions between metallic and insulating domains. Here, we extend our model by combining a correlated Random Field Ising Model with defect diffusion-segregation, enabling accurate hysteresis modeling while predicting the relationship between RRM and domain interactions. Our simulations demonstrate that the maximum RRM occurs when the turnaround temperature approaches the inflection point. This peak in RRM vs. turnaround temperature is consistent with prior transport measurements, as well as our own optical measurements reported here. Significantly, we find that increasing nearest-neighbor interactions enhances the maximum memory effect, thus providing a clear mechanism for optimizing RRM performance. Since our model employs minimal assumptions, we predict that RRM should be a widespread phenomenon in materials exhibiting patterned phase coexistence of electronic domains. This work not only advances fundamental understanding of memory behavior in TMOs but also establishes a much-needed theoretical framework for optimizing device applications.

36 MATERIALS SCIENCE↗

Evaluation of a cloudy cold-air pool in the Columbia River basin in different versions of the High-Resolution Rapid Refresh (HRRR) model

The accurate forecast of persistent orographic cold-air pools in numerical weather prediction models is essential for the optimal integration of wind energy into the electrical grid during these events. Model development efforts during the second Wind Forecast Improvement Project (WFIP2) aimed to address the challenges related to this. We evaluated three versions of the National Oceanic and Atmospheric Administration (NOAA) High-Resolution Rapid Refresh model with two different horizontal grid spacings against in situ and remote sensing observations to investigate how developments in physical parameterizations and numerical methods targeted during WFIP2 impacted the simulation of a persistent cold-air pool in the Columbia River basin. Differences amongst model versions were most apparent in simulated temperature and low-level cloud fields during the persistent phase of the cold-air pool. The model developments led to an enhanced low-level cloud cover, resulting in better agreement with the observations. This removed a diurnal cycle in the near-surface temperature bias at stations throughout the basin by reducing a cold bias during the night and a warm bias during the day. However, low-level clouds did not clear sufficiently during daytime in the newest model version, which leaves room for further model developments. The model developments also led to a better representation of the decay of the cold-air pool by slowing down its erosion.

54 ENVIRONMENTAL SCIENCES↗

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗

Quantitative relationship between structural orthorhombicity, shear modulus, and heat capacity anomaly of the nematic transition in iron-based superconductors

Electronic nematicity in iron pnictide materials has been extensively studied by various experimental techniques, yet its heat capacity anomaly at the phase transition has not been examined quantitatively. In this paper, we review the thermodynamic description of nematicity in Ba(Fe 1- x Cox) 2 As 2 using the Landau free energy, which defines the behavior of three thermodynamic quantities: the structural orthorhombicity that develops below the nematic transition, the softening shear modulus above the transition, and the discontinuous heat capacity at the transition. We derive a quantitative relationship between these three quantities, which is found to hold for a range of dopings. We report this result shows that the nematic transition is exceedingly well described by a mean-field model in the underdoped regime of the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculating Ti-6Al-4V β Transus Through a Chemistry-Based Equation Derived from Combined Element Binary Phase Diagrams

A mathematical model was developed to correlate the relationship between titanium alloying elements (Al, O, N, C, V, and Fe) in Ti-6Al-4V, and the resulting beta transus temperature (β{sub t}). The model was produced by analyzing the slope of the lines separating the alpha + beta (α + β) and beta (β) phase fields on combined element binary phase diagrams generated through the Pandat modeling software. Experimental results showed a strong agreement between the baseline data, obtained through differential thermal analysis (DTA), and the calculated β{sub t}, determined by measuring the chemistry at a location adjacent to the DTA sample. The difference between the calculated β{sub t} and the baseline DTA data was approximately ± 10 °F (± 5.56 °C). The thermodynamically ideal β{sub t} temperature was additionally modeled for each set of chemistries through the Pandat software. These results were compared to the β{sub t} temperatures determined through the chemistry-based equation and DTA methods.

36 MATERIALS SCIENCE↗