Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Real--Time Dynamic Modeling”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Condensate Accretion in Shock Tube's Expansion Fan

It has been shown that turbulence and temperature influence the droplet sizes in expansion fan induced condensation by studying the Rayleigh scattering from one port in our shock tube's test section. We have modified our set-up so as to allow, using two ports, the real time measurement of the influence of turbulence and temperature on the rate at which these droplets grow. To do this, we looked at the Rayleigh scattering from two different ports for ten Reynolds numbers at five different temperatures. We modeled the time of flight of droplets, using the equations of one-dimensional gas dynamics and the measured shock wave speed in shock tube's driven section.

Mezonlin, Ephrem-Denis↗
NASA NTRS1998

Time-Resolved O3 Chemical Chain Reaction Kinetics Via High-Resolution IR Laser Absorption Methods

Excimer laser photolysis in combination with time-resolved IR laser absorption detection of OH radicals has been used to study O3/OH(v = 0)/HO2 chain reaction kinetics at 298 K, (i.e.,(k(sub 1) is OH + 03 yields H02 + 02 and (k(sub 2) is H02 + 03 yields OH + 202). From time-resolved detection of OH radicals with high-resolution near IR laser absorption methods, the chain induction kinetics have been measured at up to an order of magnitude higher ozone concentrations ([03] less than or equal to 10(exp 17) molecules/cu cm) than accessible in previous studies. This greater dynamic range permits the full evolution of the chain induction, propagation, and termination process to be temporally isolated and measured in real time. An exact solution for time-dependent OH evolution under pseudo- first-order chain reaction conditions is presented, which correctly predicts new kinetic signatures not included in previous OH + 03 kinetic analyses. Specifically, the solutions predict an initial exponential loss (chain "induction") of the OH radical to a steady-state level ([OH](sub ss)), with this fast initial decay determined by the sum of both chain rate constants, k(sub ind) = k(sub 1) + k(sub 2). By monitoring the chain induction feature, this sum of the rate constants is determined to be k(sub ind) = 8.4(8) x 10(exp -14) cu cm/molecule/s for room temperature reagents. This is significantly higher than the values currently recommended for use in atmospheric models, but in excellent agreement with previous results from Ravishankara et al.

Kulcke, Axel↗
NASA NTRS2000

Towards a Comprehensive Dynamic-chemistry Assimilation for Eos-Chem: Plans and Status in NASA's Data Assimilation Office

In order to support the EOS-Chem project, a comprehensive assimilation package for the coupled chemical-dynamical system is being developed by the Data Assimilation Office at NASA GSFC. This involves development of a coupled chemistry/meteorology model and of data assimilation techniques for trace species and meteorology. The model is being developed using the flux-form semi-Lagrangian dynamical core of Lin and Rood, the physical parameterizations from the NCAR Community Climate Model, and atmospheric chemistry modules from the Atmospheric Chemistry and Dynamics branch at NASA GSFC. To date the following results have been obtained: (i) multi-annual simulations with the dynamics-radiation model show the credibility of the package for atmospheric simulations; (ii) initial simulations including a limited number of middle atmospheric trace gases reveal the realistic nature of transport mechanisms, although there is still a need for some improvements. Samples of these results will be shown. A meteorological assimilation system is currently being constructed using the model; this will form the basis for the proposed meteorological/chemical assimilation package. The latter part of the presentation will focus on areas targeted for development in the near and far terms, with the objective of Providing a comprehensive assimilation package for the EOS-Chem science experiment. The first stage will target ozone assimilation. The plans also encompass a reanalysis (ReSTS) for the 1991-1995 period, which includes the Mt. Pinatubo eruption and the time when a large number of UARS observations were available. One of the most challenging aspects of future developments will be to couple theoretical advances in tracer assimilation with the practical considerations of a real environment and eventually a near-real-time assimilation system.

Pawson, Steven↗
NASA NTRS2000

Nozomi Cis-Lunar Phase Orbit Determination

Japan's Institute of Space and Astronautical Science (ISAS) launched Nozomi, its first mission to the planet Mars using the newly developed M-V launch vehicle on July 3, 1998. Scientific objectives of the mission are to study the structure and dynamics of the Martian upper atmosphere and its interaction with the solar wind. Nozomi is a cooperative mission between ISAS and the National Aeronautics and Space Administration (NASA). The NASA contribution includes navigation and tracking services provided by the Jet Propulsion Laboratory (JPL). The spacecraft also serves as an engineering demonstration of basic technology for planetary exploration. One of the new technologies was a unique trajectory, developed by ISAS, which used solar gravitational perturbations at the weak stability boundary as an aid to achieve an Earth-Mars transfer orbit. This trajectory saves approximately 120 m/s of Delta V compared to direct hyperbolic insertion and is considered an enabling technology for the mission. Nozomi was the first spacecraft to employ this trajectory and provided on-orbit validation of the technique. The trajectory was achieved by initially placing the spacecraft in a highly elliptical cis-lunar phasing orbit. Six maneuvers were performed during this period to correct injection errors and target an outbound lunar swingby in September 1998. The gravity assist from the lunar swingby raised apogee to the vicinity of the weak stability boundary. After three more targeting maneuvers, Nozomi performed an inbound lunar swingby followed immediately by a powered Earth swingby in late December 1998. A 420 m/s Trans Mars Insertion (TMI) burn at the final Earth periapsis was intended to place the spacecraft on a heliocentric trajectory leading to Mars orbit insertion in October 1999. Orbit determination for Nozomi is performed in parallel by both ISAS and the Multi-Mission Navigation (MMNAV) group at JPL. This was an advantage for the mission because each group would generate solutions based on data collected from their respective tracking networks. Spacecraft events, such as sequence uplinks and maneuvers, were generally scheduled during passes at the Usuda tracking station in Japan. As a result, maneuver design and reconstruction was derived from MMNAV solutions based on JPL tracking data obtained immediately prior to or following maneuvers. Data was also exchanged between ISAS and MMNAV so orbit determination could be performed on joint data sets in support of critical targeting late in the cis-lunar phase. In this paper, information regarding the MMNAV orbit determination effort for the first six months of the mission is presented. The spacecraft trajectory is characterized first, followed by a discussion of the orbit determination estimation procedure and models. Results from selected orbit solutions are presented and compared against reconstructed trajectories. One area of emphasis in this paper is orbit determination in the vicinity of the weak stability boundary. Precise navigation was necessary to target the second lunar swingby and the powered Earth swingby. Delivery accuracy of 150 m was required for these critical encounters, but a number of factors contributed to the general degradation of orbit determination accuracy. This included the fact that the spacecraft was at apogee, at a range of 1.7 million km and moving at less than I km/sec perpendicular to the line of sight. Nozomi was also close to zero degrees declination where there are known limitations on orbit determination performance. Finally, S-band tracking data was acquired through the Nozomi backup low gain antenna. This antenna is offset from the axis of this spin stabilized spacecraft and superimposed large signatures in the Doppler and range data. These difficulties were overcome by combining long data arcs, spanning several maneuvers, with a high fidelity solar pressure model. The model included a physically accurate representation of the spacecraft structure and a high time resolution orientation model. Observation modeling included the removal of the spin induced Doppler bias, spin signature and per pass correction of range calibration errors applied for data leading up to critical events. As a result, all orbit determination goals were met. A second area of emphasis in this paper is the JPL tracking and orbit determination effort in support of the TMI maneuver. TMI occurred out of contact with ground stations and the JPL Goldstone tracking complex had the first pass following the bum. As a result, MMNAV had the responsibility to make a rapid assessment of the maneuver performance. MMNAV made the determination that a 100 m/s under bum had occurred and promptly informed ISAS via voice lines. ISAS immediately began preparations for a correction maneuver (TMIc), which had to be performed during the next Usuda pass. The near real time assessment by MMNAV provided accurate antenna frequency and pointing updates for the spacecraft acquisition at Usuda and the close coordination between the two agencies enabled the design and successful execution of the TMc maneuver. Propellant consumption during the correction burn dictated that the mission be redesigned. ISAS developed a new plan which adds 3 full solar orbits, two Earth swingbys and one lunar swingby with arrival at Mars in January 2004. The final Mars orbit will still enable the mission to achieve all of its science objectives.

Ryne, Mark↗
NASA NTRS2000

On-Board Real-Time State and Fault Identification for Rovers

For extended autonomous operation, rovers must identify potential faults to determine whether its execution needs to be halted or not. At the same time, rovers present particular challenges for state estimation techniques: they are subject to environmental influences that affect senior readings during normal and anomalous operation, and the sensors fluctuate rapidly both because of noise and because of the dynamics of the rover's interaction with its environment. This paper presents MAKSI, an on-board method for state estimation and fault diagnosis that is particularly appropriate for rovers. The method is based on a combination of continuous state estimation, wing Kalman filters, and discrete state estimation, wing a Markov-model representation.

Washington, Richard↗
NASA NTRS2002

Mars Smart Lander Simulations for Entry, Descent, and Landing

Two primary simulations have been developed and are being updated for the Mars Smart Lander Entry, Descent, and Landing (EDL). The high fidelity engineering end-to-end EDL simulation that is based on NASA Langley's Program to Optimize Simulated Trajectories (POST) and the end-to-end real-time, hardware-in-the-loop simulation testbed, which is based on NASA JPL's (Jet Propulsion Laboratory) Dynamics Simulator for Entry, Descent and Surface landing (DSENDS). This paper presents the status of these Mars Smart Lander EDL end-to-end simulations at this time. Various models, capabilities, as well as validation and verification for these simulations are discussed.

Striepe, S. A.↗
NASA NTRS2002

Adaptive Identification and Control of Flow-Induced Cavity Oscillations

Progress towards an adaptive self-tuning regulator (STR) for the cavity tone problem is discussed in this paper. Adaptive system identification algorithms were applied to an experimental cavity-flow tested as a prerequisite to control. In addition, a simple digital controller and a piezoelectric bimorph actuator were used to demonstrate multiple tone suppression. The control tests at Mach numbers of 0.275, 0.40, and 0.60 indicated approx. = 7dB tone reductions at multiple frequencies. Several different adaptive system identification algorithms were applied at a single freestream Mach number of 0.275. Adaptive finite-impulse response (FIR) filters of orders up to N = 100 were found to be unsuitable for modeling the cavity flow dynamics. Adaptive infinite-impulse response (IIR) filters of comparable order better captured the system dynamics. Two recursive algorithms, the least-mean square (LMS) and the recursive-least square (RLS), were utilized to update the adaptive filter coefficients. Given the sample-time requirements imposed by the cavity flow dynamics, the computational simplicity of the least mean squares (LMS) algorithm is advantageous for real-time control.

Kegerise, M. A.↗
NASA NTRS2003

Multiple Input Design for Real-Time Parameter Estimation in the Frequency Domain

A method for designing multiple inputs for real-time dynamic system identification in the frequency domain was developed and demonstrated. The designed inputs are mutually orthogonal in both the time and frequency domains, with reduced peak factors to provide good information content for relatively small amplitude excursions. The inputs are designed for selected frequency ranges, and therefore do not require a priori models. The experiment design approach was applied to identify linear dynamic models for the F-15 ACTIVE aircraft, which has multiple control effectors.

Morelli, Eugene↗
NASA NTRS2004

Real-Time Measurements of Aft Dome Insulation Erosion on Space Shuttle Reusable Solid Rocket Motor

Real-time erosion of aft dome internal insulation was measured with internal instrumentation on a static test of a lengthened version of the Space Shuffle Reusable Solid Rocket Motor (RSRM). This effort marks the first time that real-time aft dome insulation erosion (Le., erosion due to the combined effects of thermochemical ablation and mechanical abrasion) was measured in this kind of large motor static test [designated as Engineering Test Motor number 3 (ETM3)I. This paper presents data plots of the erosion depth versus time. The data indicates general erosion versus time behavior that is in contrast to what would be expected from earlier analyses. Engineers have long known that the thermal environment in the aft dome is severe and that the resulting aft dome insulation erosion is significant. Models of aft dome erosion involve a two-step process of computational fluid dynamics (CFD) modeling and material ablation modeling. This modeling effort is complex. The time- dependent effects are difficult to verify with only prefire and postfire insulation measurements. Nozzle vectoring, slag accumulation, and changing boundary conditions will affect the time dependence of aft dome erosion. Further study of this data and continued measurements on future motors will increase our understanding of the aft dome flow and erosion environment.

McWhorter, Bruce↗
NASA NTRS2002

Dynamical Studies of the Middle Atmosphere Using High Resolution Doppler Imager Observations

This report summarizes the activities of NASA grant NAG5-11068, "Dynamicai Studies of the Middle Atmosphere Using High Resolution Doppler Imager Observations." The High Resolution Doppler Imager (HRDI) on the Upper Atmosphere Research Satellite (UARS) has been providing direct measurements of the Earth's horizontal wind field in the stratosphere, mesosphere and lower thermosphere. Mesospheric temperatures, ozone, and O((sup 1) D) densities, and stratospheric aerosol extinctions coefficients, are also retrieved. The goal of HRDI is to measure the vector winds in the stratosphere (10-40 km), mesosphere, and lower thermosphere (approximately 50-120 km) during the day, and the lower thermosphere at night (approximately 95 km) to an accuracy of 5 m/s. The horizontal wind vector is measured by observing the Doppler shift of rotational lines of molecular oxygen along two lines of sight. In addition to winds, temperatures and volume emission rates are determined in the mesosphere and lower thermosphere, from which ozone and O((sup 1) D) concentrations can be derived, and aerosol scattering coefficients are determined in the stratosphere. UARS was launched on September 12, 1991, into a 585-km circular orbit inclined 57 degrees to the equator HRDI was activated September 28, 1991 and following a period of checkout and adjustment of the instrument parameters, scientific observations began November 2, 199 1. HRDI operated nearly continuously from launch until April 1995. At that time the UARS solar array drive failed, forcing the instruments to time-share the available power. From July 1995 to July 1996 HRDI operated approximately 50% of the time. At that point, one of the three spacecraft batteries failed and from then until September 1998 the duty cycle was less than 20% per month, At that time it was determined that HRDI could operate during each daytime pass, which increased the daytime duty cycle to close to l00%, while nighttime operations were limited to about a week per month. In the fall of 1999, the second tape recorder failed requiring a real time contact with a TRDSS satellite to retrieve that data. This resulted in about 60% data collection efficiency. Finally, in the summer of 2000, the second star sensor failed requiring the spacecraft attitude to be controlled by a three axis magnetometer and sun sensor. This resulted in a loss of attitude knowledge but operations continue with the anticipation of correcting the attitude. A new method for determining the tide and mean structure from satellite data in conjunction with a new tidal model has been devised. For brevity, it shall be referred to as the TMAT or Tide-Mean Assimilation Technique. Most previous methods of tidal analysis are based on various ways of slicing the data set.

Skinner, Wilbert↗
NASA NTRS2005

A Generic Inner-Loop Control Law Structure for Six-Degree-of-Freedom Conceptual Aircraft Design

A generic control system framework for both real-time and batch six-degree-of-freedom (6-DOF) simulations is presented. This framework uses a simplified dynamic inversion technique to allow for stabilization and control of any type of aircraft at the pilot interface level. The simulation, designed primarily for the real-time simulation environment, also can be run in a batch mode through a simple guidance interface. Direct vehicle-state acceleration feedback is required with the simplified dynamic inversion technique. The estimation of surface effectiveness within real-time simulation timing constraints also is required. The generic framework provides easily modifiable control variables, allowing flexibility in the variables that the pilot commands. A direct control allocation scheme is used to command aircraft effectors. Primary uses for this system include conceptual and preliminary design of aircraft, when vehicle models are rapidly changing and knowledge of vehicle 6-DOF performance is required. A simulated airbreathing hypersonic vehicle and simulated high-performance fighter aircraft are used to demonstrate the flexibility and utility of the control system.

Cox, Timothy H.↗
NASA NTRS2004

CDGPS-Based Relative Navigation for Multiple Spacecraft

This thesis investigates the use of Carrier-phase Differential GPS (CDGPS) in relative navigation filters for formation flying spacecraft. This work analyzes the relationship between the Extended Kalman Filter (EKF) design parameters and the resulting estimation accuracies, and in particular, the effect of the process and measurement noises on the semimajor axis error. This analysis clearly demonstrates that CDGPS-based relative navigation Kalman filters yield good estimation performance without satisfying the strong correlation property that previous work had associated with "good" navigation filters. Several examples are presented to show that the Kalman filter can be forced to create solutions with stronger correlations, but these always result in larger semimajor axis errors. These linear and nonlinear simulations also demonstrated the crucial role of the process noise in determining the semimajor axis knowledge. More sophisticated nonlinear models were included to reduce the propagation error in the estimator, but for long time steps and large separations, the EKF, which only uses a linearized covariance propagation, yielded very poor performance. In contrast, the CDGPS-based Unscented Kalman relative navigation Filter (UKF) handled the dynamic and measurement nonlinearities much better and yielded far superior performance than the EKF. The UKF produced good estimates for scenarios with long baselines and time steps for which the EKF would diverge rapidly. A hardware-in-the-loop testbed that is compatible with the Spirent Simulator at NASA GSFC was developed to provide a very flexible and robust capability for demonstrating CDGPS technologies in closed-loop. This extended previous work to implement the decentralized relative navigation algorithms in real time.

Mitchell, Megan Leigh↗
NASA NTRS2005

A Generic Inner-Loop Control Law Structure for Six-Degree-of-Freedom Conceptual Aircraft Design

A generic control system framework for both real-time and batch six-degree-of-freedom simulations is presented. This framework uses a simplified dynamic inversion technique to allow for stabilization and control of any type of aircraft at the pilot interface level. The simulation, designed primarily for the real-time simulation environment, also can be run in a batch mode through a simple guidance interface. Direct vehicle-state acceleration feedback is required with the simplified dynamic inversion technique. The estimation of surface effectiveness within real-time simulation timing constraints also is required. The generic framework provides easily modifiable control variables, allowing flexibility in the variables that the pilot commands. A direct control allocation scheme is used to command aircraft effectors. Primary uses for this system include conceptual and preliminary design of aircraft, when vehicle models are rapidly changing and knowledge of vehicle six-degree-of-freedom performance is required. A simulated airbreathing hypersonic vehicle and simulated high-performance fighter aircraft are used to demonstrate the flexibility and utility of the control system.

Cox, Timothy H.↗
NASA NTRS2005

Motion Cueing Algorithm Development: Human-Centered Linear and Nonlinear Approaches

While the performance of flight simulator motion system hardware has advanced substantially, the development of the motion cueing algorithm, the software that transforms simulated aircraft dynamics into realizable motion commands, has not kept pace. Prior research identified viable features from two algorithms: the nonlinear "adaptive algorithm", and the "optimal algorithm" that incorporates human vestibular models. A novel approach to motion cueing, the "nonlinear algorithm" is introduced that combines features from both approaches. This algorithm is formulated by optimal control, and incorporates a new integrated perception model that includes both visual and vestibular sensation and the interaction between the stimuli. Using a time-varying control law, the matrix Riccati equation is updated in real time by a neurocomputing approach. Preliminary pilot testing resulted in the optimal algorithm incorporating a new otolith model, producing improved motion cues. The nonlinear algorithm vertical mode produced a motion cue with a time-varying washout, sustaining small cues for longer durations and washing out large cues more quickly compared to the optimal algorithm. The inclusion of the integrated perception model improved the responses to longitudinal and lateral cues. False cues observed with the NASA adaptive algorithm were absent. The neurocomputing approach was crucial in that the number of presentations of an input vector could be reduced to meet the real time requirement without degrading the quality of the motion cues.

Houck, Jacob A.↗
NASA NTRS2005

Diagnostics of Carbon Nanotube Formation in a Laser Produced Plume: Spectroscopic in situ nanotube detection using spectral absorption and surface temperature measurements by black body emission

Carbon nanotubes hold great promise for material advancements in the areas of composites and electronics. The advancement of research in these areas is dependent upon the availability of carbon nanotubes to a broad spectrum of academic and industrial researchers. Although there has been much progress made in reducing the costs of carbon nanotubes and increasing the quality and purity of the products, an increase in demand for still less expensive and specific nanotubes types has also grown. This summer's work has involved two experiments that have been designed to further the understanding of the dynamics and chemical mechanisms of carbon nanotube formation. It is expected that a better understanding of the process of formation of nanotubes will aid current production designs and stimulate ideas for future production designs increasing the quantity, quality, and production control of carbon nanotubes. The first experiment involved the measurement of surface temperature of the target as a function of time with respect to the ablation lasers. A peak surface temperature of 5000 K was determined from spectral analysis of black body emission from the target surface. The surface temperature as a function of various changes in operating parameters was also obtained. This data is expected to aid the modeling of ablation and plume dynamics. The second experiment involved a time and spatial measurement of the spectrally resolved absorbance of the laser produced plume. This experiment explored the possibility of developing absorbance and fluorescence to detect carbon nanotubes during production. To attain control over the production of nanotubes with specific properties and reduce costs, a real time in situ diagnostics method would be very beneficial. Results from this summer's work indicate that detection of nanotubes during production may possibly be used for production feed back control.

DeBoer, Gary D.↗
NASA NTRS2005

On-Line Loss of Control Detection Using Wavelets

Wavelet transforms are used for on-line detection of aircraft loss of control. Wavelet transforms are compared with Fourier transform methods and shown to more rapidly detect changes in the vehicle dynamics. This faster response is due to a time window that decreases in length as the frequency increases. New wavelets are defined that further decrease the detection time by skewing the shape of the envelope. The wavelets are used for power spectrum and transfer function estimation. Smoothing is used to tradeoff the variance of the estimate with detection time. Wavelets are also used as front-end to the eigensystem reconstruction algorithm. Stability metrics are estimated from the frequency response and models, and it is these metrics that are used for loss of control detection. A Matlab toolbox was developed for post-processing simulation and flight data using the wavelet analysis methods. A subset of these methods was implemented in real time and named the Loss of Control Analysis Tool Set or LOCATS. A manual control experiment was conducted using a hardware-in-the-loop simulator for a large transport aircraft, in which the real time performance of LOCATS was demonstrated. The next step is to use these wavelet analysis tools for flight test support.

Brenner, Martin J.↗
NASA NTRS2005

Small Body GN&C Research Report: A Guidance and Control Technique for Small-Body Proximity Operations with Guaranteed Guidance Resolvability and Required Thruster Silent Time

The guidance and control (G&C) algorithms for enabling small-body proximity operations are developed by using a model predictive control approach along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. The open-loop guidance is solved ahead of time or in a resolvable, real-time manner through the use of PWG (Pseudo Way-point Generation), a technique developed in this research. The PWG scheme ensures required thruster silent times during trajectory maneuvers. The feedback control is implemented to track the PWG trajectories in a manner that guarantees the resolvability for the open-loop problem, enabling the ability to update the G&C in a model-predictive manner. The schemes incorporate gravity models and thruster ring times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure. The optimal control problem is cast as an LMI (Linear Matrix Inequality) and then solved through Semi-Definite Programming techniques in a computationally efficient manner that provides convergence and constraint guarantees.

guidance and control (G&C)↗
NASA NTRS2004

Analysis of SSEM Sensor Data Using BEAM

A report describes analysis of space shuttle main engine (SSME) sensor data using Beacon-based Exception Analysis for Multimissions (BEAM) [NASA Tech Briefs articles, the two most relevant being Beacon-Based Exception Analysis for Multimissions (NPO- 20827), Vol. 26, No.9 (September 2002), page 32 and Integrated Formulation of Beacon-Based Exception Analysis for Multimissions (NPO- 21126), Vol. 27, No. 3 (March 2003), page 74] for automated detection of anomalies. A specific implementation of BEAM, using the Dynamical Invariant Anomaly Detector (DIAD), is used to find anomalies commonly encountered during SSME ground test firings. The DIAD detects anomalies by computing coefficients of an autoregressive model and comparing them to expected values extracted from previous training data. The DIAD was trained using nominal SSME test-firing data. DIAD detected all the major anomalies including blade failures, frozen sense lines, and deactivated sensors. The DIAD was particularly sensitive to anomalies caused by faulty sensors and unexpected transients. The system offers a way to reduce SSME analysis time and cost by automatically indicating specific time periods, signals, and features contributing to each anomaly. The software described here executes on a standard workstation and delivers analyses in seconds, a computing time comparable to or faster than the test duration itself, offering potential for real-time analysis.

Zak, Michail↗