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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Additive manufacturing of nanotube-loaded thermosets via direct ink writing and radio-frequency heating and curing

Direct Ink Writing (DIW) is an extrusion-based additive manufacturing method where the print medium is a liquid-phase ‘ink’ dispensed out of nozzles and deposited along digitally defined paths. Conventional DIW of thermosetting resins relies on viscosity modifying agents, novel crosslinking chemistries, and/or long curing schedules in an oven. Here we demonstrate the use of a co-planar radio frequency applicator to generate an electric field, which can be used to rapidly heat and cure nano-filled composite resins as they are printed. This method avoids the need for an oven or post- curing step. This process consists of a sequential print-and-cure cycle which allows for printing of high-resolution, multi-layered structures. Every extruded layer is partially cured using RF before depositing the next layer; this allows the printed part to maintain structural integrity. The process enables both increased throughput and decreased touch time relative to traditional manufacturing. Commercial epoxy resin with varied nano-filler loadings were examined as DIW candidates. After printing, the thermo-mechanical properties, surface finish, and shape retention of RF-cured samples were comparable to conventionally cured samples. Furthermore, this method of manufacturing establishes RF heating as a suitable alternative to conventional methods, facilitating rapid, free-form processing of thermosetting resins without a mold.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Physics-Informed Neural Networks for Heat Transfer Problems

Abstract Physics-informed neural networks (PINNs) have gained popularity across different engineering fields due to their effectiveness in solving realistic problems with noisy data and often partially missing physics. In PINNs, automatic differentiation is leveraged to evaluate differential operators without discretization errors, and a multitask learning problem is defined in order to simultaneously fit observed data while respecting the underlying governing laws of physics. Here, we present applications of PINNs to various prototype heat transfer problems, targeting in particular realistic conditions not readily tackled with traditional computational methods. To this end, we first consider forced and mixed convection with unknown thermal boundary conditions on the heated surfaces and aim to obtain the temperature and velocity fields everywhere in the domain, including the boundaries, given some sparse temperature measurements. We also consider the prototype Stefan problem for two-phase flow, aiming to infer the moving interface, the velocity and temperature fields everywhere as well as the different conductivities of a solid and a liquid phase, given a few temperature measurements inside the domain. Finally, we present some realistic industrial applications related to power electronics to highlight the practicality of PINNs as well as the effective use of neural networks in solving general heat transfer problems of industrial complexity. Taken together, the results presented herein demonstrate that PINNs not only can solve ill-posed problems, which are beyond the reach of traditional computational methods, but they can also bridge the gap between computational and experimental heat transfer.

Engineering↗

NASA NTRS1995

Twenty years of interplanetary magnetic field variations with periods in the range of 10 days to 3 years

Twenty years of interplanetary magnetic field data collected primarily by the IMP-8 spacecraft near Earth has been analyzed by a dynamic periodogram method in search of significant periodicities in the range of 10 days to 3 years. The method has the advantage of detecting variations with time in the periodicities besides determining the power and phase of the dominant variations. It has been found that the well known periodicities near 1 year and 27 days are strongly modulated by the solar cycle. Both of these periodicities are only detected during solar minimum. During solar maximum. a number of unusual variations are observed. Special emphasis will be placed on the recently reported 1.3 year variation in solar wind parameters besides periods in the interplanetary magnetic field near 51, 73 and 154 days. Correlations with solar wind plasma and solar index variations will also be presented.

Szabo, A.↗

NASA NTRS1993

Effects of thermal-solutal convection on temperature and solutal fields under various gravitational orientations

Semiconductor crystals such as Hg(1-x)CD(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. The value of effective diffusion coefficient calculated from fitting of the data to 1D model varies with Hg(1-x)Cd(x)Te growth conditions. This indicates that the growth condition of the Hg(1-x)Cd(x)Te is not purely diffusion controlled. Because of the higher thermal conductivity in the melt than that in the crystal in the growth system, there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart, and Thorpe, Hutt and Soulsby have shown that under such conditon a fluid will become convectively unstable as a result of different diffusitivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositonal segregation in both axial and radial directions in the unidirectionally solidified crystals under various gravitational directions. To reach this goal, we start with a simplified problem to study the effects of thermal-solutal convection on the temperature and solutal fields under various gravitional orientations. We begin by reviewing model governing equations.

Wang, Jai-Ching↗

NASA NTRS1996

An Optimal Deconvolution Method for Reconstructing Pneumatically Distorted Near-Field Sonic Boom Pressure Measurements

In-flight measurements of the SR-71 near-field sonic boom were obtained by an F-16XL airplane at flightpath separation distances from 40 to 740 ft. Twenty-two signatures were obtained from Mach 1.60 to Mach 1.84 and altitudes from 47,600 to 49,150 ft. The shock wave signatures were measured by the total and static sensors on the F-16XL noseboo. These near-field signature measurements were distorted by pneumatic attenuation in the pitot-static sensors and accounting for their effects using optimal deconvolution. Measurement system magnitude and phase characteristics were determined from ground-based step-response tests and extrapolated to flight conditions using analytical models. Deconvolution was implemented using Fourier transform methods. Comparisons of the shock wave signatures reconstructed from the total and static pressure data are presented. The good agreement achieved gives confidence of the quality of the reconstruction analysis. although originally developed to reconstruct the sonic boom signatures from SR-71 sonic boom flight tests, the methods presented here generally apply to other types of highly attenuated or distorted pneumatic measurements.

Whitmore, Stephen A.↗

NASA NTRS1999

A Mode Detection Method Using the Azimuthal Directivity of a Turbofan Model

The azimuthal, far field directivity of a scale fan model was measured in high resolution. The model is a 12 inch diameter rotor with 16 blades followed by 40 stator vanes. The tests were conducted at the nominal 100% speed corresponding to a tip speed of 905 ft/sec. Measurement of the radiated sound field, forward of the fan, was made in an anechoic chamber with an inflow control device and a baffle separating the aft and forward radiated interaction noise. The acoustic field was surveyed with a circular hoop array of 16 microphones which was moved to 14 axial stations. At each axial station the hoop was rotated in half-degree increments to take 736 points in the azimuthal angle. In addition to sound pressure level, the phase angle relative to a reference microphone was measured at each point. The sound pressure level is shown to vary in patterns by 10-15 dB especially for the fundamental tone but also for the first and second harmonic. A far field mode detection method has been developed and used with the data which determines the modes generated by the fan and which then interact to form the azimuthal directivity.

Thomas, R. H.↗

DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1987

Discrete element weld model, phase 2

A numerical method was developed for analyzing the tungsten inert gas (TIG) welding process. The phenomena being modeled include melting under the arc and the flow in the melt under the action of buoyancy, surface tension, and electromagnetic forces. The latter entails the calculation of the electric potential and the computation of electric current and magnetic field therefrom. Melting may occur at a single temperature or over a temperature range, and the electrical and thermal conductivities can be a function of temperature. Results of sample calculations are presented and discussed at length. A major research contribution has been the development of numerical methodology for the calculation of phase change problems in a fixed grid framework. The model has been implemented on CHAM's general purpose computer code PHOENICS. The inputs to the computer model include: geometric parameters, material properties, and weld process parameters.

Prakash, C.↗

DOE OSTI2020

Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase) (Final Report)

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership (RCSP) Initiative. The public/private consortium, funded through the MRCSP, brings together nearly 40 industry partners in 10 contiguous states in Mid-West and North-East regions. Battelle, as the project lead, oversees research & development, operations, and coordinates activities among various partners. Since its founding, MRCSP has made significant progress towards advancing CCS and CCUS technologies for participating states in the region. Between 2008 to 2020, Phase III of the MRCSP program focused on the development of large-scale injection tests. This report is part of a series of topical reports prepared under the MRCSP Phase III. The MRCSP Final Technical Report summarizes the research conducted while monitoring large-scale injection tests at a CO 2 -EOR facility in Northern Michigan. The project successfully injected more than one million metric tons of CO 2 into depleted oil fields in the Northern Niagaran Pinnacle Reef Trend (NNPRT) for enhanced oil recovery (EOR). These commercial-scale tests provided valuable real-world knowledge that has been used to refine technologies and methods, reduce uncertainties, and demonstrate safety and effectiveness of CCUS technologies to increase public acceptance. The Final Technical Report provides an overview of Phase III efforts including advancing the understanding of the regional geology, storage capacities, deployment of many monitoring and modeling techniques and its efficacy, and technical demonstration of large-scale CO 2 injection, containment, and associated storage in the NNPRT region. Specific topical reports on the subjects covered in the Final Technical Report are also available through EDX/OSTI.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2024

Atomic resolution scanning transmission electron microscopy at liquid helium temperatures for quantum materials

Fundamental quantum phenomena in condensed matter, ranging from correlated electron systems to quantum information processors, manifest their emergent characteristics and behaviors predominantly at low temperatures. This necessitates the use of liquid helium (LHe) cooling for experimental observation. Atomic resolution scanning transmission electron microscopy combined with LHe cooling (cryo-STEM) provides a powerful characterization technique to probe local atomic structural modulations and their coupling with charge, spin and orbital degrees-of-freedom in quantum materials. However, achieving atomic resolution in cryo-STEM is exceptionally challenging, primarily due to sample drifts arising from temperature changes and noises associated with LHe bubbling, turbulent gas flow, etc. In this work, we demonstrate atomic resolution cryo-STEM imaging at LHe temperatures using a commercial side-entry LHe cooling holder. Firstly, we examine STEM imaging performance as a function of He gas flow rate, identifying two primary noise sources: He-gas pulsing and He-gas bubbling. Secondly, we propose two strategies to achieve low noise conditions for atomic resolution STEM imaging: either by temporarily suppressing He gas flow rate using the needle valve or by acquiring images during the natural warming process. Lastly, we show the applications of image acquisition methods and image processing techniques in investigating structural phase transitions in Cr 2 Ge 2 Te 6 , CuIr 2 S 4 , and CrCl 3 . In conclusion, our findings represent an advance in the field of atomic resolution electron microscopy imaging for quantum materials and devices at LHe temperatures, which can be applied to other commercial side-entry LHe cooling TEM holders.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗

NASA NTRS1973

Study of exothermic processes in shock ignited gases by the use of laser shear interferometry

Report on measurements of maximum exothermic power pulses attainable from a given chemical system. Experimental tests involved the use of a shock tube technique whereby the exothermic process of combustion was controlled by reflected shock, so that it occurred under virtually inviscid flow conditions, while the measurements were performed at a resolution commensurate with the actual rate of chemical reaction. Experimental observations were made by means of a novel method of laser shear interferometry - a cross-breed between holography and the conventional means for measuring refractive index fields, in that, on one hand, it was based on the exploitation of the phase coherence of the laser light beam, recording first a diffraction image of the wave fronts which, for the desired final result, had to be optically reconstructed, and, on the other, it yielded eventually either two-dimensional interferograms or schlieren photographs of the observed phenomena.

Meyer, J. W.↗

NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

DOE OSTI2025

Diffusion and phase formation in the γ-uranium-technetium system

Phase formation in the U-Tc binary system at 800 °C was investigated using a diffusion couple experiment. Scanning electron microscopy (SEM) and energy dispersive X-ray spectroscopy (EDS) identified four novel potential intermetallic phases - U 7 Tc 3 , U 13 Tc 12 , U 3 Tc 5 , and UTc 4 . Diffusion coefficients were calculated for the intermetallic phases using the Boltzmann-Matano method and were respectively found to be – 120, 38.2, 15.6, and 1.51 × 10 −18 m 2 /s. Tc also exhibits a solid-solution phase with high penetration into the U with a diffusion coefficient of ∼ 10 −14 m 2 /s. Furthermore, these findings expand the number of known U-Tc phases and provide the first diffusion coefficients for the U-Tc system, and contribute valuable data to the broader field of actinide metallurgy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Simple strategy for the simulation of axially symmetric large-area metasurfaces

Metalenses are composed of nanostructures for focusing light and have been widely explored in many exciting applications. However, their expanding dimensions pose simulation challenges. We propose a method to simulate metalenses in a timely manner using vectorial wave and ray tracing models. We sample the metalens’s radial phase gradient and locally approximate the phase profile by a linear phase response. Each sampling point is modeled as a binary blazed grating, employing the chosen nanostructure, to build a transfer function set. The metalens transmission or reflection is then obtained by applying the corresponding transfer function to the incoming field on the regions surrounding each sampling point. Fourier optics is used to calculate the scattered fields under arbitrary illumination for the vectorial wave method, and a Monte Carlo algorithm is used in the ray tracing formalism. We validated our method against finite-difference time domain simulations at 632 nm, and we were able to simulate metalenses larger than 3000 wavelengths in diameter on a personal computer.

Martins, Augusto (ORCID:0000000295546481)↗

NASA NTRS1977

Microscopic Lagrangian description of warm plasmas. III - Nonlinear wave-particle interaction

The averaged-Lagrangian method is applied to nonlinear wave-particle interactions in an infinite, homogeneous, magnetic-field-free plasma. The specific example of Langmuir waves is considered, and the combined effects of four-wave interactions and wave-particle interactions are treated. It is demonstrated how the latter lead to diffusion in velocity space, and the quasilinear diffusion equation is derived. The analysis is generalized to the random phase approximation. The paper concludes with a summary of the method as applied in Parts 1-3 of the paper.

Galloway, J. J.↗

DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗