Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion process”

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 1,063 records · Page 59

ION ACCELERATION AT THE QUASI-PARALLEL BOW SHOCK: DECODING THE SIGNATURE OF INJECTION

Collisionless shocks are efficient particle accelerators. At Earth, ions with energies exceeding 100 keV are seen upstream of the bow shock when the magnetic geometry is quasi-parallel, and large-scale supernova remnant shocks can accelerate ions into cosmic-ray energies. This energization is attributed to diffusive shock acceleration; however, for this process to become active, the ions must first be sufficiently energized. How and where this initial acceleration takes place has been one of the key unresolved issues in shock acceleration theory. Using Cluster spacecraft observations, we study the signatures of ion reflection events in the turbulent transition layer upstream of the terrestrial bow shock, and with the support of a hybrid simulation of the shock, we show that these reflection signatures are characteristic of the first step in the ion injection process. These reflection events develop in particular in the region where the trailing edge of large-amplitude upstream waves intercept the local shock ramp and the upstream magnetic field changes from quasi-perpendicular to quasi-parallel. The dispersed ion velocity signature observed can be attributed to a rapid succession of ion reflections at this wave boundary. After the ions’ initial interaction with the shock, they flow upstream along the quasi-parallel magnetic field. Each subsequent wavefront in the upstream region will sweep the ions back toward the shock, where they gain energy with each transition between the upstream and the shock wave frames. Within three to five gyroperiods, some ions have gained enough parallel velocity to escape upstream, thus completing the injection process.

79 ASTRONOMY AND ASTROPHYSICS↗

Finding the Remnants of the Milky Way's Last Neutron Star Mergers

The discovery of a binary neutron star merger (NSM) through both its gravitational wave and electromagnetic emission has revealed these events to be key sites of $r$-process nucleosynthesis. Here, we evaluate the prospects of finding the remnants of Galactic NSMs by detecting the gamma-ray decay lines from their radioactive $r$-process ejecta. We find that 126 Sn, which has several lines in the energy range 415-695 keV and resides close to the second $r$-process peak, is the most promising isotope, because of its half-life $t_{1/2}$ = 2:30(14) 10 5 yr being comparable to the ages of recent NSMs. Using a Monte Carlo procedure, we predict that multiple remnants are detectable as individual sources by next-generation $\gamma$-ray telescopes which achieve sub-MeV line sensitivities of ~18 -8 -10 -6 $\gamma$ cm -2 s -1 . However, given the unknown locations of the remnants, the most promising search strategy is a systematic survey of the Galactic plane and bulge extending to high Galactic latitudes. Individual known supernova remnants which may be mis-classi ed NSM remnants could also be targeted, especially those located outside the Galactic plane. Detection of a moderate sample of Galactic NSM remnants would provide important clues to unresolved issues such as the production of actinides in NSMs, properties of merging NS binaries, and even help distinguish them from rare supernovae as current Galactic $r-$process sources. We also investigate the diffuse flux from longer-lived nuclei (e.g. 182 Hf) that could in principle trace the Galactic spatial distribution of NSMs over longer timescales, but find that the detection of the diffuse flux appears challenging even with next-generation telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

Anomalous Distributions of Primary Cosmic Rays as Evidence for Time-dependent Particle Acceleration in Supernova Remnants

Recent precise measurements of cosmic-ray (CR) spectra show that the energy distribution of protons is softer than those of heavier nuclei, and there are spectral hardenings for all nuclear compositions above ∼200 GV. Models proposed for these anomalies generally assume steady-state solutions of the particle acceleration process. We show that if the diffusion coefficient has a weak dependence on the particle rigidity near shock fronts of supernova remnants (SNRs), time-dependent solutions of the linear diffusive shock acceleration at two stages of SNR evolution can naturally account for these anomalies. The high-energy component of CRs is dominated by acceleration in the free expansion and adiabatic phases with enriched heavy elements and a high shock speed. The low-energy component may be attributed to acceleration by slow shocks propagating in dense molecular clouds with low metallicity in the radiative phase. Instead of a single power-law distribution, the spectra of time-dependent solutions soften gradually with the increase of energy, which may be responsible for the “knee” of CRs.

79 ASTRONOMY AND ASTROPHYSICS↗

Incoherent Dynamics of Intense Proton Beams Under Electron Cooling

The minimum emittance of ion beams achieved using electron cooling is limited by the heating processes of Intra Beam Scattering and diffusion driven by resonance crossing of particles due to space-charge. We describe a new experiment to explore the intense space-charge regime with a transverse tune shift approaching -0.5 using 2.5 MeV protons at the Integrable Optics Test Accelerator (IOTA) at Fermilab. We also report on the results from PyORBIT simulations incorporating transverse space-charge and electron cooling with emphasis on the incoherent dynamics of the particles.

43 PARTICLE ACCELERATORS↗

A spectral analysis of the domain decomposed Monte Carlo method for linear systems

The domain decomposed behavior of the adjoint Neumann-Ulam Monte Carlo method for solving linear systems is analyzed using the spectral properties of the linear operator. Relationships for the average length of the adjoint random walks, a measure of convergence speed and serial performance, are made with respect to the eigenvalues of the linear operator. In addition, relationships for the effective optical thickness of a domain in the decomposition are presented based on the spectral analysis and diffusion theory. Using the effective optical thickness, the Wigner rational approximation and the mean chord approximation are applied to estimate the leakage fraction of stochastic histories from a domain in the decomposition as a measure of parallel performance and potential communication costs. The one-speed, two-dimensional neutron diffusion equation is used as a model problem to test the models for symmetric operators. In general, the derived approximations show good agreement with measured computational results. (authors)

97 MATHEMATICS AND COMPUTING↗

Modified Closures in Monte Carlo Algorithms for Diffusive Binary Stochastic Media Transport Problems

In a stochastic medium, the material properties at a given spatial location are known only statistically. The most common approach to solving particle transport problems involving binary stochastic media (BSM) is to use the atomic mix approximation in which the transport problem is solved using homogenized (volume-averaged) material properties. A common deterministic model developed for solving BSM particle transport problems is the Levermore-Pomraning (LP) model. Zimmerman and Adams proposed a Monte Carlo algorithm (Algorithm A) that solves the LP equations and another Monte Carlo algorithm (Algorithm B) that locally preserves the sampled material realization; we refer to these Monte Carlo algorithms as the LP and LPLRP (LP local realization preserving) algorithms, respectively. One-dimensional (1D) planar geometry benchmark studies have shown that the LPLRP algorithm is often significantly more accurate than the LP algorithm for problems with an incident angular flux as well as for problems with an interior source. The LPLRP algorithm implementation in these benchmark comparisons made explicit use of the one-dimensional nature of the problem. The LP model is derived by assuming an 'upwind' closure in the coupling term relating the two materials in a binary stochastic medium. Su and Pomraning developed a modified form of this closure by considering the small correlation length limit and requiring the modified closure to produce the correct exponential decay for a source-free half-line albedo problem in rod geometry. They concluded that the modified closure is generally not inferior to the LP closure and in some cases is significantly better. Brantley investigated the use of the Su-Pomraning (SP) closure in the Monte Carlo LP algorithm (LP-SP) for the suite of benchmark problems and concluded that 1) the LP-SP algorithm was somewhat more accurate overall than the LP algorithm and somewhat less accurate overall than the LPLRP algorithm for an incident angular flux benchmark suite, and 2) the LP-SP algorithm was generally the least accurate of the algorithms for an interior source benchmark suite. Larsen et al. performed an asymptotic analysis of the transport equation in 1D planar geometry for the situation in which the physical system is 1) a random binary stochastic medium with material macroscopic total cross sections and mean slab width values of O(1) or smaller and 2) optically thick with weak absorption and sources at each spatial point and therefore globally diffusive. The asymptotic analysis demonstrates that, under these assumptions, the transport equation limits to the conventional diffusion equation with atomically-mixed (volume-averaged) material properties. Larsen et al. also demonstrate that, under these same assumptions, the LP equations asymptotically limit to an atomically-mixed diffusion equation with a diffusion coefficient that is too large, leading to an unphysical flattening of the ensemble-averaged scalar flux distribution. Vasques and Yadav recently performed an asymptotic analysis of an adjusted Levermore-Pomraning closure in which the Markovian transition functions are rescaled such that the equations asymptotically limit to the correct atomically-mixed diffusion equation; we refer to this model as LP-VY. In this paper, we demonstrate that the adjusted LP closure proposed by Vasques and Yadav is a special case of the SP closure obtained by assuming that the absorption cross sections in the two materials are equal. We further demonstrate that the Su-Pomraning closure has the correct asymptotic limit for the diffusive physical system and is therefore an appropriate closure for these problems. We also describe how to incorporate these various closures into the Monte Carlo LP and LPLRP particle transport algorithms, and we present numerical results comparing the accuracy of the Monte Carlo LP, LP-VY, LP-SP, LPLRP, and LPLRP-SP algorithms for the set of diffusive benchmark problems. We demonstrated that the SP closure has the correct asymptotic limit for the diffusive physical system under investigation in this paper and that the VY closure is a special case of the SP closure obtained by assuming equal absorption in the two materials. Both of these closures are readily implemented in the Monte Carlo LP and LPLRP algorithms for particle transport in BSM. Through numerical comparisons to a diffusive benchmark suite, we confirmed that 1) the LP and LPLRP algorithms do not limit to the appropriate diffusion solution and 2) the LP-VY, LP-SP, and LPLRP-SP algorithms limit to the appropriate diffusion solution and are of comparable accuracy. The LPLRP algorithm has been shown to be generally more accurate than the LP algorithm for non-diffusive problems. However, the numerical results in this paper indicate that the LPLRP algorithm does not limit to the correct diffusion equation for the asymptotic limit under consideration. The work in this paper demonstrates that the SP algorithm does limit to the correct diffusion equation for the asymptotic limit considered. Future work will investigate the possibility of whether the LPLRP-SP algorithm is accurate for both diffusive and non-diffusive problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗

Three-Dimensional Tracking of Interfacial Hopping Diffusion

Theoretical predictions have suggested that molecular motion at interfaces—which influences processes including heterogeneous catalysis, (bio)chemical sensing, lubrication and adhesion, and nanomaterial self-assembly—may be dominated by hypothetical “hops” through the adjacent liquid phase, where a diffusing molecule readsorbs after a given hop according to a probabilistic “sticking coefficient.” Here, we use three-dimensional (3D) single-molecule tracking to explicitly visualize this process for human serum albumin at solid-liquid interfaces that exert varying electrostatic interactions on the biomacromolecule. Following desorption from the interface, a molecule experiences multiple unproductive surface encounters before readsorption. An average of approximately seven surface collisions is required for the repulsive surfaces, decreasing to approximately two and a half for surfaces that are more attractive. The hops themselves are also influenced by long-range interactions, with increased electrostatic repulsion causing hops of longer duration and distance. These findings explicitly demonstrate that interfacial diffusion is dominated by biased 3D Brownian motion involving bulk-surface coupling and that it can be controlled by influencing short- and long-range adsorbate-surface interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic flux and heat losses by diffusive, advective, and Nernst effects in MagLIF-like plasma

The MagLIF approach to inertial confinement fusion involves subsonic/isobaric compression and heating of a DT plasma with frozen-in magnetic flux by a heavy cylindrical liner. The losses of heat and magnetic flux from the plasma to the liner are thereby determined by plasma advection and gradient-driven transport processes, such as thermal conductivity, magnetic field diffusion and thermomagnetic effects. Theoretical analysis based on obtaining exact self-similar solutions of the classical collisional Braginskii's plasma transport equations in one dimension demonstrates that the heat loss from the hot plasma to the cold liner is dominated by the transverse heat conduction and advection, and the corresponding loss of magnetic flux is dominated by advection and the Nernst effect. For a large electron Hall parameter ω{sub e}τ{sub e} effective diffusion coefficients determining the losses of heat and magnetic flux are both shown to decrease with ω{sub e}τ{sub e} as does the Bohm diffusion coefficient, which is commonly associated with low collisionality and two-dimensional transport. This family of exact solutions can be used for verification of codes that model the MagLIF plasma dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Automated Analysis, Reporting, and Archiving for Robotic Nondestructive Assay of Holdup Deposits - 19508

To decommission deactivated gaseous diffusion enrichment facilities, miles of contaminated pipe must be measured. The current method requires thousands of manual measurements, repeated manual data transcription, and months of manual analysis. The Pipe Crawling Activity Measurement System (PCAMS), developed by Carnegie Mellon University and in commissioning for use at the DOE Portsmouth Gaseous Diffusion Enrichment Facility, uses a robot to measure Uranium-235 from inside pipes and automatically log the data. Radiation measurements, as well as imagery, geometric modeling, and precise measurement positioning data are digitally transferred to the PCAMS server. On the server, data can be automatically processed in minutes and summarized for analyst review. Measurement reports are auto-generated with the push of a button. A database specially-configured to hold heterogeneous data such as spectra, images, and robot trajectories serves as archive. This paper outlines the features and design of the PCAMS Post-Processing Software, currently in commissioning for use at the Portsmouth Gaseous Diffusion Enrichment Facility. The analysis process, the analyst interface to the system, and the content of auto-generated reports are each described. Example pipe-interior geometric surface models, illustration of how key report features apply in operational runs, and user feedback are discussed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Swelling and Diffusion during Methanol Sorption into Hydrated Nafion

Diffusion within polymer electrolyte membranes is often coincident with time-dependent processes such as swelling and polymer relaxation, which are factors that limit their ability to block molecular crossover during use. The solution-diffusion model of membrane permeation, which is the accepted theory for dense polymers, applies only to steady-state processes and does not address dynamic internal structural changes that can accompany permeation. To begin discovery of how such changes can be coupled to the permeation process, we have constructed here a stochastic multiscale reaction-diffusion model that examines time-dependent methanol uptake into and swelling of hydrated Nafion. Several potential mechanisms of diffusion and polymer response are tested. The simulation predictions are compared to real-time Fourier transform infrared attenuated total reflectance spectroscopy (FTIR-ATR) absorbance reported in the literature [ Hallinan , D. T. , Jr. ; Elabd , Y. A. J. Phys. Chem. B 2007 , 111 , 13221 - 13230 ]. Of the proposed polymer response mechanisms, only one, a reaction-limited, local response to increasing methanol concentration that takes the entire experimental time frame of 600 s, produces simulated FTIR-ATR data consistent with experiment. The simulations show that water diffusion out of the membrane is minimal during methanol sorption and that changes in the measured infrared absorbances are due primarily to the increase in methanol concentration accompanied by dilution of water during swelling. Swelling involves densification of the polymer structure even as there is an overall volume expansion of the film. Potential connections between the polymer densification and molecular-level structural changes of Nafion in methanol are discussed. These results indicate that the interaction between methanol and Nafion serves to increase Nafion's capacity to accommodate large volumes of methanol-water solutions, facilitating increased permeation across the membrane relative to pure water.

36 MATERIALS SCIENCE↗

Tunable Atomic Layer Deposition into Ultra-High-Aspect-Ratio (>60000:1) Aerogel Monoliths Enabled by Transport Modeling

Atomic layer deposition (ALD) modification of ultra-high-aspect-ratio structures (>10,000:1) is a powerful platform with applications in catalysis, filtration, and energy conversion. However, the deposition of conformal and tunable ALD coatings at these aspect ratios remains challenging, resulting in empirical tradeoffs between precursor utilization and reaction time. Here, we demonstrate tunable control of the ALD infiltration depth into an aerogel monolith (AM) and develop a reaction-diffusion model to accurately describe the coating process. Specifically, we investigate the ALD exposure time and precursor dose needed to conformally coat a silica AM with pore sizes of ~20 nm, a monolith thickness of ~2.5 mm, and aspect ratios exceeding 60,000:1. We demonstrate complete infiltration into the AM, which is quantified by elemental mapping. A reaction-diffusion model is developed, which accounts for multiple doses and precursor depletion in the ALD chamber during an exposure step. The experimentally validated model enables the prediction and tuning of infiltration depth into a tortuous, high-aspect-ratio structure such as an AM, allowing for the synthesis of rationally designed material architectures. Additionally, the model allows for co-optimization of the total deposition time and percentage of unreacted precursor, which are important for the manufacturability and sustainability of ALD processing. Lastly, we demonstrate that ultra-thin ALD Al 2 O 3 coatings can be used to stabilize silica AMs against structural degradation under high-temperature annealing conditions (700-800 °C), by limiting changes in surface area and monolith volume. As a result, this improved high-temperature stability has implications for numerous aerogel applications, including catalysis and thermal insulation.

14 SOLAR ENERGY↗

Theory and modeling of active brazing.

Active brazes have been used for many years to produce bonds between metal and ceramic objects. By including a relatively small of a reactive additive to the braze one seeks to improve the wetting and spreading behavior of the braze. The additive modifies the substrate, either by a chemical surface reaction or possibly by alloying. By its nature, the joining process with active brazes is a complex nonequilibrium non-steady state process that couples chemical reaction, reactant and product diffusion to the rheology and wetting behavior of the braze. Most of the these subprocesses are taking place in the interfacial region, most are difficult to access by experiment. To improve the control over the brazing process, one requires a better understanding of the melting of the active braze, rate of the chemical reaction, reactant and product diffusion rates, nonequilibrium composition-dependent surface tension as well as the viscosity. This report identifies ways in which modeling and theory can assist in improving our understanding.

van Swol, Frank B.↗

Aminopolymer Mobility and Support Interactions in Silica-PEI Composites for CO 2 Capture Applications: A Quasielastic Neutron Scattering Study

Composite gas sorbents, formed from an active polymer phase and a porous support, are promising materials for the separation of acid gases from a variety of gas streams. Significant changes in sorption performance (capacity, rate, stability etc.) can be achieved by tuning the properties of the polymer and the nature of interactions between polymer and support. We utilize quasielastic neutron scattering (QENS) and coarse-grained molecular dynamics (MD) simulations to characterize the dynamic behavior of the most commonly reported polymer in such materials, poly(ethylenimine) (PEI), both in bulk form and when supported in a mesoporous silica framework. The polymer chain dynamics (rotational and translational diffusion) are characterized using two neutron backscattering spectrometers that have overlapping time scales, ranging from picoseconds to a few nanoseconds. Two modes of motion are detected for the PEI molecule in QENS. At low energy transfers, a “slow process” on the time scale of ~200 ps is found and attributed to jump-mediated, center-of-mass diffusion. Second, a “fast process” at ~20 ps scale is also found and is attributed to a locally confined, jump-diffusion. Characteristic data (time scale and spectral weight) of these processes are compared to those characterized by MD, and reasonable agreement is found. For the nanopore-confined PEI, we observe a significant reduction in the time scale of polymer motion as compared to the bulk. The impacts of silica surface functionalization and of polymer fill fraction in the silica pores (controlling the portion of polymer molecules in contact with the pore walls), are both studied in detail. Hydrophobic functionalization of the silica leads to an increase of the PEI mobility above that in native silanol-terminated silica, but the dynamics are still slower than those in bulk PEI. Sorbents with faster PEI dynamics are also found to be more efficient for CO 2 capture, possibly because sorption sites are more accessible than those in systems with slower PEI dynamics. Therefore, this work supports the existence of a link between the affinity of the support for PEI and the accessibility of active sorbent functional groups.

36 MATERIALS SCIENCE↗

Reversible Tuning of the Surface Plasmon Resonance of Indium Tin Oxide Nanocrystals by Gas-Phase Oxidation and Reduction

Heavily doped oxide nanocrystals exhibit a tunable localized surface plasmon resonance (LSPR) in the infrared, a property that is promising for applications in photonics, spectroscopy, and photochemistry. Nanocrystal carrier density and, thus, spectral response are adjustable via chemical reaction; yet, the fundamental processes that govern this behavior are poorly understood. In this work, we study the time dependence of the LSPR supported by indium tin oxide (ITO) nanocrystals during O 2 and N 2 annealing with in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). We show that the LSPR red-shifts upon oxidation in O 2 and blue-shifts to its original position upon reduction in N 2 . A reaction–diffusion model allows us to rationalize the underlying physicochemical processes and quantitatively connect nanocrystal redox chemistry, solid-state diffusion, carrier density, and the LSPR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon fiber design improvements based on economic models

Circular fiber geometries are predominant in commercial carbon fiber material systems, but the use of this fiber shape has numerous limitations. Circular geometries have continuous symmetry, which is helpful for various processing considerations, but also have the largest possible maximum diffusion thickness for a given fiber area. This characteristic means that circular carbon fibers always have the highest material processing cost and lowest production throughput compared to any other fiber shape with the same area and tow count. To quantify material cost and other benefits for non-circular carbon fiber geometries, process models for polyacrylonitrile based carbon fiber production are developed in relationship to the carbon fiber shape, size, and tow count. For a given fiber shape, precursor production costs are shown to favor maximizing fiber size while conversion costs are minimized by the smallest fiber size. These competing cost trends result in a numerically optimal fiber size for a given shape and tow count, while both cost components are decreased by increasing tow count. Cost–performance tradeoffs for three lobe fiber geometries are studied by supplementing the cost trends with a numerical failure model to predict compressive strength for discrete shape variants. The shape selection is shown to be more sensitive to variations in cost than compressive strength while suboptimal shape designs can improve manufacturing robustness and achievable fiber volume fractions. Finally, an optimal three lobe carbon fiber is identified that balances the set of considerations while reducing costs and embodied energy and increasing production throughput compared to a commercial carbon fiber.

Carbon fiber↗

Synthesis, Characterization and Experimental Studies of Nano Zn–Al–Fe{sub 3}O{sub 4} Blended Alginate/Ca Beads for the Adsorption of Rhodamin B

In the present study, adsorptive removal of rhodamin-B (RMN-B) from aqueous solution was investigated by using nano Zn–Al–Fe{sub 3}O{sub 4} incorporated sodium alginate (SA) beads. SA is a natural polysaccharide, non toxic and biodegradable material. The incorporation of Zn–Al–Fe{sub 3}O{sub 4} into SA was done by ionic polymerization method using CaCl{sub 2}. A batch mode study was used for the adsorption of RMN-B onto Zn–Al–Fe{sub 3}O{sub 4}/ alginate beads. The adsorption efficiency of RMN-B was investigated by the following modeling equations like Freundlich, Langmuir, Pseudo-first order, Pseudo-second order, Intraparticle diffusion, Gibbs free energy, Enthalpy and Entropy. The predicted experimental results of various modeling equations were implied that the process follow Langmuir and pseudo-second order equations. The intraparticle diffusion parameters indicate that the diffusion between RMN-B and nano beads were predominant in throughout of process. The surface morphology of nano beads of before and after adsorption were show that RMN-B is strongly adhered on whole surface of the nano material. The FTIR spectrum indicates that the RMN-B strongly adsorbed on the surface and TGA study reveals that the alginate/Zn–Al–Fe{sub 3}O{sub 4}/Ca material has an excellent temperature withstand capacity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Microstructural Modeling of the α + β Phase in Ti-6Al-4V: A Diffusion-Based Approach

Complex heat treatment operations and advanced manufacturing processes such as laser or electron-beam welding will see the metallic workpiece experience a considerable range of temperatures and heating/cooling rates. These intrinsic conditions will have a significant bearing upon the microstructure of the material, and in turn upon the thermo-mechanical properties. In this work, a diffusion-based approach to model the growth and shrinkage of precipitates in the alpha + beta field of the common titanium alloy Ti-6Al-4V is established. Further, the numerical model is extended using a JMA-type approach to explore the dependency of the beta-transus temperature on extremely high heating rates, whereby dissolution alone is too slow to accurately describe the alpha to beta-phase transformation. Experimental heat treatments at rates of 5, 50, and 500 °C/s were performed, and metallographic analysis of the samples was used to validate the numerical modeling framework predictions for lamellar shrinkage, while data from the literature has been used to evaluate the numerical modeling framework for the growth of equiaxed microstructures. The agreement between measurements and numerical predictions was found to be good.

36 MATERIALS SCIENCE↗