Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microscopic 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Deuteration removes quantum dipolar defects from KDP crystals

Abstract Dielectric properties of the hydrogen-bonded ferroelectric crystal KH 2 PO 4 (KDP) differ significantly from those of KD 2 PO 4 (DKDP). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polarization, but a detailed microscopic model is missing. We show that all-atom path integral molecular dynamics simulations can predict the isotope effects, revealing the microscopic mechanism that differentiates KDP and DKDP. Proton tunneling generates phosphate configurations that do not contribute to the polarization. At low temperatures, these quantum dipolar defects are substantial in KDP but negligible in DKDP. These intrinsic defects explain why KDP has lower spontaneous polarization and transition entropy than DKDP. The prominent role of quantum fluctuations in KDP is related to the unusual strength of the hydrogen bonds and should be equally important in other crystals of the KDP family, which exhibit similar isotope effects.

Yang, Bingjia (ORCID:0000000340749553)↗
DOE OSTI2024

Towards Full Field-of-View Fourier Ptychography for Extreme Ultraviolet Microscope

We evaluate various Fourier ptychographic microscopy (FPM) reconstruction algorithms using both simulated and experimental data acquired from an Extreme Ultraviolet (EUV, 13.5 nm wavelength) microscope. We specifically focus on the algorithms' ability to robustly address field-dependent aberrations, which enables increased spatial resolution and quantitative phase imaging across an expanded field of view. We systematically compare the algorithms' performance under aberrations for a single zoneplate imaging system, utilizing Fourier Ring Correlation (FRC) as a systematic metric for assessing reconstruction quality. Furthermore, we explore the impact of systematic errors on the reconstruction of experimental data, aiming to increase the effective field of view by 25-fold, from the nominal 5x5 um2 diffraction-limited area. Additionally, our evaluation incorporates innovative FPM-adjacent methodologies, including the Angular Ptychographic Imaging with Closed-form method (APIC), for reconstructing EUV images.

Gu, Chaoying↗
DOE OSTI2018

Solving macroscopic and microscopic pin-fin heat sink problems by adapted spectral method

In this paper, simulation of heat transfer in a heat sink with macroscopic and microscopic scales when one pin-fin is added to the system is investigated by the proposed spectral method. In the microscopic problems, heat transfer model uses dual-phase lag formulas in contrast with macroscopic problems when Fourier law is used to formulate the governing equation. In macroscopic problem, the results are compared with COMSOL multiphysics software results and a good agreement between the results are shown. In microscopic problems, 3D Gaussian heat source is used and boundary conditions obey the Newton law. Comparisons show the efficiency of the current method, while the results are compared with existed literature. It is shown that: by adding one pin fin to the heat sink problem, temperature of the system will decrease. Numerical results are given and graphs of the temperature for various kinds of pin-fin heat sink problems are depicted.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Control over emergence and alignment of transient blisters in thermo-responsive gels using hierarchically patterned substrates

Surfaces with tunable microscale textures are vital in a large variety of technological applications, including heat transfer, antifouling and adhesion. To facilitate such broad-scale use, there is a need to create surfaces that undergo reconfigurable changes in topology and thus, enable switchable functionality. To date, there is a relative dearth of methods for engineering surfaces that can be actuated to change topography over a range of length scales, and hence, form tunable hierarchically structured layers. Combining modeling and experiments, we design a geometrically patterned, thermo-responsive poly (N-isopropylacrylamide) gel film that undergoes controllable hierarchical changes in topology with changes in temperature. At the bottom, the film is covalently bound to a solid, curved substrate; at the top, the film encompasses longitudinal rectangular ridges that are oriented perpendicular to the underlying cylindrical curves. At temperatures below lower critical solution temperature (LCST), the swollen gel exhibits 3D variations in polymer density and thickness defined by the gel's top and bottom topography. As the temperature rises above LCST, the interplay between the upper ridges and lower curves in the gel drives non-uniform, directional solvent transport, the nucleation and propagation of a phase-separated higher-density skin layer, and the resulting pressure buildup within the film. These different, interacting kinetic processes lead to an instability, which produces transient microscopic blisters in the film. Through simulations, we show how tuning the width of the ridges modifies the propagation of a skin layer and creates localized pressure build-up points, which enables control over the emergence, distribution, and alignment of the microscopic blisters. Futhermore, we provide a simple argument to predict the size of such microscopic features. Experiments confirm our predictions and further highlight how our computational model enables the rational design of topographical transitions in these tunable films. The development of actuatable, hierarchically structured films provides new routes for achieving switchable functionality in actuators, drug release systems and adhesives.

42 ENGINEERING↗
DOE OSTI2015

Atomic friction at exposed and buried graphite step edges: Experiments and simulations

The surfaces of layered materials such as graphite exhibit step edges that affect friction. Step edges can be exposed, where the step occurs at the outmost layer, or buried, where the step is underneath another layer of material. Here, we study friction at exposed and buried step edges on graphite using an atomic force microscope (AFM) and complementary molecular dynamics simulations of the AFM tip apex. Exposed and buried steps exhibit distinct friction behavior, and the friction on either step is affected by the direction of sliding, i.e., moving up or down the step, and the bluntness of the tip. These trends are analyzing in terms of the trajectory of the AFM tip as it moves over the step, which is a convolution of the topography of the surface and the tip shape.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Modeling mode competition in laser diodes

We present an advanced model for the simulation of laser dynamics based on the microscopic many-particle equations of motion and going beyond the often used rate equation model. Our model is used to simulate (Al,In)GaN based edge-emitting laser diodes showing longitudinal mode-competition phenomena which are responsible for the effect of mode hopping, i.e. the modes are switched on and off causing mode rolling from higher to lower frequencies. We show that the effect is caused by beating vibrations of the carrier density which result in an asymmetric coupling term between neighboring modes. Here, the results agree with the case where we assume that the coupling of two modes is independent of all the other modes, as long as the mode frequencies are not equidistant. If the mode frequencies are equidistant, the mode dynamics depend strongly on the initial phases and the initial noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Insights into the structure and dynamics of K + ions at the muscovite–water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. Here, in this work, we make a first step in this direction by investigating the behavior of the native K + ions at the interface between neat water and the muscovite mica (001) surface with ab-initio -based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K + ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K + diffusion between mica’s ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K + ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Ab-initio molecular dynamics↗
DOE OSTI2018

Three-dimensional model of intracellular and intercellular Ca2+ waves propagation in endothelial cells

Highlights: • Image-based 3D simulation model for intercellular and intracellular Ca{sup 2+} waves. • Metabolic reactions were divided into membrane and cytoplasm domains. • Ca{sup 2+} and IP{sub 3} diffused across the cell and gap junctions. • Ca{sup 2+} waves are propagated from focal stimulated point in numerical simulation. • Ca{sup 2+} wave speed was in agreement with those from cell culture experiments. Intracellular and intercellular Ca{sup 2+} waves play key roles in cellular functions, and focal stimulation triggers Ca{sup 2+} wave propagation from stimulation points to neighboring cells, involving localized metabolism reactions and specific diffusion processes. Among these, inositol 1,4,5-trisphosphate (IP{sub 3}) is produced at membranes and diffuses into the cytoplasm to release Ca{sup 2+} from endoplasmic reticulum (ER). In this study, we developed a three-dimensional (3D) simulation model for intercellular and intracellular Ca{sup 2+} waves in endothelial cells (ECs). 3D model of 2 cells was reconstructed from confocal microscopic images and was connected via gap junctions. Cells have membrane and cytoplasm domains, and metabolic reactions were divided into each domain. Finally, the intracellular and intercellular Ca{sup 2+} wave propagations were induced using microscopic stimulation and were compared between numerical simulations and experiments. The experiments showed that initial sharp increases in intracellular Ca{sup 2+} occurred approximately 0.3 s after application of stimuli. In addition, Ca{sup 2+} wave speeds remained constant in cells, with intracellular and intercellular speeds of approximately 35 and 15 μm/s, respectively. Simulations indicated initial increases in Ca{sup 2+} concentrations at points of stimulation, and these were then propagated across stimulated and neighboring cells. In particular, initial rapid increases in intracellular Ca{sup 2+} were delayed and subsequent intracellular and intercellular Ca{sup 2+} wave speeds were approximately 25 and 12 μm/s, respectively. Simulation results were in agreement with those from cell culture experiments, indicating the utility of our 3D model for investigations of intracellular and intercellular messaging in ECs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Interactive Web Application for Traffic Simulation Data Management and Visualization

As traffic simulation software becomes more effective for realistically simulating and analyzing traffic dynamics and vehicle interactions on the mesoscopic and microscopic level, the management, dissemination, and collaborative visualization of traffic simulation results produced by individual transportation planners presents a significant challenge. Existing online content management systems have a very limited capability in allowing users to query specific traffic simulation scenarios and geospatially visualize simulation results through shareable and interactive web interfaces. This paper presents a web-based application for promoting the archiving, sharing, and visualization of large-scale traffic simulation outputs. The application is developed to enhance cyber-physical controls, communications, and public education for collaborative transportation planning. Unique features of the web application include: (a) allowing users to upload their new traffic simulation scenarios (parameters and outputs), as well as search existing scenarios using easily accessible interfaces; (b) optimizing simulation output files with heterogeneous data formats and projected coordinate systems for web-based storage and management using a scalable and searchable data/metadata standard; (c) standardizing user-uploaded simulation outputs using web interfaces and data processing libraries with parallel computing capacity; and (d) providing shareable web visual interfaces for visualizing the traffic flow and signal information stored in simulation outputs (e.g., regional traffic patterns and individual vehicle interactions) and visually comparing multiple simulation outputs both spatially and temporally. Furthermore, the paper presents the conceptual design and implementation of this application, and demonstrates the application’s performance for sharing, comparing, and visualizing simulation outputs from VISSIM and SUMO, two commonly used traffic simulation software programs.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2021

Nature of dynamic gradients, glass formation, and collective effects in ultrathin freestanding films

Significance The origin of spatially heterogeneous alterations in the dynamics of nanoconfined liquids has been a major question in physics, chemistry, and materials science for over 30 y. We combine simulation and theory to establish a unified predictive microscopic understanding in freestanding films at low to medium molecular weight, spanning from thick to ultrathin films. The simulation-based phenomenology of altered dynamics is theoretically predicted and understood based on a combination of coupled near-interface gradients of local caging constraints and interfacial and finite size-induced changes in a collective elastic barrier to motion. No invocation of any realized or incipient critical-like phenomenon is necessary, which may have broad implications for other confined systems and the mechanism of the bulk glass transition.

nanoconfinement↗
DOE OSTI2022

Measurement and Simulation of the Magnetic Fields from a 555 Timer Integrated Circuit Using a Quantum Diamond Microscope and Finite-Element Analysis

Quantum diamond microscope (QDM) magnetic field imaging is an emerging interrogation and diagnostic technique for integrated circuits (ICs). To date, the ICs measured with a QDM have been either too complex for us to predict the expected magnetic fields and benchmark the QDM performance or too simple to be relevant to the IC community. In this paper, we establish a 555 timer IC as a “model system” to optimize QDM measurement implementation, benchmark performance, and assess IC device functionality. To validate the magnetic field images taken with a QDM, we use a spice electronic circuit simulator and finite-element analysis (FEA) to model the magnetic fields from the 555 die for two functional states. Furthermore, we compare the advantages and the results of three IC-diamond measurement methods, confirm that the measured and simulated magnetic images are consistent, identify the magnetic signatures of current paths within the device, and discuss using this model system to advance QDM magnetic imaging as an IC diagnostic tool.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Electric field control of chirality

Polar textures have attracted substantial attention in recent years as a promising analog to spin-based textures in ferromagnets. Here, using optical second-harmonic generation–based circular dichroism, we demonstrate deterministic and reversible control of chirality over mesoscale regions in ferroelectric vortices using an applied electric field. The microscopic origins of the chirality, the pathway during the switching, and the mechanism for electric field control are described theoretically via phase-field modeling and second-principles simulations, and experimentally by examination of the microscopic response of the vortices under an applied field. The emergence of chirality from the combination of nonchiral materials and subsequent control of the handedness with an electric field has far-reaching implications for new electronics based on chirality as a field-controllable order parameter.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng–Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young–Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young–Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young–Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests are carried out to verify the effectiveness of the proposed multi-scale method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Free-Energy Landscape of the Dissolution of Gibbsite at High pH

The individual elementary reactions involved in the dissolution of a solid into solution remain mostly speculative due to a lack of direct experimental probes. In this regard, in this paper we have applied atomistic simulations to map the free-energy landscape of the dissolution of gibbsite from a step edge as a model of metal hydroxide dissolution. The overall reaction combines kink formation and kink propagation. Two individual reactions were found to be rate-limiting for kink formation, that is, the displacement of Al from a step site to a ledge adatom site and its detachment from ledge/terrace adatom sites into the solution. Finally, as a result, a pool of mobile and labile adsorbed species, or adatoms, exists before the release of Al into solution. Because of the quasi-hexagonal symmetry of gibbsite, kink site propagation can occur in multiple directions. Overall, our results will enable the development of microscopic mechanistic models of metal oxide dissolution.

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

Molecular origins of the zeta potential

The zeta potential (ZP) is an oft-reported measure of the macroscopic charge state of solid surfaces and colloidal particles in contact with solvents. However, the origin of this readily measurable parameter has remained divorced from the molecular-level processes governing the underlying electrokinetic phenomena, which limits its usefulness. Here, we connect the macroscopic measure to the microscopic realm through nonequilibrium molecular dynamics simulations of electroosmotic flow between parallel slabs of the hydroxylated (110) rutile (TiO 2 ) surface. These simulations provided streaming mobilities, which were converted to ZP via the commonly used Helmholtz-Smoluchowski equation. A range of rutile surface charge densities (0.1 to –0.4 C/m 2 ), corresponding to pH values between about 2.8 and 9.4, in RbCl, NaCl, and SrCl 2 aqueous solutions, were modeled and compared to experimental ZPs for TiO 2 particle suspensions. Simulated ZPs qualitatively agree with experiment and show that “anomalous” ZP values and inequalities between the point of zero charge derived from electrokinetic versus pH titration measurements both arise from differing co- and counterion sorption affinities. We show that at the molecular level the ZP arises from the delicate interplay of spatially varying dynamics, structure, and electrostatics in a narrow interfacial region within about 15 Å of the surface, even in dilute salt solutions. This contrasts fundamentally with continuum descriptions of such interfaces, which predict the ZP response region to be inversely related to ionic strength. In reality the properties of this interfacial region are dominated by relatively immobile and structured water. Furthermore, viscosity values are substantially greater than in the bulk, and electrostatic potential profiles are oscillatory in nature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗