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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗
NASA NTRS1991

Neutral line chaos and phase space structure

Phase space structure and chaos near a neutral line are studied with numerical surface-of-section (SOS) techniques and analytic methods. Results are presented for a linear neutral line model with zero crosstail electric field. It was found that particle motion can be divided into three regimes dependening on the value of the conserved canonical momentum, Py, and the conserved Hamiltonian, h. The phase space structure, using Poincare SOS plots, is highly sensitive to bn = Bn/B0 variations, but not to h variations. It is verified that the slow motion preserves the action, Jz, as evaluated by Sonnerup (1971), when the period of the fast motion is smaller than the time scale of the slow motion. Results show that the phase space structure and particle chaos depend sensitively upon Py and bn, but are independent of h.

Burkhart, Grant R.↗
DOE OSTI2024

Unveiling mesoscopic structures in distorted lamellar phases through deep learning-based small angle neutron scattering analysis

Hypothesis: The formation of distorted lamellar phases, distinguished by their arrangement of crumpled, stacked layers, is frequently accompanied by the disruption of long-range order, leading to the formation of interconnected network structures commonly observed in the sponge phase. Nevertheless, traditional scattering functions grounded in deterministic modeling fall short of fully representing these intricate structural characteristics. Our hypothesis posits that a deep learning method, in conjunction with the generalized leveled wave approach used for describing structural features of distorted lamellar phases, can quantitatively unveil the inherent spatial correlations within these phases. Experiments and Simulations: This report outlines a novel strategy that integrates convolutional neural networks and variational autoencoders, supported by stochastically generated density fluctuations, into a regression analysis framework for extracting structural features of distorted lamellar phases from small angle neutron scattering data. To evaluate the efficacy of our proposed approach, we conducted computational accuracy assessments and applied it to the analysis of experimentally measured small angle neutron scattering spectra of AOT surfactant solutions, a frequently studied lamellar system. Findings: The findings unambiguously demonstrate that deep learning provides a dependable and quantitative approach for investigating the morphology of wide variations of distorted lamellar phases. It is adaptable for deciphering structures from the lamellar to sponge phase including intermediate structures exhibiting fused topological features. In conclusion, this research highlights the effectiveness of deep learning methods in tackling complex issues in the field of soft matter structural analysis and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dispersive optical systems for scalable Raman driving of hyperfine qubits

Hyperfine atomic states are among the most promising candidates for qubit encoding in quantum information processing. In atomic systems, hyperfine transitions are typically driven through a two-photon Raman process by a laser field which is amplitude modulated at the hyperfine qubit frequency. Here we introduce a method for generating amplitude modulation by phase modulating a laser and reflecting it from a highly dispersive optical element known as a chirped Bragg grating. This approach is passively stable, offers high efficiency, and is compatible with high-power laser sources, enabling large Rabi frequencies and improved quantum coherence. We benchmark this approach by globally driving an array of approximately 300 neutral 87 Rb atomic qubits trapped in optical tweezers and obtain Rabi frequencies of 2 MHz with photon-scattering error rates of less than 2×10 -4 per π pulse. This robust approach can be directly integrated with local addressing optics in both neutral atom and trapped ion systems to facilitate high-fidelity single-qubit operations for quantum information processing.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1991

Development of N/P AlGaAs free-standing top solar cells for tandem applications

The combination of a free standing AlGaAs top solar cell and an existing bottom solar cell is the highest performance, lowest risk approach to implementing the tandem cell concept. The solar cell consists of an AlGaAs substrate layer, an AlGaAs base layer, an AlGaAs emitter, and an ultra-thin AlGaAs window layer. The window layer is compositionally graded which minimizes reflection at the window layer/emitter interface and creates a built-in electric field to improve quantum response in the blue region of the spectrum. Liquid phase epitaxy (LPE) is the only viable method to produce this free standing top solar cell. Small (0.125 sq cm), transparent p/n AlGaAs top solar cells were demonstrated with optimum bandgap for combination with a silicon bottom solar cell. The efficiency of an AlGaAs/Si stack using the free standing AlGaAs device upon an existing silicon bottom solar cell is 24 pct. (1X, Air Mass Zero (AM0). The n/p AlGaAs top solar cell is being developed in order to facilitate the wiring configuration. The two terminal tandem stack will retain fit, form, and function of existing silicon solar cells. Progress in the development of large area (8 and 16 sq cm), free standing AlGaAs top solar cells is discussed.

Negley, Gerald H.↗
NASA NTRS1990

Heat transfer in space systems; Proceedings of the Symposium, AIAA/ASME Thermophysics and Heat Transfer Conference, Seattle, WA, June 18-20, 1990

Theoretical and experimental studies of heat-tranfer in a space environment are discussed in reviews and reports. Topics addressed include a small-scale two-phase thermosiphon to cool high-power electronics, a low-pressure-drop heat exchanger with integral heat pipe, an analysis of the thermal performance of heat-pipe radiators, measurements of temperature and concentration fields in a rectangular heat pipe, and a simplified aerothermal heating method for axisymmetric blunt bodies. Consideration is given to entropy production in a shock wave, bubble-slug transition in a two-phase liquid-gas flow under microgravity, plasma arc welding under normal and zero gravity, the Microgravity Thaw Experiment, the flow of a thin film on stationary and rotating disks, an advanced ceramic fabric body-mounted radiator for Space Station Freedom phase 0 design, and lunar radiators with specular reflectors.

Chan, S. H.↗
DOE OSTI2026

Adaptively coupled phase retrieval in multi-peak Bragg coherent diffraction imaging

Recent advances in Bragg coherent diffraction imaging (BCDI) experimental techniques permit routine measurement of multiple Bragg peaks from a single crystalline grain. The resulting images contain the full lattice distortion vector field which can be differentiated to provide lattice strain and rotation. With the advent of fourth-generation synchrotron light sources, such multi-peak datasets are produced at high rates, facilitating the need for rapid phase retrieval of the multiple peaks and subsequent image analysis. Here we describe and demonstrate a new implementation of a coupled phase retrieval technique for multi-peak BCDI which simultaneously treats each Bragg peak of the dataset and produces a three-dimensional image of the crystal's morphology and lattice distortion field. In addition, this method uses the redundant information contained in the various Bragg diffraction patterns to detect and suppress spurious signal appearing on the detector in a subset of the measurements. Compared with manual data editing, adaptive coupling produces a more consistent phase profile in reciprocal space and sharper surfaces in direct space, with no significant difference in computational cost. These improvements reduce the need for manual preprocessing and enable robust high-throughput analysis of multi-peak BCDI data, supporting near-real-time strain microscopy at modern synchrotron facilities.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase Transitions in Organic and Organic/Inorganic Aerosol Particles

The phase state of aerosol particles can impact numerous atmospheric processes, including new particle growth, heterogeneous chemistry, cloud condensation nucleus formation, and ice nucleation. In this article, the phase transitions of inorganic, organic, and organic/inorganic aerosol particles are discussed, with particular focus on liquid-liquid phase separation (LLPS). The physical chemistry that determines whether LLPS occurs, at what relative humidity it occurs, and the resultant particle morphology is explained using both theoretical and experimental methods. The known impacts of LLPS on aerosol processes in the atmosphere are discussed. Finally, potential evidence for LLPS from field and chamber studies is presented. By understanding the physical chemistry of the phase transitions of aerosol particles, we will acquire a better understanding of aerosol processes, which in turn impact human health and climate.

Chemistry↗
NASA NTRS2000

Directionally Solidified Ceramics Produced

Produced Multiphase, interpenetrating structures are an alternative route to obtaining structural ceramic materials with adequate strength, toughness, and stability for high-temperature aerospace applications. The eutectic architecture, a continuous-reinforcing phase within a higher volume phase or matrix, can be described as a naturally occurring, in situ composite. The phases of a eutectic are thermodynamically compatible at high homologous temperatures. Strong and stable materials have been produced. Toughness, however, remains a technical obstacle. The potential for producing materials with enhanced toughness along with adequate strength and stability was demonstrated using the laser-heated float zone (LHFZ) growth method at the NASA Glenn Research Center at Lewis Field. LHFZ growth at Glenn provides a means to efficiently produce and record the underlying growth phenomena associated with two-phase structures. To initiate directional solidification, a seed of single-crystal sapphire (<0001> direction) was lowered onto the molten liquid until wetting occurred and then withdrawn at a constant rate. Neither the crystal nor the source rod was rotated. The materials produced were tested mechanically in tension, and the resulting microstructure was examined with a scanning electron microscope. Both the inherent properties of the constituent phases and the properties of the interface between them affect the mechanical behavior and the fracture surfaces. The following scanning electron micrographs show the microstructures of two different materials that were tested to failure in tension. In the left micrograph, the flat fracture surface is typical of a material that is strong but has low toughness. In the right micrograph, the crack is effectively deflected at the interface between the two phases, achieving higher toughness at moderately lower strength levels. Conducting mechanical tests to determine the high temperature properties of these materials is the next step in determining their eventual suitability.

Farmer, Serene C.↗
NASA NTRS1990

Linear-phase approximation in the triangular facet near-field physical optics computer program

Analyses of reflector antenna surfaces use a computer program based on a discrete approximation of the radiation integral. The calculation replaces the actual surface with a triangular facet representation; the physical optics current is assumed to be constant over each facet. Described here is a method of calculation using linear-phase approximation of the surface currents of parabolas, ellipses, and shaped subreflectors and compares results with a previous program that used a constant-phase approximation of the triangular facets. The results show that the linear-phase approximation is a significant improvement over the constant-phase approximation, and enables computation of 100 to 1,000 lambda reflectors within a reasonable time on a Cray computer.

Imbriale, W. A.↗
NASA NTRS1995

Electronic Compensation For Distortion Of Antenna Reflector

Proposed method of obtaining approximately desired radiation or reception pattern from antenna that includes reflector based on concept of superposition of electromagnetic fields generated by multiple feedhorns or feed antenna elements arrayed at various positions near reflector and excited at electronically adjustable magnitudes and phases. In intended application, reflector nominally paraboloidal, feed elements N feedhorns in hexagonal array, and method used to compensate for deviations of real reflector surface from nominal paraboloidal shape. Method and concept also applicable to electronic beam steering and electronic antenna compensation in other situations.

Bailey, M. C.↗
DOE OSTI2025

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

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

Aluminum alloy↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2009

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗
DOE OSTI2021

Inter-well connectivity detection in CO 2 WAG projects using statistical recurrent unit models

Routine well-wise injection and production measurements contain significant information on subsurface structure and properties. Data-driven technology that interprets surface data into subsurface structure or properties can assist operators in making informed decisions by providing a better understanding of field assets. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO 2 EOR projects utilizing the water-alternating-gas (WAG) process. SRU is a special type of recurrent neural network (RNN) that allows for better characterization of temporal trends, by learning various statistics of the input at different time scales. In our application, the complete states (injection rate, pressure and cumulative injection) at injectors and pressure states at producers are fed to SRU as the input and the phase rates at producers are treated as the output. Once the SRU is trained and validated, it is then used to assess the connectivity of each injector to any producer using permutation variable importance method, wherein inputs corresponding to an injector are shuffled and the increase in prediction error at a given producer is recorded as the importance (connectivity metric) of the injector to the producer. This method is tested in both synthetic and field-scale cases. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. This significantly improves confidence in our data-driven procedure. The novelty of this work is that it is purely data-driven method and can directly interpret routine surface measurements to intuitive subsurface knowledge. Furthermore, the streamline-based validation procedure provides physics-based backing to the results obtained from data analytics. This study results in a reliable and efficient data analytics framework that is well-suited for large field applications.

42 ENGINEERING↗
DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗