Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 667 records · Page 37

Isentropic phase changes in dissociating gases and the method of sound dispersion for the investigation of homogeneous gas reactions with very high speed : Conclusion

The analytical results of Part I are also applied to sound dispersion by friction and heat conduction, An irreversible change of momentum, energy, and type of particle corresponding to friction, heat conduction, and diffusion effects can appear both in the direction of the sound field and traverse to it. Longitudinal damping, the coupling of longitudinal damping and that due to chemical and physical changes, and coupling of diffusion and compositional changes are treated for a plane sound wave of infinite extent. The same principles are also applied to sound effects in cylindrical tubes. The limitations of the method are discussed in some detail.

SOUND VELOCITY - GASES↗

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗

Development and testing of coarse-grained models for ultrasonic simulations of cast austenitic stainless steel

Ultrasonic inspection of cast austenitic stainless steel (CASS) in the nuclear industry is particularly challenging because of sound field scatter and attenuation caused by the coarse-grained microstructure. Modeling and simulation are important tools in ultrasonic testing, as they can be used to help address key aspects of inspections, such as developing new probe designs, predicting inspection reliability, and testing phased-array focal laws. However, developing a useful and reliable CASS model is challenging due to the many grain interfaces and crystalline orientations that must be captured. We demonstrate a method of creating a realistic CASS model that is usable in CIVA, a commercially available modeling and simulation software platform. Using polished and chemically etched sections, we generate models of a coarse-grained equiaxed specimen and a columnar specimen. We also test an alternative method of generating a coarse-grained model using Voronoi regions. We qualitatively compare sound field scatter and quantitatively compare sound field attenuation and beam partitioning in simulated sound fields to those of laboratory-measured sound fields. Results show that the Voronoi models perform as well as or better than the models based on actual grain morphology. Here we also show that model-to-model randomness in Voronoi grain structure can impact the magnitude of a simulated echo response by a factor of two or more. Although CASS models are potentially a good depiction of reality for a given scenario, they should not be considered representative since CASS morphology can change significantly from specimen to specimen or within the same specimen.

36 MATERIALS SCIENCE↗

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Flare-induced MHD disturbances in the corona - Moreton waves and type 2 shocks

The propagation in the corona of the magnetohydrodynamic (MHD) disturbance possibly emitted at the explosive stage in the initial phase of a flare is considered. The behavior of the MHD fast-mode wavefront, whose source is located at the flare, is calculated by using eiconal-characteristic method in the High Altitude Observatory (HAO) realistic models of coronal magnetic field and density for the days of some particular flare events. It is shown as the result that the peculiar behavior of Moreton's surface wave and the peculiar appearance in the shape and position of the type 2 burst sources can be consistently understood by considering the refraction, focussing, and formation of shocks of MHD fast-mode disturbance in the actual distribution of Alfven velocity in the corona. Moreton waves seem to appear only when the flare explosion happens to occur at the edge of an active region and faces a low-Alfven-velocity region lying on the surface. The wave, which is initially emitted isotropically is refracted into a direction in which the condition for down-refraction holds to allow chromospheric reentrance of disturbance.

Uchida, Y.↗

Remote sensing of vegetation and soil using microwave ellipsometry

A method is described of determining vegetation height and water content of vegetation from the intensity and state of elliptical polarization of a reflected train of microwaves. The method comprises the steps of reflecting a circularly polarized train of microwaves from vegetation at a predetermined angle of incidence and detecting the reflected train of microwaves. The ratio of the intensities of the electric field vector components is determined, the phase difference of the components is measured, and the refractive index and thickness of the layer of vegetation are computed from a formula. The refractive index is given essentially by the water content of the vegetation.

Auer, S. O.↗

VIRGO: Experiment for helioseismology and solar irradiance monitoring

The scientific objectives of the variability of solar irradiance and gravity oscillations (VIRGO) experiment are as follows: to determine the characteristics of pressure and internal gravity oscillations by observing irradiance and radiance variations; to measure the solar total and spectral irradiance, and to quantify their variability. Helioseismological methods can be applied to these data in order to probe the solar interior. Certain convection characteristics and their interaction with magnetic fields will be studied from the results of the irradiance monitoring and from the comparison of the amplitudes and phases of the oscillations as observed from the brightness by VIRGO and from velocity by the global oscillations at low frequency (GOLF) experiment. The VIRGO experiment contains two active-cavity radiometers that monitor the solar constant, two three-channel sunphotometers that measure the spectral irradiance, and a low resolution imager with 12 pixels that measures the radiance distribution over the solar disk at 500 nm. The scientific objectives of VIRGO are presented, the instruments and the data acquisition and control system are described, and their measured performances are given.

Froehlich, Claus↗

Background Oriented Schlieren Implementation in a Jet-Surface Interaction Test

Many current and future aircraft designs rely on the wing or other aircraft surfaces to shield the engine noise from observers on the ground. However the available data regarding how a planar surface interacts with a jet to shield and/or enhance the jet noise are currently limited. Therefore, the Jet-Surface Interaction Tests supported by NASA's Fundamental Aeronautics Program's Fixed Wing Project were undertaken to supply experimental data covering a wide range of surface geometries and positions interacting with high-speed jet flows in order to support the development of noise prediction methods. Phase 1 of the Test was conducted in the Aero-Acoustic Propulsion Laboratory at NASA Glenn Research Center and consisted of validating noise prediction schemes for a round nozzle interacting with a planar surface. Phased array data and far-field acoustic data were collected for both the shielded and reflected sides of the surface. Phase 1 results showed that the broadband shock noise was greatly reduced by the surface when the jet was operated at the over-expanded condition, however, it was unclear whether this reduction was due a change in the shock cell structure by the surface. In the present study, Background Oriented Schlieren is implemented in Phase 2 of the Jet-Surface Interaction Tests to investigate whether the planar surface interacts with the high-speed jet ow to change the shock cell structure. Background Oriented Schlieren data are acquired for under-expanded, ideally-expanded, and over-expanded ow regimes for multiple axial and radial positions of the surface at three different plate lengths. These data are analyzed with far-field noise measurements to relate the shock cell structure to the broadband shock noise produced by a jet near a surface.

Clem, Michelle M.↗

Data Assimilation and Uncertainties in Early Solar Cycle Predictions

Stochastic nature of solar activity variations together with our limited knowledge of the dynamo mechanism and subsurface dynamics causes uncertainty in predictions of the solar cycle. For improving the physics-­‐based predictions we can take advantage of the mathematical data assimilation approach that allows us to take into account both, observational errors and model uncertainties, and provide estimates of the next solar cycle along with prediction uncertainties. In this study we use the Parker's migratory dynamo model together with the equation of magnetic helicity balance, which reproduces main properties of the sunspot cycles and allow us to minimize discrepancies between the observed global activity variations and the model solution. The test simulation runs show that a reliable prediction can be obtained for two phases of preceding solar cycle: 1) if the polar field reversals shortly after the solar maxima (strong toroidal field and weak poloidal field), and 2) during the solar minima (strongest poloidal and weak toroidal fields). The early estimate of Cycle 25 obtained by this method shows that this cycle will start in 2019 – 2020, reach the maximum in 2023 – 2024, and that the mean sunspot number at the maximum will be about 90 (for the v2.0 sunspot number series).

Data Assimilation↗

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗

Group-equivariant autoencoder for identifying spontaneously broken symmetries

We introduce the group-equivariant autoencoder (GE autoencoder), a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE autoencoder such that the encoder learns an order parameter invariant to these “never-broken” symmetries. This procedure produces a dramatic reduction in the number of free parameters such that the GE-autoencoder size is independent of the system size. We include symmetry regularization terms in the loss function of the GE autoencoder so that the learned order parameter is also equivariant to the remaining symmetries of the system. By examining the group representation by which the learned order parameter transforms, we are then able to extract information about the associated spontaneous symmetry breaking. We test the GE autoencoder on the 2D classical ferromagnetic and antiferromagnetic Ising models, finding that the GE autoencoder (1) accurately determines which symmetries have spontaneously broken at each temperature; (2) estimates the critical temperature in the thermodynamic limit with greater accuracy, robustness, and time efficiency than a symmetry-agnostic baseline autoencoder; and (3) detects the presence of an external symmetry-breaking magnetic field with greater sensitivity than the baseline method. Lastly, we describe various key implementation details, including a quadratic-programming-based method for extracting the critical temperature estimate from trained autoencoders and calculations of the DNN initialization and learning rate settings required for fair model comparisons.

42 ENGINEERING↗

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗

Multiple scattering in cloud layers; some results

Theoretical methods are discussed for calculating radiative effects of aerosols. Experimental determination is emphasized for relevant aerosol parameters on a global basis to arrive at realistic estimates of heating and cooling. Internal radiation fields in very thin and very thick slabs are reviewed. Phase functions, polarization, emission by internal sources, and path length distribution are also considered.

Vandehulst, H. C.↗

Ignition Transients of Large Segmented Solid Rocket Boosters

A model is described which provides a means for analyzing the complexities of ignition transients and pressure peaks of large, high performance, segmented solid rocket boosters. The method accounts for: (1) temporal and spatial development of the flow field set up by the head end igniter discharge, (2) ignition and flame spreading coupled to chamber flow, (3) the steep velocity, pressure, and temperature gradients that occur during the early phases of ignition, and (4) the interactions that produce ignition spikes (i.e., compression of chamber gases during pressurization, erosive burning, and mass added effect of igniter discharge). The technique differs from earlier models in that the flow interactions between the slots and main chamber are accounted for, and the original computer program for monolithic motors is improved. The procedures were used to predict the ignition transients of the current design for the space shuttle booster.

Caveny, L. H.↗

High altitude chemically reacting gas particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

The overall contractual effort and the theory and numerical solution for the Reacting and Multi-Phase (RAMP2) computer code are described. The code can be used to model the dominant phenomena which affect the prediction of liquid and solid rocket nozzle and orbital plume flow fields. Fundamental equations for steady flow of reacting gas-particle mixtures, method of characteristics, mesh point construction, and numerical integration of the conservation equations are considered herein.

Smith, S. D.↗

Encapsulation materials research

Encapsulation materials for solar cells were investigated. The different phases consisted of: (1) identification and development of low cost module encapsulation materials; (2) materials reliability examination; and (3) process sensitivity and process development. It is found that outdoor photothermal aging devices (OPT) are the best accelerated aging methods, simulate worst case field conditions, evaluate formulation and module performance and have a possibility for life assessment. Outdoor metallic copper exposure should be avoided, self priming formulations have good storage stability, stabilizers enhance performance, and soil resistance treatment is still effective.

Willis, P. B.↗