Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed point 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Loop current fluctuations and quantum critical transport

We study electrical transport at quantum critical points (QCPs) associated with loop current ordering in a metal, focusing specifically on models of the “Hertz-Millis” type. At the infrared (IR) fixed point and in the absence of disorder, the simplest such models have infinite DC conductivity and zero incoherent conductivity at nonzero frequencies. However, we find that a particular deformation, involving N N species of bosons and fermions with random couplings in flavor space, admits a finite incoherent, frequency-dependent conductivity at the IR fixed point, \sigma(\omega>0)\sim\omega^{-2/z} σ ( ω > 0 ) ∼ ω − 2 / z , where z z is the boson dynamical exponent. Leveraging the non-perturbative structure of quantum anomalies, we develop a powerful calculational method for transport. The resulting "anomaly-assisted large N N expansion" allows us to extract the conductivity systematically. Although our results imply that such random-flavor models are problematic as a description of the physical N = 1 N = 1 system, they serve to illustrate some general conditions for quantum critical transport as well as the anomaly-assisted calculational methods. In addition, we revisit an old result that irrelevant operators generate a frequency-dependent conductivity, \sigma(\omega>0) \sim \omega^{-2(z-2)/z} σ ( ω > 0 ) ∼ ω − 2 ( z − 2 ) / z , in problems of this kind. We show explicitly, within the scope of the original calculation, that this result does not hold for any order parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

(U) Using the MCNP6 Perturbation Capability for Source Nuclide Density Sensitivities

The Taylor series (differential operator) perturbation method, as implemented in the PERT capability in MCNP6 (Ref. 1), can be used for first- and second-order sensitivity analyses in fixed-source problems. Unfortunately, the PERT capability does not yet work with point-detector (F5) or pulse-height (F8) tallies. Also, the second-order Taylor series term lacks the cross terms that would allow the calculation of mixed sensitivities. Finally, the PERT capability does not compute any sensitivities involving the source in a fixed-source problem. When neutrons or gamma rays are emitted due to radioactive decay, the source rate density depends intrinsically on the nuclide densities in the source material. When the PERT capability in MCNP6 is used to compute sensitivities of a response to nuclide and material densities, those sensitivities do not include the effect of the source.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Robust Automatic EXAFS First-Shell Fits

Extended X-ray absorption fine structure (EXAFS) is a widely used technique for atomic structure determination. Fourier transformation connects EXAFS in k space and R space. However, determining the appropriate k-range for the transformation can be challenging, but critical for the first-shell fit. In this study, we present an automatic method to determine the k-range using the Larch package and a Python program. The first step is to estimate spectral noise across a series of k-ranges with a fixed minimum value and identify the optimal maximum value in the k-range (k max ). The k max is determined by an empirical noise threshold that marks the point where the noise level in the Fourier transformed spectrum changes dramatically. Using the obtained k max value, the first shell is modeled to determine the minimum k value (k min ) by optimizing the background function through alignment of the spectrum with theory. The optimal k min corresponds to the point of the minimum R-factor, which quantifies the difference between the experimental and fitted spectrum. Our method was tested on various typical datasets and yielded suitable k-ranges for Fourier transformation and accurate first-shell fits. This approach helps avoid unreliable, irreproducible data analysis, especially for noisy data from diluted samples, and enables robust automatic first-shell EXAFS fitting.

EXAFS analysis↗
DOE OSTI2023

Discordant Clinical and Microbiological Outcomes Are Associated With Late Clinical Relapse in Clinical Trials for Complicated Urinary Tract Infections

Abstract Background Current US Food and Drug Administration guidance recommends that the primary endpoint for complicated urinary tract infection clinical trials be a composite of the clinical and microbiological responses, assessed at a fixed point after therapy. Although some participants meet the criteria for clinical success, they do not meet the criteria for microbiological eradication and are classified as treatment failures. These discordant outcomes have raised questions about the utility of the microbiological endpoint. Methods We analyzed participant data from 13 phase 3 clinical trials submitted to the US Food and Drug Administration (N = 4842). Outcomes were determined at the test of cure (TOC) visit, recommended to occur at least 5 days after therapy and at the late follow-up visit, recommended to occur 21 to 28 days after randomization. Clinical and microbiological success were defined as the resolution of complicated urinary tract infection symptoms present at study entry, with no new symptoms (clinical cure), and a reduction in density of the original pathogen to <103 CFU/mL on urine culture (microbiological eradication). Results Among included participants, 70.7% were concordant successes at the TOC visit, 18.0% were discordant failures (clinical cure/microbiological persistence), and 6.7% were concordant failures (clinical failure/microbiological persistence). Discordant participants were at an increased risk for clinical failure at the late follow-up visit, and the risk of late clinical failure increased with time. Conclusions Discordant clinical and microbiological outcomes at the TOC visit were associated with an increased risk of late clinical failure. Microbiological outcomes appear to be an important component of the endpoint.

Immunology↗
DOE OSTI2020

Linear point and sound horizon as purely geometric standard rulers

The baryon acoustic oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the purely geometric BAO methods. To date two different purely geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable (e.g., the BAO). We define the linear point and sound horizon and we fully characterize and compare the two rulers’ cosmological-parameter dependence. We find that they are both geometrical (i.e., independent of the primordial cosmological parameters) within the required accuracy, and that they have the same parameter dependence for a wide range of parameter values. We estimate the rulers’ best-fit values and errors, given the cosmological constraints obtained by the Planck Satellite team from their measurements of the cosmic microwave background temperature and polarization anisotropies. We do this for three different cosmological models encompassed by the purely geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by 0.5σ when we do not fix the spatial geometry to be flat in Λ CDM. Finally, this points toward a sensitivity of the rulers to different cosmological models when they are estimated from the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2020

Bounded motion design in the Earth zonal problem using differential algebra based normal form methods

Establishing long-term relative bounded motion between orbits in perturbed dynamics is a key challenge in astrodynamics to enable cluster flight with minimum propellant expenditure. Herein, we present an approach that allows for the design of long-term relative bounded motion considering a zonal gravitational model. Entire sets of orbits are obtained via high-order Taylor expansions of Poincarè return maps about reference fixed points. The high-order normal form algorithm is used to determine a change in expansion variables of the map into normal form space, in which the phase space behavior is circular and can be easily parameterized by action–angle coordinates. The action–angle representation of the normal form coordinates is then used to parameterize the original Poincarè return map and average it over a full phase space revolution by a path integral along the angle parameterization. As a result, the averaged nodal period and drift in the ascending node are obtained, for which the bounded motion conditions are straightforwardly imposed. Sets of highly accurate bounded orbits are obtained, extending over several thousand kilometers, and valid for decades.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Reconstructing the reflectivity of liquid surfaces from grazing incidence X-ray off-specular scattering data

The capillary wave model of a liquid surface predicts both the X-ray specular reflection and the diffuse scattering around it. A quantitative method is presented to obtain the X-ray reflectivity (XRR) from a liquid surface through the diffuse scattering data around the specular reflection measured using a grazing incidence X-ray off-specular scattering (GIXOS) geometry at a fixed horizontal offset angle with respect to the plane of incidence. With this approach the entire Q z -dependent reflectivity profile can be obtained at a single, fixed incident angle. This permits a much faster acquisition of the profile than with conventional reflectometry, where the incident angle must be scanned point by point to obtain a Q z -dependent profile. The XRR derived from the GIXOS-measured diffuse scattering, referred to in this paper as pseudo-reflectivity, provides a larger Q z range compared with the reflectivity measured by conventional reflectometry. Transforming the GIXOS-measured diffuse scattering profile to pseudo-XRR opens up the GIXOS method to widely available specular XRR analysis software tools. Here the GIXOS-derived pseudo-XRR is compared with the XRR measured by specular reflectometry from two simple vapor–liquid interfaces at different surface tension, and from a hexadecyltrimethylammonium bromide monolayer on a water surface. For the simple liquids, excellent agreement (beyond 11 orders of magnitude in signal) is found between the two methods, supporting the approach of using GIXOS-measured diffuse scattering to derive reflectivities. Pseudo-XRR obtained at different horizontal offset angles with respect to the plane of incidence yields indistinguishable results, and this supports the robustness of the GIXOS-XRR approach. The pseudo-XRR method can be extended to soft thin films on a liquid surface, and criteria are established for the applicability of the approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Efficient Quantum Gibbs Samplers with Kubo–Martin–Schwinger Detailed Balance Condition

Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Convex Optimization for Nonequilibrium Steady States on a Hybrid Quantum Processor

Finding the transient and steady state properties of open quantum systems is a central problem in various fields of quantum technologies. Here, in this work, we present a quantum-assisted algorithm to determine the steady states of open system dynamics. By reformulating the problem of finding the fixed point of Lindblad dynamics as a feasibility semidefinite program, we bypass several well-known issues with variational quantum approaches to solving for steady states. We demonstrate that our hybrid approach allows us to estimate the steady states of higher dimensional open quantum systems and discuss how our method can find multiple steady states for systems with symmetries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Frequency Domain Reflectometry (FDR) Simulation Techniques for Digital Twin Representation of an Electrical Cable

Simulation of cable system response to frequency domain reflectometry (FDR) tests can be instrumental to understanding these tests and the nature and influence of various cable anomalies on test signatures. Reflectometry simulations are based upon a finite element representation of cable conductors and insulation to produce an S-parameter at each evaluated frequency. The aggregate collection of cable model S-parameters can simulate the influence on a test signal injected into a physical cable. Such an approach was undertaken in this work to produce a digital twin simulation of a low-voltage electrical cable. The electrical cable digital twin examined the influence of test simulation parameters and the relative influence of cable anomalies, including thermal aging, water or moisture exposure, water or moisture ingress, and other anomalies. The digital twin in this work included modeling of the conductors, insulation, jacket, and surrounding environment (air, water, etc.). The digital twin could be expanded to include cable bends, junctions and splices, branch or T systems, and termination impedances of motors or instruments. Observations and conclusions of this work include: 1. Fully 3D digital twin simulation of an electrical cable using an FDR approach is possible. However, there are tradeoffs between simulation fidelity and solution time, which must be balanced to ensure the simulation solves in an adequate amount of time (e.g., less than 20 minutes). Simulation parameters to balance include frequency bandwidth, number of frequencies, mesh density, connection impedance, and permittivity tolerance. 2. The digital twin simulation can explain FDR sensitivity to various cable anomalies, including entry and exit from an oven or water bath. 3. The digital twin simulation FDR response attenuates with distance along the cable and is further affected by the frequency bandwidth, which is similar to that observed with physical measurements. 4. The resolution of the digital twin FDR peaks increased with increasing bandwidth and with increasing number of frequencies, again similar to physical measurements. 5. The presence of multiple anomalies in the digital twin does not substantially attenuate the FDR response to anomalies located beyond the first encountered anomaly and impedance mismatch. 6. Spectral variation of the permittivity did not have a significant effect on the FDR response compared to a fixed nominal value. 7. Extension of the digital twin to 1000 ft still allowed for detection of distal anomalies near the far end of the electrical cable from the instrument connection point. 8. The ARENA test bed facilitates efficient NDE evaluations of well understood cable anomalies with various NDE methods without risking actual plant damage.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Spectral function of fermions in a highly occupied non-Abelian plasma

We develop a method to obtain fermion spectral functions non-perturbatively in a non-Abelian gauge theory with high occupation numbers of gauge fields. After recovering the free field case, we extract the spectral function of fermions in a highly occupied non-Abelian plasma close to its non-thermal fixed point, i.e., in a self-similar regime of the non-equilibrium dynamics. We find good agreement with hard loop perturbation theory for medium-induced masses, dispersion relations and quasiparticle residues. We also extract the full momentum dependence of the damping rate of the collective excitations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Polarization Selectivity of Aloof-Beam Electron Energy-Loss Spectroscopy in One-Dimensional ZnO Nanorods

Orientation-dependent electronic properties of wurtzite zinc oxide nanorods are characterized by aloof-beam electron energy-loss spectroscopy (EELS) carried out in a scanning transmission electron microscope (STEM). The two key crystal orientation differentiating transitions specific to the in-plane (13.0 eV) and out-of-plane (11.2 eV) directions with respect to the wurtzite structure are examined by first-principles density-functional theory calculations. We note some degree of orientation dependence at the onset of direct band gap transition near 3.4 eV. We demonstrate that good polarization selectivity can be achieved by placing the electron probe at different locations around the specimen with increasing impact parameter while keeping the beam-specimen orientation fixed. The observed results are qualitatively elucidated in terms of the perpendicular electric fields generated by the fast electron (60 kV) used in the microscope. The fact that good polarization selectivity can be achieved by aloof-beam EELS without the requirement of sample reorientation is an attractive aspect from the characterization method point of view in the STEM-EELS community.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Real-Time Wave Energy Converter Control Using Instantaneous Frequency

Wave Energy Converters (WECs) rely on effective Power Take-Off (PTO) control strategies to maximize energy absorption under dynamic sea conditions. Traditional hydrodynamic modeling techniques may require computationally intensive convolution calculations, making real-time control implementation challenging. This paper presents an alternative approach by leveraging instantaneous frequency estimation to dynamically adjust PTO damping in response to varying wave frequencies. Two real-time frequency estimation methods are explored: the Hilbert Transform (HT) and Phase-Locked Loop (PLL). The Hilbert Transform method provides accurate frequency tracking but introduces a delayed response due to its dependence on causal data. Conversely, the PLL approach demonstrates strong potential in frequency tracking but requires careful gain tuning, particularly in complex sea states. Comparative evaluations across multiple test cases—including sinusoidal variations, amplitude steps, frequency step changes, and real-world JONSWAP spectrum waves—highlight the strengths and limitations of each method. The two different PTO control techniques across the various frequency estimation methods were tested under real-sea states using a state-space model of a point-absorbing Wave Energy Converter. The Capture Width Ratio (CWR) is used as a performance metric, with results showing that the HT achieves a 10.6% improvement, while the PLL estimation yields a 0.9% improvement relative to the fixed parameter control baseline. These results highlight the effectiveness of real-time frequency estimation in improving energy absorption compared to static control parameters.

WEC control↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2020

Taming the fixed-node error in diffusion Monte Carlo via range separation

By combining density-functional theory (DFT) and wave function theory via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction known as configuration interaction using a perturbative selection made iteratively (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of mu = 0.5 bohr(-1) provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.

Scemama, Anthony↗