Search NASA⌕ Search

SEARCH · Search NASA

Results for “ENERGY METHOD STRESS CALCULATION”

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.

69 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A numerical study of two-dimensional vortex shedding from rectangular cylinders

An efficient time-marching, non-iterative calculation method is used to analyze time-dependent flows around rectangular cylinders. The turbulent flow in the wake region of a square section cylinder is analyzed using an anisotropic k-epsilon model. Initiation and subsequent development of the vortex shedding phenomenon is naturally captured once a perturbation is introduced in the flow. Transient calculations using standard eddy-viscosity and an anisotropic k-epsilon model averaged over an integral number of cycles to get the fluctuating energy (organized and turbulent) are compared with experimental data. It is shown that the anisotropic k-epsilon model resolves the anisotropy of the Reynolds stresses and gives mean energy distribution closer to the experiment than the standard k-epsilon model.

Hadid, A. H.↗
NASA NTRS1973

Diagnostics of sources of disturbances and distribution of vibrations over the width of a tape in tape-feed mechanisms

Disturbances created by certain assemblies and components of tape-feed mechanisms (TFM) and acting on a moving magnetic tape are studied. The method, based on elements of digital logic, is established by stress-strain diagrams of the longitudinal deformations and vibrations across the width of a magnetic tape. Experimental studies were carried out for determination of the functional relationships of longitudinal deformations in a section of magnetic tape to the magnitude of roller play, tension vibrations rate of movement, and elasticity of magnetic tapes. A block diagram of the measurements is shown. Appropriate digital computer algorithms and programs were proposed for statistical analysis of the data obtained. Estimates of the mathematical expectation, dispersion, intercorrelation function, energy spectral density, and distribution pattern of the random process values were calculated.

Kenstavichyus, A. B. B.↗
NASA NTRS2004

Fracture Mechanics Analysis of an Annular Crack in a Three-concentric-cylinder Composite Model

A boundary-value problem governing a three-phase concentric-cylinder model was analytically modeled to analyze annular interfacial crack problems with Love s strain functions in order to find the stress intensity factors (SIFs) and strain energy release rates (SERRs) at the tips of an interface crack in a nonhomogeneous medium. The complex form of a singular integral equation (SIE) of the second kind was formulated using Bessel s functions in the Fourier domain, and the SIF and total SERR were calculated using Jacoby polynomials. For the validity of the SIF equations to be established, the SIE of the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the results were compared with the previous results of Erdogan. A preliminary set of parametric studies was carried out to show the effect of interphase properties on the SERR. The method presented here provides insight about the effect of interphase properties on the crack driving force.

Kuguoglu, Latife H.↗
NASA NTRS1995

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗
NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗
NASA NTRS1999

Fracture Mechanics of Thin, Cracked Plates Under Tension, Bending and Out-of-Plane Shear Loading

Cracks in the skin of aircraft fuselages or other shell structures can be subjected to very complex stress states, resulting in mixed-mode fracture conditions. For example, a crack running along a stringer in a pressurized fuselage will be subject to the usual in-plane tension stresses (Mode-I) along with out-of-plane tearing stresses (Mode-III like). Crack growth and initiation in this case is correlated not only with the tensile or Mode-I stress intensity factor, K(sub I), but depends on a combination of parameters and on the history of crack growth. The stresses at the tip of a crack in a plate or shell are typically described in terms of either the small deflection Kirchhoff plate theory. However, real applications involve large deflections. We show, using the von-Karman theory, that the crack tip stress field derived on the basis of the small deflection theory is still valid for large deflections. We then give examples demonstrating the exact calculation of energy release rates and stress intensity factors for cracked plates loaded to large deflections. The crack tip fields calculated using the plate theories are an approximation to the actual three dimensional fields. Using three dimensional finite element analyses we have explored the relationship between the three dimensional elasticity theory and two dimensional plate theory results. The results show that for out-of-plane shear loading the three dimensional and Kirchhoff theory results coincide at distance greater than h/2 from the crack tip, where h/2 is the plate thickness. Inside this region, the distribution of stresses through the thickness can be very different from the plate theory predictions. We have also explored how the energy release rate varies as a function of crack length to plate thickness using the different theories. This is important in the implementation of fracture prediction methods using finite element analysis. Our experiments show that under certain conditions, during fatigue crack growth, the presence of out-of-plane shear loads induces a great deal of contact and friction on the crack surfaces, dramatically reducing crack growth rate. A series of experiments and a proposed computational approach for accounting for the friction is discussed.

Zehnder, Alan T.↗
NASA NTRS1988

Tensile and creep rupture behavior of P/M processed Nb-base alloy, WC-3009

Due to its high strength at temperatures up to 1600 K, fabrication of niobium base alloy WC-3009 (Nb30Hf9W) by traditional methods is difficult. Powder metallurgy (P/M) processing offers an attractive fabrication alternative for this high strength alloy. Spherical powders of WC-3009 produced by electron beam atomizing (EBA) process were successfully consolidated into a one inch diameter rod by vacuum hot pressing and swaging techniques. Tensile strength of the fully dense P/M material at 300-1590 K were similar to the arc-melted material. Creep rupture tests in vacuum indicated that WC-3009 exhibits a class 1 solid solution (glide controlled) creep behavior in the 1480 to 1590 K temperature range and stress range of 14 to 70 MPa. The creep behavior was correlated with temperature and stress using a power law relationship. The calculated stress exponent n, was about 3.2 and the apparent activation energy, Q, was about 270 kJ/mol. The large creep ductility exhibited by WC-3009 was attributed to its high strain rate sensitivity.

Hebsur, Mohan G.↗
NASA NTRS1986

Analysis of a resistance-energy balance method for estimating daily evaporation from wheat plots using one-time-of-day infrared temperature observations

Accurate estimates of evaporation over field-scale or larger areas are needed in hydrologic studies, irrigation scheduling, and meteorology. Remotely sensed surface temperature might be used in a model to calculate evaporation. A resistance-energy balance model, which combines an energy balance equation, the Penman-Monteith (1981) evaporation equation, and van den Honert's (1948) equation for water extraction by plant roots, is analyzed for estimating daily evaporation from wheat using postnoon canopy temperature measurements. Additional data requirements are half-hourly averages of solar radiation, air and dew point temperatures, and wind speed, along with reasonable estimates of canopy emissivity, albedo, height, and leaf area index. Evaporation fluxes were measured in the field by precision weighing lysimeters for well-watered and water-stressed wheat. Errors in computed daily evaporation were generally less than 10 percent, while errors in cumulative evaporation for 10 clear sky days were less than 5 percent for both well-watered and water-stressed wheat. Some results from sensitivity analysis of the model are also given.

Choudhury, B. J.↗
NASA NTRSDate not supplied

Certifying Acoustic Emission for the Crack Inspection of Integrated Pressure Components

To address the limitations of other nondestructive evaluation (NDE) methods on integrated components due to access constraints, a novel Acoustic Emission (AE) pressure vessel health test was developed in accordance with ASME BPVC.V Article 12. This test was then validated using notched Ti-6Al-4V panels taken to a stress level equivalent to 1.1 times the Maximum Allowable Working Pressure (MAWP) of a flight pressure vessel. The assessment involved utilizing JSC small sample Probability of Detection (POD) calculations, which represents the initial step towards obtaining flight certification for an NDE technique. A robust filtering technique that specifically identifies crack growth was formulated with the assistance of industry partners. This detection method exhibited superior detection ability compared to other NDE methods or microscopy. Within industry, the presence of a group of triangulated energetic crack-like indications typically leads to supplemental inspection or the failure of a tank. The experimental results confirm that this approach, commonly applied to more robust tanks, can be viably extended to integrated aerospace pressure components in scenarios where components can be subjected to pressures up to 1.1 times the MAWP. The testing process encountered delays due to sensor delivery and panel fabrication. Thus far, samples with crack initiators have shown energy levels 10 times higher than blunted controls. This suggests a trend where high individual and total cluster energy could be employed as acceptance criteria. Moreover, the cumulative energy of clusters might serve as a foundation for estimating crack sizes; however, this aspect requires further assessment and localized supplemental inspections are preferred.

certifying↗
NASA NTRS2008

Assessment of Proper Bonding Methods and Mechanical Characterization FPGA CQFPs

This presentation discusses fractured leads on field-programmable gate array (FPGA) during flight vibration. Actions taken to determine root cause and resolution of the failure include finite element analysis (FEA) and vibration testing and scanning electron microscopy (with X-ray microanalysis) and energy dispersive spectrometry (SEM/EDS) failure assessment. Bonding methods for surface mount parts is assessed, including critical analysis and assessment of random fatigue damage. Regarding ceramic quad flat pack (CQFP) lead fracture, after disassembling the attitude control electronics (ACE) configuration, photographs showed six leads cracked on FPGA RTSX72SU-1 CQ208B package located on the RWIC card. An identical package (FPGA RTSX32SU-1 CQ208B) mounted on the RWIC did not results in cracked pins due to vibration. FPGA lead failure theories include workmanship issues in the lead-forming, material defect in the leads of the FPGA packages, and the insecure mounting of the board in the card guides, among other theories. Studies were conducted using simple calculations to determine the response and fatigue life of the package. Shorter packages exhibited more response when loaded by out-of-plane displacement of PCB while taller packages exhibit more response when loaded by in-plane acceleration of PCB. Additionally, under-fill did not contribute to reducing stress in leads due to out-of-plane PCB loading or from component twisting, as much as corner bonding. The combination of corner bond and under-fill is best to address mechanical and thermal S/C environment. Test results of bonded parts showed reduced (dampened) amplitude and slightly shifted peaks at the un-bonded natural frequency and an additional response at the bonded frequency. Stress due to PCBB out-of-plane loading was decreased on in the corners when only a corner bond was used. Future work may address CQFP fatigue assessment, including the investigation of discrepancy in predicted fatigue damage, as well as comparing fatigue life and fatigue damage cycle ration computed using FEA and Miner's rule to results from a fatigue assessment software program.

Davis, Milton C.↗
NASA NTRS1996

Structure Identification Within a Transitioning Swept-Wing Boundary Layer

Extensive measurements are made in a transitioning swept-wing boundary layer using hot-film, hot-wire and cross-wire anemometry. The crossflow-dominated flow contains stationary vortices that breakdown near mid-chord. The most amplified vortex wavelength is forced by the use of artificial roughness elements near the leading edge. Two-component velocity and spanwise surface shear-stress correlation measurements are made at two constant chord locations, before and after transition. Streamwise surface shear stresses are also measured through the entire transition region. Correlation techniques are used to identify stationary structures in the laminar regime and coherent structures in the turbulent regime. Basic techniques include observation of the spatial correlations and the spatially distributed auto-spectra. The primary and secondary instability mechanisms are identified in the spectra in all measured fields. The primary mechanism is seen to grow, cause transition and produce large-scale turbulence. The secondary mechanism grows through the entire transition region and produces the small-scale turbulence. Advanced techniques use Linear Stochastic Estimation (LSE) and Proper Orthogonal Decomposition (POD) to identify the spatio-temporal evolutions of structures in the boundary layer. LSE is used to estimate the instantaneous velocity fields using temporal data from just two spatial locations and the spatial correlations. Reference locations are selected using maximum RMS values to provide the best available estimates. POD is used to objectively determine modes characteristic of the measured flow based on energy. The stationary vortices are identified in the first laminar modes of each velocity component and shear component. Experimental evidence suggests that neighboring vortices interact and produce large coherent structures with spanwise periodicity at double the stationary vortex wavelength. An objective transition region detection method is developed using streamwise spatial POD solutions which isolate the growth of the primary and secondary instability mechanisms in the first and second modes, respectively. Temporal evolutions of dominant POD modes in all measured fields are calculated. These scalar POD coefficients contain the integrated characteristics of the entire field, greatly reducing the amount of data to characterize the instantaneous field. These modes may then be used to train future flow control algorithms based on neural networks.

Chapman, Keith↗
NASA NTRS2002

Spatial and temporal traction response in human airway smooth muscle cells

Tractions that cells exert on their substrates are essential in cell spreading, migration, and contraction. These tractions can be determined by plating the cells on a flexible gel and measuring the deformation of the gel by using fluorescent beads embedded just below the surface of the gel. In this article we describe the image correlation method (ICM) optimized for determining the displacement field of the gel under a contracting cell. For the calculation of the traction field from the displacement field we use the recently developed method of Fourier transform traction cytometry (FTTC). The ICM and FTTC methods are applied to human airway smooth muscle cells during stimulation with the contractile agonist histamine or the relaxing agonist isoproterenol. The overall intensity of the cell contraction (the median traction magnitude, the energy transferred from the cell to the gel, and the net contractile moment) increased after activation with histamine, and decreased after treatment with isoproterenol. Cells exhibited regional differences in the time course of traction during the treatment. Both temporal evolution and magnitude of traction increase induced by histamine varied markedly among different cell protrusions, whereas the nuclear region showed the smallest response. These results suggest that intracellular mediators of cell adhesion and contraction respond to contractile stimuli with different rates and intensities in different regions of the cell.

NASA Discipline Cell Biology↗
NASA NTRS1995

Aerodynamic sound of flow past an airfoil

The long term objective of this project is to develop a computational method for predicting the noise of turbulence-airfoil interactions, particularly at the trailing edge. We seek to obtain the energy-containing features of the turbulent boundary layers and the near-wake using Navier-Stokes Simulation (LES or DNS), and then to calculate the far-field acoustic characteristics by means of acoustic analogy theories, using the simulation data as acoustic source functions. Two distinct types of noise can be emitted from airfoil trailing edges. The first, a tonal or narrowband sound caused by vortex shedding, is normally associated with blunt trailing edges, high angles of attack, or laminar flow airfoils. The second source is of broadband nature arising from the aeroacoustic scattering of turbulent eddies by the trailing edge. Due to its importance to airframe noise, rotor and propeller noise, etc., trailing edge noise has been the subject of extensive theoretical (e.g. Crighton & Leppington 1971; Howe 1978) as well as experimental investigations (e.g. Brooks & Hodgson 1981; Blake & Gershfeld 1988). A number of challenges exist concerning acoustic analogy based noise computations. These include the elimination of spurious sound caused by vortices crossing permeable computational boundaries in the wake, the treatment of noncompact source regions, and the accurate description of wave reflection by the solid surface and scattering near the edge. In addition, accurate turbulence statistics in the flow field are required for the evaluation of acoustic source functions. Major efforts to date have been focused on the first two challenges. To this end, a paradigm problem of laminar vortex shedding, generated by a two dimensional, uniform stream past a NACA0012 airfoil, is used to address the relevant numerical issues. Under the low Mach number approximation, the near-field flow quantities are obtained by solving the incompressible Navier-Stokes equations numerically at chord Reynolds number of 104. The far-field noise is computed using Curle's extension to the Lighthill analogy (Curle 1955). An effective method for separating the physical noise source from spurious boundary contributions is developed. This allows an accurate evaluation of the Reynolds stress volume quadrupoles, in addition to the more readily computable surface dipoles due to the unsteady lift and drag. The effect of noncompact source distribution on the far-field sound is assessed using an efficient integration scheme for the Curle integral, with full account of retarded-time variations. The numerical results confirm in quantitative terms that the far-field sound is dominated by the surface pressure dipoles at low Mach number. The techniques developed are applicable to a wide range of flows, including jets and mixing layers, where the Reynolds stress quadrupoles play a prominent or even dominant role in the overall sound generation.

Wang, Meng↗
NASA NTRS2004

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗
NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗