Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time scale separation”

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 703 records · Page 39

Structure of the Jovian Magnetodisk Current Sheet: Initial Galileo Observations

The ten-degree tilt of the Jovian magnetic dipole causes the magnetic equator to move back and forth across Jupiter's rotational equator and tile Galileo orbit that lies therein. Beyond about 24 Jovian radii, the equatorial current sheet thins and tile magnetic structure changes from quasi-dipolar into magnetodisk-like with two regions of nearly radial but antiparallel magnetic field separated by a strong current layer. The magnetic field at the center of the current sheet is very weak in this region. Herein we examine tile current sheet at radial distances from 24 55 Jovian radii. We find that the magnetic structure very much resembles tile structure seen at planetary magnetopause and tail current sheet crossings. Tile magnetic field variation is mainly linear with little rotation of the field direction, At times there is almost no small-scale structure present and the normal component of the magnetic field is almost constant through the current sheet. At other times there are strong small-scale structures present in both the southward and northward directions. This small-scale structure appears to grow with radial distance and may provide the seeds for tile explosive reconnection observed at even greater radial distances oil tile nightside. Beyond about 40 Jovian radii, the thin current sheet also appears to be almost constantly in oscillatory motion with periods of about 10 min. The amplitude of these oscillations also appears to grow with radial distance. The source of these fluctuations may be dynamical events in tile more distant magnetodisk.

Russell, C. T.↗

Space Technology 5 Multipoint Observations of Temporal and Spatial Variability of Field-Aligned Currents

Space Technology 5 (ST5) is a constellation mission consisting of three microsatellites. It provides the first multipoint magnetic field measurements in low Earth orbit, which enables us to separate spatial and temporal variations. In this paper, we present a study of the temporal variability of field-aligned currents using the ST5 data. We examine the field-aligned current observations during and after a geomagnetic storm and compare the magnetic field profiles at the three spacecraft. The multipoint data demonstrate that mesoscale current structures, commonly embedded within large-scale current sheets, are very dynamic with highly variable current density and/or polarity in approx.10 min time scales. On the other hand, the data also show that the time scales for the currents to be relatively stable are approx.1 min for mesoscale currents and approx.10 min for large-scale currents. These temporal features are very likely associated with dynamic variations of their charge carriers (mainly electrons) as they respond to the variations of the parallel electric field in auroral acceleration region. The characteristic time scales for the temporal variability of mesoscale field-aligned currents are found to be consistent with those of auroral parallel electric field.

Le, G.↗

A Suzaku X-ray Observation of One Orbit of the Supergiant Fast X-ray Transient IGR J16479-4514

We report on a 250 ks long X-ray observation of the supergiant fast X-ray transient (SFXT) IGR J16479-4514 performed with Suzaku in 2012 February. During this observation, about 80% of the short orbital period (P(sub orb) approximates 3.32 days) was covered as continuously as possible for the first time. The source light curve displays variability of more than two orders of magnitude, starting with a very low emission state (10(exp −13) erg / sq cm/s; 1-10 keV) lasting the first 46 ks, consistent with being due to the X-ray eclipse by the supergiant companion. The transition to the uneclipsed X-ray emission is energy dependent. Outside the eclipse, the source spends most of the time at a level of 6-7X10)(exp−12) erg/sq. cm/s) punctuated by two structured faint flares with a duration of about 10 and 15 ks, respectively, reaching a peak flux of 3-4X10(exp −11) erg/sq. cm./S, separated by about 0.2 in orbital phase. Remarkably, the first faint flare occurs at a similar orbital phase of the bright flares previously observed in the system. This indicates the presence of a phase-locked large scale structure in the supergiant wind, driving a higher accretion rate onto the compact object. The average X-ray spectrum is hard and highly absorbed, with a column density, NH, of 10*exp 23)/sq cm, clearly in excess of the interstellar absorption. There is no evidence for variability of the absorbing column density, except that during the eclipse, where a less absorbed X-ray spectrum is observed. A narrow Fe K-alpha emission line at 6.4 keV is viewed along the whole orbit, with an intensity which correlates with the continuum emission above 7 keV. The scattered component visible during the X-ray eclipse allowed us to directly probe the wind density at the orbital separation, resulting in rho(sub w)=7X10(exp −14) g/cubic cm. Assuming a spherical geometry for the supergiant wind, the derived wind density translates into a ratio M(sub w)/v(sub infinity) = 7X10(exp −17) Solar M/km which, assuming terminal velocities in a large range 500-3000 km/s, implies an accretion luminosity two orders of magnitude higher than that observed. As a consequence, a mechanism should be at work reducing the mass accretion rate. Different possibilities are discussed.

J164794514↗

UAS Air Traffic Controller Acceptability Study: Evaluating Detect and Avoid Technology and Communication Delays in Simulation - 2

This study evaluated the effects of communications delays and winds on air traffic controller ratings of acceptability of horizontal miss distances (HMDs) for encounters between Unmanned Aircraft Systems (UAS) and manned aircraft in a simulation of the Dallas-Ft. Worth (DFW) airspace. Fourteen encounters per hour were staged in the presence of moderate background traffic. Seven recently retired controllers with experience at DFW served as subjects. Guidance provided to the UAS pilots for maintaining a given HMD was provided by information from Detect and Avoid (DAA) self-separation algorithms (Stratway+) displayed on the Multi-Aircraft Control System. This guidance consisted of amber "bands" on the heading scale of the UAS navigation display indicating headings that would result in a loss of well clear between the UAS and nearby traffic. Winds tested were successfully handled by the DAA algorithms and did not affect the controller acceptability ratings of the HMDs. Voice communications delays for the UAS were also tested and included one-way delay times of 0, 400, 1200, and 1800 msec. For longer communications delays, there were changes in strategy and communications flow that were observed and reported by the controllers. The aim of this work is to provide useful information for guiding future rules and regulations applicable to flying UAS in the NAS. Information from this study will also be of value to the Radio Technical Commission for Aeronautics (RTCA) Special Committee 228 - Minimum Performance Standards for UAS.

Comstock, James R., Jr.↗

Measurement and Analysis of Terminal Shock Oscillation and Buffet Forcing Functions on a Launch Vehicle Payload Fairing

The buffet loads on a launch vehicle payload shroud can be impacted by the unsteadiness associated with a terminal shock at high subsonic speeds. At these conditions, flow accelerates to supersonic speeds on the nose of the payload fairing and is terminated by a normal shock on the cylindrical section downstream of the nose cone/cylinder shoulder. The location of the terminal shock and associated separated boundary layer is affected by the freestream Mach number, Reynolds number, and the pitch/yaw of the launch vehicle. Furthermore, even when the freestream conditions and vehicle attitude are constant, this terminal shock oscillates on the surface of the vehicle. The time-varying surface pressure associated with the terminal shock results in unsteady aerodynamic loads that may interact with vehicle structural dynamic modes and the guidance and control of the vehicle. Buffet testing of a 3-percent scale rigid buffet model of a launch vehicle cargo configuration with a tangent-ogive payload shroud was conducted in 2012 and in 2016. Initial buffet forcing functions (BFFs) utilized a coarse pressure sensor distribution on the vehicle surface in which a single longitudinal station with eight sensors observed the terminal shock environment at Mach 0.90. An examination of these circumferential pressures reveal large impulse-like pressure fluctuations and an asymmetry in pressure when the vehicle is at a nonzeroangle of attack that result in high BFFs. Revisions to the shock integration region were made based on computational fluid dynamics and shadowgraph video of shock motion to better represent the BFFs and reduce the high loads resulting from this environment. To more clearly understand this terminal shock environment, a second wind tunnel test was conducted with a dense distribution of 256 sensors at the terminal shock location. These sensor arrays presents a unique opportunity to observe the unsteady terminal shock environment and to characterize the impact of various integration schemes on the BFFs. This paper presents a summary of the development of BFFs for this terminal shock and a detailed analyses of shock region pressure coefficients, coherence, BFFs, shock location time histories, and power spectral density to help guide development of BFFs for other launch vehicle test and analysis programs.

Piatak, David J.↗

Observations of, and sources of the spatial and temporal variability of ozone in the middle atmosphere on climatological time scales (OZMAP) and equatorial dynamics: Seasonal variations of ozone trends

The long term trends (least square linear regression with time) of ozone content at seven European, seven North American, three Japanese and two tropical stations during 21 years (1964 to 1984) are analyzed. In all regions negative trends are observed during the 1970s, but are partly compensated by limited periods of positive trends during the late 1960s and late 1970s. Solely the North American ozone data show negative trends in all 10 year periods. When the long term ozone trends are evaluated for each month of the year separately, a seasonal variation is revealed, which in Europe and North America has largest negative trends in late winter and spring. While in Europe the negative trends in winter/spring are partly compensated by positive trends in summer, in North America the summer values reach only zero, retaining the significant negative trend in annual mean values. In contrast to the antarctic ozone hole, the spring reduction of ozone in Europe and in North America is associated with stratospheric temperatures increasing in the analyzed period and therefore is consistent with the major natural ozone production and loss processes.

Entzian, G.↗

Unveiling the underlying interactions in Ta 2 NiSe 5 from photoinduced lifetime change

We present a generic procedure for quantifying the interplay of electronic and lattice degrees of freedom in photodoped insulators through a comparative analysis of theoretical many-body simulations and time- and angle-resolved photoemission spectroscopy (TR-ARPES) of the transient response of the candidate excitonic insulator Ta 2 NiSe 5 . Our analysis demonstrates that the electron-electron interactions dominate the electron-phonon ones. In particular, a detailed analysis of the TR-ARPES spectrum enables a clear separation of the dominant broadening (electronic lifetime) effects from the much smaller band-gap renormalization. Theoretical calculations show that the observed strong spectral broadening arises from the electronic scattering of the photoexcited particle-hole pairs and cannot be accounted for in a model in which electron-phonon interactions are dominant. The competing interactions were quantified using the scaling analysis in the weak fluence regime. We demonstrate that the magnitude of the weaker subdominant band-gap renormalization sensitively depends on the distance from the semiconductor/semimetal transition in the high-temperature state, which could explain the apparent contradictions between various TR-ARPES experiments. The analysis presented here indicates that electron-electron interactions play a vital role (albeit not the sole one) in stabilizing the insulating state, and establishes the comparison between lifetime and gap evolution as an important probe of correlated insulators.

36 MATERIALS SCIENCE↗

Coupling Waste Feedstocks to Microbial Protein Production in a Circular Food System

Global food production is a major contributor to greenhouse gas emissions, water consumption, and land use. As an alternative to conventional agriculture, the production of waste-derived microbial protein (MP) holds promise for reducing environmental impacts. MP can be mass-produced in volumetrically scalable fermentation processes on short time scales, enabling facile scale-up with lower greenhouse gas emissions, land use, and water impacts than animal and, in some cases, plant protein. MP can also be produced from waste feedstocks, diverting waste from landfills or the natural environment. This Perspective explores the availability and suitability of waste feedstocks for MP production, suggesting that MP generated from waste feedstocks in the United States could fulfill twice the current national protein demand. Here, we also discuss the biotechnological and separations processes required to produce food-grade MP for human consumption from waste. Key challenges include MP consistency, consumer and regulatory acceptance, and the process utilities (electricity, heat, and nutrients) that account for up to 85% of MP costs and most environmental impacts, all of which present opportunities for innovation in the microbiology and process design spaces. Overall, this work highlights the potential of MP to contribute to a more circular, resilient, and sustainable food system.

09 BIOMASS FUELS↗

Influence of Intraseasonal Variability on the Development of Monsoon Depressions

Previous studies showed that the activity of monsoon low-pressure systems (LPS), which produce a large fraction of the South Asian monsoon's total rainfall, is modulated by intraseasonal variability. Using satellite-derived products and atmospheric reanalyzes, this study examines how the boreal summer intraseasonal oscillation (ISO) separately modulates the occurrence of weaker LPS (lows) and stronger LPS (depressions). It is found that the genesis of lows is insensitive to ISO phase, while depressions exhibit a strong preference for genesis during the phase that is convectively active over the northern Indian Ocean. Essentially, development of LPS into depressions depends upon the timing of genesis of the initial disturbance. Evidence is presented supporting the hypothesis that the development of lows into depressions is fostered by large-scale atmospheric conditions governed by the ISO. Furthermore, results also suggest that while lows have no preference for forming over ocean compared to land, depressions mostly form over ocean.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. Herein, a simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The spatially nonuniform convergence of the numerical solution of flows

The spatial distribution of the numerical disturbances that are generated during the numerical solution of a flow is examined. It is shown that the distribution of the disturbances is not uniform. In regions where the structure of a flow is simple, the magnitude of the generated disturbances is small and their decay is fast. However, in complex flow regions, as in separation and vortical areas, large magnitude disturbances appear and their decay may be very slow. The observed nonuniformity of the numerical disturbances makes possible the reduction of the calculation time by application of what may be called the partial-grid calculation technique, in which a major part of the calculation procedure is applied in selective subregions, where the velocity disturbances are large, and not within the whole grid. This technique is expected to prove beneficial in large-scale calculations such as the flow about complete aircraft configurations at high angle of attack. Also, it has been shown that if the Navier-Stokes equations are written in a generalized coordinate system, then in regions in which the grid is fine, such as near solid boundaries, the norms become infinitesimally small, because in these regions the Jacobian has very large values. Thus, the norms, unless they are unscaled by the Jacobians, reflect only the changes that happen at the outer boundaries of the computation domain, where the value of the Jacobian approaches unity, and not in the whole flow field.

Panaras, Argyris G.↗

Inflow measurement made with a laser velocimeter on a helicopter model in forward flight. Volume 3: Rectangular planform blades at an advance ratio of 0.30

An experimental investigation was conducted in the 14- by 22-Foot Subsonic Tunnel at NASA Langley Research Center to measure the inflow into a scale model helicopter rotor in forward flight (micron sub infinity = 0.30). The measurements were made with a two component Laser Velocimeter (LV) one chord above the plane formed by the path of the rotor tips (tip path plane). A conditional sampling technique was employed to determine the azimuthal position of the rotor at the time that each velocity measurement was made so that the azimuthal fluctuations in velocity could be determined. Measurements were made at a total of 180 separate locations in order to clearly define the inflow character. These data are presented without analysis.

Elliott, Joe W.↗

Inflow measurements made with a laser velocimeter on a helicopter model in forward flight. Volume 4: Tapered planform blades at an advance ratio of 0.15

An experimental investigation was conducted in the 14- by 22-Foot Subsonic Tunnel at NASA Langley Research Center to measure the inflow into the scale model helicopter rotor in forward flight (mu sub infinity = 0.15). The measurements were made with a two-component Laser Velocimeter (LV) one chord above the plane formed by the path of the rotor tips (tip path plane). A conditional sampling technique was employed to determine the position of the rotor at the time that each velocity measurement was made so that the azimuthal fluctuations in velocity could be determined. Measurements were made at a total of 146 separate locations in order to clearly define the inflow character. This data is presented herein without analysis. In order to increase the availability of the resulting data, both the mean and azimuthally dependent values are included as part of this report on two 5.25 inch floppy disks in MS-DOS format.

Althoff, Susan L.↗

The spatially non-uniform convergence of the numerical solutions of flows

The spatial distribution of the numerical disturbances that are generated during the numerical solution of a flow is examined. It is shown that the distribution of the disturbances is not uniform. In regions where the structure of a flow is simple, the magnitudes of the generated disturbances is small and their decay is fast. However, in complex flow regions, as in separation and vortical areas, large magnitude disturbances appear and their decay may be very slow. The observed nonuniformity of the numerical disturbances makes possible the reduction of the calculation time by application of what may be called the partial-grid calculation technique, in which a major part of the calculation procedure is applied in selective subregions, where the velocity disturbances are large, and not within the whole grid. This technique is expected to prove beneficial in large-scale calculations such as the flow about complete aircraft configurations at high angle of attack. Also, it has been shown that if the Navier-Stokes equations are written in a generalized coordinate system, then in regions in which the grid is fine, such as near solid boundaries, the norms become infinitesimally small, because in these regions the Jacobian has very large values. Thus, the norms, unless they are unscaled by the Jacobians, reflect only the changes that happen at the outer boundaries of the computation domain, where the value of the Jacobian approaches unity, and not in the whole flow field.

Panaras, Argyris G.↗

Automated Processing of ISIS Topside Ionograms into Electron Density Profiles

Modeling of the topside ionosphere has for the most part relied on just a few years of data from topside sounder satellites. The widely used Bent et al. (1972) model, for example, is based on only 50,000 Alouette 1 profiles. The International Reference Ionosphere (IRI) (Bilitza, 1990, 2001) uses an analytical description of the graphs and tables provided by Bent et al. (1972). The Alouette 1, 2 and ISIS 1, 2 topside sounder satellites of the sixties and seventies were ahead of their times in terms of the sheer volume of data obtained and in terms of the computer and software requirements for data analysis. As a result, only a small percentage of the collected topside ionograms was converted into electron density profiles. Recently, a NASA-funded data restoration project has undertaken and is continuing the process of digitizing the Alouette/ISIS ionograms from the analog 7-track tapes. Our project involves the automated processing of these digital ionograms into electron density profiles. The project accomplished a set of important goals that will have a major impact on understanding and modeling of the topside ionosphere: (1) The TOPside Ionogram Scaling and True height inversion (TOPIST) software was developed for the automated scaling and inversion of topside ionograms. (2) The TOPIST software was applied to the over 300,000 ISIS-2 topside ionograms that had been digitized in the fkamework of a separate AISRP project (PI: R.F. Benson). (3) The new TOPIST-produced database of global electron density profiles for the topside ionosphere were made publicly available through NASA s National Space Science Data Center (NSSDC) ftp archive at . (4) Earlier Alouette 1,2 and ISIS 1, 2 data sets of electron density profiles from manual scaling of selected sets of ionograms were converted fiom a highly-compressed binary format into a user-friendly ASCII format and made publicly available through nssdcftp.gsfc.nasa.gov. The new database for the topside ionosphere established as a result of this project, has stimulated a multitude of new studies directed towards a better description and prediction of the topside ionosphere. Marinov et al. (2004) developed a new model for the upper ion transition height (Oxygen to Hydrogen and Helium) and Bilitza (2004) deduced a correction term for the I N topside electron density model. Kutiev et al. (2005) used this data to develop a new model for the topside ionosphere scale height (TISH) as a function of month, local time, latitude, longitude and solar flux F10.7. Comparisons by Belehaki et al. (2005) show that TISH is in general agreement with scale heights deduced from ground ionosondes but the model predicts post-midnight and afternoon maxima whereas the ionosonde data show a noon maximum. Webb and Benson (2005) reported on their effort to deduce changes in the plasma temperature and ion composition from changes in the topside electron density profile as recorded by topside sounders. Limitations and possible improvements of the IRI topside model were discussed by Coisson et al. (2005) including also the possible use of the NeQuick model, Our project progressed in close collaboration and coordination with the GSFC team involved in the ISIS digitization effort. The digitization project was highly successful producing a large amount of digital topside ionograms. Several no-cost extensions of the TOPIST project were necessary to keep up with the pace and volume of the digitization effort.

Reinisch, bodo W.↗

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Event-by-event direction reconstruction of solar neutrinos in a high light-yield liquid scintillator

The direction of individual B 8 solar neutrinos has been reconstructed using the SNO + liquid scintillator detector. Prompt, directional Cherenkov light was separated from the slower, isotropic scintillation light using time information, and a maximum likelihood method was used to reconstruct the direction of individual scattered electrons. A clear directional signal was observed, correlated with the solar angle. The observation was aided by a period of low primary fluor concentration that resulted in a slower scintillator decay time. This is the first time that event-by-event direction reconstruction in high light-yield liquid scintillator has been demonstrated in a large-scale detector. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Preliminary comparisons of VHF radar maps of F-region irregularities with scintillations in the equatorial region

Multiantenna 50 MHz radar backscatter maps of echo power from night-time F-region equatiorial irregularities obtained at Jicamarca, Peru were compared with simultaneous VHF scintillation observations from Huancayo at 137 and 254 MHz during the period 20 November to 12 December 1975. Saturation of VHF scintillations in excess of 20 dB was observed at both these frequencies during times when radar maps showed large intense plume structures rising into the topside ionosphere. On nights when only thin layers of bottomside irregularities were observed, moderate to weak scintillations were recorded at VHF. Preliminary values of east-west horizontal irregularity drift velocities were obtained and compared with scintillation rate observations. Using the 1.5-deg and 4.5-deg longitudinal separation between the Jicamarca radar and ionospheric observation points of the two satellites from Huancayo, information was derived regarding large-scale east-west structure during the development phase of the irregularities.

Basu, S.↗