Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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 775 records · Page 43

MATRIX-VBS (v1.0): Implementing an Evolving Organic Aerosol Volatility in an Aerosol Microphysics Model

The gas-particle partitioning and chemical aging of semi-volatile organic aerosol are presented in a newly developed box model scheme, where its effect on the growth, composition, and mixing state of particles is examined. The volatility-basis set (VBS) framework is implemented into the aerosol microphysical scheme MATRIX (Multiconfiguration Aerosol TRacker of mIXing state), which resolves mass and number aerosol concentrations and in multiple mixing-state classes. The new scheme, MATRIX-VBS, has the potential to significantly advance the representation of organic aerosols in Earth system models by improving upon the conventional representation as non-volatile particulate organic matter, often also with an assumed fixed size distribution. We present results from idealized cases representing Beijing, Mexico City, a Finnish forest, and a southeastern US forest, and investigate the evolution of mass concentrations and volatility distributions for organic species across the gas and particle phases, as well as assessing their mixing state among aerosol populations. Emitted semi-volatile primary organic aerosols evaporate almost completely in the intermediate-volatility range, while they remain in the particle phase in the low-volatility range. Their volatility distribution at any point in time depends on the applied emission factors, oxidation by OH radicals, and temperature. We also compare against parallel simulations with the original scheme, which represented only the particulate and non-volatile component of the organic aerosol, examining how differently the condensed-phase organic matter is distributed across the mixing states in the model. The results demonstrate the importance of representing organic aerosol as a semi-volatile aerosol, and explicitly calculating the partitioning of organic species between the gas and particulate phases.

volatility-basis set↗

Human motion planning based on recursive dynamics and optimal control techniques

This paper presents an efficient optimal control and recursive dynamics-based computer animation system for simulating and controlling the motion of articulated figures. A quasi-Newton nonlinear programming technique (super-linear convergence) is implemented to solve minimum torque-based human motion-planning problems. The explicit analytical gradients needed in the dynamics are derived using a matrix exponential formulation and Lie algebra. Cubic spline functions are used to make the search space for an optimal solution finite. Based on our formulations, our method is well conditioned and robust, in addition to being computationally efficient. To better illustrate the efficiency of our method, we present results of natural looking and physically correct human motions for a variety of human motion tasks involving open and closed loop kinematic chains.

Non-NASA Center↗

Improving the chi-squared approximation for bivariate normal tolerance regions

Let X be a two-dimensional random variable distributed according to N2(mu,Sigma) and let bar-X and S be the respective sample mean and covariance matrix calculated from N observations of X. Given a containment probability beta and a level of confidence gamma, we seek a number c, depending only on N, beta, and gamma such that the ellipsoid R = (x: (x - bar-X)'S(exp -1) (x - bar-X) less than or = c) is a tolerance region of content beta and level gamma; i.e., R has probability gamma of containing at least 100 beta percent of the distribution of X. Various approximations for c exist in the literature, but one of the simplest to compute -- a multiple of the ratio of certain chi-squared percentage points -- is badly biased for small N. For the bivariate normal case, most of the bias can be removed by simple adjustment using a factor A which depends on beta and gamma. This paper provides values of A for various beta and gamma so that the simple approximation for c can be made viable for any reasonable sample size. The methodology provides an illustrative example of how a combination of Monte-Carlo simulation and simple regression modelling can be used to improve an existing approximation.

Feiveson, Alan H.↗

Echo mapping of active galactic nuclei broad-line regions: Fundamental algorithms

We formulate and test a series of algorithms for echo mapping the emission-line regions near active galactic nuclei from measurements of correlated variability in their line and continuum light curves. The linear regularization method (LRM) employs a direct inversion of evenly spaced light-curve data, with a regularization parameter that can be used to control the trade-off between noise and resolution. Matrix formulas express the formal solution as well as its variance and covariance in terms of uncertainties in the measurements. Unlike the maximum-entropy method (MEM), LRM applies to kernels with both positive and negative values, but the results are somewhat limited by ringing effects. A positivity constraint proves effective in controlling the ringing. MEM combines regularization and positivity in a natural way, but similar results are also found using positivity constraints with nonentropic regularization functions. Direct inversions of unevenly sampled light curves require interpolating the noisy data. In this case better results are found by solving for both the continuum light curve and kernel function in a simultaneous fit to the data. Our conclusion is that while echo mapping currently gives ambiguous results, the algorithms are not the limiting factor. Progress depends on efforts to increase the accuracy and completeness of sampling of the observed light curves.

Vio, Roberto↗

Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.

Swanson, Roy C., Jr.↗

Assimilation of satellite microwave observations over the rainbands of tropical cyclones

A novel Bayesian Monte Carlo integration (BMCI) technique was developed to retrieve geophysical variables from satellite microwave radiometer data in the presence of tropical cyclones. The BMCI technique includes three steps: generating a stochastic database, simulating satellite brightness temperatures using a radiative transfer model, and retrieving geophysical variables such as profiles of temperature, relative humidity, and cloud liquid and ice water content from real observations. The technique also provides uncertainty estimates for each retrieval and can output the error covariance matrix of selected parameters. The measurements from the Advanced Technology Microwave Sounder (ATMS) on board Suomi National Polar-Orbiting Partnership (Suomi NPP) and the Global Precipitation Measurement (GPM) Microwave Imager (GMI) were used as input. A new technique was developed to correct the ATMS and GMI observations for the beam-filling effect, which is due to small-scale variability of precipitation and clouds when compared with the instrument footprint and also the nonlinear relation between the brightness temperature and precipitation. In addition, the assimilation of the BMCI retrievals into the NASA GEOS model is discussed for Hurricane Maria. The results show that assimilating the BMCI retrievals can influence the dynamical features of the cyclone, including a stronger warm core, a symmetric eye, and vertically aligned wind columns. Two possible factors that may limit the impact of the BMCI retrievals include 1) the resolution of the model (about 25 km), which was too coarse to show the potential of the BMCI data in improving the representation of tropical storms in the model forecast, and 2) the data assimilation system not being able to consider vertically correlated observation errors.

Isaac Moradi↗

The effects of gravity darkening on the ultraviolet continuum polarization produced by circumstellar disks

We investigate the effects of gravity darkening on the UV continuum polarization produced by an axisymmetric disk that surrounds a rapidly rotating star. Although the model is a single scattering approximation, we do include the effects of attenuation (electron scattering plus hydrogen bound-free absorption) by the disk, using an approach similar to that of Sobolev (1963). Because of the gravity darkening of the star and the attenuation within the disk, the radiation field is not axially symmetric about the radius vector. This implies that the polarization source functions are no longer provided by the finite disk depolarization factors of Cassinelli, Nordsieck, & Murison (1987), which are functions of the intensity moments in a spherically symmetric atmosphere. We reformulate the polarization source functions using generalized intensity moment tensors (J, H(sub i), K(sub ij)) that are valid for an arbitrary radiation field and envelope geometry. We find that the polarization source functions are simplest when using intensity moments in the observer's reference frame. On the other hand, the intensity moments are most easily evaluated in the stellar reference frame. Using the rotation transformation properties of the generalized intensity moments, we relate the observer's moments to those evaluated in the stellar reference frame. Our procedure for determining the polarization source functions thus merely involves choosing a set of Euler angles for the coordinate rotations, and then evaluating the associated rotation matrix. The geometrical complications of polarization transfer are thus reduced to obtaining a coordinate rotation matrix.

Bjorkman, J. E.↗

Interactions between integrin receptors and fibronectin are required for calvarial osteoblast differentiation in vitro

We previously showed that anti-fibronectin antibodies or soluble fibronectin fragments containing the central cell-binding domain inhibit formation of mineralized nodules by fetal calvarial osteoblasts in vitro. These findings suggest a critical role for fibronectin in osteoblast differentiation and morphogenesis. In this study we tested the hypothesis that fibronectin's effects on osteogenesis are mediated via direct interactions with integrin receptors for fibronectin on osteoblasts. Immunocytochemical analysis identified the integrin fibronectin receptor alpha5ss1 in fetal rat calvarial tissue and in cultured osteoblasts at all stages of differentiation. Three other integrins, alpha3ss1, alpha8ss1 and alphavss3, which can bind fibronectin, as well as other matrix components, were also identified in tissue and at all stages of cell culture. Immunoprecipitation data showed that alpha5ss1 levels are constant throughout osteoblast differentiation whereas levels of alpha3ss1 and alpha8ss1 decline in mature mineralized cultures. To determine whether integrin fibronectin receptors are required for osteoblast formation of mineralized nodules, we examined the extent of nodule formation in the presence and absence of function-perturbing anti-integrin antibodies. The antibodies were present continuously in cultures beginning at confluence (day 3), and nodule formation was measured at days 10 and 20. An anti-alpha5 integrin subunit antibody reduced nodule formation to less than 5% of control values at both time points. Inhibition of nodule formation was reversible and did not affect cell attachment and viability. Function-perturbing antibodies against alpha3ss1 and alpha8ss1 also reduced nodule formation, to less than 20% of control values. In contrast, function-perturbing antibodies to alphavss3 and alphavss5 did not affect nodule formation, indicating that the inhibitions noted were indeed specific. To determine the effect of antibody treatment on gene expression, steady-state mRNA expression was examined and found to be suppressed for osteoblast markers alkaline phosphatase and osteocalcin. Together, these results indicate that direct osteoblast interactions with the extracellular matrix are mediated by a select group of integrin receptors that includes alpha5ss1, alpha3ss1 and alpha8ss1. We further conclude that the specific alpha5ss1 fibronectin receptor mediates critical interactions between osteoblasts and fibronectin required for both bone morphogenesis and osteoblast differentiation.

NASA Center ARC↗

Influence of Finite Element Size in Residual Strength Prediction of Composite Structures

The sensitivity of failure load to the element size used in a progressive failure analysis (PFA) of carbon composite center notched laminates is evaluated. The sensitivity study employs a PFA methodology previously developed by the authors consisting of Hashin-Rotem intra-laminar fiber and matrix failure criteria and a complete stress degradation scheme for damage simulation. The approach is implemented with a user defined subroutine in the ABAQUS/Explicit finite element package. The effect of element size near the notch tips on residual strength predictions was assessed for a brittle failure mode with a parametric study that included three laminates of varying material system, thickness and stacking sequence. The study resulted in the selection of an element size of 0.09 in. X 0.09 in., which was later used for predicting crack paths and failure loads in sandwich panels and monolithic laminated panels. Comparison of predicted crack paths and failure loads for these panels agreed well with experimental observations. Additionally, the element size vs. normalized failure load relationship, determined in the parametric study, was used to evaluate strength-scaling factors for three different element sizes. The failure loads predicted with all three element sizes provided converged failure loads with respect to that corresponding with the 0.09 in. X 0.09 in. element size. Though preliminary in nature, the strength-scaling concept has the potential to greatly reduce the computational time required for PFA and can enable the analysis of large scale structural components where failure is dominated by fiber failure in tension.

Satyanarayana, Arunkumar↗

Stress failure of pulmonary capillaries: role in lung and heart disease

Pulmonary capillaries have extremely thin walls to allow rapid exchange of respiratory gases across them. Recently it has been shown that the wall stresses become very large when the capillary pressure is raised, and in anaesthetised rabbits, ultrastructural damage to the walls is seen at pressures of 40 mm Hg and above. The changes include breaks in the capillary endothelial layer, alveolar epithelial layer, and sometimes all layers of the wall. The strength of the thin part of the capillary wall can be attributed to the type IV collagen in the extracellular matrix. Stress failure of pulmonary capillaries results in a high-permeability form of oedema, or even frank haemorrhage, and is apparently the mechanism of neurogenic pulmonary oedema and high-altitude pulmonary oedema. It also explains the exercise-induced pulmonary haemorrhage that occurs in all racehorses. Several features of mitral stenosis are consistent with stress failure. Overinflation of the lung also leads to stress failure, a common cause of increased capillary permeability in the intensive care environment. Stress failure also occurs if the type IV collagen of the capillary wall is weakened by autoantibodies as in Goodpasture's syndrome. Neutrophil elastase degrades type IV collagen and this may be the starting point of the breakdown of alveolar walls that is characteristic of emphysema. Stress failure of pulmonary capillaries is a hitherto overlooked and potentially important factor in lung and heart disease.

Review↗

Leveraging metal matrix composites to reduce costs in space mechanisms

Advanced metal matrix composites may be one of the most promising technologies for reducing cost in structural components without compromise to strength or stiffness. A microlight 12.50 N (2.81 lb), two-axis, solar array drive assembly (SADA) was made for the Advanced Materials Applications to Space Structures (AMASS) Program flight experiment. The SADA had both its inner and outer axis housings fabricated from silicon carbide particulate reinforced alumimun. Two versions of the housings were made. The first was machined from a solid billet of material. The second was plaster cast to a near net shape that required minimal finish machining. Both manufacturing methods were compared upon completion. Results showed a cost savings with the cast housing was possible for quantities greater than one and probable for quantities greater than two. For quantities approaching ten, casting resulted in a reduction factor of almost three in the cost per part.

Nye, Ted↗

Convergence Acceleration for Multistage Time-Stepping Schemes

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 could be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. Numerical dissipation operators (based on the Roe scheme, a matrix formulation, and the CUSP scheme) as well as the number of RK stages are considered in evaluating the RK/implicit scheme. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. In two dimensions, turbulent flows over an airfoil at subsonic and transonic conditions are computed. The effects of mesh cell aspect ratio on convergence are investigated for Reynolds numbers between 5.7 x 10(exp 6) and 100.0 x 10(exp 6). Results are also obtained for a transonic wing flow. For both 2-D and 3-D problems, the computational time of a well-tuned standard RK scheme is reduced at least a factor of four.

Swanson, R. C.↗

Closed-form analysis of fiber-matrix interface stresses under thermo-mechanical loadings

Closed form techniques for calculating fiber matrix (FM) interface stresses, using repeating square and diamond regular arrays, were presented for a unidirectional composite under thermo-mechanical loadings. An Airy's stress function micromechanics approach from the literature, developed for calculating overall composite moduli, was extended in the present study to compute FM interface stresses for a unidirectional graphite/epoxy (AS4/3501-6) composite under thermal, longitudinal, transverse, transverse shear, and longitudinal shear loadings. Comparison with finite element results indicate excellent agreement of the FM interface stresses for the square array. Under thermal and longitudinal loading, the square array has the same FM peak stresses as the diamond array. The square array predicted higher stress concentrations under transverse normal and longitudinal shear loadings than the diamond array. Under transverse shear loading, the square array had a higher stress concentration while the diamond array had a higher radial stress concentration. Stress concentration factors under transverse shear and longitudinal shear loadings were very sensitive to fiber volume fraction. The present analysis provides a simple way to calculate accurate FM interface stresses for both the square and diamond array configurations.

Naik, Rajiv A.↗

Supporting Crew Autonomy in Deep Space Exploration: Preliminary Onboard Capability Requirements and Proposed Research Questions. Technical Report of the Autonomous Crew Operations Technical Interchange Meeting

Communication delays are a critical challenge posed by long duration deep space exploration. Space missions historically have relied on an ever-present Mission Control Center (MCC) to direct operations in near real-time. As unanticipated anomalies that defeat fault detection and resolution systems do arise, the lack of real-time communication will significantly weaken what the MCC support represents: a reliable safety net for the flight crew through its deep and diverse areas of expertise and investigative resources. As a consequence, future space vehicles and habitats need to be equipped with capabilities to support the flight crew to operate with little or no ground support. Considerations must be given to vehicle and mission designs that will fortify the traditionally ground-centered safety net and forge new support systems, when communication delays exist. In August 2018, NASA’s Human Research Program, through its Human Factors and Behavioral Performance Element, convened a Technical Interchange Meeting (TIM) on Autonomous Crew Operations at NASA Ames Research Center. The goal of the meeting was to gather input from NASA centers, industry, academia, and branches of the Department of Defense (DoD) to address how intelligent technologies can be applied to augment onboard capabilities to support crew anomaly response. The TIM featured 24 presentations by 29 speakers and hosted a total of 59 attendees, including 43 from 5 NASA centers (Ames, Johnson, Langley, Marshall, and Jet Propulsion Lab) and 4 from the DoD (3 from Army Research Lab and 1 from Naval Postgraduate School), with remaining attendees from academia (e.g., UC Davis, CMU) and industry (e.g., IBM, Siemens). Discussions were centered around three themes: standards and guidelines, lessons learned in analog environments, and technologies. To help provide a framework for discussion, a concept matrix describing anomaly response processes was created prior to the TIM (Figure 1, page 6). The matrix captures the steps involved (monitoring and detection, diagnosis, solution development and evaluation, solution implementation and verification, resolution documentation) as well as the resources and capabilities required to support these steps (data, knowledge, analysis, synthesis, resource management). A wallpaper size printout of the matrix was utilized at the TIM to solicit attendee inputs along the three themes; the activity garnered 108 submissions of ideas. Overall, what emerged from TIM discussions was a picture of mismatch between crew anomaly response needs and support that can be provided by existing intelligent technologies. The needs are broad, spanning multiple steps and processes/resources, with many of which lacking support from existing technologies, such as knowledge management throughout the steps of problem solving (especially in resolution documentation) and manpower management. The solutions provided by existing intelligent technologies are specific to the steps/processes that they are designed to support and constrained to solving only problems similar to those that have occurred before. What is lacking from technologies is typically made up by humans, specifically their complex critical thinking, creative problem solving, and domain expertise. In the end, the TIM highlighted the pressing need to support responses to onboard anomalies during autonomous crew operations, particularly those that have eluded the system tests, inspection, and other assurance processes. Such anomalies can potentially threaten crew and vehicle safety, as well as significantly impact overall operations with additional workload. These fairly rare events are difficult to anticipate and prepare for, given the state-of-the-art in intelligent technologies. This is true even for anomalies that stem from “unknown knowns”—cases in which there is sufficient external information to characterize the problem but the overall pattern fails to be recognized by the problem solver, or in which the internal knowledge needed to solve a problem is held tacitly and potentially accessible by the problem solver but not articulated. It follows that the ability to tackle anomalies lies not only with the availability of relevant information and knowledge but also their accessibility in times of need. To that end, we propose research questions along the following three broad themes: • How intelligent technologies can help make relevant knowledge and information available? • How intelligent technologies can help make relevant knowledge and information accessible? • How intelligent technologies can help support the crew operating as a team in anomaly response processes?

autonomous crew operations↗

Strength Scaling in Fiber Composites

A research program was initiated to study and isolate the factors responsible for scale effects in the tensile strength of graphite/epoxy composite laminates. Four layups were chosen with appropriate stacking sequences so as to highlight individual and interacting failure modes. Four scale sizes were selected for investigation including full scale size, 3/4, 2/4, and 1/4, with n = to 4, 3, 2, and 1, respectively. The full scale specimen sizes was 32 piles thick as compared to 24, 16, and 8 piles for the 3/4, 2/4, and 1/4 specimen sizes respectively. Results were obtained in the form of tensile strength, stress-strain curves and damage development. Problems associated with strength degradation with increasing specimen sizes are isolated and discussed. Inconsistencies associated with strain measurements were also identified. Enhanced x ray radiography was employed for damage evaluation, following step loading. It was shown that fiber dominated layups were less sensitive to scaling effects compared to the matrix dominated layups.

Kellas, Sotiris↗

Precision of Sensitivity in the Design Optimization of Indeterminate Structures

Design sensitivity is central to most optimization methods. The analytical sensitivity expression for an indeterminate structural design optimization problem can be factored into a simple determinate term and a complicated indeterminate component. Sensitivity can be approximated by retaining only the determinate term and setting the indeterminate factor to zero. The optimum solution is reached with the approximate sensitivity. The central processing unit (CPU) time to solution is substantially reduced. The benefit that accrues from using the approximate sensitivity is quantified by solving a set of problems in a controlled environment. Each problem is solved twice: first using the closed-form sensitivity expression, then using the approximation. The problem solutions use the CometBoards testbed as the optimization tool with the integrated force method as the analyzer. The modification that may be required, to use the stiffener method as the analysis tool in optimization, is discussed. The design optimization problem of an indeterminate structure contains many dependent constraints because of the implicit relationship between stresses, as well as the relationship between the stresses and displacements. The design optimization process can become problematic because the implicit relationship reduces the rank of the sensitivity matrix. The proposed approximation restores the full rank and enhances the robustness of the design optimization method.

Patnaik, Surya N.↗

Strength scaling in fiber composites

A research program was initiated to study and isolate the factors responsible for scale effects in the tensile strength of graphite/epoxy composite laminates. Four layups were chosen with appropriate stacking sequences so as to highlight individual and interacting failure modes. Four scale sizes were selected for investigation including full scale size, 3/4, 2/4, and 1/4, with n = to 4, 3, 2, and 1, respectively. The full scale specimen sizes was 32 piles thick as compared to 24, 16, and 8 piles for the 3/4, 2/4, and 1/4 specimen sizes respectively. Results were obtained in the form of tensile strength, stress-strain curves and damage development. Problems associated with strength degradation with increasing specimen sizes are isolated and discussed. Inconsistencies associated with strain measurements were also identified. Enchanced X-ray radiography was employed for damage evaluation, following step loading. It was shown that fiber dominated layups were less sensitive to scaling effects compared to the matrix dominated layups.

Kellas, Sotiris↗

Peru Health and Air Quality II: Leveraging Earth Observations and Health Data to Map Outbreak Risk and Inform Public Health Interventions for Zoonotic Disease Prevention

Peru's Madre de Dios region is a hotspot for dengue fever and leishmaniasis due to its tropical Amazonian climate. Though treatable, these zoonotic diseases are debilitating for under-resourced communities whose already close proximity to mosquito and sandfly vectors continues to increase via rapid urbanization and deforestation. Peru’s Ministries of Health (MINSA) and Environment (MINAM) are working to better understand the environmental factors amplifying the risk of dengue fever and leishmaniasis transmission. The first term of this project classified the land use and land cover of Madre de Dios’ 11 districts for 2010, 2015, and 2020 and identified a correlation between both diseases and urbanization. Our team expanded this analysis by creating urban-forest edge maps and incorporating climatic and topographic variables with data from Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI), the Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals for GPM (IMERG), and the Shuttle Radar Topography Mission (SRTM). We determined these variables’ impacts on disease incidence by assessing existing literature and running regression models. Dengue fever correlated with urban-forest edge, urban area, slope, temperature, and precipitation. Leishmaniasis primarily correlated with forest-edge area and elevation, but lacking additional statistical significance prevented further work, a decision supported by the literature. Thus, the risk matrix and risk map which we scripted in R to visualize the risk of disease posed to districts alongside health post locations addresses only dengue fever. The results and products will inform MINSA and MINAM in public health interventions, resource distribution, and policy initiatives.

Jennifer Rogers↗