Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Electro-Mechanical Simulation of a Large Aperture MOEMS Fabry-Perot Tunable Filter

We are developing a micro-machined electrostatically actuated Fabry-Perot tunable filter with a large clear aperture for application in high through-put wide-field imaging spectroscopy and lidar systems. In the first phase of this effort, we are developing key components based on coupled electro-mechanical simulations. In particular, the movable etalon plate design leverages high coating stresses to yield a flat surface in drum-head tension over a large diameter (12.5 mm). In this approach, the cylindrical silicon movable plate is back etched, resulting in an optically coated membrane that is suspended from a thick silicon support ring. Understanding the interaction between the support ring, suspended membrane, and coating is critical to developing surfaces that are flat to within stringent etalon requirements. In this work, we present the simulations used to develop the movable plate, spring suspension system, and electrostatic actuation mechanism. We also present results from tests of fabricated proof of concept components.

Kuhn, Jonathan L.↗
DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗
NASA NTRS2015

Next-Generation Ion Thruster Design Tool

Computational tools that accurately predict the performance of electric propulsion devices are highly desirable and beneficial to NASA and the broader electric propulsion community. The current state of the art in electric propulsion modeling relies heavily on empirical data and numerous computational "knobs." In Phase I of this project, Tech-X Corporation developed the most detailed ion engine discharge chamber model that currently exists. This kinetic model simulates all particles in the discharge chamber along with a physically correct simulation of the electric fields. In addition, kinetic erosion models are included for modeling the ion-impingement effects on thruster component erosion. In Phase II, Tech-X developed a user-friendly computer program for NASA and other governmental and industry customers. Tech-X has implemented a number of advanced numerical routines to bring the computational time down to a commercially acceptable level. NASA now has a highly sophisticated, user-friendly ion engine discharge chamber modeling tool.

Stolz, Peter↗
DOE OSTI2022

Simulating the non-monotonic strain response of nanoporous multiferroic composites under electric field control

In this work, we simulate and analyze the mechanical response of a class of multiferroic materials consisting of a templated porous nanostructure made out of cobalt ferrite (CFO) partially filled by atomic layer deposition (ALD) with a ferroelectric phase of lead zirconate titanate (PZT). The strain in the device is measured when an electric field is applied for varying ALD thicknesses, displaying a non-monotonic dependence with a maximum strain achieved for a coating thickness of 3 nm. To understand this behavior, we apply finite element modeling to the smallest repeatable unit of the nanoporous template and simulate the mechanical response as a function of PZT coating thickness. We find that this non-monotonic response is caused by the interplay between two driving forces opposing one another. First, increased porosity works toward increasing the strain due to a reduced system stiffness. Second, decreased porosity involves a larger mass fraction of PZT, which drives the electro-mechanical response of the structure, thus leading to a larger strain. The balance between these two driving forces is controlled by the shear coupling at the CFO/PZT interface and the effective PZT cross section along the direction of the applied electric field. Here, our numerical results show that considering a nonlinear piezoelectric response for PZT leads to an improved agreement with the experimental data, consistent with ex situ poling of the nanostructure prior to magnetic measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantum phase transitions in long-range interacting hyperuniform spin chains in a transverse field

Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of the interesting cases of disordered hyperuniformity are provided by complex many-body systems such as liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. Such spin-chain systems are ideal models for exploring the effects of quantum mechanics on hyperuniformity. It is well-known that the transverse field Ising model shows a quantum phase transition (QPT) at zero temperature. Under the quantum effects of a transverse magnetic field, classical hyperuniform spin chains are expected to lose their hyperuniformity. High-precision simulations of these cases are complicated because of the presence of highly nontrivial long-range interactions. We perform an extensive analysis of these systems using density matrix renormalization group simulations to study the possibilities of phase transitions and the mechanism by which they lose hyperuniformity. Even for a spin chain of length 30, we see discontinuous changes in properties like the “τ order metric” of the ground state, the measure of hyperuniformity, and the second cumulant of the total magnetization along the x-direction, all suggestive of first-order QPTs. An interesting feature of the phase transitions in these disordered hyperuniform spin chains is that, depending on the parameter values, the presence of a transverse magnetic field may lead remarkably to an increase in the order of the ground state as measured by the “τ order metric,” even if hyperuniformity is lost. Therefore, it would be possible to design materials to target specific novel quantum behaviors in the presence of a transverse magnetic field. Our numerical investigations suggest that these spin chains can show no more than two QPTs. We further analyze the long-range interacting spin chains via the Jordan-Wigner mapping onto a system of spinless fermions, showing that under the pairwise-interaction approximation and a mean-field treatment, there can be at most two quantum phase transitions. Here, based on these numerical and theoretical explorations, we conjecture that for spin chains with long-range pair interactions that have convergent cosine transforms, there can be a maximum of two quantum phase transitions at zero temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the “many-body energy” (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. In this study, MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

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

Multidisciplinary design optimization in Architecture, Engineering, and Construction: a detailed review and call for collaboration

The design of buildings has become a complex and multidisciplinary problem involving multiple conflicting objectives as architects and designers address competing technical, economic, environmental, and societal concerns. This has been driving research in Architecture, Engineering, and Construction (AEC) toward rigorous multidisciplinary decision-making frameworks that generate and evaluate numerous design alternatives using multi-objective optimization in concert with simulation and analysis models of varying fidelity and computational expense. While such frameworks are well known and widely employed in the aerospace and systems engineering domains, efforts by design professionals and researchers in the AEC field are scattered at best. In this paper, we provide a detailed review of recent developments in optimization frameworks in the AEC field and subsequently highlight how such developments are largely compartmentalized into separate domains such as structural, energy, daylighting, and other performance factors. Herein, we further discuss the technical challenges involved in concurrent coupled multidisciplinary design optimization (MDO) in the AEC field such as interoperability issues between Building Information Modeling (BIM) environments, analysis/simulation tools, and optimization frameworks. We conclude by outlining research needed for more unified modeling and simulation-based optimization frameworks to aid in complex and multidisciplinary building design processes. We also highlight the need for the identification and development of multi-fidelity simulation tools for use across multiple design phases. As such, this paper contributes a novel roadmap to leverage aerospace and systems engineering research in MDO into the field of AEC, along with a call for researchers in the MDO community to seek collaborations in AEC field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Emergent electric field induced by dissipative sliding dynamics of domain walls in a Weyl magnet

The dynamic motion of topological defects in magnets induces an emergent electric field, as exemplified by the continuous flow of skyrmion vortices. However, the electrodynamics underlying this emergent field remains poorly understood. In this context, magnetic domain walls—one-dimensional topological defects with two collective modes, sliding and spin-tilt—offer a promising platform for exploration. Here we demonstrate that the dissipative motion of domain walls under oscillatory current excitation generates an emergent electric field. We image domain patterns and quantify the domain-wall length under applied magnetic fields in mesoscopic devices based on the magnetic Weyl semimetal NdAlSi. These devices exhibit exceptionally strong domain-wall scattering and a pronounced emergent electric field, as observed in the imaginary component of the complex impedance. Spin dynamics simulations reveal that domain-wall sliding dominates over spin-tilting, in which the phase delay of the domain-wall motion with respect to the driving force impacts the emergent electric field. Our findings establish domain-wall dynamics as a platform for studying emergent electromagnetic fields and motivate further investigations of the coupled motion of magnetic solitons and conduction electrons.

Yamada, Rinsuke [The University of Tokyo, Japan]↗
DOE OSTI2025

MCCCS-MN

The MCCCS‒MN (Monte Carlo for Complex Chemical Systems‒Minnesota) software is developed by the Siepmann research group at the University of Minnesota. MCCCS‒MN allows for the simulation of multi-component molecular systems in the canonical, isobaric-isothermal (including constant stress for solids), grand-canonical, semi-grand, and Gibbs (NVT, NPT, and more than two simulation boxes) ensembles. It uses the configurational-bias Monte Carlo method to efficiently sample phase space for linear, branched and cyclic chain molecules, the adiabatic nuclear and electronic sampling Monte Carlo method to treat many-body polarization effects, and the aggregation-volume-bias Monte Carlo algorithm to efficiently sample the spatial distribution of associating molecules. MCCCS-MN employs a molecular representation of the system where force fields contain bonded and non-bonded terms. Funding for the development of MCCCS-MN through grants from the National Science Foundation (simulation of fluid phase equilibria and chromatography) and the Department of Energy (simulation of adsorption equilibria) is gratefully acknowledged.

Siepmann, J.Ilja [University of Minnesota - Twin C↗
DOE OSTI2020

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: Presentation to DOE-NE [Slides]

The research goal of this project is to develop at Idaho National Laboratory (INL) a first-of-its-kind Multiphysics Object-Oriented Simulation Environment (MOOSE)-based, multiscale, multiphysics spark plasma sintering (SPS) modeling and simulation code application, termed “Freya.” Freya will simulate the thermo-mechanical-electrical aspects of the SPS fabrication process and will be paired with lower length scale sub-models, such as phase-field, to predict the resulting microstructure. SPS is an advanced manufacturing process that can be used to solve a variety of material manufacturing challenges; however, this process is an extremely challenging problem for modeling and simulation. The SPS process is inherently multiphysics and multi-scale, with the macroscale electro-thermo-mechanical behavior linked intricately to the microstructure evolution of the part being sintered. Accurate modeling and simulation tools, specifically geared towards the SPS process, are needed to predict the influence of the multiple variables involved in the manufacturing process. Modeling and simulation accuracy is achieved and demonstrated through comparison to multiple validation experiments. The validation efforts for Freya include both separate effects and complete multiphysics SPS process experiments. One of the key benefits this Laboratory Directed Research & Development (LDRD) project offers stems from the emphasis placed on experimental validation of the Freya models, both on the individual length scales and of the final coupled multiscale multiphysics simulations. Experimental validation of Freya’s multiscale coupling capability provides the technical credibility necessary for potential future industry and research partners to accept the simulation predictions.

36 MATERIALS SCIENCE↗
NASA NTRS1973

Heating and flow-field studies on a straight-wing hypersonic reentry vehicle at angles of attack from 20 to 80 deg with simulation of real-gas trends

Data are presented from a series of phase-change heat transfer and flow visualization tests at Mach 7.4, 8, and 10.3 in air, Mach 19.5 in nitrogen, Mach 20.3 in helium, and Mach 6 in tetrafluoromethane (CF4) on the windward surface of a straight wing hypersonic reentry configuration for angles of attack from 20 deg to 80 deg. The results indicate that: (1) for hypersonic stream Mach numbers, the flow field over the straight-wing configuration is essentially independent of Mach number, (2) transition Reynolds number decreases with increasing angle of attack, (3) at some critical angle of attack, the wing-shock standoff distance is greatly increased and the stagnation line moves downstream from the wing leading edge, (4) value of the critical angle of attack is very sensitive to the flow shock density ratio or effective gamma, and (5) at angles of attack above the critical value for all gases, the nondimensional level of heat transfer to the wing is higher for the higher shock density ratio flows.

Hunt, J. L.↗
NASA NTRS2011

Modeling the Superstorm in November 2003

The superstorm on 20.21 November 2003 was the largest geomagnetic storm in solar cycle 23 as measured by Dst, which attained a minimum value of .422 nT. We have simulated this storm to understand how particles originating from the solar wind and ionosphere get access to the magnetosphere and how the subsequent transport and energization processes contribute to the buildup of the ring current. The global electromagnetic configuration and the solar wind H+ distribution are specified by the Lyon-Fedder-Mobarry (LFM) magnetohydrodynamics model. The outflow of H+ and O+ ions from the ionosphere are also considered. Their trajectories in the magnetosphere are followed by a test-particle code. The particle distributions at the inner plasma sheet established by the LFM model and test-particle calculations are then used as boundary conditions for a ring current model. Our simulations reproduce the rapid decrease of Dst during the storm main phase and the fast initial phase of recovery. Shielding in the inner magnetosphere is established at early main phase. This shielding field lasts several hours and then breaks down at late main phase. At the peak of the storm, strong penetration of ions earthward to L shell of 1.5 is revealed in the simulation. It is surprising that O+ is significant but not the dominant species in the ring current in our calculation for this major storm. It is very likely that substorm effects are not well represented in the models and O+ energization is underestimated. Ring current simulation with O+ energy density at the boundary set comparable to Geotail observations produces excellent agreement with the observed symH. As expected in superstorms, ring current O+ is the dominant species over H+ during the main to mid-recovery phase of the storm.

Fok, Mei-Ching↗
DOE OSTI2024

Evidence for field induced quantum spin liquid behavior in a spin-1/2 honeycomb magnet

One of the most important issues in modern condensed matter physics is the realization of fractionalized excitations, such as the Majorana excitations in the Kitaev quantum spin liquid. To this aim, the 3d-based Kitaev material Na 2 Co 2 TeO 6 is a promising candidate whose magnetic phase diagram of B // a* contains a field-induced intermediate magnetically disordered phase within 7.5 T < |B| < 10 T. The experimental observations, including the restoration of the crystalline point group symmetry in the angle-dependent torque and the coexisting magnon excitations and spinon-continuum in the inelastic neutron scattering spectrum, provide strong evidence that this disordered phase is a field induced quantum spin liquid with partially polarized spins. Our variational Monte Carlo simulation with the effective K-J 1 -Γ-Γ'-J 3 model reproduces the experimental data and further supports this conclusion.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Noise transmission characteristics of a large scale composite fuselage model

Results from an experimental test undertaken to study the basic noise transmission characteristics of a realistic, built-up composite fuselage model are presented. The floor-equipped stiffened composite cylinder was exposed to a number of different exterior noise source configurations in a large anechoic chamber. These exterior source configurations included two point sources located in the same plane on opposite sides of the cylinder, a single point source and a propeller simulator. The results indicate that the interior source field is affected strongly by exterior noise source phasing. Sidewall treatment is seen to reduce the overall interior sound pressure levels and dampen dominant acoustic resonances so that other acoustic modes can affect interior noise distribution.

Beyer, Todd B.↗
NASA NTRS1992

STARDUST: A simulation experiment of cosmic dust analogues production in microgravity conditions

The aim, activity, and some preliminary results of the STARDUST program are presented. The condensation of solid materials from the vapor phase is important in several scientific fields such as chemical vapor deposition, air pollution, and the formation of refractory cosmic dust around stars. Conventional studies of refractory grain formation, using high temperature furnace and shock tube techniques, are restricted to short time scales and suffer from buoyancy induced convection that limit their accuracy. In order to simulate more accurately the condensation of refractory grains near stars and to investigate the advantages of performing condensation studies in microgravity conditions, an experimental investigation of vapor phase condensation in microgravity was undertaken. The experimental equipment currently used is reported. The results from the first flight series and particle aggregation modeling efforts are presented.

Ferguson, Frank↗
DOE OSTI2020

Phase-shift speckle-shearing interferometry

We have numerically simulated the process of measuring stress – strain states by the method of speckle-shearing interferometry using the phase-shift technique. A computer model with the possibility of setting its strain and roughness is developed, which includes a model of a diffusely reflecting test object corresponding to the characteristics of a real membrane, as well as a speckle interferometer model that allows speckle interferograms to be obtained for different speckle sizes and angles between interfering beams. The process of reconstructing the object surface topogram from model speckle interferograms by the phase-shift technique is implemented. Using the developed models, a two-dimensional shearogram are obtained, which is a derivative of the strain field of a circular membrane. Comparison of the results of numerical simulation with experimental data shows that the differences (rms deviations) do not exceed 0.02 μm. It is also shown that the error of interferogram reconstruction by the phase-shift technique increases significantly when the test object strains exceed 12 μm. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗