Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “higher form”

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 361 records · Page 20

Switched steerable multiple beam antenna system

A steerable multibeam five element cross-feed cluster antenna system is described. The feed power is divided into five branches. Each branch includes a switching network comprised of a plurality of time delay elements each individually controlled by a respective electromagnetic latching switch. Frequency independent individual two-dimensional beam steering at intermediate (IF) scanning frequencies is thereby provided wherein discrete incremental time delays are introduced by the switching networks into each branch and the signals recombined thereafter to form each beam. The electromagnetic latched switching reduces power consumption and permits higher power switching and reciprocal coincident tranmsit and receive operation. Frequency independence due to incremental time delay switching permits coincident reciprocal operation and steering for transmit-receive signal paths carrying different transmit-receive frequencies. Diagonal quarter wave plates in the waveguides alter polarization from the circular to orthogonal linear to provide transmitter-receiver isolation.

Richard S Iwasaki↗

Fracture-tough, corrosion-resistant bearing steels

The fundamental principles allowing design of stainless bearing steels with enhanced toughness and stress corrosion resistance has involved both investigation of basic phenomena in model alloys and evaluation of a prototype bearing steel based on a conceptual design exercise. Progress in model studies has included a scanning Auger microprobe (SAM) study of the kinetics of interfacial segregation of embrittling impurities which compete with the kinetics of alloy carbide precipitation in secondary hardening steels. These results can define minimum allowable carbide precipitation rates and/or maximum allowable free impurity contents in these ultrahigh strength steels. Characterization of the prototype bearing steel designed to combine precipitated austenite transformation toughening with secondary hardening shows good agreement between predicted and observed solution treatment response including the nature of the high temperature carbides. An approximate equilibrium constraint applied in the preliminary design calculations to maintain a high martensitic temperature proved inadequate, and the solution treated alloy remained fully austenitic down to liquid nitrogen temperature rather than transforming above 200 C. The alloy can be martensitically transformed by cryogenic deformation, and material so processed will be studied further to test predicted carbide and austenite precipitation behavior. A mechanistically-based martensitic kinetic model was developed and parameters are being evaluated from available kinetic data to allow precise control of martensitic temperatures of high alloy steels in future designs. Preliminary calculations incorporating the prototype stability results suggest that the transformation-toughened secondary-hardening martensitic-stainless design concept is still viable, but may require lowering Cr content to 9 wt. pct. and adding 0.5 to 1.0 wt. pct. Al. An alternative design approach based on strain-induced martensitic transformation during cryogenic forming, thus removing the high martensitic constraint, may permit alloy compositions offering higher fracture roughness.

Olson, Gregory B.↗

Molten (Mg0.88Fe0.12)2SiO4 at lower mantle conditions - Melting products and structure of quenched glasses

Infrared spectra of quenched magnesium silicate glasses synthesized by fusing olivine at pressures in excess of 50 GPa and temperatures greater than 2500 K demonstrate that silicon is dominantly present in four-fold coordination with respect to oxygen within these quenched glasses. This low coordination is attributed, by analogy with the structural behavior of glasses compressed at 300 K, to the instability of higher coordinations in glasses of these compositions on decompression. Spectra of glasses formed in a hydrous environment document that water is extensively soluble in melts at these high pressures and temperatures. Also, these results are consistent with the melting of (Mg0.88Fe0.12)2SiO4 compositions to liquids near pyroxene in stoichiometry under these conditions, with iron-rich magnesiowuestite being the liquidus phase.

Williams, Quentin↗

Discreteness noise versus force errors in N-body simulations

A low accuracy in the force calculation per time step of a few percent for each particle pair is sufficient for collisionless N-body simulations. Higher accuracy is made meaningless by the dominant discreteness noise in the form of two-body relaxation, which can be reduced only by increasing the number of particles. Since an N-body simulation is a Monte Carlo procedure in which each particle-particle force is essentially random, i.e., carries an error of about 1000 percent, the only requirement is a systematic averaging-out of these intrinsic errors. We illustrate these assertions with two specific examples in which individual pairwise forces are deliberately allowed to carry significant errors: tree-codes on supercomputers and algorithms on special-purpose machines with low-precision hardware.

Hernquist, Lars↗

Are there two disk star cluster systems in the LMC?

The surface distribution of star clusters in the LMC has revealed the existence of two elliptical systems superimposed. A small surface density outer ellipse outlining the large system and a higher density elliptical inner system. The major axis of the two subsystems form an angle of 60 deg. The central subsystem contains all youngest populous globular star clusters from both the stellar and dynamical point of view. The large outer system seems to contain all old globular clusters, showing a dumpy distribution and outlining two arms at the northeast and southwest. The small LMC clusters were found to have masses less than 1000 solar masses and are distributed all over the large elliptical system.

Kontizas, Mary↗

Poly(Arylene Ether)s Containing Pendent Ethynyl Groups

Poly(arylene ether)s containing pendent ethynyl groups synthesized. Offer advantage over linear poly(arylene ether)s in upon heating, pendent ethynyl groups react to form cross-linked molecular structures exhibiting greater resistance to solvents and higher glass-transition temperatures and tensile moduli. Polymers useful as adhesives, moldings, films, and matrices of composite materials.

Hergenrother, Paul M.↗

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗

Mechanisms for the Formations of the Thymine Under Astrophysical Conditions and Implications for the Origin of Life

Nucleobases are the carriers of the genetic information in ribonucleic acid and deoxyribonucleic acid (DNA) for all life on Earth. Their presence in meteorites clearly indicates that compounds of biological importance can form via non-biological processes in extraterrestrial environments. Recent experimental studies have shown that the pyrimidine-based nucleobases uracil and cytosine can be easily formed from the ultraviolet irradiation of pyrimidine in H2O-rich ice mixtures that simulate astrophysical processes. In contrast, thymine, which is found only in DNA, is more difficult to form under the same experimental conditions, as its formation usually requires a higher photon dose. Earlier quantum chemical studies confirmed that the reaction pathways were favorable provided that several H2O molecules surrounded the reactants. However, the present quantum chemical study shows that the formation of thymine is limited because of the inefficiency of the methylation of pyrimidine and its oxidized derivatives in an H2O ice, as supported by the laboratory studies. Our results constrain the formation of thymine in astrophysical environments and thus the inventory of organic molecules delivered to the early Earth and have implications for the role of thymine and DNA in the origin of life.

Origin↗

Optical beam forming techniques for phased array antennas

Conventional phased array antennas using waveguide or coax for signal distribution are impractical for large scale implementation on satellites or spacecraft because they exhibit prohibitively large system size, heavy weight, high attenuation loss, limited bandwidth, sensitivity to electromagnetic interference (EMI) temperature drifts and phase instability. However, optical beam forming systems are smaller, lighter, and more flexible. Three optical beam forming techniques are identified as applicable to large spaceborne phased array antennas. They are (1) the optical fiber replacement of conventional RF phased array distribution and control components, (2) spatial beam forming, and (3) optical beam splitting with integrated quasi-optical components. The optical fiber replacement and the spatial beam forming approaches were pursued by many organizations. Two new optical beam forming architectures are presented. Both architectures involve monolithic integration of the antenna radiating elements with quasi-optical grid detector arrays. The advantages of the grid detector array in the optical process are the higher power handling capability and the dynamic range. One architecture involves a modified version of the original spatial beam forming approach. The basic difference is the spatial light modulator (SLM) device for controlling the aperture field distribution. The original liquid crystal light valve SLM is replaced by an optical shuffling SLM, which was demonstrated for the 'smart pixel' technology. The advantages are the capability of generating the agile beams of a phased array antenna and to provide simultaneous transmit and receive functions. The second architecture considered is the optical beam splitting approach. This architecture involves an alternative amplitude control for each antenna element with an optical beam power divider comprised of mirrors and beam splitters. It also implements the quasi-optical grid phase shifter for phase control and grid amplifier for RF power. The advantages are no SLM is required for this approach, and the complete antenna system is capable of full monolithic integration.

Wu, Te-Kao↗

Early Energetic Particle Irradiation of the HED Parent Body Regolith

Previous studies have shown that many individual grains within the dark phase of the Kapoeta howardite were irradiated with energetic particles while residing on the surface of the early HED regolith. Particle tracks in these grains vary in density by more than an order of magnitude and undoubtedly were formed by energetic heavy (Fe) ions associated with early solar flares. Early Irradiation of HED Regolith: Concentrations of excess Ne alone are not sufficient to decide between competing galactic and solar irradiation models. However, from recent studies of depth samples of oriented lunar rocks, we have shown that the cosmogenic 21-Ne/22-Ne ratio produced in feldspar differs substantially between Galactic Cosmic Radiation (GCR) and solar protons, and that this difference is exactly that predicted from cross-section data. Using Ne literature data and new isotopic data we obtained on acid-etched, separated feldspar from both the light and dark phases of Kapoeta, we derive 21-Ne/22-Ne = 0.80 for the recent GCR irradiation and 21-Ne/22-Ne = 0.68 for the early regolith irradiation. This derived ratio indicates that the early Ne production in the regolith occurred by both galactic and solar protons. If we adopt a likely one-component regolith model in which all grains were exposed to galactic protons but individual grains had variable exposure to solar protons, we estimate that this early GCR irradiation lasted for about 3-6 m.y. More complex two-component regolith models involving separate solar and galactic irradiation would permit this GCR age to be longer. Higher-energy solar protons would permit the GCR to be longer. Higher-energy solar protons would permit the GCR age to be shorter. Further, cosmogenic 126(Xe) in Kapoeta dark is no more than a factor of about 2 higher than that observed in Kapoeta light. Because 126(Xe) can only be formed by galactic protons and not solar protons, these data support a short GCR irradiation for the HED regolith. This would also be the maximum time peRiod for the solar irradiation. Various asteroidal regolith models, based on Monte Carlo modeling of impact rates as a function of size and on irradiation features of meteorites, predict surface exposure times of about 0.1 to 10 m.y., and depend on such factors as gravity, rock mechanical properties, and micrometeoroid flux. Because the depth at which solar Fe tracks are produced (is much less than 1 micrometer) is much less than the depth at which Solar Cosmic Rays (SCR) Ne is produced (about 1 cm), for a reasonably well-stirred HED regolith the "surface exposure time" for SCR 21-Ne production should be significantly longer than that for solar tracks and some other surface irradiation features. Enhanced Solar Proton Irradiation: For bulk samples of Kapoeta dark feldspar and a one-component regolith model, the derived ratio of 21-Ne/22-Ne = 0.68 implies that the early production ratio of SCR 21-Ne to GCR 21-Ne was about 0.5-1.5. This ratio is independent of any assumptions about the fraction of dark grains that are irradiated or of the variability in the degree of solar irradiation among grains. The 21-Ne SCR/GCR ratio indirectly derived from bulk Kapoeta pyroxene is somewhat larger, as is the ratio derived for simple two-component regolith models. Individual feldspar grains that were extensively solar irradiated would require even larger 21-Ne SCR/GCR production ratios. In contrast, the theoretical SCR/GCR production ratio for lunar feldspar with 0 g/CM2 shield ing is is less than or equal to 2, and the lowest ratio observed in near-surface samples of lunar anorthosites is less than or equal to 1. Considering the greater solar distance of Vesta (compared to the Moon), the likelihood that SCR 21-Ne was acquired under some shielding where production rates are lower, and the likelihood that the exposure time to galactic protons exceeded the exposure time to solar protons because of their very different penetration depths, the 21-Ne SCR/GCR production ratio on the HED parent body was probably < 0.1. The relatively large difference between the derived 21-Ne SCR/GCR ratio in Kapoeta dark feldspar and the estimated production ratio strongly indicates that the early solar irradiation involved a flux -20-50x the recent solar flux. This enhanced proton flux was probably associated with an overall greater solar activity in the first approximately 10(exp 7) to 10(exp 8) years of solar history.

Bogard, D. D.↗

Proton affinity determinations and proton-bound dimer structure indications in C2 to C15, (alpha),(omega)-alkyldiamines

The 'kinetic method' was used to determine the proton affinity (PA) of a,coalkyldiamines from collision induced dissociation (CID) studies of protonated heterodimers. These PA values were consistently lower than those reported in the proton affinity scale. The apparent discrepancy was rationalized in terms of differences in the conformation of the protonated diamine monomers. The minimum energy species, formed by equilibrium proton transfer processes, have a cyclic conformation and the ion charge is shared by both amino-groups which are bridged by the proton. On the other hand, the species formed through dissociation of protonated dimers have a linear structure and the charge is localized on one of the amino-groups. Thus, the difference in the PA values obtained by both methods is a measure of the additional stability acquired by the protonated diamines through cyclization and charge delocalization. The major collision dissociation pathway of the protonated diamine monomers involved elimination of an ammonia moiety. Other reactions observed included loss of the second amino-group and several other bond cleavages. CID of the protonated dimers involved primarily formation of a protonated monomer through cleavage of the weaker hydrogen bond and subsequently loss of ammonia at higher collision energies. As observed from the CID studies, doubly charged ions were also formed from the diamines under conditions of the electrospray ionization.

Karpas, Z.↗

Silicon crystals: Process for manufacturing wafer-like silicon crystals with a columnar structure

Wafer-like crystals suitable for making solar cells are formed by pouring molten Si containing suitable dopants into a mold of the desired shape and allowing it to solidify in a temperature gradient, whereby the large surface of the melt in contact with the mold is kept at less than 200 D and the free surface is kept at a temperature of 200-1000 D higher, but below the melting point of Si. The mold can also be made in the form of a slit, whereby the 2 sides of the mold are kept at different temperatures. A mold was milled in the surface of a cylindrical graphite block 200 mm in diameter. The granite block was induction heated and the bottom of the mold was cooled by means of a water-cooled Cu plate, so that the surface of the mold in contact with one of the largest surfaces of the melt was held at approximately 800 D. The free surface of the melt was subjected to thermal radiation from a graphite plate located 2 mm from the surface and heated to 1500 D. The Si crystal formed after slow cooling to room temperature had a columnar structure and was cut with a diamond saw into wafers approximately 500 mm thick. Solar cells prepared from these wafers had efficiencies of 10 to 11%.

Authier, B.↗

Single-Chip T/R Module for 1.2 GHz

A single-chip CMOS-based (complementary-metal-oxide-semiconductorbased) transmit/receive (T/R) module is being developed for L-band radar systems. Previous T/R module implementations required multiple chips employing different technologies (GaAs, Si, and others) combined with off-chip transmission lines and discrete components including circulators. The new design eliminates the bulky circulator, significantly reducing the size and mass of the T/R module. Compared to multi-chip designs, the single-chip CMOS can be implemented with lower cost. These innovations enable cost-effective realization of advanced phased array and synthetic aperture radar systems that require integration of thousands of T/R modules. The circulator is a ferromagnetic device that directs the flow of the RF (radio frequency) power during transmission and reception. During transmission, the circulator delivers the transmitted power from the amplifier to the antenna, while preventing it from damaging the sensitive receiver circuitry. During reception, the circulator directs the energy from the antenna to the low-noise amplifier (LNA) while isolating the output of the power amplifier (PA). In principle, a circulator could be replaced by series transistors acting as electronic switches. However, in practice, the integration of conventional series transistors into a T/R chip introduces significant losses and noise. The prototype single-chip T/R module contains integrated transistor switches, but not connected in series; instead, they are connected in a shunt configuration with resonant circuits (see figure). The shunt/resonant circuit topology not only reduces the losses associated with conventional semiconductor switches but also provides beneficial transformation of impedances for the PA and the LNA. It provides full singlepole/ double-throw switching for the antenna, isolating the LNA from the transmitted signal and isolating the PA from the received signal. During reception, the voltage on control line RX/TX (raised bar) is high, causing the field-effect transistor (FET) switch S1 to be closed, forming a parallel resonant tank circuit L1||C1. This circuit presents high impedance to the left of the antenna, so that the received signal is coupled to the LNA. At the same time, FET switches S2 and S3 are open, so that C2 is removed from the circuit (except for a small parasitic capacitance). The combination of L2 and C3 forms a matching network that transforms the antenna impedance of 50 ohms to a higher value from the perspective of the LNA input terminal. This transformation of impedance improves LNA noise figure by increasing the received voltage delivered to the input transistor. This allows lower transconductance and therefore a smaller transistor, which makes it possible to design the CMOS LNA for low power consumption. During transmission, the voltage on control line RX/TX (raised bar) is low, causing switch S1 to be open. In this configuration, the combination of L1 and C1 transforms the antenna impedance to a lower value from the perspective of the PA. This low impedance is helpful in producing a relatively high output power compatible with the low CMOS operating potential. At the same time, switches S2 and S3 are closed, forming the parallel resonant tank circuit L2||C2. This circuit presents high impedance to the right of the antenna, directing the PA output signal to the antenna and away from the LNA. During this time, S3 presents a short circuit across the LNA input terminals to guarantee that the voltage seen by the LNA is small enough to prevent damage.

Moussessian, Alina↗

Full-Scale Wind-Tunnel Tests of a Series of Metal Propellers on a VE-7 Airplane

An adjustable blade metal propeller was tested at five different angle settings, forming a series varying in pitch. The propeller was mounted on a VE-7 airplane in the twenty-foot propeller research tunnel of the National Advisory Committee for Aeronautics. The efficiencies were found to be from 4 to 7 per cent higher than those of standard wood propellers operating under the same conditions. The results are given in convenient form for use in selecting propellers for aircraft.

Weick, Fred E↗

Ion Beam Heating in the Auroral Zone

Recent satellite observations at high altitudes (greater than 5000 km) in the auroral zone have shown the existence of hybrid or bimodal ion beam distributions that are evidence of both parallel and perpendicular ion acceleration. Acceleration parallel to the magnetic field is most likely due to quasi-static electric fields (double layers) which can create outflowing ion beams; since ions of different mass will have different drift speeds due to this acceleration, a plasma configuration unstable to the ion-ion two-stream acoustic mode develops. When the net drift velocity (U) between the two ion species is greater than the sound speed (C(sub 0)), the ion-ion instability has maximum growth at oblique wave propagation. To study the nonlinear effects of the ion-ion instability in terms of plasma heating, a numerical simulation parametric study has been performed. It was found that the parallel acceleration that forms the ion beams occurs on a time scale faster than ion- ion wave growth at low drifts; thus ion-ion wave growth is expected to occur primarily for higher drift speeds (U greater than C(sub 0)) which results in strong oblique heating of the ions (both hydrogen and oxygen) forming elevated ion conics (sometimes called 'bowl' distributions). Also, strong parallel electron heating in the direction of the ion beams can occur, and electrons near the top of the acceleration region may attain a net upward drift along with the elevated ion conics. Variation of the oxygen density greatly affects the ion heating due to the ion-ion instability; as the oxygen density decreases, oxygen heating increases, in agreement with observations (Collin et al., 1987). Ion-ion electrostatic wave properties and the plasma heating that results over a wide range of auroral zone parameters are included.

Schriver, David↗

The crystallization of tough thermoplastic resins in the presence of carbon fibers

The presence of carbon fibers increased the crystallization rates of both PEEK and PPS thermoplastic polymers. The effect was most pronounced at higher crystallization temperatures. Isothermal crystallization rates were analyzed by applying classical phenomenological nucleation theory. Unusually high values of the so-called Avrami exponent were found for neat PEEK. Isothermal crystallization of PEEK and PPS polymers produced crystalline samples having a wide variety of melting temperatures. The melting as observed by differential scanning calorimetry occurred as dual endotherms which were called primary (higher temperature) and secondary melting peaks. Each primary peak accounted for most of the crystallinity present. The secondary peaks represented the melting of crystallites formed later than those attributable to the primary endotherms. The presence of carbon fibers increased the thermal stability of both PEEK and PPS crystallites as manifested by higher temperatures for the primary melting peaks. This may be attributable to increased crystallite size, greater crystallite perfection, and/or favorable modification of the crystallite interface. Over the range studied, crystallization temperature strongly influenced the positions of the secondary peaks but not the primary peaks.

Theil, Michael H.↗

Time Series Explorer: Bayesian Blocks with Generalized Profiles and in Higher Dimensions

The Time Series Explorer (TSE) is a project aimed at provided new and advanced time series analysis algorithms in two forms: a tool kit and an automated pipeline applying selected tools in machine learning settings. I will present a sketch of TSE with emphasis on time-domain modeling in general and recent improvements of the Bayesian Block (BB) algorithm in particular. This includes generalizing the shape of the elementary blocks from the current constant-rate model to general shapes, such as two-sided exponentials. Related topics will include extension of BB to higher dimensions and a novel way to detect and characterize short time-scale bursts in time-tagged event data. The Fermi Gamma Ray Space Telescope light curve for the Crab Nebula will be used as an example for all of the algorithms discussed. This work is in collaboration with Tom Loredo.

Bayesian Block (BB) algorithm↗