SEARCH · Search NASA

Results for “discrete phase model”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Colloidal Clusters and Networks Formed by Oppositely Charged Nanoparticles with Varying Stiffnesses

Colloidal clusters and gels are ubiquitous in science and technology. Particle softness has a strong effect on interparticle interactions; however, our understanding of the role of this factor in the formation of colloidal clusters and gels is only beginning to evolve. Here, we report the results of experimental and simulation studies of the impact of particle softness on the assembly of clusters and networks from mixtures of oppositely charged polymer nanoparticles (NPs). Experiments were performed below or above the polymer glass transition temperature, at which the interaction potential and adhesive forces between the NPs were significantly varied. Hard NPs assembled in fractal clusters that subsequently organized in a kinetically arrested colloidal gel, while soft NPs formed dense precipitating aggregates, due to the NP deformation and the decreased interparticle distance. Importantly, interactions of hard and soft NPs led to the formation of discrete precipitating NP aggregates at a relatively low volume fraction of soft NPs. A phenomenological model was developed for interactions of oppositely charged NPs with varying softnesses. Furthermore, the experimental results were in agreement with molecular dynamics simulations based on the model. This work provides insight on interparticle interactions before, during, and after the formation of hard–hard, hard–soft, and soft–soft contacts and has impact for numerous applications of reversible colloidal gels, including their use as inks for additive manufacturing.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

DOE OSTI2022

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

constrained wave↗

DOE OSTI2022

Perturbation theory under the truncated Wigner approximation: How system-environment entanglement formation drives quantum decoherence

In discrete quantum systems, quantum decoherence is the disappearance of simple phase relations as a result of interactions with an environment. For many applications, the question is not necessarily how to avoid (inevitable) system-environment interactions, but rather how to design environments that optimally preserve a system's phase relations in spite of such interactions. The formation of system-environment entanglement is a major driving mechanism for decoherence, and a detailed understanding of this process could inform strategies for conserving coherence optimally. This requires scalable, flexible, and systematically improvable quantum dynamical methods that retain detailed information about the entanglement properties of the environment, yet very few current methods offer this combination of features. Here, in this work, we address this need by introducing a theoretical framework wherein we combine the truncated Wigner approximation with standard time-dependent perturbation theory, allowing for computing expectation values of operators in the combined system-environment Hilbert space. We demonstrate the utility of this framework by applying it to the spin-boson model, representative of qubits and simple donor-acceptor systems. For this model, our framework provides an analytical description of perturbative contributions to expectation values. We monitor how quantum decoherence at zero temperature is accompanied by entanglement formation with individual environmental degrees of freedom. Based on this entanglement behavior, we find that the selective suppression of low-frequency environmental modes is particularly effective for mitigating quantum decoherence.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Direct recovery of interfacial topography from coherent X-ray reflectivity: model calculations for a 1D interface

The use of coherent X-ray reflectivity to recover interfacial topography is described using model calculations for a one-dimensional interface. The results reveal that the illuminated topography can be recovered directly from the measured reflected intensities. This is achieved through an analysis of the Patterson function, the Fourier transform of the scattering intensity (as a function of lateral momentum transfer, Q // , at fixed vertical momentum transfer, Q z ). Specifically, a second-order Patterson function is defined that reveals the discrete set of separations and contrast factors (i.e., the product of changes in the effective scattering factor) associated with discontinuities in the effective interfacial topography. Finally, it is shown that the topography is significantly over-determined by the measurements, and an algorithm is described that recovers the actual topography through a deterministic sorting of this information.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Physical, numerical, and computational challenges of modeling neutrino transport in core-collapse supernovae

The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributions or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Realistic simulations of spin squeezing and cooperative coupling effects in large ensembles of interacting two-level systems

We describe an efficient numerical method for simulating the dynamics of interacting spin ensembles in the presence of dephasing and decay. The method builds on the discrete truncated Wigner approximation for isolated systems, which combines the mean-field dynamics of a spin ensemble with a Monte Carlo sampling of discrete initial spin values to account for quantum correlations. Here we show how this approach can be generalized for dissipative spin systems by replacing the deterministic mean-field evolution by a stochastic process, which describes the decay of coherences and populations while preserving the length of each spin. We demonstrate the application of this technique for simulating nonclassical spin-squeezing effects or the dynamics and steady states of cavity QED models with 10 5 interacting two-level systems. This opens up the possibility to perform accurate real-scale simulations of a diverse range of experiments in quantum optics or with solid-state spin ensembles under realistic laboratory conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

A Fast Butterfly-Compressed Hadamard–Babich Integrator for High-Frequency Helmholtz Equations in Inhomogeneous Media with Arbitrary Sources

Here we present a butterfly-compressed representation of the Hadamard-Babich (HB) ansatz for the Green's function of the high-frequency Helmholtz equation in smooth inhomogeneous media. For a computational domain discretized with Nv discretization cells, the proposed algorithm first solves and tabulates the phase and HB coefficients via eikonal and transport equations with observation points and point sources located at the Chebyshev nodes using a set of much coarser computation grids, and then butterfly compresses the resulting HB interactions from all Nv cell centers to each other. The overall CPU time and memory requirement scale as O(Nv log2 Nv) for any bounded two-dimensional (2D) domains with arbitrary excitation sources. A direct extension of this scheme to bounded 3D domains yields an O(Nv4/3) CPU complexity, which can be further reduced to quasi-linear complexities with proposed remedies. The scheme can also efficiently handle scattering problems involving inclusions in inhomogeneous media. Although the current construction of our HB integrator does not accommodate caustics, the resulting HB integrator itself can be applied to certain sources, such as concave-shaped sources, to produce caustic effects. Compared to finite-difference frequency domain methods, the proposed HB integrator is free of numerical dispersion and requires fewer discretization points per wavelength. As a result, it can solve wave propagation problems well beyond the capability of existing solvers. Remarkably, the proposed scheme can accurately model wave propagation in 2D domains with 640 wavelengths per direction and in 3D domains with 54 wavelengths per direction on a state-of-the-art supercomputer at Lawrence Berkeley National Laboratory.

Hadamard--Babich ansatz↗

DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗

DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗

DOE OSTI2021

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Wavelength Scaling of Widely-Tunable Terahertz Quantum-Cascade Metasurface Lasers

Terahertz (THz) external-cavity lasers based on quantum-cascade (QC) metasurfaces are emerging as widely-tunable, single-mode sources with the potential to cover the 1--6 THz range in discrete bands with milliwatt-level output power. By operating on an ultra-short cavity with a length on the order of the wavelength, the QC vertical-external-cavity surface-emitting-laser (VECSEL) architecture enables continuous, broadband tuning while producing high quality beam patterns and scalable power output. The methods and challenges for designing the metasurface at different frequencies are discussed. As the QC-VECSEL is scaled below 2 THz, the primary challenges are reduced gain from the QC active region, increased metasurface quality factor and its effect on tunable bandwidth, and larger power consumption due to a correspondingly scaled metasurface area. At frequencies above 4.5 THz, challenges arise from a reduced metasurface quality factor and the excess absorption that occurs from proximity to the Reststrahlen band. The results of four different devices — with center frequencies 1.8 THz, 2.8 THz, 3.5 THz, and 4.5 THz — are reported. Each device demonstrated at least 200 GHz of continuous single-mode tuning, with the largest being 650 GHz around 3.5 THz. The limitations of the tuning range are well modeled by a Fabry-Pérot cavity which accounts for the reflection phase of the metasurface and the effect of the metasurface quality factor on laser threshold. Lastly, the effect of different output couplers on device performance is studied, demonstrating a significant trade-off between the slope efficiency and tuning bandwidth.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Profusion of symmetry-protected qubits from stable ergodicity breaking

We show how combining a discrete symmetry with topological Hilbert space fragmentation can give rise to exponentially many topologically stable qubits protected by a single discrete symmetry. We illustrate this explicitly with the example of the CZ𝑝 model, where the encoded qubits are prethermally stable to arbitrary symmetry-respecting perturbations for parametrically long times, substantially enhancing the robustness of a recently proposed construction based on nontopological fragmentation. In this model, the encoded qubits naturally come in pairs for which a universal set of transversal logical gates can be performed, ruling out (by the Eastin-Knill theorem) the possibility of using them for quantum error correction. We also comment on the combination of symmetry enrichment and topological fragmentation more generally, and the implications for use of systems exhibiting Hilbert space fragmentation as quantum memories.

kinetically constrained models↗

DOE OSTI2022

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗

DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗

DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗