Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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

Characterization of Epoxy Functionalized Graphite Nanoparticles and the Physical Properties of Epoxy Matrix Nanocomposites

This work presents a novel approach to the functionalization of graphite nanoparticles. The technique provides a mechanism for covalent bonding between the filler and matrix, with minimal disruption to the sp2 hybridization of the pristine graphene sheet. Functionalization proceeded by covalently bonding an epoxy monomer to the surface of expanded graphite, via a coupling agent, such that the epoxy concentration was measured as approximately 4 wt.%. The impact of dispersing this material into an epoxy resin was evaluated with respect to the mechanical properties and electrical conductivity of the graphite-epoxy nanocomposite. At a loading as low as 0.5 wt.%, the electrical conductivity was increased by five orders of magnitude relative to the base resin. The material yield strength was increased by 30% and Young s modulus by 50%. These results were realized without compromise to the resin toughness.

Miller, Sandi G.↗
NASA NTRS2011

Drop Size Distribution Measurements Supporting the NASA Global Precipitation Measurement Mission: Infrastructure and Preliminary Results

Global Precipitation Measurement Mission (GPM) retrieval algorithm validation requires datasets that characterize the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (5 -- 50 km). Key to this process is the combined use of disdrometer and polarimetric radar platforms. Here the disdrometer measurements serve as a reference for up-scaling dual-polarimetric radar observations of the PSD to the much larger volumetric sampling domain of the radar. The PSD observations thus derived provide a much larger data set for assessing DSD variability, and satellite-based precipitation retrieval algorithm assumptions, in all three spatial dimensions for a range of storm types and seasons. As one component of this effort, the GPM Ground Validation program recently acquired five 3rd generation 2D Video disdrometers as part of its Disdrometer and Radar Observations of Precipitation Facility (DROP), currently hosted in northern Alabama by the NASA Marshall Space Flight Center and the University of Alabama in Huntsville. These next-generation 2DVDs were operated and evaluated in different phases of data collection under the scanning domain of the UAH ARMOR C-band dual-polarimetric radar. During this period approximately 7500 minutes of PSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSDs were then used as input to T-matrix code for computation of polarimetric radar moments at C-band. The combined dataset was then analyzed with two basic objectives in mind: 1) the investigation of seasonal variability in the rain PSD parameters as observed by the 2DVDs; 2) the use of combined polarimetric moments and observed gamma distribution parameters in a functional form to retrieve PSD parameters in 4-D using the ARMOR radar for precipitation occurring in different seasons and for different rain system types. Preliminary results suggest that seasonal variations in the DSD parameters do occur, but are most pronounced when comparing tropical PSDs to either winter or summer convective precipitation. For example the previously documented shift to relatively smaller drop diameters in higher number concentrations for equivalent rain rate bins was observed in tropical storm rainbands occurring over Huntsville. On a more inter seasonal basis empirical fits between parameters such as D0 and ZDR do not appear to exhibit robust seasonal biases- i.e., one fit seems to work for all seasons within acceptable standard error (O[10%]) for estimates of D0. In polarimetric retrievals of the vertical variability in PSD (rain layer) for a tropical rainband we find that the Do varies with height when partitioned by specified precipitation categories (e.g., convective or stratiform, heavy and light stratiform etc.) but this variation is of order 10-20% and is smaller than the difference in D0 observed between the basic delineation of convective and stratiform precipitation types. Currently we are expanding our analysis of the vertical structure of the PSD to include several seasonally and/or dynamically-different storm system types (e.g., winter convection and stratiform events; summer mid-latitude convective etc.) sampled by ARMOR. The study will present the results of our combined analyses.

Petersen, Walter A.↗
NASA NTRS2012

Investigation of Effects of Material Architecture on the Elastic Response of a Woven Ceramic Matrix Composite

To develop methods for quantifying the effects of the microstructural variations of woven ceramic matrix composites on the effective properties and response of the material, a research program has been undertaken which is described in this paper. In order to characterize and quantify the variations in the microstructure of a five harness satin weave, CVI SiC/SiC, composite material, specimens were serially sectioned and polished to capture images that detailed the fiber tows, matrix, and porosity. Open source quantitative image analysis tools were then used to isolate the constituents and collect relevant statistics such as within ply tow spacing. This information was then used to build two dimensional finite element models that approximated the observed section geometry. With the aid of geometrical models generated by the microstructural characterization process, finite element models were generated and analyses were performed to quantify the effects of the microstructure and its variation on the effective stiffness and areas of stress concentration of the material. The results indicated that the geometry and distribution of the porosity appear to have significant effects on the through-thickness modulus. Similarly, stress concentrations on the outer surface of the composite appear to correlate to regions where the transverse tows are separated by a critical amount.

Goldberg, Robert K.↗
NASA NTRS1977

Estimation of the geophysical properties of the ocean surface using aircraft microwave measurements

An improved model of the effects of sea state on microwave signature has been developed which incorporates the different effects of whitecaps and streaks to define the response of microwave channels to wind speed. This model has been demonstrated to agree with recent measurements. An approximation model has also been incorporated to describe the effects of precipitation on microwave radiation through a computationally rapid routine. The use of these models and a new technique to allow the selection of the most climatologically appropriate D-matrix is demonstrated in the inversion of data collected over the bering Sea. Surface wind speed agrees very well with observations while good results are obtained for integrated water vapor, and liquid water.

Fowler, M. G.↗
NASA NTRS1984

Analysis of hourglass instabilities and control in underintegrated finite element methods

Belytschko et al. (1981, 1984) has developed stabilization methods for the treatment of underintegrated FEM problems; these methods involve the computation of an underintegrated stiffness matrix, which is rank-deficient, and the addition of a stabilization matrix which effectively eliminates the spurious modes. An attempt is presently made to give this a priori stabilization method a mathematical means of support. Attention is also given to an a posteriority stabilization method for hourglass control, in which an approximate solution of the underintegrated system is obtained and then subjected to a special projection in order to eliminate the hourglass modes. A proof is obtained for the convergence of this stabilized underintegrated approximation to the exact solution of a model problem at almost the same rate (as the mesh is refined) as the fully integrated solutions.

Jacquotte, O.-P.↗
NASA NTRS1991

Measurements and mechanisms of localized aqueous corrosion in aluminum-lithium-copper alloys

Summary information is included for electrochemical aspects of stress corrosion cracking in alloy 2090 and an introduction to the work to be initiated on the new X2095 (Weldalite) alloy system. Stress corrosion cracking (SCC) was studied in both S-T and L-T orientations in alloy 2090. A constant load TTF test was performed in several environments with a potentiostatically applied potential. In the same environments the electrochemical behavior of phases found along subgrain boundaries was assessed. It was found that rapid failure due to SCC occurred when the following criteria was met: E(sub BR,T1) is less than E(sub applied) is less than E(sub Br, matrix phase). Although the L-T orientation is usually considered more resistant to SCC, failures in this orientation occurred when the stated criteria was met. This may be due to the relatively isotropic geometry of the subgrains which measure approximately 12 to 25 microns in diameters. Initial studies of alloy X2095 includes electrochemical characterization of three compositional variations each at three temperatures. The role of T(sub 1) dissolution in SCC behavior is addressed using techniques similar to those used in the research of 2090 described. SCC susceptibility is also studied using alternate immersion facilities at Reynolds Metals Corporation. Pitting is investigated in terms of stability, role of precipitate phases and constituent particles, and as initiation sites for SCC. In all research endeavors, attempts are made to link electrochemistry to microstructure. Previous work on 2090 provides a convenient basis for comparison since both alloys contain T(sub 1) precipitates but with different distributions. In 2090 T(sub 1) forms preferentially on subgrain boundaries whereas in X2095 the microstructure appears to be more homogeneous with finer T(sub 1) particles. Another point for comparison is the delta prime strengthening phase found in 2090 but absent in X2095.

Wall, Douglas↗
NASA NTRS1990

Real-time simulation of supersonic inlets

A previously published real-time simulation algorithm, the matrix stability region placement (MSRP) method, is used to simulate a small perturbation model of the NASA Lewis Mach 2.5 40-60 mixed compression inlet. The model is representative of high-speed internal flow propulsion systems which can be approximated as quasi-one-dimensional flows. The resulting system of equations, which is stiff, is also simulated by the second-order Adam-Bashforth method. It is shown that the MSRP method can be used to simulate small perturbation models of high-speed internal flow propulsion systems in real time.

Mossayebi, F.↗
NASA NTRS2002

Simulation of Energy Response of the ATIC Calorimeter

ATIC (Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CR particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2001. The ATIC flight collected approximately 25 million events. For reconstruction of primary spectra from spectra of energy deposits measured in the experiment, correlations between kinetic energy of a primary particle E(sub kin) and energy deposit in the calorimeter E(sub d) should be known. For this purpose, simulations of energy response of the calorimeter on energy spectra of different nuclei were done. The simulations were performed by GEANT-3.21 code with QGSM generator for nucleus - nucleus interactions. The incident flux was taken as isotropic in the ATIC aperture. Primary spectra power-law by momentum were used as inputs according to standard models of cosmic ray acceleration. These spectra become power-law by kinetic energy at E(sub kin) higher than approx.20Mc(sup 2), where M is primary nucleus mass. It should be noted that energy deposit spectra measured by ATIC illustrate similar behavior. Distributions of ratio E(sub kin)/E(sub d) are presented for different energy deposits and for a set of primaries. For power-law regions of energy spectra at E(sub d)> or equal to 20Mc(sup 2) the obtained mean value of E(sub kin)/E(sub d) increases from approx.2.4 for protons to approx.3.1 for iron, while rms/ decreases from 50% for protons to about 15% for iron. These values were obtained for the spectral index gamma=1.6

Batkov, K. E.↗
NASA NTRS2004

The Origin and Impact History of Lunar Meteorite Yamato 86032

Yamato (Y) 86032 is a feldspathic lunar highland breccia having some characteristics of regolith breccia. The absence of KREEP components in the matrix in Y86032 indicates that these meteorites came from a long distance from Mare Imbrium, perhaps from the far-side of the moon. One ferroan anorthosite (FAN) clast in Y86032 has a very old Ar-Ar age of approximately 4.35-4.4 Ga. The negative Nd of this clast may suggest a direct link with the primordial magma ocean. The facts indicate that Y86032 contains components derived from a protolith of the original lunar crust. Detailed petrologic characterization of each component in this breccia is essential to understand the early impact history and origin of the lunar highland crust. We made a large slab (5.2 x 3.6 cm x 3-5 mm) of Y86032 to better understand the relationship of various lithologies and their petrologic origin.

Yamaguchi, A.↗
NASA NTRS2015

An Upper Bound on Orbital Debris Collision Probability When Only One Object has Position Uncertainty Information

Upper bounds on high speed satellite collision probability, P (sub c), have been investigated. Previous methods assume an individual position error covariance matrix is available for each object. The two matrices being combined into a single, relative position error covariance matrix. Components of the combined error covariance are then varied to obtain a maximum P (sub c). If error covariance information for only one of the two objects was available, either some default shape has been used or nothing could be done. An alternative is presented that uses the known covariance information along with a critical value of the missing covariance to obtain an approximate but useful P (sub c) upper bound. There are various avenues along which an upper bound on the high speed satellite collision probability has been pursued. Typically, for the collision plane representation of the high speed collision probability problem, the predicted miss position in the collision plane is assumed fixed. Then the shape (aspect ratio of ellipse), the size (scaling of standard deviations) or the orientation (rotation of ellipse principal axes) of the combined position error ellipse is varied to obtain a maximum P (sub c). Regardless as to the exact details of the approach, previously presented methods all assume that an individual position error covariance matrix is available for each object and the two are combined into a single, relative position error covariance matrix. This combined position error covariance matrix is then modified according to the chosen scheme to arrive at a maximum P (sub c). But what if error covariance information for one of the two objects is not available? When error covariance information for one of the objects is not available the analyst has commonly defaulted to the situation in which only the relative miss position and velocity are known without any corresponding state error covariance information. The various usual methods of finding a maximum P (sub c) do no good because the analyst defaults to no knowledge of the combined, relative position error covariance matrix. It is reasonable to think, given an assumption of no covariance information, an analyst might still attempt to determine the error covariance matrix that results in an upper bound on the P (sub c). Without some guidance on limits to the shape, size and orientation of the unknown covariance matrix, the limiting case is a degenerate ellipse lying along the relative miss vector in the collision plane. Unless the miss position is exceptionally large or the at-risk object is exceptionally small, this method results in a maximum P (sub c) too large to be of practical use. For example, assuming that the miss distance is equal to the current ISS alert volume along-track (+ or -) distance of 25 kilometers and that the at-risk area has a 70 meter radius. The maximum (degenerate ellipse) P (sub c) is about 0.00136. At 40 kilometers, the maximum P (sub c) would be 0.00085 which is still almost an order of magnitude larger than the ISS maneuver threshold of 0.0001. In fact, a miss distance of almost 340 kilometers is necessary to reduce the maximum P (sub c) associated with this degenerate ellipse to the ISS maneuver threshold value. Such a result is frequently of no practical value to the analyst. Some improvement may be made with respect to this problem by realizing that while the position error covariance matrix of one of the objects (usually the debris object) may not be known the position error covariance matrix of the other object (usually the asset) is almost always available. Making use of the position error covariance information for the one object provides an improvement in finding a maximum P (sub c) which, in some cases, may offer real utility. The equations to be used are presented and their use discussed.

Frisbee, Joseph H., Jr.↗
NASA NTRS1980

Collisional redistribution of radiation. I - The density matrix

As a first step in deriving accurate, calculable expressions for collisional redistribution of radiation, valid outside the impact regime a method for calculating an accurate binary-collision density operator for an atom in the presence of a driving field and perturbers is presented. Projection-operator techniques to establish first the precise validity of the Markoff approximation for radiative relaxation are used. Using the same techniques, it is shown how the collisional-relaxation problem may be analyzed in a practical manner outside the Markoff approximation. The effect of correlations between radiative and collisional events is included in a consistent way and the physical implications of these correlation effects are demonstrated for a simple example, the two-level atom with nondegenerate levels.

Burnett, K.↗
NASA NTRS2000

Radioprotective Thiolamines WR-1065 and WR-33278 Selectively Denature Nonhistone Nuclear Proteins

Differential scanning calorimetry was used to study the interactions of nuclei isolated from Chinese hamster V79 cells with the radioprotector WR-1065, other thiol compounds, and polyamines. Differential scanning calorimetry monitors denaturation of macromolecules and resolves the major nuclear components (e.g. constrained and relaxed DNA, nucleosome core, and nuclear matrix) of intact nuclei on the basis of thermal stability. WR-1065 treatment (0.5-10 mM) of isolated nuclei led to the irreversible denaturation of nuclear proteins, a fraction of which are nuclear matrix proteins. Denaturation of 50% of the total nonhistone nuclear protein content of isolated nuclei occurred after exposure to 4.7 mM WR-1065 for 20 min at 23 C. In addition, a 22% increase in the insoluble protein content of nuclei isolated from V79 cells that had been treated with 4 mM WR-1065 for 30 min at 37 C was observed, indicating that WR-1065-induced protein denaturation occurs not only in isolated nuclei but also in the nuclei of intact cells. From the extent of the increase in insoluble protein in the nucleus, protein denaturation by WR-1065 is expected to contribute to drug toxicity at concentrations greater than approximately 4 mM. WR-33278, the disulfide form of WR1065, was approximately twice as effective as the free thiol at denaturing nuclear proteins. The proposed mechanism for nucleoprotein denaturation is through direct interactions with protein cysteine groups with the formation of destabilizing protein-WR-1065 disulfides. In comparison to its effect on nuclear proteins in isolated nuclei, WR-1065 had only a very small effect on non-nuclear proteins of whole cells, isolated nuclear matrix, or the thiol-rich Ca (2+) ATPase of sarcoplasmic reticulum, indicating that WR-1065 can effectively denature protein only inside an intact nucleus, probably due to the increased concentration of the positively charged drug in the vicinity of DNA.

Booth, Valerie K.↗
NASA NTRS1976

An investigation of flap-lag stability of wind turbine rotors in the presence of velocity gradients and helicopter rotors in forward flight

The flap-lag equations of motion of a torsionally rigid and centrally hinged spring-restrained rigid blade are developed using quasi-steady blade-element aerodynamic theory. These equations have periodic coefficients and are applicable to both wind-turbine rotors with velocity gradients and helicopter rotors in forward flight. By solving these equations both by the Floquet-Liapunov method and by an approximate method, the range of applicability of the latter method is established. Flap-lag and pure flapping stability boundaries illustrating the effects of flow and rotor parameters are presented. Finally, some discussion on the techniques for generating the Floquet transition matrix and on the relative merits of the coordinate axes system is presented.

Kaza, K. R. V.↗
NASA NTRS1990

R-matrix calculations for electron-impact excitation of C(+), N(2+), and O(3+) including fine structure

The new R-matrix package for comprehensive close-coupling calculations for electron scattering with the first three ions in the boron isoelectronic sequence, the astrophysically significant C(+), N(2+), and O(3+), is presented. The collision strengths are calculated in the LS coupling approximation, as well as in pair-coupling scheme, for the transitions among the fine-structure sublevels. Calculations are carried out at a large number of energies in order to study the detailed effects of autoionizing resonances.

Luo, D.↗
NASA NTRS1990

Application of global sensitivity equations in multidisciplinary aircraft synthesis

The present paper investigates the applicability of the Global Sensitivity Equation (GSE) method in the multidisciplinary synthesis of aeronautical vehicles. The GSE method provides an efficient approach for representing a large coupled system by smaller subsystems and accounts for the subsystem interactions by means of first-order behavior sensitivities. This approach was applied in an aircraft synthesis problem with performance constraints stemming from the disciplines of structures, aerodynamics, and flight mechanics. Approximation methods were considered in an attempt to reduce problem dimensionality and to improve the efficiency of the optimization process. The influence of efficient constraint representations, the choice of design variables, and design variable scaling on the conditioning of the system matrix was also investigated.

Hajela, P.↗
NASA NTRS1994

A simple suboptimal least-squares algorithm for attitude determination with multiple sensors

Three-axis attitude determination is equivalent to finding a coordinate transformation matrix which transforms a set of reference vectors fixed in inertial space to a set of measurement vectors fixed in the spacecraft. The attitude determination problem can be expressed as a constrained optimization problem. The constraint is that a coordinate transformation matrix must be proper, real, and orthogonal. A transformation matrix can be thought of as optimal in the least-squares sense if it maps the measurement vectors to the reference vectors with minimal 2-norm errors and meets the above constraint. This constrained optimization problem is known as Wahba's problem. Several algorithms which solve Wahba's problem exactly have been developed and used. These algorithms, while steadily improving, are all rather complicated. Furthermore, they involve such numerically unstable or sensitive operations as matrix determinant, matrix adjoint, and Newton-Raphson iterations. This paper describes an algorithm which minimizes Wahba's loss function, but without the constraint. When the constraint is ignored, the problem can be solved by a straightforward, numerically stable least-squares algorithm such as QR decomposition. Even though the algorithm does not explicitly take the constraint into account, it still yields a nearly orthogonal matrix for most practical cases; orthogonality only becomes corrupted when the sensor measurements are very noisy, on the same order of magnitude as the attitude rotations. The algorithm can be simplified if the attitude rotations are small enough so that the approximation sin(theta) approximately equals theta holds. We then compare the computational requirements for several well-known algorithms. For the general large-angle case, the QR least-squares algorithm is competitive with all other know algorithms and faster than most. If attitude rotations are small, the least-squares algorithm can be modified to run faster, and this modified algorithm is faster than all but a similarly specialized version of the QUEST algorithm. We also introduce a novel measurement averaging technique which reduces the n-measurement case to the two measurement case for our particular application, a star tracker and earth sensor mounted on an earth-pointed geosynchronous communications satellite. Using this technique, many n-measurement problems reduce to less than or equal to 3 measurements; this reduces the amount of required calculation without significant degradation in accuracy. Finally, we present the results of some tests which compare the least-squares algorithm with the QUEST and FOAM algorithms in the two-measurement case. For our example case, all three algorithms performed with similar accuracy.

Brozenec, Thomas F.↗
NASA NTRS1992

On adaptive weighted polynomial preconditioning for Hermitian positive definite matrices

The conjugate gradient algorithm for solving Hermitian positive definite linear systems is usually combined with preconditioning in order to speed up convergence. In recent years, there has been a revival of polynomial preconditioning, motivated by the attractive features of the method on modern architectures. Standard techniques for choosing the preconditioning polynomial are based only on bounds for the extreme eigenvalues. Here a different approach is proposed, which aims at adapting the preconditioner to the eigenvalue distribution of the coefficient matrix. The technique is based on the observation that good estimates for the eigenvalue distribution can be derived after only a few steps of the Lanczos process. This information is then used to construct a weight function for a suitable Chebyshev approximation problem. The solution of this problem yields the polynomial preconditioner. In particular, we investigate the use of Bernstein-Szego weights.

Fischer, Bernd↗
NASA NTRS2011

Sensitivity of C-Band Polarimetric Radar-Based Drop Size Distribution Measurements to Maximum Diameter Assumptions

The estimation of rain drop size distribution (DSD) parameters from polarimetric radar observations is accomplished by first establishing a relationship between differential reflectivity (Z(sub dr)) and the central tendency of the rain DSD such as the median volume diameter (D0). Since Z(sub dr) does not provide a direct measurement of DSD central tendency, the relationship is typically derived empirically from rain drop and radar scattering models (e.g., D0 = F[Z (sub dr)] ). Past studies have explored the general sensitivity of these models to temperature, radar wavelength, the drop shape vs. size relation, and DSD variability. Much progress has been made in recent years in measuring the drop shape and DSD variability using surface-based disdrometers, such as the 2D Video disdrometer (2DVD), and documenting their impact on polarimetric radar techniques. In addition to measuring drop shape, another advantage of the 2DVD over earlier impact type disdrometers is its ability to resolve drop diameters in excess of 5 mm. Despite this improvement, the sampling limitations of a disdrometer, including the 2DVD, make it very difficult to adequately measure the maximum drop diameter (D(sub max)) present in a typical radar resolution volume. As a result, D(sub max) must still be assumed in the drop and radar models from which D0 = F[Z(sub dr)] is derived. Since scattering resonance at C-band wavelengths begins to occur in drop diameters larger than about 5 mm, modeled C-band radar parameters, particularly Z(sub dr), can be sensitive to D(sub max) assumptions. In past C-band radar studies, a variety of D(sub max) assumptions have been made, including the actual disdrometer estimate of D(sub max) during a typical sampling period (e.g., 1-3 minutes), D(sub max) = C (where C is constant at values from 5 to 8 mm), and D(sub max) = M*D0 (where the constant multiple, M, is fixed at values ranging from 2.5 to 3.5). The overall objective of this NASA Global Precipitation Measurement Mission (GPM/PMM Science Team)-funded study is to document the sensitivity of DSD measurements, including estimates of D0, from C-band Z(sub dr) and reflectivity to this range of D(sub max) assumptions. For this study, GPM Ground Validation 2DVD's were operated under the scanning domain of the UAHuntsville ARMOR C-band dual-polarimetric radar. Approximately 7500 minutes of DSD data were collected and processed to create gamma size distribution parameters using a truncated method of moments approach. After creating the gamma parameter datasets the DSD's were then used as input to a T-matrix model for computation of polarimetric radar moments at C-band. All necessary model parameterizations, such as temperature, drop shape, and drop fall mode, were fixed at typically accepted values while the D(sub max) assumption was allowed to vary in sensitivity tests. By hypothesizing a DSD model with D(sub max) (fit) from which the empirical fit to D0 = F[Z(sub dr)] was derived via non-linear least squares regression and a separate reference DSD model with D(sub max) (truth), bias and standard error in D0 retrievals were estimated in the presence of Z(sub dr) measurement error and hypothesized mismatch in D(sub max) assumptions. Although the normalized standard error for D0 = F[Z(sub dr)r] can increase slightly (as much as from 11% to 16% for all 7500 DSDs) when the D(sub max) (fit) does not match D(sub max) (truth), the primary impact of uncertainty in D(sub max) is a potential increase in normalized bias error in D0 (from 0% to as much as 10% over all 7500 DSDs, depending on the extent of the mismatch between D(sub max) (fit) and D(sub max) (truth)). For DSDs characterized by large Z(sub dr) (Z(sub dr) > 1.5 to 2.0 dB), the normalized bias error for D0 estimation at C-band is sometimes unacceptably large (> 10%), again depending on the extent of the hypothesized D(sub max) mismatch. Modeled errors in D0 retrievals from Z(sub dr) at C-band are demonstrated in detail and comparedo similar modeled retrieval errors at S-band and X-band where the sensitivity to D(sub max) is expected to be less. The impact of D(sub max) assumptions to the retrieval of other DSD parameters such as Nw, the liquid water content normalized intercept parameter, are also explored. Likely implications for DSD retrievals using C-band polarimetric radar for GPM are assessed by considering current community knowledge regarding D(sub max) and quantifying the statistical distribution of Z(sub dr) from ARMOR over a large variety of meteorological conditions. Based on these results and the prevalence of C-band polarimetric radars worldwide, a call for more emphasis on constraining our observational estimate of D(sub max) within a typical radar resolution volume is made

Carey, Lawrence D.↗