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 1,027 records · Page 57

Phase Transition in Cylindrical Ising Nanowires and Nanotubes: Approximation of the Molecular Field Theory

Within the approximation of the molecular field theory, the critical behavior of cylindrical Ising nanowires and nanotubes is studied. The model considered in this work consists of ferromagnetic spins S{sub c} situated in the core of the system and ferromagnetic spins S{sub s} situated on the surface wall; they are interconnected by the exchange coupling J{sub 1}. The transition temperature T{sub c} for such systems has been calculated as a function of exchange interaction parameters. The influence of the surface exchange interaction and surface coupling on the magnetodynamic behavior of the system have been studied. Some characteristic properties obtained on the phase diagram are related to the ratio of characteristic physical parameters on the surface and core of the abovementioned nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modelling the complex evaporated gas flow and its impact on particle spattering during laser powder bed fusion

The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.

36 MATERIALS SCIENCE↗

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗

Historical total ozone radiative forcing derived from CMIP6 simulations

Radiative forcing (RF) time series for total ozone from 1850 up to the present day are calculated based on historical simulations of ozone from 10 climate models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). In addition, RF is calculated for ozone fields prepared as an input for CMIP6 models without chemistry schemes and from a chemical transport model simulation. A radiative kernel for ozone is constructed and used to derive the RF. The ozone RF in 2010 (2005–2014) relative to 1850 is 0.35 W m -2 [0.08–0.61] (5–95% uncertainty range) based on models with both tropospheric and stratospheric chemistry. One of these models has a negative present-day total ozone RF. Excluding this model, the present-day ozone RF increases to 0.39 W m -2 [0.27–0.51] (5–95% uncertainty range). The rest of the models have RF close to or stronger than the RF time series assessed by the Intergovernmental Panel on Climate Change in the fifth assessment report with the primary driver likely being the new precursor emissions used in CMIP6. The rapid adjustments beyond stratospheric temperature are estimated to be weak and thus the RF is a good measure of effective radiative forcing.

54 ENVIRONMENTAL SCIENCES↗

RANS Simulation ADM of the NREL Phase VI wind turbine modeled as MHK Turbine

Attached are the .cas and .dat files for the Reynolds Averaged Navier-Stokes (RANS) simulation of a single lab-scaled DOE RM1 turbine implemented in ANSYS FLUENT CFD-package. In this case study the flow field around and in the wake of the NREL Phase VI wind turbine, modeled is MHK turbine, is simulated using Actuator Disk Model (ADM) (a.k.a Porous Media) by solving RANS equations coupled with a turbulence closure model. It should be highlighted that in this simulation the actual geometry of the rotor blade is not modeled. The effect of turbine rotating blades are modeled using the Actuator Disk Theory (see the stated section of attached M.Sc. thesis for more details).

16 TIDAL AND WAVE POWER↗

Charged nontopological soliton of the SU(2) Multiplication-Sign U(1) gauge model

The SU(2) Multiplication-Sign U(1) gauge model that is the bosonic sector of the Standard Model of electroweak interactions is considered. The existence of electrically charged nontopological solitons is shown to be possible in this model. Some properties of a charged nontopological soliton are investigated. Asymptotic expressions are derived for the soliton radius, energy, and phase frequency in the thin-wall regime by the method of trial functions. Numerical solutions of the model field equations corresponding to electrically charged non-topological solitions have been obtained. The dependences of the soliton energy and charge on phase frequency are given for several model parameters. It follows from the data obtained that there exists a domain of parameters in which a charged nontopological soliton is stable to the transition to a plane-wave field configuration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Alpha phase precipitation from phase-separated beta phase in a model Ti-Mo-Al alloy studied by direct coupling of transmission electron microscopy and atom probe tomography

The benefit of direct coupling of APT with TEM dark field imaging to investigate early stages of phase transformation in multicomponent alloys is demonstrated by analyzing alpha phase precipitated in a model Ti-10 at% Mo-10 at% Al alloy during annealing at 400oC. Through such a direct coupling approach a thermodynamically unexpected solute partitioning trend between beta matrix and alpha precipitate is observed in the early stages of precipitation, which is explained based on possible nucleation of alpha phase in the Ti rich (Mo and Al depleted regions) created as a result of phase separation in beta matrix. On further higher temperature annealing at 600oC for 1 hour, the alpha precipitates were shown to grow and get enriched in Al and further depleted in Mo reaching the thermodynamic equilibrium.

Environmental Molecular Sciences Laboratory↗

Local measurement of thermal conductivity and diffusivity

Simultaneous measurement of local thermal diffusivity and conductivity is demonstrated on a range of ceramic samples. This was accomplished by measuring the temperature field spatial profile of samples excited by an amplitude modulated continuous wave laser beam. A thin gold film is applied to the samples to ensure strong optical absorption and to establish a second boundary condition that introduces an expression containing the substrate thermal conductivity. The diffusivity and conductivity are obtained by comparing the measured phase profile of the temperature field to a continuum based model. A sensitivity analysis is used to identify the optimal film thickness for extracting the both substrate conductivity and diffusivity. Proof of principle studies were conducted on a range of samples having thermal properties that are representatives of current and advanced accident tolerant nuclear fuels. It is shown that by including the Kapitza resistance as an additional fitting parameter, the measured conductivity and diffusivity of all the samples considered agreed closely with the literature values. A distinguishing feature of this technique is that it does not require a priori knowledge of the optical spot size which greatly increases measurement reliability and reproducibility.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Topological excitonic insulator with tunable momentum order

Correlated topological materials often maintain a delicate balance among physical symmetries. Many topological orders are symmetry protected, whereas most correlated phenomena arise from spontaneous symmetry breaking. Cases where symmetry breaking induces a non-trivial topological phase are rare. Here we demonstrate the presence of two such phases in Ta 2 Pd 3 Te 5 , where Coulomb interactions form excitons that condense below 100 K, one with zero and the other with finite momentum. We observed a full spectral bulk gap, which stems from exciton condensation. This topological excitonic insulator state spontaneously breaks mirror symmetries but involves a weak structural coupling. Scanning tunnelling microscopy shows gapless boundary modes in the bulk insulating phase. Their magnetic field response, together with theoretical modelling, indicates a topological origin. These observations establish Ta 2 Pd 3 Te 5 as a topological excitonic insulator in a three-dimensional crystal. Finally, our results manifest a unique sequence of topological exciton condensations in a bulk crystal, offering exciting opportunities to study critical behaviour and excitations.

36 MATERIALS SCIENCE↗

Mesoscale modeling of intergranular bubble percolation in nuclear fuels

Phase-field simulations are used to examine the variability of intergranular fission gas bubble growth and percolation on uranium dioxide grain boundaries on a mesoscopic length scale. Three key parameters are systematically varied in this study: the contact angle (or dihedral angle) defining the bubble shape, the initial bubble density on the grain boundary plane, and the ratio of the gas diffusivity on the grain boundary versus the grain interiors. The simulation results agree well with previous experimental data obtained for bubble densities and average bubble areas during coalescence events. Interestingly, the rate of percolation is found to be highly variable, with a large dependency on the contact angle and the initial bubble density and little-to-no dependency on the grain boundary gas diffusivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Critical number of flavors in QED

We demonstrate that in unquenched quantum electrodynamics (QED), chiral symmetry breaking ceases to exist above a critical number of fermion flavors N{sub f}. This is a necessary and sufficient consequence of the fact that there exists a critical value of electromagnetic coupling {alpha} beyond which dynamical mass generation gets triggered. We employ a multiplicatively renormalizable photon propagator involving leading logarithms to all orders in {alpha} to illustrate this. We study the flavor and coupling dependence of the dynamically generated mass analytically as well as numerically. We also derive the scaling laws for the dynamical mass as a function of {alpha} and N{sub f}. Up to a multiplicative constant, these scaling laws are related through ({alpha},{alpha}{sub c}){r_reversible}(1/N{sub f},1/N{sub f}{sup c}). Calculation of the mass anomalous dimension {gamma}{sub m} shows that it is always greater than its value in the quenched case. We also evaluate the {beta} function. The criticality plane is drawn in the ({alpha},N{sub f}) phase space which clearly depicts how larger N{sub f} is required to restore chiral symmetry for an increasing interaction strength.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Urgent need for warming experiments in tropical forests

Abstract Although tropical forests account for only a fraction of the planet's terrestrial surface, they exchange more carbon dioxide with the atmosphere than any other biome on Earth, and thus play a disproportionate role in the global climate. In the next 20 years, the tropics will experience unprecedented warming, yet there is exceedingly high uncertainty about their potential responses to this imminent climatic change. Here, we prioritize research approaches given both funding and logistical constraints in order to resolve major uncertainties about how tropical forests function and also to improve predictive capacity of earth system models. We investigate overall model uncertainty of tropical latitudes and explore the scientific benefits and inevitable trade‐offs inherent in large‐scale manipulative field experiments. With a Coupled Model Intercomparison Project Phase 5 analysis, we found that model variability in projected net ecosystem production was nearly 3 times greater in the tropics than for any other latitude. Through a review of the most current literature, we concluded that manipulative warming experiments are vital to accurately predict future tropical forest carbon balance, and we further recommend the establishment of a network of comparable studies spanning gradients of precipitation, edaphic qualities, plant types, and/or land use change. We provide arguments for long‐term, single‐factor warming experiments that incorporate warming of the most biogeochemically active ecosystem components (i.e. leaves, roots, soil microbes). Hypothesis testing of underlying mechanisms should be a priority, along with improving model parameterization and constraints. No single tropical forest is representative of all tropical forests; therefore logistical feasibility should be the most important consideration for locating large‐scale manipulative experiments. Above all, we advocate for multi‐faceted research programs, and we offer arguments for what we consider the most powerful and urgent way forward in order to improve our understanding of tropical forest responses to climate change.

Cavaleri, Molly A.↗

Gravity field determination for Mars Observer

Mars Observer will be the first near circular, low altitude and short-periodic orbiter of Mars. From previous Mars orbiters, such as Mariner 9 and Vikings 1 and 2, a large gravitational oblateness and in general a gravity field ten times stronger than the earth's field have been determined. Because of these, the gravity field will dominate the evolution of the spacecraft's orbital motion and is the principal error source for reconstructing the spacecraft's motion. Thus a gravity calibration (GC) period has been established prior to the mapping phase in order to refine the gravitational field model. This paper describes the plan being made by Navigation to implement the GC strategy. By analyzing simulated Doppler data acquired during the seven day GC period, we give (1) the expected improvement in the gravity field model and (2) the resulting improvement in reconstructed spacecraft orbital motion. In addition, we have developed a strategy for converging the GC gravity model.

Esposito, Pasquale B.↗

Phase diagram of chirally imbalanced QCD matter

We compute the QCD phase diagram in the plane of the chiral chemical potential and temperature using the linear sigma model coupled to quarks and to the Polyakov loop. The chiral chemical potential accounts for effects of imbalanced chirality due to QCD sphaleron transitions which may emerge in heavy-ion collisions. We found three effects caused by the chiral chemical potential: the imbalanced chirality (i) tightens the link between deconfinement and chiral phase transitions; (ii) lowers the common critical temperature; (iii) strengthens the order of the phase transition by converting the crossover into the strong first order phase transition passing via the second order end point. Since the fermionic determinant with the chiral chemical potential has no sign problem, the chirally imbalanced QCD matter can be studied in numerical lattice simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unified description of the dynamics of quintessential scalar fields

Using the dynamical system approach, we describe the general dynamics of cosmological scalar fields in terms of critical points and heteroclinic lines. It is found that critical points describe the initial and final states of the scalar field dynamics, but that heteroclinic lines give a more complete description of the evolution in between the critical points. In particular, the heteroclinic line that departs from the (saddle) critical point of perfect fluid-domination is the representative path in phase space of quintessence fields that may be viable dark energy candidates. We also discuss the attractor properties of the heteroclinic lines, and their importance for the description of thawing and freezing fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Diagnostics for Equilibrium Reconstructions in the Presence of Nonaxisymmetric Eddy Current Distributions in Tokamaks

The lithium tokamak experiment LTX is a modest-sized spherical tokamak R0=0.4 m and a =0.26 m designed to investigate the low-recycling lithium wall operating regime for magnetically confined plasmas. LTX will reach this regime through a lithium-coated shell internal to the vacuum vessel, conformal to the plasma last-closed-flux surface, and heated to 300-400 oC. This structure is highly conductive and not axisymmetric. The three-dimensional nature of the shell causes the eddy currents and magnetic fields to be three-dimensional as well. In order to analyze the plasma equilibrium in the presence of three-dimensional eddy currents, an extensive array of unique magnetic diagnostics has been implemented. Sensors are designed to survive high temperatures and incidental contact with lithium and provide data on toroidal asymmetries as well as full coverage of the poloidal cross-section. The magnetic array has been utilized to determine the effects of nonaxisymmetric eddy currents and to model the start-up phase of LTX. Measurements from the magnetic array, coupled with two-dimensional field component modeling, have allowed a suitable field null and initial plasma current to be produced. For full magnetic reconstructions, a three-dimensional electromagnetic model of the vacuum vessel and shell is under development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase diagram and critical properties in the Polyakov-Nambu-Jona-Lasinio model

We investigate the phase diagram of the so-called Polyakov-Nambu-Jona-Lasinio model at finite temperature and nonzero chemical potential with three quark flavours. Chiral and deconfinement phase transitions are discussed, and the relevant order-like parameters are analyzed. The results are compared with simple thermodynamic expectations and lattice data. A special attention is payed to the critical end point: as the strength of the flavour-mixing interaction becomes weaker, the critical end point moves to low temperatures and can even disappear.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗