Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling discrepancy”

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 613 records · Page 34

Interference by rain scatter

The data from Japan and the U.S. (the Virginia Precipitation Scatter Experiment) show excellent agreement between the two-component rain scatter model predictions and bistatic scatter measurements. In employing the model, all the scattering geometries should be classified as backscattering as defined by Crane (1974). The forward scatter model should only be used for great circle paths with both antennas pointed at the horizon and at each other in a typical troposcatter communication system geometry. The forward scatter model can also be used for main-lobe, side-lobe coupling when one antenna is pointed toward the other along the great circle path. The forward scatter observations made over the Prospect Hill - Mt Tug path show that the two-component model is incomplete. Much stronger signals were observed at Ku-band than expected based on simultaneous C-band measurements. The discrepancies may be due to: (1) scattering by ice/snow at height (posssible in April) at the 1 km height of the scattering volume), (2) the coherent effects of turbulent fluctuations in the hydrometeor number densities and (3) errors in the modeling of the statistical relationship between attenuation along the path and scattering in the common volume.

Crane, Robert K.↗

Analysis of Gaspra lightcurves using Galileo shape and photometric models

Galileo-based models for the shape of 951 Gaspra and the global-average photometric behavior of its surface have been used to model a representative subset of the asteroid's telescopic lightcurves. Fitting the synthetic lightcurves to the observed timing of lightcurve extrema, and knowing the orientation of Gaspra's axes at the time of the Galileo flyby, leads to a sidereal rotation period for the asteroid of 7.042024 +/- 0.000020 hr, a slight change from the period reported by Magnusson et al. (1992). Initially, the shapes, amplitudes, and absolute photometry of the synthetic and observed lightcurves agree with each other to within 0.05-0.1 mag. Small modifications to the Gaspra shape model on sides of the asteroid poorly imaged by Galileo (changes of 700 m or less in the southern hemisphere at longitudes 90 deg-270 deg W) reduce the typical discrepancies to approximately 0.05 mag in lightcurve shape and less than 0.03 mag in absolute photometry. The result demonstrates that Earth-based lightcurves can be used to refine the shape of a spacecraft-imaged irregular object in areas that are poorly constrained by the spacecraft observations. The consistency and phase-angle dependence of the Galileo-based model for Gaspra photometry, supports the accuracy of the absolute calibration of the Galileo SSI camera, and confirms the Earth-based determination of the V-filter geometric albedo of the asteroid (0.22 +/- 0.03; Tholen et al., submitted for publication). Remaining discrepancies between the synthetic and observed lightcurves show no indication of systematic latitudinal variations in albedo and also cannot be explained entirely by isolated albedo spots. These discrepancies are most likely caused by (1) small, remaining, hard-to-constrain errors in the Gaspra shape model and/or (2) moderate variations in macroscopic roughness across the asteroid's surface, in particular making longitudes 130 deg to 300 deg W moderately rougher than the opposite hemisphere.

Simonelli, Damon P.↗

Zonally averaged thermal balance and stability models for nitrogen polar caps on Triton

Voyager four-color imaging data of Triton are analyzed to calculate the bolometric hemispheric albedo as a function of latitude and longitude. Zonal averages of these data have been incorporated into a thermal balance model involving insolation, reradiation, and latent heat of sublimation of N2 ice for the surface. The current average bolometric albedo of Triton's polar caps is 0.8, implying an effective temperature of 34.2 K and a surface pressure of N2 of 1.6 microbar for unit emissivity. This pressure is an order of magnitude lower than the surface pressure of 18 microbar inferred from Voyager data (Broadfoot et al., 1989; Conrath et al., 1989), a discrepancy that can be reconciled if the emissivity of the N2 on Triton's surface is 0.66. The model predicts that Triton's surface north of 15 deg N latitude is experiencing deposition of N2 frosts, as are the bright portions of the south polar cap near the equator. This result explains why the south cap covers nearly the entire southern hemisphere of Triton.

Stansberry, John A.↗

Modeling of the Dust and Gas Outflows from OH 26.5+0.6: The Superwind

We have observed the extreme OH/IR star, OH 26.5+0.6, in the infrared dust continuum and in the sub- millimeter rotational lines of CO. Mid-infrared images reveal the compact nature of the circumstellar shell (less than 0.5 sec). A deep 9.7 microns absorption feature and an absorption at 18 microns show that the dust mass-loss rate is very high. However, the low antenna temperatures of CO J = 1-0 and 2-1 lines suggest that the outer part of the circumstellar shell is much more tenuous. In order to resolve this discrepancy, we have observed the J = 3-2 and 4-3 CO rotational transitions. We have developed a model for the circumstellar shell for OH 26.5 + 0.6 which is consistent with the infrared and submillimeter observations. The dust and gas data are well fitted by a two-shell model, consisting of a dense shell surrounded by a more tenuous shell. The former we identify with the superwind (M = 5.5 x 10(exp -4) solar mass/ yr), and the latter we identify with mass loss on the asymptotic giant branch (AGB) (M = 10(exp -6) solar mass/ yr). The transition between the two mass-loss phases is shown to be rather abrupt ((Delta)t less than 150 yr). Depending on the mass of the progenitor, this superwind phase may be the last thermal pulse (for M(sub *) less than 1.5 solar mass), or the first of a series of the superwind phases (for up to 8 solar mass), punctuated by a period of low mass-loss rates, before the star evolves off the AGB.

Justtanont, K.↗

Future Changes in Major Stratospheric Warmings in CCMI Models

Major stratospheric warmings (MSWs) are one of the most important phenomena of wintertime Arctic stratospheric variability. They consist of a warming of the Arctic stratosphere and a deceleration of the polar night jet, triggered by an anomalously high injection of tropospheric wave activity into the stratosphere. Due to the relevance and the impact of MSWs on the tropospheric circulation, several model studies have investigated their potential responses to climate change. However, a wide range of results has been obtained, extending from a future increase in the frequency of MSWs to a decrease. These discrepancies might be explained by different factors such as a competition of radiative and dynamical contributors with opposite effects on the Arctic polar vortex, biases of models to reproduce the related processes, or the metric chosen for the identification of MSWs. In this study, future changes in wintertime Arctic stratospheric variability are examined in order to obtaina more precise picture of future changes in the occurrence of MSWs. In particular, transient REFC2 simulations of different CCMs involved in the Chemistry Climate Model Initiative (CCMI) are used. These simulations extend from 1960 to 2100 and include forcings by halogens and greenhouse gases following the specifications of the CCMI-REF-C2 scenario. Sea surface temperatures (SSTs) and sea-ice distributions are either prescribed from coupled climate model integrations or calculated internally in the case of fully coupled atmosphere-ocean CCMs. Potential changes in the frequency and main characteristics of MSWs in the future are investigated with special focus on the dependence of the results on the criterion for the identification of MSWs and the tropospheric forcing of these phenomena.

Chemistry-Climate Modeling↗

Distributed reacceleration of cosmic rays

We develop a model in which cosmic rays, in addition to their initial acceleration by a strong shock, are continuously reaccelerated (e.g., by weak shocks) while propagating through the galaxy. The equations describing this acceleration scheme are solved analytically (approximating ionization losses by a cutoff) and numerically. Solutions for the spectra of primary and secondary cosmic rays are given in a closed analytic form, and they allow a rapid search in parameter space for viable propagation models with distributed reacceleration included. The observed boron-to-carbon ratio can be reproduced by the reacceleration theory over a range of escape parameters, some of them quite different from the standard "leaky box" model. It is also shown that even a very modest amount of reacceleration by strong shocks causes the boron-to carbon ratio to level off at sufficiently high energies, and this effect may be observed in the CRNE data. Several other curiosities in the data may be explained naturally if a modest amount of distributed reacceleration is invoked, including (a) the apparent truncation at low energy in the otherwise exponential pathlength distribution associated with the leaky box model, (b) the sub-iron isotopic anomalies and other effects noted by Silberberg et al., and (c) the discrepancy between the reported 10Be lifetime and the lifetime of cosmic rays in the dense strata of the galactic disk.

NASA Discipline Radiation Health↗

Simulation-Aided Design and Interpretation of Results of Microgravity-Based Solidification Experiments

The preparation and interpretation of microgravity experiments conducted on the International Space Station (ISS) can be challenging because of the remote location of the equipment. Although a “ground unit” is typically available on Earth, it is not located in the laboratory of the principal investigator. Furthermore, test experiments in the ground unit cannot fully duplicate heat transfer conditions of the weightless environment, and thus be a reliable reference for experiment design. Under such conditions, numerical simulations are a valuable tool both for preparing the experiments and interpreting the results obtained under terrestrial and micro-gravity conditions. In this work, we developed a detailed 3D model of the “SUBSA” furnace located at the ISS, and demonstrate its use in our present investigation which includes melt growth and vapor growth of InI crystals on the ISS. In 2002, the SUBSA furnace was designed and built for the growth of InSb crystals in the Microgravity Science Glovebox (MSG) at the ISS. The SUBSA furnace features a transparent section without heating elements, where the axial temperature gradient is high. In our present investigation, we have used the same furnace to grow four InI crystals from the melt, and two InI crystals from the vapor phase (2018-2021). InI crystals have ~20 times lower thermal conductivity than InSb crystals. As a result, during melt growth, the InI crystals were exposed to extremely high temperature gradients. To reduce the axial temperature gradient, we added external jackets made from high thermal conductivity material to the growth ampoules. The SUBSA furnace could not be modeled as asymmetric because of its rectangular shape and right-angled transparent section, which contains a convoluted geometry. The rectangular shape of the transparent section further complicates modelling of radiative heat transfer, which dominates in the zone. At present, numerical simulations are being conducted to resolve an apparent discrepancy between the readings of the 5 thermocouples located in the SUBSA furnace, and the observed position of the solid and molten material. The 3D simulation will be presented, along with the calibration and growth experiment results conducted in microgravity and in the SUBSA ground unit. The developed model should be useful for upcoming microgravity investigations planned to be conducted in the SUBSA facility.

V Riabov↗

A hydrodynamic approach to cosmology: The mixed dark matter cosmological scenario

We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.

Cen, Renyue↗

Studies on effects of boundary conditions in confined turbulent flow predictions

The differences in k epsilon model predictions of plane and axisymmetric expansion flows is investigated. The prediction of the coaxial jet for different velocity ratios of the annular to central jet is presented. The effects of inlet kinetic energy and the energy dissipation rate profiles are investigated for swirling and nonswirling flows. The effects of expansion ration and Reynolds number on the reattachment length are also presented. The results show that the inlet k and epsilon profiles have the most significant effect on the reattachment length and flow redevelopment for the case of coaxial jet of high velocity ratio. A comparison of k epsilon model predictions for the pipe expansion flow by the PHOENICS and TEACH codes reveals some discrepancies in the predicted results. TEACH prediction seems to produce unrealistic kinetic energy profiles in some regions of the flow. PHOENICS code produces a long tail in the recirculation region under certain conditions.

Nallasamy, M.↗

Experimental evaluation of damping of graphite fiber composites

An investigation is described that was conducted to measure the damping of graphite fiber composites such as are to be used on the Hubble Space Telescope. The measurements were made on a composite cylinder under normal atmospheric conditions and in a vacuum chamber with varying thermal conditions. Such composites can outgas under vacuum conditions and this emission of trapped gases may have some influence on the material damping. Some results are presented for the damping coefficient measured under normal atmospheric conditions for a graphite fiber composite. The experimental damping data thus obtained were compared with the least square estimate obtained through a statistical analysis. The statistical analysis was based on a linear mathematical model consisting of a spring in parallel with a dashpot. It is observed that the discrepancy in the damping ratio between logarithmic decrement result and the least square fit is less than 6 percent for all the different tests made on the composite cylinder. Signal-to-noise ratio was maintained to be greater than 10.

Rao, M. D.↗

Spectrophotometry of comets Giacobini-Zinner and Halley

Optical window spectrophotometry was performed on comets Giacobini-Zinner and Halley over the interval 300-1000 nm. Band and band-sequence fluxes were obtained for the brightest features of OH, CN, NH, and C2, special care having been given to determinations of extinction, instrumental sensitivities, and corrections for Fraunhofer lines. C2 Swan band-sequence flux ratios were determined with unprecedented accuracy and compared with the predictions of the detailed equilibrium models of Krishna Swamy et al. (1977, 1979, 1981, and 1987). It is found that these band sequences do not agree with the predictions, which calls into question the assumptions made in deriving the model, namely resonance fluorescence statistical equilibrium. Suggestions are made as to how to resolve this discrepancy.

Tegler, Stephen C.↗

On estimating the basin-scale ocean circulation from satellite altimetry. Part 1: Straightforward spherical harmonic expansion

Direct estimation of the absolute dynamic topography from satellite altimetry has been confined to the largest scales (basically the basin-scale) owing to the fact that the signal-to-noise ratio is more unfavorable everywhere else. But even for the largest scales, the results are contaminated by the orbit error and geoid uncertainties. Recently a more accurate Earth gravity model (GEM-T1) became available, providing the opportunity to examine the whole question of direct estimation under a more critical limelight. It is found that our knowledge of the Earth's gravity field has indeed improved a great deal. However, it is not yet possible to claim definitively that our knowledge of the ocean circulation has improved through direct estimation. Yet, the improvement in the gravity model has come to the point that it is no longer possible to attribute the discrepancy at the basin scales between altimetric and hydrographic results as mostly due to geoid uncertainties. A substantial part of the difference must be due to other factors; i.e., the orbit error, or the uncertainty of the hydrographically derived dynamic topography.

Tai, Chang-Kou↗

The ice band in IRAS 09371 + 1212

Spectra from 2 to 4 microns of the unusual IRAS source 09371 + 1212, a bipolar nebula, reveal a deep H2O ice absorption band near 3.1 microns. No long-wavelength wing is present, implying that the absorbing grain mantles are of high purity. The structure of the absorption feature is reasonably well modeled by a mixture of amorphous and crystalline H2O ice mantles on silicate grains, but discrepancies remain that could indicate other absorbers also.

Geballe, T. R.↗

Growth and characterization of high current density, high-speed InAs/AlSb resonant tunneling diodes

InAs/AlSb double-barrier resonant tunneling diodes with peak current densities up to 370,000 A/sq cm and high peak-to-valley current ratios of 3.2 at room temperature have been fabricated. The peak current density is well-explained by a stationary-state transport model with the two-band envelope function approximation. The valley current density predicted by this model is less than the experimental value by a factor that is typical of the discrepancy found in other double-barrier structures. It is concluded that threading dislocations are largely inactive in the resonant tunneling process.

Soderstrom, J. R.↗

Plate tectonics from VLBI and SLR global data

This study is based on data derived from fifteen years of observations of the SLR (side-looking radar) network and six years of the VLBI (very long baseline interferometry) network. In order to use all available information VLBI and SLR global data sets were combined in a least squares fashion to calculate station horizontal velocities. All significant data pertaining to a single site contribute to the station horizontal motion. The only constraint on the solution is that no vertical motion is allowed. This restriction does not greatly affect the precision of the overall solution given the fact that the expected vertical motion for most stations, even those experiencing post glacial uplift, is well under 1 cm/yr. Since the average baseline is under 4,000 km, only a small fraction of the station vertical velocity is translated into baseline rates so that the error introduced in the solution by restricting up-down station movement is minimal. As a reference, station velocities were then compared to the ones predicted by the NUVEL-1 geological model of DeMets et al. (1990). The focus of the study is on analyzing these discrepancies for global plate tectonics as well as regional tectonic settings. The method used also allows us not only to derive horizontal motion for individual stations but also to calculate Euler vectors for those plates that have enough stations located on the stable interior like North America, Pacific, Eurasia, and Australia.

Harrison, Christopher G. A.↗

Long-Term Trends in the Concentrations of SF6, CHClF2, and COF2 in the Lower Stratosphere from Analysis of High-Resolution Infrared Solar Occultation Spectra

Long-term trends in the concentrations of SF6, CHClF2 (CFC-22), and COF2 in the lower stratosphere have been derived from analysis of ca. 1980 and more recent infrared solar occultation spectra recorded near 32 deg N latitude at approx. 0.02/ cm resolution. Consistent sets of line parameters and spectral calibration methods have been used in the retrievals to minimize systematic error effects. Quoted error limits are 1 sigma estimated precisions. The SF6 and CHClF2 results are based on spectra recorded by balloon-borne interferometers in March 1981 and June 1988 and a comparison of these results with the Atmospheric Trace Molecule Spectroscopy (ATMOS) Experiment/Spacelab 3 measurements obtained in May 1985 near 30 deg N latitude. In the 13-18 km altitude range the mean measured SF6 mixing ratio in parts per trillion by volume (pptv) increased from 1.17 +/- 0.21 in March 1981 to 2.02 +/- 0.20 pptv in June 1988, and the CHClF2 mixing ratio below 15 km altitude increased from 51 +/- 8 pptv in March 1981 to 102 +/- 10 pptv in June 1988. The CHClF2 retrievals used new empirical CHClF2 line parameters derived from 0.03/cm resolution laboratory spectra recorded at six temperatures between 203 and 293 K; the derived mixing ratios are approx. 30% higher than obtained with earlier sets of line parameters, thereby removing a large discrepancy noted previously between IR and in situ measurements of CHClF2. Assuming an exponential growth model for fitting the trends, SF6 and CHClF2 mean increase rates of 7.4% +/- 1.9% and 9.4% +/- 1.3% /year, are obtained, respectively, which correspond to cumulative increases by factors of approx. 1.7 and -2.0 in the concentrations of these gases over the 7.2-year measurement period. Analysis of spectra recorded in October 1979 and April 1989 yields COF2 volume mixing ratios that are respectively 0.44 +/- 0.17 and 1.21 +/- 0.24 times the ATMOS/Spacelab 3 values, from which an average COF2 increase rate of 10.3 +/- 1.8%/ year over this time period has been estimated. The present results are compared with previously reported observations and trends and with one-dimensional model calculations. The model calculated trends are in reasonably good agreement with the observations.

Rinsland, C. P.↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗