Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear correlation”

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 577 records · Page 32

Optimal Initial Perturbations for Ensemble Prediction of the Madden-Julian Oscillation during Boreal Winter

An initialization strategy, tailored to the prediction of the Madden-Julian oscillation (MJO), is evaluated using the Goddard Earth Observing System Model, version 5 (GEOS-5), coupled general circulation model (CGCM). The approach is based on the empirical singular vectors (ESVs) of a reduced-space statistically determined linear approximation of the full nonlinear CGCM. The initial ESV, extracted using 10 years (1990-99) of boreal winter hindcast data, has zonal wind anomalies over the western Indian Ocean, while the final ESV (at a forecast lead time of 10 days) reflects a propagation of the zonal wind anomalies to the east over the Maritime Continent an evolution that is characteristic of the MJO. A new set of ensemble hindcasts are produced for the boreal winter season from 1990 to 1999 in which the leading ESV provides the initial perturbations. The results are compared with those from a set of control hindcasts generated using random perturbations. It is shown that the ESV-based predictions have a systematically higher bivariate correlation skill in predicting the MJO compared to those using the random perturbations. Furthermore, the improvement in the skill depends on the phase of the MJO. The ESV is particularly effective in increasing the forecast skill during those phases of the MJO in which the control has low skill (with correlations increasing by as much as 0.2 at 20 25-day lead times), as well as during those times in which the MJO is weak.

Ham, Yoo-Geun↗

Parameter identification of material constants in a composite shell structure

One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Martinez, David R.↗

Progressive Failure Analysis Methodology for Laminated Composite Structures

A progressive failure analysis method has been developed for predicting the failure of laminated composite structures under geometrically nonlinear deformations. The progressive failure analysis uses C(exp 1) shell elements based on classical lamination theory to calculate the in-plane stresses. Several failure criteria, including the maximum strain criterion, Hashin's criterion, and Christensen's criterion, are used to predict the failure mechanisms and several options are available to degrade the material properties after failures. The progressive failure analysis method is implemented in the COMET finite element analysis code and can predict the damage and response of laminated composite structures from initial loading to final failure. The different failure criteria and material degradation methods are compared and assessed by performing analyses of several laminated composite structures. Results from the progressive failure method indicate good correlation with the existing test data except in structural applications where interlaminar stresses are important which may cause failure mechanisms such as debonding or delaminations.

Sleight, David W.↗

Nonlinear interactions of electromagnetic waves with a thermally generated cesium plasma.

The interaction of a large‐amplitude disturbing wave with a thermally generated, highly ionized cesium plasma has been investigated using a small‐amplitude sensing wave as a probe. The subsequent interaction between the disturbing and sensing waves is related to the plasma electron temperature. The plasma was contained in a magnetic field of 500–1000 G and had an electron plasma frequency near 10 GHz and an electron gyromagnetic frequency around 2 GHz. The sensing wave was a 10 GHz plane wave of approximately 1.5 mW input power and the disturbing wave was 7.5 GHz plane wave of about 1.2 W peak power. Measurements of the sensing‐wave transmitted power were made using standard microwave techniques. A notch filter in the sensing‐wave receiving circuit reduced the scattered disturbing‐wave interference to an immeasurably low level. Langmuir probes were used to detect disturbing‐wave modulation and to measure plasma temperatures and densities. Measurements of the interaction were correlated with a solution of the Boltzmann transport equation for a perturbed Maxwellian distribution function with binary Coulomb collisions.

Jack G Hewitt, Jr↗

A full potential flow analysis with realistic wake influence for helicopter rotor airload prediction

A 3-D, quasi-steady, full potential flow solver was adapted to include realistic wake influence for the aerodynamic analysis of helicopter rotors. The method is based on a finite difference solution of the full potential equation, using an inner and outer domain procedure for the blade flowfield to accommodate wake effects. The nonlinear flow is computed in the inner domain region using a finite difference solution method. The wake is modeled by a vortex lattice using prescribed geometry techniques to allow for the inclusion of realistic rotor wakes. The key feature of the analysis is that vortices contained within the finite difference mesh (inner domain) were treated with a vortex embedding technique while the influence of the remaining portion of the wake (in the outer domain) is impressed as a boundary condition on the outer surface of the finite difference mesh. The solution procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. The method has been applied to both hover and forward flight conditions. Correlation with subsonic and transonic hover airload data is shown which demonstrates the merits of the approach.

Egolf, T. Alan↗

Traveling-Wave Tube Amplifier Model to Predict High-Order Modulation Intersymbol Interference

Demands for increased data rates in satellite communications necessitate higher order modulation schemes, larger system bandwidth, and minimum distortion of the modulated signal as it is passed through the traveling wave tube amplifier (TWTA). One type of distortion that the TWTA contributes to is intersymbol interference (ISI), and this becomes particularly disruptive with wide-band, complex modulation schemes. It is suspected that in addition to the dispersion of the TWT, frequency dependent reflections due to mismatches within the TWT are a significant contributor to ISI. To experimentally investigate the effect of these mismatches within the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers in addition to the acquisition of the required digital hardware. In an attempt to develop a more accurate model to correlate IS1 with the TWTA and the operational signal, a fully three-dimensional (3D), time-dependent, TWT interaction model has been developed using the electromagnetic particle-in-cell (PIC) code MAFIA (solution of Maxwell's equations by the Finite-Integration-Algorithm). The model includes a user defined slow-wave circuit with a spatially tapered region of loss to implement a sever, and spatially varied geometry (such as helical pitch) to implement a phase velocity taper. The model also includes user defined input/output coupling and an electron beam contained by solenoidal, electrostatic, or periodic permanent magnet (PPM) focusing allowing standard or novel TWTs to be investigated. This model comprehensively takes into account the effects of frequency dependent nonlinear distortions (MAM and AMPM); gain ripple due to frequency dependent reflections at the input/output coupling, severs, and mismatches from dynamic pitch variations; drive induced oscillations; harmonic generation; intermodulation products; and backward waves.

Kory, Carol L.↗

Assessing the Tangent Linear Behaviour of Common Tracer Transport Schemes and Their Use in a Linearised Atmospheric General Circulation Model

The linearity of a selection of common advection schemes is tested and examined with a view to their use in the tangent linear and adjoint versions of an atmospheric general circulation model. The schemes are tested within a simple offline one-dimensional periodic domain as well as using a simplified and complete configuration of the linearised version of NASA's Goddard Earth Observing System version 5 (GEOS-5). All schemes which prevent the development of negative values and preserve the shape of the solution are confirmed to have nonlinear behaviour. The piecewise parabolic method (PPM) with certain flux limiters, including that used by default in GEOS-5, is found to support linear growth near the shocks. This property can cause the rapid development of unrealistically large perturbations within the tangent linear and adjoint models. It is shown that these schemes with flux limiters should not be used within the linearised version of a transport scheme. The results from tests using GEOS-5 show that the current default scheme (a version of PPM) is not suitable for the tangent linear and adjoint model, and that using a linear third-order scheme for the linearised model produces better behaviour. Using the third-order scheme for the linearised model improves the correlations between the linear and non-linear perturbation trajectories for cloud liquid water and cloud liquid ice in GEOS-5.

Tangent Linear Model↗

Soft Soil Impact Testing and Simulation of Aerospace Structures

In June 2007, a 38-ft/s vertical drop test of a 5-ft-diameter, 5-ft-long composite fuselage section that was retrofitted with a novel composite honeycomb Deployable Energy Absorber (DEA) was conducted onto unpacked sand. This test was one of a series of tests to evaluate the multi-terrain capabilities of the DEA and to generate test data for model validation. During the test, the DEA crushed approximately 6-in. and left craters in the sand of depths ranging from 7.5- to 9-in. A finite element model of the fuselage section with DEA was developed for execution in LS-DYNA, a commercial nonlinear explicit transient dynamic code. Pre-test predictions were generated in which the sand was represented initially as a crushable foam material MAT_CRUSHABLE_FOAM (Mat 63). Following the drop test, a series of hemispherical penetrometer tests were conducted to assist in soil characterization. The penetrometer weighed 20-lb and was instrumented with a tri-axial accelerometer. Drop tests were performed at 16-ft/s and crater depths were measured. The penetrometer drop tests were simulated as a means for developing a more representative soil model based on a soil and foam material definition MAT_SOIL_AND FOAM (Mat 5) in LS-DYNA. The model of the fuselage with DEA was reexecuted using the updated soil model and test-analysis correlations are presented.

Fasanella, Edwin L.↗

Experimental Evaluation of Fatigue Damage Progression in Postbuckled Single Stringer Composite Specimens

The durability and damage tolerance of postbuckled composite structures are not yet completely understood, and remain difficult to predict due to the nonlinearity of the geometric response and its interaction with local damage modes. A research effort was conducted to investigate experimentally the quasi-static and fatigue damage progression in a single-stringer compression (SSC) specimen. Three specimens were manufactured with a hat-stiffener, and an initial defect was introduced with a Teflon film embedded between one flange of the stringer and the skin. One of the specimens was tested under quasi-static compressive loading, while the remaining two specimens were tested by cycling in postbuckling. The tests were performed at the NASA Langley Research Center under controlled conditions and with instrumentation that allows a precise evaluation of the postbuckling response and of the damage modes. Three-dimensional digital image correlation VIC-3D systems were used to provide full field displacements and strains on the skin and the stringer. Passive thermal monitoring was conducted during the fatigue tests using an infrared camera that showed the location of the delamination front while the specimen was being cycled. The live information from the thermography was used to stop the fatigue tests at critical stages of the damage evolution to allow detailed ultrasonic scans.

Bisagni, Chiara↗

Electron Scattering by High-Frequency Whistler Waves at Earth's Bow Shock

Electrons are accelerated to non-thermal energies at shocks in space and astrophysical environments. While different mechanisms of electron acceleration have been proposed, it remains unclear how non-thermal electrons are produced out of the thermal plasma pool. Here, we report in situ evidence of pitch-angle scattering of non-thermal electrons by whistler waves at Earths bow shock. On 2015 November 4, the Magnetospheric Multiscale (MMS) mission crossed the bow shock with an Alfvn Mach number is approximately 11 and a shock angle of approximately 84deg. In the ramp and overshoot regions, MMS revealed bursty enhancements of non-thermal (0.52 keV) electron flux, correlated with high-frequency (0.2 - 0.4 Omega(sub ce), where Omega(sub ce) is the cyclotron frequency) parallel-propagating whistler waves. The electron velocity distribution (measured at 30 ms cadence) showed an enhanced gradient of phase-space density at and around the region where the electron velocity component parallel to the magnetic field matched the resonant energy inferred from the wave frequency range. The flux of 0.5 keV electrons (measured at 1ms cadence) showed fluctuations with the same frequency. These features indicate that non-thermal electrons were pitch-angle scattered by cyclotron resonance with the high-frequency whistler waves. However, the precise role of the pitch-angle scattering by the higher-frequency whistler waves and possible nonlinear effects in the electron acceleration process remains unclear.

Oka, M.↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar Rao↗

Large scale structure in universes dominated by cold dark matter

The theory of Gaussian random density field peaks is applied to a numerical study of the large-scale structure developing from adiabatic fluctuations in models of biased galaxy formation in universes with Omega = 1, h = 0.5 dominated by cold dark matter (CDM). The angular anisotropy of the cross-correlation function demonstrates that the far-field regions of cluster-scale peaks are asymmetric, as recent observations indicate. These regions will generate pancakes or filaments upon collapse. One-dimensional singularities in the large-scale bulk flow should arise in these CDM models, appearing as pancakes in position space. They are too rare to explain the CfA bubble walls, but pancakes that are just turning around now are sufficiently abundant and would appear to be thin walls normal to the line of sight in redshift space. Large scale streaming velocities are significantly smaller than recent observations indicate. To explain the reported 700 km/s coherent motions, mass must be significantly more clustered than galaxies with a biasing factor of less than 0.4 and a nonlinear redshift at cluster scales greater than one for both massive neutrino and cold models.

Bond, J. Richard↗

Nature and the nonlinear evolution of electrostatic waves associated with the AMPTE solar wind releases

The nonlinear evolution of the electrostatic waves observed during the AMPTE (Active Magnetosphere Particle Tracer Explorers) solar wind releases is investigated. Previous linear studies indicate two distinct sets of instabilities may be responsible for the generation. One set consists of ion-acoustic type instabilities which are insensitive to the presence of a background magnetic field, while the other group corresponds to the modified two stream instabilities and requires the solar wind flow to be across the ambient magnetic field. To establish which set of instabilities is more viable a detailed linear Vlasov theory has been conducted by numerically solving the full electromagnetic dispersion relation. In addition both the plasma wave as well as the magnetic field measurements by the IRM (Ion Release Module) spacecraft were used to correlate the frequency and the power of the observed waves with the magnitude and the direction of the solar wind magnetic field. The results of these analyses indicate that the ion-acoustic type instabilities have growth rates that are an order of magnitude or more larger than those of the modified two stream instabilities. Results show that both the solar wind protons and the released ions may be heated and accelerated in the directions oblique to the solar wind flow velocity.

Omidi, N.↗

2007 Research and Engineering Annual Report

Selected research and technology activities at NASA Dryden Flight Research Center are summarized. These following activities exemplify the Center's varied and productive research efforts: Developing a Requirements Development Guide for an Automatic Ground Collision Avoidance System; Digital Terrain Data Compression and Rendering for Automatic Ground Collision Avoidance Systems; Nonlinear Flutter/Limit Cycle Oscillations Prediction Tool; Nonlinear System Identification Using Orthonormal Bases: Application to Aeroelastic/Aeroservoelastic Systems; Critical Aerodynamic Flow Feature Indicators: Towards Application with the Aerostructures Test Wing; Multidisciplinary Design, Analysis, and Optimization Tool Development Using a Genetic Algorithm; Structural Model Tuning Capability in an Object-Oriented Multidisciplinary Design, Analysis, and Optimization Tool; Extension of Ko Straight-Beam Displacement Theory to the Deformed Shape Predictions of Curved Structures; F-15B with Phoenix Missile and Pylon Assembly--Drag Force Estimation; Mass Property Testing of Phoenix Missile Hypersonic Testbed Hardware; ARMD Hypersonics Project Materials and Structures: Testing of Scramjet Thermal Protection System Concepts; High-Temperature Modal Survey of the Ruddervator Subcomponent Test Article; ARMD Hypersonics Project Materials and Structures: C/SiC Ruddervator Subcomponent Test and Analysis Task; Ground Vibration Testing and Model Correlation of the Phoenix Missile Hypersonic Testbed; Phoenix Missile Hypersonic Testbed: Performance Design and Analysis; Crew Exploration Vehicle Launch Abort System-Pad Abort-1 (PA-1) Flight Test; Testing the Orion (Crew Exploration Vehicle) Launch Abort System-Ascent Abort-1 (AA-1) Flight Test; SOFIA Flight-Test Flutter Prediction Methodology; SOFIA Closed-Door Aerodynamic Analyses; SOFIA Handling Qualities Evaluation for Closed-Door Operations; C-17 Support of IRAC Engine Model Development; Current Capabilities and Future Upgrade Plans of the C-17 Data Rack; Intelligent Data Mining Capabilities as Applied to Integrated Vehicle Health Management; STARS Flight Demonstration No. 2 IP Data Formatter; Space-Based Telemetry and Range Safety (STARS) Flight Demonstration No. 2 Range User Flight Test Results; Aerodynamic Effects of the Quiet Spike(tm) on an F-15B Aircraft; F-15 Intelligent Flight Controls-Increased Destabilization Failure; F-15 Integrated Resilient Aircraft Control (IRAC) Improved Adaptive Controller; Aeroelastic Analysis of the Ikhana/Fire Pod System; Ikhana: Western States Fire Missions Utilizing the Ames Research Center Fire Sensor; Ikhana: Fiber-Optic Wing Shape Sensors; Ikhana: ARTS III; SOFIA Closed-Door Flutter Envelope Flight Testing; F-15B Quiet Spike(TM) Aeroservoelastic Flight Test Data Analysis; and UAVSAR Platform Precision Autopilot Flight Results.

Stoliker, Patrick↗

Synthesis of Two-Photon Materials and Two-Photon Liquid Crystals

The duration of the grant was interrupted by two major accidents that the PI met with-- an auto accident in Pasadena, CA during her second summer at JPL which took almost eight months for recovery and a second accident during Fall 2000 that left her in crutches for the entire semester. Further, the time released agreed by the University was not given in a timely fashion. The candidate has been given post-grant expire time off. In spite of all these problems, the PI synthesized a number of new two-photon materials and studied the structure-activity correlation to arrive at the best-optimized structure. The PI's design proved to be one of the best in the sense that these materials has a hitherto unreported two-photon absorption cross section. Many materials based on PI's design was later made by the NASA colleague. This is Phase 1. Phase II of this grant is to orate liquid crystalline nature into this potentially useful materials and is currently in progress. Recent observations of nano- and pico-second response time of homeotropically aligned liquid crystals suggest their inherent potentials to act as laser hardening materials, i.e., as protective devices against short laser pulses. The objective of the current project is to exploit this potential by the synthesis of liquid crystals with high optical nonlinearity and optimizing their performance. The PI is trying structural variations to bring in liquid crystalline nature without losing the high two-photon cross section. Both Phase I and Phase II led to many invited presentations and publications in reputed journals like 'Science' and 'Molecular Crystals'. The list of presentations and reprints are enclosed. Another important and satisfying outcome of this grant is the opportunity that this grant offered to the budding undergraduate scientists to get involved in a visible research of international importance. All the students had a chance to learn a lot during research, had the opportunity to present their work at the National level conferences. They continue to retain their interest in their research and went on to accomplish further laurels.

Subramaniam, Girija↗

Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series

The healthy heartbeat is traditionally thought to be regulated according to the classical principle of homeostasis whereby physiologic systems operate to reduce variability and achieve an equilibrium-like state [Physiol. Rev. 9, 399-431 (1929)]. However, recent studies [Phys. Rev. Lett. 70, 1343-1346 (1993); Fractals in Biology and Medicine (Birkhauser-Verlag, Basel, 1994), pp. 55-65] reveal that under normal conditions, beat-to-beat fluctuations in heart rate display the kind of long-range correlations typically exhibited by dynamical systems far from equilibrium [Phys. Rev. Lett. 59, 381-384 (1987)]. In contrast, heart rate time series from patients with severe congestive heart failure show a breakdown of this long-range correlation behavior. We describe a new method--detrended fluctuation analysis (DFA)--for quantifying this correlation property in non-stationary physiological time series. Application of this technique shows evidence for a crossover phenomenon associated with a change in short and long-range scaling exponents. This method may be of use in distinguishing healthy from pathologic data sets based on differences in these scaling properties.

NASA Discipline Number 14-10↗

The evolution of voids in the adhesion approximation

We apply the adhesion approximation to study the formation and evolution of voids in the universe. Our simulations-carried out using 128(exp 3) particles in a cubical box with side 128 Mpc-indicate that the void spectrum evolves with time and that the mean void size in the standard Cosmic Background Explorer Satellite (COBE)-normalized cold dark matter (CDM) model with H(sub 50) = 1 scals approximately as bar D(z) = bar D(sub zero)/(1+2)(exp 1/2), where bar D(sub zero) approximately = 10.5 Mpc. Interestingly, we find a strong correlation between the sizes of voids and the value of the primordial gravitational potential at void centers. This observation could in principle, pave the way toward reconstructing the form of the primordialpotential from a knowledge of the observed void spectrum. Studying the void spectrum at different cosmological epochs, for spectra with a built in k-space cutoff we find that the number of voids in a representative volume evolves with time. The mean number of voids first increases until a maximum value is reached (indicating that the formation of cellular structure is complete), and then begins to decrease as clumps and filaments erge leading to hierarchical clustering and the subsequent elimination of small voids. The cosmological epoch characterizing the completion of cellular structure occurs when the length scale going nonlinear approaches the mean distance between peaks of the gravitaional potential. A central result of this paper is that voids can be populated by substructure such as mini-sheets and filaments, which run through voids. The number of such mini-pancakes that pass through a given void can be measured by the genus characteristic of an individual void which is an indicator of the topology of a given void in intial (Lagrangian) space. Large voids have on an average a larger measure than smaller voids indicating more substructure within larger voids relative to smaller ones. We find that the topology of individual voids is strongly epoch dependent, with void topologies generally simplifying with time. This means that as voids grow older they become progressively more empty and have less structure within them. We evaluate the genus measure both for individual voids as well as for the entire ensemble of voids predicted by CDM model. As a result we find that the topology of voids when taken together with the void spectrum is a very useful statistical indicator of the evolution of the structure of the universe on large scales.

Sahni, Varun↗