Search NASA⌕ Search

SEARCH · Search NASA

Results for “Signed Distance”

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 19 records

Exact signed distance fields using parallel Fast Sweeping Method

Signed distance fields are often used in multiphysics simulations to track material interfaces. We present a simple methodology based on the fast sweeping method to generate the exact signed distance from triangular meshes and linear paths on Cartesian grids. The methodology propagates the closest primitive to the boundary to the rest of the domain following the characteristics. A local upwind criterion is used to decide between the new and existing closest primitive at each grid point while capturing the correct sign of the global function. The methodology has optimal computational complexity and runs efficiently in distributed-memory architectures. We include 2D and 3D test cases along with a resolution study up to 0.512 trillion zones and 1,000 computer cores. The solution strategy can also be applied to other types of meshes or collections of primitives.

97 MATHEMATICS AND COMPUTING↗

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗

"Coulombic Viscosity" In Granular Materials: Planetary and Astrophysical Implications

The term "Coulombic viscosity" is introduced here to define an empirically observed phenomenon from experiments conducted in both microgravity, and in ground-based 1-g conditions. In the latter case, a sand attrition device was employed to test the longevity of aeolian materials by creating two intersecting grain-circulation paths or cells that would lead to most of the grain energy being expended on grain-to-grain collisions (simulating dune systems). In the areas in the device where gravitationally-driven grain-slurries recycled the sand, the slurries moved with a boundary-layer impeded motion down the chamber walls. Excessive electrostatic charging of the grains during these experiments was prevented by the use of an a.c. corona (created by a Tesla coil) through which the grains passed on every cycle. This created both positive and negative ions which neutralized the triboelectrically-generated grain charges. When the corona was switched on, the velocity of the wall-attached slurries increased by a factor of two as approximately determined by direct observation. What appeared to be a freely-flowing slurry of grains impeded only by intergranular mechanical friction, had obviously been significantly retarded in its motion by electrostatic forces between the grains; with the charging reduced, the grains were able to move past one another without a flow "viscosity" imposed by the Coulombic intergranular forces. A similar phenomenon was observed during microgravity experiments aboard Space Shuttle in USML-1 & USML-2 spacelabs where freely-suspended clouds of sand were being investigated for their potential to for-m aggregates. In this environment, the grains were also charged electrostatically (by natural processes prior to flight), but were free from the intervention of gravity in their interactions. The grains were dispersed into dense clouds by bursts of air turbulence and allowed to form aggregates as the ballistic and turbulent motions damped out. During this very brief (30-60 sec) damping period, motion of the grains was observed to be retarded by the electrostatic interactions. The fact that the grains almost instantly formed aggregates was evidence that their ballistic motions had been constrained and redirected by the dipole-dipole interactions that led to filamentary aggregate development. Undoubtedly, the "Coulombic viscosity" of the cloud assisted in damping grain motion so rapidly. The electrostatically-induced grain-cloud viscosity or drag exerted on grain motion, is a complex function of three major parameters: charge magnitude, charge sign, and mean intergranular distance. The above experiments illustrate one particular type of granular behavior. The discussion here will therefore be restricted to drag relationships: (a) between grains that are naturally charged triboelectrically and thus exhibit dipole-dipole attractions between one another even if there are slight net charges present (which can be overwhelmed by dipole coupling at short distances), and (b) between grains that are densely spaced where the intergranular distance varies between zero and some value (usually tens or hundreds of grain diameters) that permits each grain to detect the dipole moment of another grain -- the distance is not so great that other grains appears as neutral electrical "singularities. I. Aeolian transport: During motion of grains in a saltation cloud (on Earth, Mars, or Venus), triboelectric charging must occur as a result of multiple grain contacts, and by friction with the entraining air. A situation might develop that is similar to the one described above in the attrition device: grain motion becoming significantly retarded (reduced flux) as grains find it increasingly difficult to either separate from the surface, or to pass one another without Coulombic retarding forces. A "Coulombic drag" will exist at flux initiation and increase with time to work in direct opposition to the aerodynamic drag that drives the grain motion. It is predicted that this will lead to an increase with time of both the aerodynamic and bed-dilatancy thresholds (3). Because of Paschen discharge effects in the martian atmosphere, the electrostatic charging in a saltation cloud may be partially abated, but this will lead to greater grain mobility, more charging, and thus to a charge-discharge steady state mediated by mechanical interactions. II. Dry colluvial systems: Sand avalanches on dunes, dry debris flows, talus flows, avalanches, and pyroclastic surges are examples of gravity-driven, dense granular flows where rock/grain fragmentation and grain-to-grain interactions cause triboelectrification (sometimes augmented by other electrical charging processes), and where the grain densities of the systems are such that strong dipole-dipole interactions between grains might be expected to be present. Because it is expected that the Coulombic forces between grains will cause a sluggishness or enhanced granular-flow viscosity, the motion of a grain mass will be retarded or damped so that this will assist, ultimately, in terminating the flow. The greatest Coulombic viscosity will be created in the most highly charged systems, which will also be the most energetic. Thus, grain flows have some tendency to be self-limiting by internal energy partitioning; gravitational potential is converted to Coulombic potential, which manifests itself as a drag force between the grains. III. Volcanic eruption plumes and impact ejecta curtains: The violence of these systems leads to powerful electrical charging of particulates. Lightning storms emanating from volcanic plumes are a testimony to the levels of charging. As pyroclastic grains interact forcefully and frequently within eruption plumes, it is reasonable to predict that the internal turbulent motions of the plume will be significantly damped by the Coulombic viscosity exerted by grain charges. Additional information is contained in the original.

Marshall, J. R.↗

Substituent and Heteroatom Effects on π–π Interactions: Evidence That Parallel-Displaced π-Stacking is Not Driven by Quadrupolar Electrostatics

Stacking interactions are a recurring motif in supramolecular chemistry and biochemistry, where a persistent theme is a preference for parallel-displaced aromatic rings rather than face-to-face π-stacking. This is usually explained in terms of quadrupole–quadrupole interactions between the arene moieties but that interpretation is inconsistent with accurate calculations, which reveal that the quadrupolar picture is qualitatively wrong. At typical π-stacking distances, quadrupolar electrostatics may differ in sign from an exact calculation based on charge densities of the interacting arenes. We apply symmetry-adapted perturbation theory to dimers composed of substituted benzene and various aromatic heterocycles, which display a wide range of electrostatic interactions, and we investigate the interplay of Pauli repulsion, dispersion, and electrostatics as it pertains to parallel-displaced π-stacking. Profiles of energy components along cofacial slip-stacking coordinates support a prominent role for the “van der Waals model” (dispersion in competition with Pauli repulsion), even for polar monomers where electrostatic interactions are significant. While electrostatic interactions are necessary to explain the optimal face-to-face π-stacking distance and to account for the relative orientation of one polar arene with respect to another, we find no evidence to support continued invocation of quadrupolar electrostatics as a basis for π-stacking. Our results suggest that a driving force for offset-stacking exists even in the absence of electrostatic interactions. Consequently, tuning electrostatics via functionalization does not guarantee that slip-stacking can be avoided. This has implications for rational design of soft materials and other supramolecular architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The very slow nova He 3-558

The emission-line object He 3-558 has been found to be a very slow nova which rose to maximum brightness in 1948 and is probably still fading. Its spectrum and light curve are described, and a distance of 3.5 kpc is derived. No sign of an earlier outburst was found in the Harvard collection of photographs.

Henize, K. G.↗

The non-existence of the Oort cometary shell

The procedures adopted as theory for a shell of comets are shown to be invalid. Any plot of numbers of long-period comets against 1/a will automatically exhibit a peak at small values of this parameter, and cannot be inverted to demonstrate a high volume-density of aphelia in space. The positions of actual aphelion-points show no sign of any concentration at any range. Further, the aphelion-distance undergoes large almost random changes owing to planetary perturbations at each return, and present values can yield no indication of original positions. A recent attempt to save the shell-theory by making use solely of comets with large perihelion-distance is shown to rest on exactly the same errors as have done all the earlier presentations. The plain conclusions emerge that the shell-theory is devoid of any support by facts, and that the alleged shell of comets is non-existent.

Lyttleton, R. A.↗

Errors in Hadamard spectroscopy or imaging caused by imperfect masks

An analysis is given of the errors in Hadamard spectroscopy that are caused by transparent slits in the mask being systematically wider or else narrower than they should be. It is shown that if the input spectrum consists of a single line, the distorted spectrum that is actually calculated consists of this line, plus four small blips. When the transparent slits are too wide, these blips are of equal height and the same sign, one pair surrounding the line, and another pair displaced a certain distance from it. When the slits are too narrow, the displaced blips have the same amplitude but are negative.

Tai, M. H.↗

Mean sea surface and geoid gradient comparisons with TOPEX altimeter data

Cycles 4 to 54 of TOPEX data have been analyzed through comparisons with the mean sea surface given on the disturbed geophysical data record (GDR). Two inverted barometer correction procedures were considered for the data reduction. One used a constant atmospheric pressure for all data while the one adopted for use, for most computations, introduced a cycle average pressure. The maximum difference between the two estimates was 3.0 cm with a clear annual signal. With the modified correction the TOPEX sea surface was compared to The Ohio State University (OSU) mean sea surface, given on the GDR, to estimate three translations ( delta x = -2.3 cm; delta y = 25.0 cm; delta z = -0.3 cm) and a bias (43.3 cm) between the two surfaces. The only significant translation is delta y which indicates the reference frame of the TOPEX system differs from that used in the OSU mean sea surface system. The bias between the TOPEX mean sea surface and the OSU mean sea surface was used to estimate an equatorial radius of 6,378,136.55 m based on an 18-cm biased estimate of the TOPEX altimeter. Examination of the average difference, by cycle, between the TOPEX sea surface and the OSU mean sea surface suggested a bias change of 3.1 +/- 2.2 mm/yr with a positive sign indicating the average ocean surface is rising or the altimeter measured distance is decreasing. Models were implemented that solved directly for a bias, bias rate annual/semiannual, and tide correction terms. The computations indicated that a simultaneous solution for this bias, bias rate, and annual/semiannual terms gave the most accurate results. Nonsimultaneous solutions led to slightly different bias rate values. The root mean square difference between the TOPEX sea surface and OSU sea surface, after translation and bias correction, was +/- 17 cm for a typical cycle. Some locations were indentified where the difference could reach 2.3 cm and were repeated over several cycles indicating errors in the mean sea surface. Most of the large differences occur in regions lacking altimeter data prior to the TOPEX/POSEIDON mission and/or areas of significant bathymetric signature. Geoid gradients are needed for the reduction of individual track data to a reference track. The accuracy of the determination of such gradients was determined through the comparison of predicted along-track gradients to the observed gradients. Among four mean sea surfaces tested the best agreement was found with the OSU mean sea surface placed on the TOPEX geophysical data record.

Rapp, Richard H.↗

Errors in Hadamard spectroscopy or imaging caused by imperfect masks

An analysis is given of the errors in Hadamard spectroscopy that are caused by transparent slits in the mask being systematically wider or narrower than they should be. It is shown that if the input spectrum consists of a single line, the distorted spectrum that is actually calculated consists of this line plus four small blips. When the transparent slits are too wide, these blips are of equal height and the same sign with one pair surrounding the line and another pair 5isplaced a certain distance from it. When the slits are too narrow, the displaced blips have the same amplitude, but are negative. The response to an arbitrary input spectrum is then determined from this. The same method of analysis may also be used to handle other types of errors.

Tai, M. H.↗

Scale Dependence of Magnetic Helicity in the Solar Wind

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k approximately equal to 2AU(sup -1) (or frequency nu approximately equal to 2 microHz) at distances below 2.8AU and at k approximately equal to 30AU(sup -1) (or nu approximately equal to 25 microHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10(sup 45) Mx(sup 2) cycle(sup -1) at large scales and to a three times lower value at smaller scales.

Brandenburg, Axel↗

Autonomous Formation Flight: Project Overview

Objectives: a) Map the vortex effects; b) Formation Auto-Pilot Requirements. Two NASA F/A-18 aircraft in formation: a) NASA 845 Systems Research Aircraft; b) NASA 847 Support Aircraft. Flight Conditions: M = 0.56, 25000 feet (Subsonic condition); b) M = 0.86, 36000 feet (Transonic condition). Nose-To-Tail (N2T) Distances: 20, 55, 110 and 190 feet. Lessons learned: a) Controllable flight in vortex is possible with pilot feedback (displays); b) Position hold at best C(sub D), is attainable; c) Best drag location is close to max rolling moment; e) Drag reductions demonstrated up to 22% (WFE up to 20%); f) Induced drag results compare favorably with simple prediction model; g) "Sweet Spot" (lateral & vertical area > 25%) is larger than predicted; h) Larger wing overlaps result in sign reversals in roll, yaw; i) As predicted, favorable effects degrade gradually with increased nose-to-tail distances after peaking at 3 span lengths aft; and j) Demonstrated - over 100 N mi (>15%) range improvement and 650 lbs (14%) fuel savings on actual simulated F/A-18 cruise mission.

Cole, Jennifer↗

The effect of opposing unsteady vorticity on turbulent wall flow

A cylinder and a thin plate were placed close together in a flow adjacent to a wall to study the effects on the turbulent boundary layer. Different spacings of the cylinder and plate within the shear flow were investigated to assess the possibility of lowering the production of fluctuating vorticity in the boundary layer by generating a fluctuating vorticity of opposite sign. Streakline photographs visualized changes in the flow induced by alterations in the cylinder/plate separation distance, the flow velocity and the angle of attack of the thin plate. Drag data were also acquired with varying thicknesses of the thin plate and diameters of the cylinder. Downstream skin friction reductions were obtained with the production of unsteady control vortices with a low turbulence boundary layer. Up to 4 percent drag reduction was also observed when the cylinder was sufficiently far from the wall.

Goodman, W. L.↗

East–West Proton Flux Anisotropy Observed with the PAMELA Mission

We present a study of the east–west anisotropy of trapped-proton fluxes in low-Earth orbit based on the measurements of the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics(PAMELA)experiment. The differential intensities of eastward- and westward-traveling protons detected in the South Atlantic Anomaly region were estimated as a function of equatorial pitch angle and drift shell, for six energy bins between80 MeV and 2 GeV. We found that, as a consequence of the strong atmospheric gradient coupled with the large gyro radius in this energy range, the intensities of eastward fluxes exceed those of westward fluxes by a factor of∼10–20. However, the reported directional asymmetry also depends on the sign of the local flux gradient, resulting in more intense westward fluxes beyond the radial distances where the inner belt peaks. PAMELA observations can be used to improve the description of the near-Earth radiation environment at lowest altitudes and highest trapping energies, where current theoretical and empirical models are affected by the largest uncertainties. Unified Astronomy Thesaurus concepts: Cosmic rays(329);Van Allen radiation belts(1758)1. Introduction Low-altitude inner-belt protons are strongly influenced by the density distribution of Earth’s atmosphere, mostly through interactions with its neutral constituents, which induce significant flux anisotropies. In prim is, the atmospheric loss cone results in a steep pitch-angle distribution, which becomes narrower for lower drift shells. A further, azimuthal anisotropy originates from finite gyro radius effects at proton energies in excess of a few tens of MeV(Haerendel1962; Lenchek & Singer1962). In fact, for a given spacecraft position, protons with the same pitch angle but different gyro phase—the azimuth angle associated with the gyration motion—have their guiding centers on different drift shells. In particular, protons from the west and from the east gyrate around magnetic field lines located at higher and lower altitudes, respectively. The guiding-center separationΔhincreaseswith increasing energy, so that protons moving eastward will encounter progressively lower drift-averaged densities, thus experiencing less atmospheric absorption; the opposite situation will occur for protons traveling westward, resulting in an east–west asymmetry of flux intensities whenΔhbecomes comparable to or larger than the flux scale height(Garmire1963; Heckman &Nakano1963).The trapped-flux anisotropy is a relevant aspect of the modeling of the low Earth orbit(LEO)radiation environment, given the significant engineering implications, especially for The Astrophysical Journal,919:114(6pp), 2021 October 1https://doi.org/10.3847/1538-4357/ac1677© 2021. The American Astronomical Society.

A. Bruno↗

Vortex interactions and decay in aircraft wakes

The dynamic interactions of aircraft wake vortices are investigated using both inviscid and viscous models. The phenomenon of vortex merging resulting in the rapid aging of a vortex wake is examined in detail. It is shown that the redistribution of vorticity from convection and diffusion during merging is a mechanism effective in reducing the hazard of a wake. Inviscid computations show that the merging phenomenon may be sensitive to small changes in spanwise load distribution and that the fuselage vortex shed from the wing-fuselage junction can play a significant role in promoting merging of wing-tip and flap vortices. Vortex-wake-merging computations using a second-order closure model of turbulent transport indicate that a low-hazard wake occurs when the generating aircraft flap and wing-tip vortices of the same strength and sign. This optimum is achieved when the flap vortex is located outboard approximately 40% of the distance to the tip vortex.

Bilanin, A. J.↗

Clipping polygon faces through a polyhedron of vision

A flight simulator combines flight data and polygon face terrain data to provide a CRT display at each window of the simulated aircraft. The data base specifies the relative position of each vertex of each polygon face therein. Only those terrain faces currently appearing within the pyramid of vision defined by the pilots eye and the edges of the pilots window need be displayed at any given time. As the orientation of the pyramid of vision changes in response to flight data, the displayed faces are correspondingly displaced, eventually moving out of the pyramid of vision. Faces which are currently not visible (outside the pyramid of vision) are clipped from the data flow. In addition, faces which are only partially outside of pyramid of vision are reconstructed to eliminate the outside portion. Window coordinates are generated defining the distance between each vertex and each of the boundary planes forming the pyramid of vision. The sign bit of each window coordinate indicates whether the vertex is on the pyramid of vision side of the associated boundary panel (positive), or on the other side thereof (negative). The set of sign bits accompanying each vertex constitute the outcode of that vertex. The outcodes (O.C.) are systematically processed and examined to determine which faces are completely inside the pyramid of vision (Case A--all signs positive), which faces are completely outside (Case C--All signs negative) and which faces must be reconstructed (Case B--both positive and negative signs).

Rohner, Michel A.↗