Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “beam propagation method”

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 127 records · Page 7

Effects of various cavity designs on the performance of a CO2 TEA laser with an unstable resonator

Unstable resonator modeling has been carried out for an injection-seeded CO2 transversely excited atmosphere (TEA) laser in the NOAA/ERL/Wave Propagation Laboratory (WPL) Doppler lidar to examine the effects of various cavity designs on the quality of the output beam. The results show the effects of an injection pinhole, electrode spacing, mirror tilt, and radial reflectivity function of the output coupler. The electrode spacing in this laser has negligible effect. The injection pinhole, however, produces complicated structures in the output patterns. If the pinhole is removed, the output pattern is much smoother, and the frequency jitter is smaller. Misalignment sensitivity is very closely related to the radial reflectivity function. The superparabolic function provides the highest coupling efficiency, largest beam size, and good collimation, but produces a slightly higher misalignment sensitivity compared with a parabolic function. The Gaussian function provides the lowest misalignment sensitivity, but it produces the smallest beam size and the largest beam divergence. Also, the coupling coefficient is 50 percent lower than the optimum value. Methods for using a flat diffraction grating in unstable resonators are also investigated. The best way is to use a flat grating/positive lens combination to replace the back concave mirror.

Zhao, Yanzeng↗

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

Electromagnetic Properties of Filled, Pixellated Arrays

Filled arrays of bolometers are currently being employed for use in astronomy from the far-infrared through millimeter parts of the electromagnetic spectrum. Because of the large range of wavelengths for which this is applicable, the number of modes supported by a bolometric pixel will vary according to specific application. Withington et al. (2003) have built a formalism for treating the electromagnetic properties of such bolometers by propagating the second order statistical properties of the radiation through a canonical optical system. In this work, we use this formalism to construct beam pattern images of square pixels for various ratios of p/lambda where p is the pixel size and lambda is the wavelength. In the low mode limit, the diffraction effects cause the beam pattern to be circular with a quadrupole dependence of Stokes Q and U. High mode cases approach the geometric limit. The polarization in these cases can be seen to trace the pixel edges. The effective size of the beam has a direct impact on the inter-pixel coupling and sets the number of independent detectors in an astronomical focal plane. This technique illustrates and quantifies the relationship between pixel size and angular resolution limits for a given wavelength and telescope. This is especially true in the limit of low p/lambda. In this case, the diffraction due to pixelization is non-negligible for the calculation of the overall angular resolution of the telescope. In addition, for instruments that are polarization sensitive, this method also provides a quantitative method for determining the contribution of the instrument to the measured polarization.

Chuss, David T.↗

Endothelium Preserving Microwave Treatment for Atherosclerosis

Method and apparatus are provided to treat atherosclerosis wherein the artery is partially closed by dilating the artery while preserving the vital and sensitive endothelial layer thereof Microwave energy having a frequency from 3 GHz to 300 GHz is propagated into the arterial wall to produce a desired temperature profile therein at tissue depths sufficient for thermally necrosing connective tissue and softening fatty and waxy plaque while limiting heating of surrounding tissues including the endothelial laser and/or other healthy tissue, organs, and blood. The heating period for raising the temperature a potentially desired amount, about 20 C., within the atherosclerotic lesion may be less than about one second. In one embodiment of the invention, a radically beveled waveguide antenna is used to deliver microwave energy at frequencies from 25 GHz or 30 GHz to about 300 GHz and is focused towards a particular radial sector of the artery. Because the atherosclerotic lesions are often asymmetrically disposed, directable of focussed heating preserves healthy sectors of the artery and applies energy to the asymmetrically positioned lesion faster than a non-directed beam. A computer simulation predicts isothermic temperature profiles for the given conditions and man be used in selecting power, pulse duration, beam width, and frequency of operation to maximize energy deposition and control heat rise within the atherosclerotic lesion without harming healthy tissues or the sensitive endothelium cells.

Carl, James R.↗

Endothelium Preserving Microwave Treatment for Atherosclerois

Method and apparatus are provided to treat atherosclerosis wherein the artery is partially closed by dilating the artery while preserving the vital and sensitive endothelial layer thereof. Microwave energy having a frequency from 3 GHz to 300 GHz is propagated into the arterial wall to produce a desired temperature profile therein at tissue depths sufficient for thermally necrosing connective tissue and softening fatty and waxy plaque while limiting heating of surrounding tissues including the endothelial layer and/or other healthy tissue, organs, and blood. The heating period for raising the temperature a potentially desired amount about 20 C, with the atherosclerotic lesion may be less than about one second. In one embodiment of the invention, a radically beveled waveguide antenna is used to deliver microwave energy at frequencies from 25 GHz or 30 GHz to about 300 GHz and is focused towards a particular radical sector of the artery. Because the atherosclerotic lesions are often asymmetrically disposed, directable or focussed heating preserves healthy sectors of the artery and applies energy to the asymmetrically positioned lesion faster than a non-directed beam. A computer simulation predicts isothermic temperature profiles for the given conditions and may be used in selecting power, pulse duration, beam width, and frequency of operation to maximize energy deposition and control heat rise within the atherosclerotic lesion without harming healthy tissues or the sensitive endothelium cells.

Carl, James R.↗

Endothelium Preserving Microwave Treatment for Atherosclerosis

Method and apparatus are provided to treat atherosclerosis wherein the artery is partially closed by dilating the artery while preserving the vital and sensitive endothelial layer thereof. Microwave energy having a frequency from 3 GHz to 300 GHz is propagated into the arterial wall to produce a desired temperature profile therein at tissue depths sufficient for thermally necrosing connective tissue and softening fatty and waxy plaque while limiting heating of surrounding tissues including the endothelial layer and/or other healthy tissue, organs, and blood. The heating period for raising the temperature a potentially desired amount, about 20 C. within the atherosclerotic lesion may be less than about one second. In one embodiment of the invention, a radically beveled waveguide antenna is used to deliver microwave energy at frequencies from 25 GHz or 30 GHz to about 300 GHz and is focused towards a particular radial sector of the artery. Because the atherosclerotic lesions are often asymmetrically disposed directable or focussed heating preserves healthy sectors of the artery and applies energy to the asymmetrically positioned lesion faster than a non-directed beam. A computer simulation predicts isothermic temperature profiles for the given conditions and may be used in selecting power, pulse duration, beam width, and frequency of operation to maximize energy deposition and control heat rise within the atherosclerotic lesion without harming healthy tissues or the sensitive endothelium cells.

Carl, James R.↗

Unusually Long Path Length for a Nearly Scatter-Free Solar Particle Event Observed By Solar Orbiter at 0.43 au

Context. After their acceleration and release at the Sun, solar energetic particles (SEPs) are injected into the interplanetary medium and are bound to the interplanetary magnetic field (IMF) by the Lorentz force. The expansion of the IMF close to the Sun focuses the particle pitch-angle distribution, and scattering counteracts this focusing. Solar Orbiter observed an unusual solar particle event on 9 April 2022 when it was at 0.43 astronomical units (au) from the Sun. Aims. We show that the inferred IMF along which the SEPs traveled was about three times longer than the nominal length of the Parker spiral and provide an explanation for this apparently long path. Methods. We used velocity dispersion analysis (VDA) information to infer the spiral length along which the electrons and ions traveled and infer their solar release times and arrival direction. Results. The path length inferred from VDA is approximately three times longer than the nominal Parker spiral. Nevertheless, the pitch-angle distribution of the particles of this event is highly anisotropic, and the electrons and ions appear to be streaming along the same IMF structures. The angular width of the streaming population is estimated to be approximately 30 degrees. The highly anisotropic ion beam was observed for more than 12 h. This may be due to the low level of fluctuations in the IMF, which in turn is very probably due to this event being inside an interplanetary coronal mass ejection The slow and small rotation in the IMF suggests a flux-rope structure. Small flux dropouts are associated with very small changes in pitch angle, which may be explained by different flux tubes connecting to different locations in the flare region. Conclusions. The unusually long path length along which the electrons and ions have propagated virtually scatter-free together with the short-term flux dropouts offer excellent opportunities to study the transport of SEPs within interplanetary structures. The 9 April 2022 solar particle event offers an especially rich number of unique observations that can be used to limit SEP transport models.

Robert F. Wimmer-Schweingruber↗

Endothelium Preserving Microwave Treatment for Atherosclerosis

Method and apparatus are provided to treat atherosclerosis wherein the artery is partially closed by dilating the artery while preserving the vital and sensitive endothelial layer thereof. Microwave energy having a frequency from 3 GHz to 300 GHz is propagated into the arterial wall to produce a desired temperature profile therein at tissue depths sufficient for thermally necrosing connective tissue and softening fatty and waxy plaque while limiting heating of surrounding tissues including the endothelial layer and/or other healthy tissue, organs, and blood. The heating period for raising the temperature a potentially desired amount, about 20 C. within the atherosclerotic lesion may be less than about one second. In one embodiment of the invention, a radically beveled waveguide antenna is used to deliver microwave energy at frequencies from 25 GHz or 30 GHz to about 300 GHz and is focused towards a particular radial sector of the artery. Because the atherosclerotic lesions are often asymmetrically disposed, directable or focussed heating preserves healthy sectors or the artery and applies energy to the asymmetrically positioned lesion faster than a non-directed bean. A computer simulation predicts isothermic temperature profiles for the given conditions and may be used in selecting power, pulse duration, beam width, and frequency of operation to maximize energy deposition and control heat rise within the atherosclerotic lesion without harming healthy tissues or the sensitive endothelium cells.

James R Carl↗

Recent Trends in the Analysis of Quasioptical Systems

The recent trend in microwave instruments is the use of multiple millimeter and submillimeter wavelength bands. These systems are typically analyzed by using physical optics, Gaussian beams or ray tracing techniques. Physical optics offers high accuracy at the expense of computation time. This trade-off becomes particularly apparent in the analysis of multiple reflector antennas, such as beam waveguide antennas, where physical optics is used to compute the current on each reflector from the current on the previous reflector. At the other end of the spectrum is ray tracing approaches that ignore diffraction effects entirely. These methods are fast but sacrifice the ability to predict some effects accurately. An intermediate approach is to use an appropriate set of expansion functions to model the field between the reflectors. If the set is chosen wisely only a few coefficients need to be determined from each reflector current. The field is then computed at the next reflector through the use of the expansion functions and their coefficients rather than by using the previous reflector current. For a beam waveguide system with no enclosing tubes an excellent set of expansion functions is the Gaussian beam mode set. In many cases a preliminary design which includes the effects on diffraction may be obtained by considering only the fundamental mode and a thin lens model for the reflectors. Higher-order modes are included to model the effects of the curved reflector, which include asymmetric distortion of the beam, cross polarization, and beam truncation. This paper describes a computer code implementing higher-order Gaussian beam scattering by multiple reflector systems. There are four essential steps in the algorithm. (1) Compute the current on the first reflector using physical optics using either a feed model or by an incident set of Gaussian beam modes. (2) Find the direction of propagation for the reflected Gaussian beam-set using ray tracing. (3) Determine the waist size and location for the output beam set by examining the amplitude and phase distribution of the current on the reflector. (4) Compute the amplitudes of the individual modes in the output mode set. These steps are then repeated for each addition reflector in the chain. In each of these cases the previous Gaussian beam set provides the input field for the current calculation. Details of the four steps discussed above will be discussed. Examples will compare results from the Gaussian beam approach to pure physical optics, illustrating both its merits and limitations. Hybrid approaches capable of eliminating some of the limitations will also be discussed.

Hoppe, D. J.↗

Simulated Service and Stress Corrosion Cracking Testing for Friction Stir Welded Spun Form Domes

Damage tolerance testing development was required to help qualify a new spin forming dome fabrication process for the Ares 1 program at Marshall Space Flight Center (MSFC). One challenge of the testing was due to the compound curvature of the dome. The testing was developed on a sub-scale dome with a diameter of approximately 40 inches. The simulated service testing performed was based on the EQTP1102 Rev L 2195 Aluminum Lot Acceptance Simulated Service Test and Analysis Procedure generated by Lockheed Martin for the Space Shuttle External Fuel Tank. This testing is performed on a specimen with an induced flaw of elliptical shape generated by Electrical Discharge Machining (EDM) and subsequent fatigue cycling for crack propagation to a predetermined length and depth. The specimen is then loaded in tension at a constant rate of displacement at room temperature until fracture occurs while recording load and strain. An identical specimen with a similar flaw is then proof tested at room temperature to imminent failure based on the critical offset strain achieved by the previous fracture test. If the specimen survives the proof, it is then subjected to cryogenic cycling with loads that are a percentage of the proof load performed at room temperature. If all cryogenic cycles are successful, the specimen is loaded in tension to failure at the end of the test. This standard was generated for flat plate, so a method of translating this to a specimen of compound curvature was required. This was accomplished by fabricating a fixture that maintained the curvature of the specimen rigidly with the exception of approximately one-half inch in the center of the specimen containing the induced flaw. This in conjunction with placing the center of the specimen in the center of the load train allowed for successful testing with a minimal amount of bending introduced into the system. Stress corrosion cracking (SCC) tests were performed using the typical double beam assembly and with 4-point loaded specimens under alternate immersion conditions in a 3.5% NaCl environment for 90 days. In addition, experiments were conducted to determine the threshold stress intensity factor for SCC (K1SCC) of Al-Li 2195 which to our knowledge has not been determined previously. The successful simulated service and stress corrosion testing helped to provide confidence to continue to Ares 1 scale dome fabrication.

Stewart, Thomas J.↗

First Year Wilkinson Microwave Anisotropy Probe(WMAP) Observations: Data Processing Methods and Systematic Errors Limits

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of l/f noise. We infer and correct for a small (approx. 1 %) transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

Hinshaw, G.↗

Aligning Three Off-Axis Mirrors with Help of a DOE

A proposed method based on the use of a special-purpose diffractive optical element (DOE) would simplify (relative to prior methods) the alignment of three off-axis mirrors that constitute an imaging optical system. The method would exploit the fact that a DOE can be fabricated lithographically with high accuracy by electron-beam lithography in a thin film of poly(methyl methacrylate). The method would effectively transfer much of the problem of obtaining the needed accuracy from the mechanical-mirror-alignment domain to the lithographic domain. Unlike other methods that depend on specific symmetries (e.g., sphericity and/or concentricity), this method is expected to apply with equal ease and accuracy to mirrors of any configuration including aspherical, decentered mirrors. Assuming that one of the mirrors of a general three-mirror imaging optical system can serve as a reference for the alignment of the other two mirrors, such a system has 12 degrees of freedom in alignment. In the proposed method, one would use an interferometer in combination with a DOE to effect precise and relatively rapid and easy alignment of two of the mirrors with respect to each other, thus reducing the alignment task to that of the six degrees of freedom of the remaining mirror. The figure depicts a representative three-mirror off-axis imaging system, wherein the primary and tertiary mirrors (M1 and M3, respectively) are concave and the secondary mirror (M2) is convex. The DOE for aligning this system would be fabricated on the right surface of an optical flat and could be made to have either negative or positive focusing power, depending on the requirements of the specific application. The DOE could be designed to be placed at any convenient distance from M1 and M3 - again, depending on the application. The DOE would be illuminated with light coming from the left, generated by an interferometer. First, assuming the optical flat is of high quality, the plane of the DOE would be aligned perpendicular to the collimated beam by use of light reflected from the left face of the optical flat. The DOE would comprise two independent areas: one dedicated to M1, the other to M3. The portions of the collimated beam passing through those areas would be diffracted towards the corresponding mirrors. A mask, not shown in the figure, could be used to prevent light from passing through the rest of the area of the optical flat. Light rays reflected from M1 and M3 would retrace their paths through the DOE and would propagate leftward to the interferometer.

Mouroulis, Pantazis↗

Remote Sensing of Hail from Space: Retrievals, Climatologies, and Challenges Going Forward

In addition to the myriad threats that severe hailstorms pose to society, infrastructure and agriculture, severe hail is difficult to measure in situ, and surface-based hail reporting and detection methods are inconsistent and subject to geographical or societal biases. This motivates the use of spaceborne remote-sensing platforms to retrieve hail and construct climatologies in the most globally consistent way. Passive-microwave algorithms leverage the sensitivity of spaceborne passive-microwave radiometers to scattering by hail, particularly in the channels from 10 to 89 GHz. These retrievals are used to construct global climatologies of severe hail. The Bang and Cecil (2019) retrievals has been applied to several different spaceborne sensors: the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), Global Precipitation Measurement Mission (GPM) Microwave Imager, Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and Advanced Microwave Scanning Radiometer 2 (AMSR2) sensors and used to construct [near] global climatologies of severe hail. This retrieval, and others, have been tested using Global Precipitation Measurement (GPM) Ku-band precipitation radar, to assess their effectiveness and regional variability. A successful retrieval and climatology are those that correspond tightly to radar reflectivity and give minimal appearance of regional biases, especially with latitude. Satellite platforms offer consistent observations, even in remote, data-sparse, and oceanic regions that ground-based networks exclude. There are, however, potential disconnects between the processes identified aloft by the satellite and the resultant weather at the ground, leading to uncertainties in the retrievals that may propagate into satellite-based climatologies, particularly in the Tropics, where there are abundant strong - but not necessarily hailing - storms that are strongly represented in the current satellite climatologies. There are ongoing efforts to assess and mitigate the contributing factors to these uncertainties, chiefly among them the effects of non-uniform beam filling in the passive-microwave footprint, and the relationships between the size distributions of hailstones aloft and the dynamic processes and environments with which they interact throughout their trajectories.

Sarah D Bang↗

Spaceborne Remote Sensing of Hail: Retrievals, Climatologies, and Challenges Going Forward

In addition to the myriad threats that severe hailstorms pose to society, infrastructure and agriculture, severe hail is difficult to measure in situ, and surface-based hail reporting and detection methods are inconsistent and subject to geographical or societal biases. This motivates the use of spaceborne remote-sensing platforms to retrieve hail and construct climatologies in the most globally consistent way. Passive-microwave algorithms leverage the sensitivity of spaceborne passive-microwave radiometers to scattering by hail, particularly in the channels from 10 to 89 GHz. These retrievals are used to construct global climatologies of severe hail. The Bang and Cecil (2019) retrievals has been applied to several different spaceborne sensors: the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), Global Precipitation Measurement Mission (GPM) Microwave Imager, Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and Advanced Microwave Scanning Radiometer 2 (AMSR2) sensors and used to construct [near] global climatologies of severe hail. This retrieval, and others, have been tested using Global Precipitation Measurement (GPM) Ku-band precipitation radar, to assess their effectiveness and regional variability. A successful retrieval and climatology are those that correspond tightly to radar reflectivity and give minimal appearance of regional biases, especially with latitude. Satellite platforms offer consistent observations, even in remote, data-sparse, and oceanic regions that ground-based networks exclude. There are, however, potential disconnects between the processes identified aloft by the satellite and the resultant weather at the ground, leading to uncertainties in the retrievals that may propagate into satellite-based climatologies, particularly in the Tropics, where there are abundant strong - but not necessarily hailing - storms that are strongly represented in the current satellite climatologies. There are ongoing efforts to assess and mitigate the contributing factors to these uncertainties, chiefly among them the effects of non-uniform beam filling in the passive-microwave footprint, and the relationships between the size distributions of hailstones aloft and the dynamic processes and environments with which they interact throughout their trajectories.

Sarah D Bang↗

Detecting Hail from Space: Algorithms, Climatologies, and Challenges Going Forward

In addition to the myriad threats that severe hailstorms pose to society, infrastructure and agriculture, severe hail is difficult to measure in situ, and surface-based hail reporting and detection methods are inconsistent and subject to geographical or societal biases. This motivates the use of spaceborne remote-sensing platforms to retrieve hail and construct climatologies in a globally uniform way. We have developed a hail detection algorithm that leverages the sensitivity of spaceborne passive-microwave radiometers to scattering by hail, particularly in the channels from 10 to 89 GHz. We use this retrieval to construct global climatologies of severe hail using several different spaceborne sensors: the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), Global Precipitation Measurement Mission (GPM) Microwave Imager, Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and Advanced Microwave Scanning Radiometer 2 (AMSR2) sensors and are working to extend into the late 1980’s using the Special Sensor Microwave Imager/Sounder (SSMI/(S)) data. Using coincident Global Precipitation Measurement (GPM) Ku-band precipitation radar, we assessed this retrieval and several others in the literature for their effectiveness and regional variability. We developed a passive-microwave algorithm that corresponds tightly to radar reflectivity and gives the least appearance of regional biases compared to other passive-microwave approaches in the literature. Satellite platforms offer consistent observations, even in remote, data-sparse, and oceanic regions that ground-based networks exclude. There are, however, potential disconnects between the processes identified aloft by the satellite and the resultant weather at the ground, leading to uncertainties in the retrievals that may propagate into satellite-based climatologies, particularly in the Tropics, where there are abundant strong - but not necessarily hailing - storms that are strongly represented in the current satellite climatologies. We will discuss ongoing efforts to assess and mitigate the contributing factors to these uncertainties, chiefly among them the effects of non-uniform beam filling in the passive-microwave footprint, and the relationships between the size distributions of hailstones aloft and the dynamic processes and environments with which they interact throughout their trajectories.

Sarah D Bang↗

Spaceborne Passive-Microwave Hail Detection: Global Climatologies, Validation, and Challenges Going Forward

In addition to the myriad threats that severe hailstorms pose to society, infrastructure and agriculture, severe hail is difficult to measure in situ, and surface-based hail reporting and detection methods are inconsistent and vulnerable to nonmeteorological biases. This motivates the use of spaceborne remote-sensing platforms to retrieve hail and construct climatologies in a globally uniform way. We have developed a hail detection algorithm that leverages the sensitivity of spaceborne passive-microwave radiometers to scattering by hail, particularly in the channels from 10 to 89 GHz. We use this retrieval to construct global passive-microwave climatologies of severe hail using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), Global Precipitation Measurement Mission (GPM) Microwave Imager, Advanced Microwave Scanning Radiometer for EOS (AMSR-E), and Advanced Microwave Scanning Radiometer 2 (AMSR2) sensors and are working to extend into the pre-TRMM era to the Special Sensor Microwave Imager/Sounder (SSMI/(S)) data. Using coincident Global Precipitation Measurement (GPM) Ku-band precipitation radar, we assessed this retrieval and several others in the literature for their effectiveness and regional variability. We find that this retrieval, which leverages a signature in the Minimum 19-GHz polarization corrected temperature (PCT) combined with the 37-GHz PCT depression normalized by tropopause height constrains the radar reflectivity most tightly, and gives the least appearance of regional biases compared to other passive-microwave approaches in the literature. Satellite platforms offer consistent observations, even in remote, data-sparse, and oceanic regions that ground-based networks exclude. There are, however, potential disconnects between the processes identified aloft by the satellite and the resultant weather at the ground, leading to uncertainties in the retrievals that may propagate into satellite-based climatologies, particularly in the Tropics, where there are abundant strong - but not necessarily hailing - storms that are strongly represented in the current satellite climatologies. We will discuss ongoing efforts to assess and mitigate the contributing factors to these uncertainties, chiefly among them the effects of non-uniform beam filling in the passive-microwave footprint, and the relationships between the size distributions of hailstones aloft and the dynamic processes and environments with which they interact throughout their trajectories.

Sarah D. Bang↗

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗