Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Efficient solution methods for the Navier-Stokes equations

Implicit finite difference schemes for solving two-dimensional and three-dimensional Euler and thin layer Navier-Stokes equations are addressed. The methods are demonstrated in fully vectorized codes for a Cray type architecture. The Beam and Warming implicit approximate factorization algorithm in generalized coordinates is used. The methods are either time accurate or accelerated non-time accurate steady state schemes. Acceleration and efficiency modifications such as matrix reduction, diagonalization, and flux split schemes are presented. Two dimensional inviscid and viscous calculations (e.g., airfoils with a deflected spoiler, circulation control airfoils, and unsteady buffeting) and of three dimensional viscous elliptical bodies, exhausting boattails, and generic oblique wing computations are discussed.

Pulliam, T. H.↗
NASA NTRS1987

Electrodeposition and codeposition under low gravity/nonconvecting conditions

An experimental electrodeposition system was developed for modeling the behavior of inert particles codepositing in an electroplating matrix under low-gravity conditions. The device consists of a Co-electrodeposition cell operating in a convectionless mode (cathode over anode) and containing polystyrene particles with density approximating that of the electroplating solution. Data were obtained in shielded cells at 1 g, and the experiment was duplicated at 0.01 g on a KC-135 flight. No difference was found between convection-free bench experiments and the 0.01 g KC-135 experiments with 0.0900 M CuSO4. Codeposition experiments using 11.8-micron matched-density polystyrene spheres in 1 M CoSO4 have revealed that a noticeable particle gradient is created as the solution density mismatches in the vicinity of the cathode; gentle stirring was required to maintain a homogeneous particle suspension. Cr3C2 dust, which readily disperses at 1 g, tended to coagulate into spherical globules at 0.01 g, when stirred.

Riley, Clyde↗
NASA NTRS1991

Atomic oxygen erosion of polymer matrix composites

Results are presented on atomic oxygen erosion of circular tubes and flat plates fabricated from different epoxy matrix composites containing, boron and Kevlar fiber reinforcements. The UTIAS experiment was located at station D-12, at 90 deg to the LDEF leading edge which was yawed about approximately 8 to 10 deg relative to the velocity vector. This resulted in enhanced atomic oxygen exposure at D-12, providing a total fluence of about 4 x 10(exp 20) atoms/sq cm. Thickness loss measurements were obtained for both flat plates and cylindrical tubes. The effect of reflected atomic oxygen on erosion yield and surface morphology on back surfaces (not directly exposed to AO) are graphically demonstrated. Erosion yields are also compared with other space flight data for similar materials, and ground-based tests in the UTIAS AO simulator. The effect of angle of incidence on surface morphology, erosion profiles and thickness losses are presented and the results correlated with simulator tests as well.

Tennyson, R. C.↗
NASA NTRS1992

Pulsed single-blow regenerator testing

A pulsed single-blow method has been developed for testing of Stirling regenerator materials performance. The method uses a tubular flow arrangement with a steady gas flow passing through a regenerator matrix sample that packs the flow channel for a short distance. A wire grid heater spanning the gas flow channel is used to heat a plug of gas by approximately 2 K for approximately 350 ms. Foil thermocouples monitor the gas temperature entering and leaving the sample. Data analysis based on a 1D incompressible-flow thermal model allows the extraction of Stanton number. A figure of merit involving heat transfer and pressure drop is used to present results for steel screens and steel felt. The observations show a lower figure of merit for the materials tested than is expected based on correlations obtained by other methods.

Oldson, J. C.↗
NASA NTRS1994

A mathematical model for simulating noise suppression of lined ejectors

A mathematical model containing the essential features embodied in the noise suppression of lined ejectors is presented. Although some simplification of the physics is necessary to render the model mathematically tractable, the current model is the most versatile and technologically advanced at the current time. A system of linearized equations and the boundary conditions governing the sound field are derived starting from the equations of fluid dynamics. A nonreflecting boundary condition is developed. In view of the complex nature of the equations, a parametric study requires the use of numerical techniques and modern computers. A finite element algorithm that solves the differential equations coupled with the boundary condition is then introduced. The numerical method results in a matrix equation with several hundred thousand degrees of freedom that is solved efficiently on a supercomputer. The model is validated by comparing results either with exact solutions or with approximate solutions from other works. In each case, excellent correlations are obtained. The usefulness of the model as an optimization tool and the importance of variable impedance liners as a mechanism for achieving broadband suppression within a lined ejector are demonstrated.

Watson, Willie R.↗
NASA NTRS1995

Thermographic imaging for high-temperature composite materials: A defect detection study

The ability of a thermographic imaging technique for detecting flat-bottom hole defects of various diameters and depths was evaluated in four composite systems (two types of ceramic matrix composites, one metal matrix composite, and one polymer matrix composite) of interest as high-temperature structural materials. The holes ranged from 1 to 13 mm in diameter and 0.1 to 2.5 mm in depth in samples approximately 2-3 mm thick. The thermographic imaging system utilized a scanning mirror optical system and infrared (IR) focusing lens in conjunction with a mercury cadmium telluride infrared detector element to obtain high resolution infrared images. High intensity flash lamps located on the same side as the infrared camera were used to heat the samples. After heating, up to 30 images were sequentially acquired at 70-150 msec intervals. Limits of detectability based on depth and diameter of the flat-bottom holes were defined for each composite material. Ultrasonic and radiographic images of the samples were obtained and compared with the thermographic images.

Roth, Don J.↗
NASA NTRS1994

Mineralogy and Ar-39 - Ar-40 of an old pristine basalt: Thermal history of the HED parent body

Previous investigations of mineral chemistry and Rb-Sr and Sm-Nd ages indicated that clast,84 from eucrite Yamato 75011 had preserved the pristine nature of its initial crystallization during an early stage of the HED parent body. Microscale mineralogy and Ar-39-Ar-40 ages of this clast, however, revealed local disturbance of microtextures and partially reset ages. This evidence suggests that, in addition to initial crystallization and rapid cooling, the Y75011,84 clast experienced shock deformation, reheating of short duration at higher temperature, and brecciation. These characteristics suggest two or more impact events. Fe-rich olivine filling fractures in pyroxene may have been introduced during the accompanying shock fracturing. The inferred Ar-39-Ar-40 degassing ages for Y75011 matrix and clast, 84 are 3.94 +/- 0.04 Ga and 3.98 +/- 0.03 Ga, respectively. The suggested degassing age for a clast from Y790020, believed to be paired with Y75011, is approximately 4.03 Ga, but could be younger. We consider it likely that all three samples experienced a common degassing event 3.95 +/- 0.05 Ga ago, but we cannot rule out two or more events spaced over a approximately 0.1 Ga interval. Higher temperature extractions of the two clast samples show significantly older apparent ages up to approximately 4.5 Ga and suggest that the time/temperature regime of this event was not sufficient to degas Ar totally. Most likely, the K-Ar ages were reset by thermal metamorphism associated with one or more impact events associated with shock fracturing, formation of Fe-rich olivine veins, and/or meteorite brecciation. The pyroxene annealing that commonly occurs in many eucrites is likely to be a much earlier process than the impact-produced textural changes and reset K-Ar ages observed in these meteorites. The existence of mineralogical and chronological evidence for metamorphism in an otherwise pristine eucrite suggests that the HED parent body experienced an extensive degree of early cratering.

Takeda, Hiroshi↗
NASA NTRS1993

On the synthesis of resonance lines in dynamical models of structured hot-star winds

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles are then derived from formal solution integration using this source function. Two more approximate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10% or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Puls, J.↗
NASA NTRS2010

Aero-Heating of Shallow Cavities in Hypersonic Freestream Flow

The purpose of these experiments and analysis was to augment the heating database and tools used for assessment of impact-induced shallow-cavity damage to the thermal protection system of the Space Shuttle Orbiter. The effect of length and depth on the local heating disturbance of rectangular cavities tested at hypersonic freestream conditions has been globally assessed using the two-color phosphor thermography method. These rapid-response experiments were conducted in the Langley 31-Inch Mach 10 Tunnel and were initiated immediately prior to the launch of STS-114, the initial flight in the Space Shuttle Return-To-Flight Program, and continued during the first week of the mission. Previously-designed and numerically-characterized blunted-nose baseline flat plates were used as the test surfaces. Three-dimensional computational predictions of the entire model geometry were used as a check on the design process and the two-dimensional flow assumptions used for the data analysis. The experimental boundary layer state conditions were inferred using the measured heating distributions on a no-cavity test article. Two test plates were developed, each containing 4 equally-spaced spanwise-distributed cavities. The first test plate contained cavities with a constant length-to-depth ratio of 8 with design point depth-to-boundary-layer-thickness ratios of 0.1, 0.2, 0.35, and 0.5. The second test plate contained cavities with a constant design point depth-to-boundary-layer-thickness ratio of 0.35 with length-to-depth ratios of 8, 12, 16, and 20. Cavity design parameters and the test condition matrix were established using the computational predictions. Preliminary results indicate that the floor-averaged Bump Factor (local heating rate nondimensionalized by upstream reference) at the tested conditions is approximately 0.3 with a standard deviation of 0.04 for laminar-in/laminar-out conditions when the cavity length-to-boundary-layer thickness is between 2.5 and 10 and for cavities in the depth-to-boundary-layer-thickness range of 0.3 to 0.8. Over this same range of conditions and parameters, preliminary results also indicate that the maximum Bump Factor on the cavity centerline falls between 2.0 and 2.75, as long as the cavity-exit conditions remain laminar. Cavities with length-to-boundary-layer-thickness ratio less than 2.5 can not be easily classified with this approach and require further analysis.

Everhart, Joel L.↗
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

Benchmark neutrinoless double- β decay matrix elements in a light nucleus

We compute nuclear matrix elements of neutrinoless double-β decay mediated by light Majorana-neutrino exchange in the $\textit{A}$ = 6 system. The goal for this work is to benchmark two many-body approaches, the no-core shell model and the multireference in-medium similarity renormalization group (SRG). We use a SRG-evolved next-to-next-to-next-to-leading chiral order two-body potential for the nuclear interaction, and make the approximation that the isospin is conserved. We compare the results of the two approaches as a function of the cutoff on the many-body basis space. Although differences are seen in the predicted nuclear radii, the ground-state energies and neutrinoless double-β decay matrix elements produced by the two approaches show significant agreement. We discuss the implications for calculations in heavier nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

Thermal Conductivity and Erosion Durability of Composite Two-Phase Air Plasma Sprayed Thermal Barrier Coatings

To enhance efficiency of gas turbines, new thermal barrier coatings (TBCs) must be designed which improve upon the thermal stability limit of 7 wt% yttria stabilized zirconia (7YSZ), approximately 1200 C. This tenant has led to the development of new TBC materials and microstructures capable of improved high temperature performance. This study focused on increasing the erosion durability of cubic zirconia based TBCs, traditionally less durable than the metastable t' zirconia based TBCs. Composite TBC microstructures composed of a low thermal conductivity/high temperature stable cubic Low-k matrix phase and a durable t' Low-k secondary phase were deposited via APS. Monolithic coatings composed of cubic Low-k and t' Low-k were also deposited, in addition to a 7YSZ benchmark. The thermal conductivity and erosion durability were then measured and it was found that both of the Low-k materials have significantly reduced thermal conductivities, with monolithic t' Low-k and cubic Low-k improving upon 7YSZ by approximately 13 and approximately 25%, respectively. The 40 wt% t' Low-k composite (40 wt% t' Low-k - 60 wt% cubic Low-k) showed a approximately 22% reduction in thermal conductivity over 7YSZ, indicating even at high levels, the t' Low-k secondary phase had a minimal impact on thermal in the composite coating. It was observed that a mere 20 wt% t' Low-k phase addition can reduce the erosion of a cubic Low-k matrix phase composite coating by over 37%. Various mixing rules were then investigated to assess this non-linear composite behavior and suggestions were made to further improve erosion durability.

Rare earth zirconia↗
NASA NTRS2004

Modified Coulomb-Dipole Theory for 2e Photoionization

In the light of recent experiment on 2e photoionization of Li near threshold, we have considered a modification of the Coulomb-dipole theory, retaining the basic assumption that the threshold is dominated by asymmetric events in phase space [implies r(sub 1), k(sub 1)) greater than or equal to 2(r(sub 2), k(sub )]. In this region [in a collinear model, 2/r(sub 12) approached + 2/(r(sub 1)+r(sub 2)] the interaction reduces to V(rIsub 1) is greater than or equal to 2r(sub 2) is identically equal to [(-Z/r(sub 2)-(A-1)/r(sub 1)] + [(-2r(sub 2)/r(sub 1 exp 2)] is identically equal to V(sub c)+[V(sub d)]. For two electron emission Z = 2, thus both electrons see a Coulomb potential (V(sub c)) asymptotically, albeit each seeing a different charge. The residual potential (V(sub d)) is dipole in character. Writing the total psi = psi (sub c) + psi(sub d) = delta psi, and noting that. (T+V(sub c)-E)psy(sub c) = 0 and (T+V(sub c))psi(sub d) = 0 can be solved exactly, we find, substituting psi into the complete Schrod. Eq., that delta psi = -(H-E)(exp -1)(V(sub d) psi(sub 0)+V(sub c psi (sub 1). Using the fact that the absolute value of V(sub c) is much more than the absolute value of V(sub d) in almost all of configuration space, we can replace H by H(sub 0) in 9H-E)(exp -1) to obtain an improved approximation psi (improved) = psi(sub c) + psi(sub d) -(H(sub 0)-E)(exp -1) (V(sub c) psi (sub 0) + V(sub c) psi(sub 1). Here's the Green's function (H(sub 0)-E)(exp -1), can be exhibited explicitly, but the last term in psi (improved) is small, compared to the first two terms. Inserting them into the transition matrix element, which one handles in the usual way, we obtain in the limit E approaches 0, the threshold law: Q(E) alpha E + M(E)E(exp 5/4) + higher order (Eq. 1a). The modulation function, M(E), is a well-defined (but very non-trivial integral, but it is expected to be well approximated by a sinusoidal function containing a dipole phase term (M(E) = c sin[alpha log (E) + micron] (Eq. 1b). Experimental results show definite modulations, and are well fitted by Eqs (1).

Source record↗
DOE OSTI2024

Investigating sol–gel matrix loading capacity toward producing surrogate nuclear explosive debris with realistic composition

Post detonation nuclear forensic materials which resemble the size, color, elemental composition, and radionuclide content of real nuclear debris would be valuable for developing and validating new nuclear forensic techniques. As nuclear fallout types vary significantly, the ability to tailor each of these parameters accurately is desired to produce materials capable of testing analytical methods under a wide array of forensic scenarios. Sol–gel synthesis techniques can provide tunability of size, shape and composition for producing a wide variety of solid nuclear forensics benchmarking materials. Further, the sol–gel process consists of forming a metal oxide material, often silica, through polymerization of a metal-alkoxy precursor. In this work, we characterize the ability to load sol–gel particles with secondary elemental components such as iron, aluminum, and calcium toward producing benchmarking materials approximating the elemental composition of historic nuclear debris from the Nevada National Security Site. We also demonstrate quantitative radionuclide encapsulation toward producing benchmarking materials with controllable radionuclide content. Finally, we employ these techniques to produce nuclear debris benchmarking materials with controllable elemental matrix composition and radionuclide content and compare these samples with the composition of a historic fallout sample previously reported from the Nevada National Security Site.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Properties of silicon carbide fiber-reinforced silicon nitride matrix composites

The mechanical properties of NASA Lewis developed SiC/RBSN composites and their thermal and environmental stability havd been studied. The composites consist of nearly 30 vol pct of aligned 142 micron diameter chemically vapor-deposited SiC fibers in a relatively porous silicon nitride matrix. In the as-fabricated condition, the unidirectional and 2-D composites exhibited metal-like stress-strain behavior, graceful failure, and showed improved properties when compared with unreinforced matrix of comparable density. Furthermore, the measured room temperature tensile properties were relativley independent of tested volume and were unaffected by artifical notches normal to the loading direction or by thermal shocking from temperatures up to 800 C. The four-point bend strength data measured as a function of temperature to 1400 C in air showed that as-fabricated strength was maintained to 1200 C. At 1400 C, however, nearly 15 pct loss in strength was observed. Measurement of room temperature tensile strength after 100 hr exposure at temperatures to 1400 C in a nitrogen environment indicated no loss from the as-fabricated composite strength. On the other hand, after 100 hr exposure in flowing oxygen at 1200 and 1400 C, the composites showed approximately 40 pct loss from their as-fabricated ultimate tensile strength. Those exposed between 400 to 1200 C showed nearly 60 pct strength loss. Oxidation of the fiber/matrix interface as well as internal oxidation of the porous Si3N4 matrix are likely mechanisms for strength degradation. The excellent strength reproducibility, notch insensitivity, and high temperature strength of the composite makes it an ideal candidate for advanced heat engine applications provided coating or densification methods are developed to avoid internal oxidation attack.

Bhatt, Ramakrishna T.↗
DOE OSTI2024

ZeroCAL: Eliminating Carbon Dioxide Emissions from Limestone’s Decomposition to Decarbonize Cement Production

Limestone (calcite, CaCO 3 ) is an abundant and cost-effective source of calcium oxide (CaO) for cement and lime production. However, the thermochemical decomposition of limestone (~800 °C, 1 bar) to produce lime (CaO) results in substantial carbon dioxide (CO 2(g) ) emissions and energy use, i.e., ~1 tonne [t] of CO 2 and ~1.4 MWh per t of CaO produced. Here, we describe a new pathway to use CaCO 3 as a Ca source to make hydrated lime (portlandite, Ca(OH) 2 ) at ambient conditions (p, T) while nearly eliminating process CO 2(g) emissions (as low as 1.5 mol. % of the CO 2 in the precursor CaCO 3 , equivalent to 9 kg of CO 2(g) per t of Ca(OH) 2 ) within an aqueous flowelectrolysis/ pH-swing process that coproduces hydrogen (H 2(g) ) and oxygen (O 2(g) ). Because Ca(OH) 2 is a zero-carbon precursor for cement and lime production, this approach represents a significant advancement in the production of zero-carbon cement. The Zero CArbon Lime (ZeroCAL) process includes dissolution, separation/recovery, and electrolysis stages according to the following steps: (Step 1) chelator (e.g., ethylenediaminetetraacetic acid, EDTA)-promoted dissolution of CaCO 3 and complexation of Ca 2+ under basic (>pH 9) conditions, (Step 2a) Ca enrichment and separation using nanofiltration (NF), which allows separation of the Ca-EDTA complex from the accompanying bicarbonate (HCO 3 – ) species, (Step 2b) acidity-promoted decomplexation of Ca from EDTA, which allows near-complete chelator recovery and the formation of a Ca-enriched stream, and (Step 3) rapid precipitation of Ca(OH) 2 from the Ca-enriched stream using electrolytically produced alkalinity. These reactions can be conducted in a seawater matrix yielding coproducts including hydrochloric acid (HCl) and sodium bicarbonate (NaHCO 3 ), resulting from electrolysis and limestone dissolution, respectively. Careful analysis of the reaction stoichiometries and energy balances indicates that approximately 1.35 t of CaCO 3 , 1.09 t of water, 0.79 t of sodium chloride (NaCl), and ~2 MWh of electrical energy are required to produce 1 t of Ca(OH) 2 , with significant opportunity for process intensification. This approach has major implications for decarbonizing cement production within a paradigm that emphasizes the use of existing cement plants and electrification of industrial operations, while also creating approaches for alkalinity production that enable cost-effective and scalable CO 2 mineralization via Ca(OH) 2 carbonation.

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

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗