Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 199 records · Page 11

Anthropogenic enhancements to production of highly oxygenated molecules from autoxidation

As particulate matter abatement strategies reduce anthropogenic nitrogen oxide (NOx) emissions, the chemical regime of an urban atmosphere may change to one in which autoxidation becomes increasingly important. Developed countries such as the United States have already experienced large declines in nitrogen oxides because of controls on emissions, and the possibility for increased aerosol exists as autoxidation results in highly oxidized molecules (HOM) which efficiently form particulate matter. Using a-pinene reaction with hydroxyl radicals as an archetype system, a molecular representation of organic aerosol was created and evaluated against laboratory experiments. The applica-tion of parameters constrained by these experiments to the ambient atmosphere indicated autoxidation is becoming more competitive in urban areas of the United States while at the same time HOM production rates are declining due to decreases in oxidant abundance. For locations such as Guangzhou, China, where emissions of NOx generally increased from the 1990s to late 2000s, HOM production was also predicted to increase. These results demonstrate that anthro-pogenic emissions can enhance HOM and the resulting biogenic organic aerosol formation through oxidant abundance. Ambient particulate matter is responsible for a significant fraction of the global burden of disease (Lim et al., 2012) and affects climate (Myhre et al., 2013). The organic portion of ambient aerosol largely forms due to the condensation of low vapor pressure (Jimenez et al., 2009) or highly soluble (Saxena and Hildemann, 1996) gases in the atmosphere . Oxidation of monoterpenes from vegetation, results in significant formation of particle mass (Pandis et al., 1991), and is likely the dominant source of organic aerosol in the southeastern United States (Zhang et al., 2018; Xu et al., 2018). While monoterpenes are included as a secondary organic aerosol (SOA) source in most chemical transport models, many current regional to global scale model parameterizations (Pye et al., 2010, 2017; Koo et al., 2014) lack a mechanistic dependence of monoterpene SOA on NOx and oxidant identity (e.g. OH vs. ozone). This limits the ability of models to predict how changes in emissions will affect ambient concentrations and new particle formation events. Building a mechanistic understanding of formation pathways is critical for being able to determine historic drivers of ambient particulate matter as well as inform predictions of future levels

Pye, Havala O.↗

Sandia Unstructured Triangle Table Generator

The software generates data tables for thermodynamic and transport properties of materials as described by a set of input models. For each input model parameterization, an associated table is created on an unstructured triangular grid. These grids all conform to the same topology. A statistical accuracy guarantee is provided for the tabular representation of the model. Details of the model and table specification are given in a XML input deck.

Carpenter, JohnH.↗

Application of aerosol-ice nucleating particle closure to establish the leading parameters governing ice crystal number concentration under commonly observed mixed-phase cloud conditions

This is the final technical report for DOE DE-SC0021034 award entitled “Application of aerosol-ice nucleating particle closure to establish the leading parameters governing ice crystal number concentration under commonly observed mixed-phase cloud conditions”. The overarching objective of this project was to gain a predictive understanding of the processes that lead from aerosols to ice-nucleating particles (INPs) to ice crystal number concentrations (ICNCs) in stratiform mixed-phase clouds where supercooled water droplets and ice crystals coexist, and to inform climate model parameterizations. During this four-year project we addressed the three main research goals of this objective: 1) Study the aerosol-INP link; 2) Develop an aerosol-cloud column model (AC-1D); and 3) Study the INP-ICNC link. Those activities significantly advanced our predictive understanding of how to represent INPs in could models and its implications on ice crystal number concentrations. The results have been published in five peer-reviewed publications with three more in review and preparation. It allowed for collaboration resulting in additional fourteen peer-reviewed publications that are concerned with the objective of this research project. This research has been communicated to the scientific community at national and international meetings and conferences by numerous talks, seminars, and posters including fifteen invited presentations. This project allowed training of four graduate students and one postdoctoral research associate resulting in one M.S. thesis and advancement to Ph.D. candidacy.

54 ENVIRONMENTAL SCIENCES↗

Incorporating Valve Switching Losses Into a Static Optimal Control Algorithm for the Hybrid Hydraulic-Electric Architecture (HHEA)

The Hybrid Hydraulic-Electric Architecture (HHEA) has, in recent years, been proposed as an energy efficient alternative to conventional load-sensing architectures in mobile machines such as excavators and wheel-loaders. HHEA leverages the advantages of hydraulic power and electric power to eliminate throttling valves while also improving the energy and control performance of the system. The architecture utilizes a set of common pressure rails to provide a majority of power and and a small electric motor driven pump to modulate this power to meet the exact demand. Previous work has developed a computationally efficient Lagrange Multiplier approach for determining the optimal pressure rail selections that minimizes the energy losses in the system. The static model used considers only the energy use for each pressure rail selection but not the losses associated with the valves during the transition. This paper presents an approach to include the switching losses in the model and in the optimization procedure. To capture the switching losses, switching events between different rails and at various input and output flow rates were simulated with consideration of valve spool dynamics. A parameterized model that summarizes the losses is then obtained, allowing switching losses to be added to the previous energy analysis. The performance of the switching loss model was compared with reference data obtained from a high-fidelity simulation model. To incorporate the switching losses into optimal control algorithm, an efficient dynamic programming approach that prevents frequent switching is adopted in place of the Lagrange multiplier approach. The overall effect of switching losses on the energy consumption and optimal control decisions is presented. In general, switching losses contribute to about 9–10% of input energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Understanding the Effect of Aerosol Properties on Cloud Droplet Formation during TCAP Field Campaign Report

The formation of clouds is an essential element in understanding the Earth’s radiative budget. Liquid water clouds form when the relative humidity exceeds saturation and condensedphase water nucleates on atmospheric particulate matter. The effect of aerosol properties such as size, morphology, and composition on cloud droplet formation has been studied theoretically as well as in the laboratory and field. Almost without exception these studies have been limited to parallel measurements of aerosol properties and cloud formation or collection of material after the cloud has formed, at which point nucleation information has been lost. Studies of this sort are adequate when a large fraction of the aerosol activates, but correlations and resulting model parameterizations are much more uncertain at lower supersaturations and activated fractions.

54 ENVIRONMENTAL SCIENCES↗

Polynomial chaos expansions for dependent random variables

Polynomial chaos expansions (PCE) are well-suited to quantifying uncertainty in models parameterized by independent random variables. The assumption of independence leads to simple strategies for building multivariate orthonormal bases and for sampling strategies to evaluate PCE coefficients. In contrast, the application of PCE to models of dependent variables is much more challenging. Three approaches can be used to construct PCE of models of dependent variables. The first approach uses mapping methods where measure transformations, such as the Nataf and Rosenblatt transformation, can be used to map dependent random variables to independent ones; however we show that this can significantly degrade performance since the Jacobian of the map must be approximated. A second strategy is the class of dominating support methods. In these approaches a PCE is built using independent random variables whose distributional support dominates the support of the true dependent joint density; we provide evidence that this approach appears to produce approximations with suboptimal accuracy. A third approach, the novel method proposed here, uses Gram-Schmidt orthogonalization (GSO) to numerically compute orthonormal polynomials for the dependent random variables. This approach has been used successfully when solving differential equations using the intrusive stochastic Galerkin method, and in this paper we use GSO to build PCE using a non-intrusive stochastic collocation method. The stochastic collocation method treats the model as a black box and builds approximations of the input-output map from a set of samples. Building PCE from samples can introduce ill-conditioning which does not plague stochastic Galerkin methods. To mitigate this ill-conditioning we generate weighted Leja sequences, which are nested sample sets, to build accurate polynomial interpolants. We show that our proposed approach, GSO with weighted Leja sequences, produces PCE which are orders of magnitude more accurate than PCE constructed using mapping or dominating support methods.

97 MATHEMATICS AND COMPUTING↗

A fully Bayesian method for jointly fitting instrumental calibration and X-ray spectral models

Owing to a lack of robust principled methods, systematic instrumental uncertainties have generally been ignored in astrophysical data analysis despite wide recognition of the importance of including them. Ignoring calibration uncertainty can cause bias in the estimation of source model parameters and can lead to underestimation of the variance of these estimates. We previously introduced a pragmatic Bayesian method to address this problem. The method is 'pragmatic' in that it introduced an ad hoc technique that simplified computation by neglecting the potential information in the data for narrowing the uncertainty for the calibration product. Following that work, we use a principal component analysis to efficiently represent the uncertainty of the effective area of an X-ray (or γ-ray) telescope. Here, however, we leverage this representation to enable a principled, fully Bayesian method that coherently accounts for the calibration uncertainty in high-energy spectral analysis. In this setting, the method is compared with standard analysis techniques and the pragmatic Bayesian method. The advantage of the fully Bayesian method is that it allows the data to provide information not only for estimation of the source parameters but also for the calibration product—here the effective area, conditional on the adopted spectral model. In this way, it can yield more accurate and efficient estimates of the source parameters along with valid estimates of their uncertainty. Provided that the source spectrum can be accurately described by a parameterized model, this method allows rigorous inference about the effective area by quantifying which possible curves are most consistent with the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Interior-Point Methods for Estimating Seasonal Parameters in Discrete-Time Infectious Disease Models

Infectious diseases remain a significant health concern around the world. Mathematical modeling of these diseases can help us understand their dynamics and develop more effective control strategies. In this work, we show the capabilities of interior-point methods and nonlinear programming (NLP) formulations to efficiently estimate parameters in multiple discrete-time disease models using measles case count data from three cities. These models include multiplicative measurement noise and incorporate seasonality into multiple model parameters. Our results show that nearly identical patterns are estimated even when assuming seasonality in different model parameters, and that these patterns show strong correlation to school term holidays across very different social settings and holiday schedules. We show that interior-point methods provide a fast and flexible approach to parameterizing models that can be an alternative to more computationally intensive methods.

60 APPLIED LIFE SCIENCES↗

Number of sites-based solver for determining coverages from steady-state mean-field micro-kinetic models

Kinetic models parameterized by ab-initio calculations have led to significant improvements in understanding chemical reactions in heterogeneous catalysis. These studies have been facilitated by implementations which determine steady-state coverages and rates of mean-field micro-kinetic models. As implemented in the open-source kinetic modeling program, CatMAP, the conventional solution strategy is to use a root-finding algorithm to determine the coverage of all intermediates through the steady-state expressions, constraining all coverages to be non-negative and to properly sum to unity. Though intuitive, this root-finding strategy causes issues with convergence to solution due to these imposed constraints. In this work, we avoid explicitly imposing these constraints, solving the mean-field steady-state micro-kinetic model in the space of number of sites instead of solving it in the space of coverages. Here, we transform the constrained root-finding problem to an unconstrained least-squares minimization problem, leading to significantly improved convergence in solving micro-kinetic models and thus enabling the efficient study of more complex catalytic reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterization of Exoplanet Atmospheres with the Optical Coronagraph on WFIRST

Wide-Field Infrared Survey Telescope (WFIRST)-CGI is a NASA technology demonstration mission that is charged with demonstrating key technologies for future exo-Earth imaging missions in space. In the process, it will obtain images and low-resolution spectra of a handful to a dozen extrasolar planets and possibly protoplanetary disks. Its unprecedented contrast levels in the optical will provide astronomers’ with their first direct look at mature, Jupiter-sized planets at moderate separations. This paper addresses the question: what science can be done with such data? An analytic noise model, which is informed by the ongoing engineering developments, is used to compute maximum achievable signal-to-noise ratios and scientifically viable integration times for hypothetical star–planet systems, as well as to investigate the constraining power of various combinations of WFIRST-CGI photometric and spectral observations. This work introduces two simple models for planetary geometric albedos, which are inspired largely by the solar system’s gas giants. The first planet model is a hybrid Jupiter–Neptune model, which separately treats the short and long wavelengths where chromophores and methane dominate absorption, respectively. The second planet model fixes cloud and haze properties in CoolTLusty to match Jupiter’s albedo spectrum, it then perturbs only the metallicity. MCMC retrievals performed on simulated observations are used to assess the precision with which planet model parameters can be measured subject to different exposure times and observing cases. Planet radius is recovered within ±15% for all observing cases with both the hybrid model and the CoolTLusty metallicity grid. Fit results for both models’ parameterizations of geometric albedo spectra demonstrate that a rough indication of the metallicity or methane content should be possible for some WFIRST-CGI targets. We conclude that real observations will likely be able to differentiate between extreme cases using these models, but will lack the precision necessary to uncover subtle trends.

79 ASTRONOMY AND ASTROPHYSICS↗

Population of anatomically variable 4D XCAT adult phantoms for imaging research and optimization

Purpose: The authors previously developed the 4D extended cardiac-torso (XCAT) phantom for multimodality imaging research. The XCAT consisted of highly detailed whole-body models for the standard male and female adult, including the cardiac and respiratory motions. In this work, the authors extend the XCAT beyond these reference anatomies by developing a series of anatomically variable 4D XCAT adult phantoms for imaging research, the first library of 4D computational phantoms. Methods: The initial anatomy of each phantom was based on chest-abdomen-pelvis computed tomography data from normal patients obtained from the Duke University database. The major organs and structures for each phantom were segmented from the corresponding data and defined using nonuniform rational B-spline surfaces. To complete the body, the authors manually added on the head, arms, and legs using the original XCAT adult male and female anatomies. The structures were scaled to best match the age and anatomy of the patient. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from the template XCAT phantom (male or female) to the target patient model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. Each new phantom was refined by checking for anatomical accuracy via inspection of the models. Results: Using these methods, the authors created a series of computerized phantoms with thousands of anatomical structures and modeling cardiac and respiratory motions. The database consists of 58 (35 male and 23 female) anatomically variable phantoms in total. Like the original XCAT, these phantoms can be combined with existing simulation packages to simulate realistic imaging data. Each new phantom contains parameterized models for the anatomy and the cardiac and respiratory motions and can, therefore, serve as a jumping point from which to create an unlimited number of 3D and 4D variations for imaging research. Conclusions: A population of phantoms that includes a range of anatomical variations representative of the public at large is needed to more closely mimic a clinical study or trial. The series of anatomically variable phantoms developed in this work provide a valuable resource for investigating 3D and 4D imaging devices and the effects of anatomy and motion in imaging. Combined with Monte Carlo simulation programs, the phantoms also provide a valuable tool to investigate patient-specific dose and image quality, and optimization for adults undergoing imaging procedures.

60 APPLIED LIFE SCIENCES↗

A set of 4D pediatric XCAT reference phantoms for multimodality research

Purpose: The authors previously developed an adult population of 4D extended cardiac-torso (XCAT) phantoms for multimodality imaging research. In this work, the authors develop a reference set of 4D pediatric XCAT phantoms consisting of male and female anatomies at ages of newborn, 1, 5, 10, and 15 years. These models will serve as the foundation from which the authors will create a vast population of pediatric phantoms for optimizing pediatric CT imaging protocols. Methods: Each phantom was based on a unique set of CT data from a normal patient obtained from the Duke University database. The datasets were selected to best match the reference values for height and weight for the different ages and genders according to ICRP Publication 89. The major organs and structures were segmented from the CT data and used to create an initial pediatric model defined using nonuniform rational B-spline surfaces. The CT data covered the entire torso and part of the head. To complete the body, the authors manually added on the top of the head and the arms and legs using scaled versions of the XCAT adult models or additional models created from cadaver data. A multichannel large deformation diffeomorphic metric mapping algorithm was then used to calculate the transform from a template XCAT phantom (male or female 50th percentile adult) to the target pediatric model. The transform was applied to the template XCAT to fill in any unsegmented structures within the target phantom and to implement the 4D cardiac and respiratory models in the new anatomy. The masses of the organs in each phantom were matched to the reference values given in ICRP Publication 89. The new reference models were checked for anatomical accuracy via visual inspection. Results: The authors created a set of ten pediatric reference phantoms that have the same level of detail and functionality as the original XCAT phantom adults. Each consists of thousands of anatomical structures and includes parameterized models for the cardiac and respiratory motions. Based on patient data, the phantoms capture the anatomic variations of childhood, such as the development of bone in the skull, pelvis, and long bones, and the growth of the vertebrae and organs. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. Conclusions: The development of patient-derived pediatric computational phantoms is useful in providing variable anatomies for simulation. Future work will expand this ten-phantom base to a host of pediatric phantoms representative of the public at large. This can provide a means to evaluate and improve pediatric imaging devices and to optimize CT protocols in terms of image quality and radiation dose.

61 RADIATION PROTECTION AND DOSIMETRY↗

Polar primary aerosols across the ocean-sea ice-snow-atmosphere interface: From sources to impacts

Primary aerosols play a critical role in polar climate systems, influencing cloud formation, precipitation, radiative balance, and surface energy budgets. This paper provides a comprehensive synthesis of primary aerosol sources, transformation and removal processes, and broader atmospheric impacts in polar regions, emphasizing their links to ocean and sea ice biogeochemistry. These aerosols (including sea salt, primary organic aerosol, and primary biological aerosol particles) originate from marine and cryospheric environments and are emitted through physical processes, such as wave breaking, bubble bursting, and blowing snow. Emission sources include seawater, sea ice, snow, and freshwater from river discharge and glacial runoff. Once airborne, these particles can serve as a chemical reservoir, influencing atmospheric composition and reactivity, and as seeds for cloud droplet and ice crystal formation, influencing cloud microphysics and polar climate. Despite their importance, many of the processes governing primary aerosol emissions and transformations remain poorly constrained. The most pressing knowledge gaps pertain to emission processes, limited spatiotemporal observational coverage, instrumentation constraints, parameterization development, and the integration of interdisciplinary expertise. To improve our understanding of primary aerosol drivers and their response to climate, future research efforts should prioritize strategically coordinated and cross-disciplinary process studies, advancements in measurement technologies and coverage, and close collaboration between modelers and observational scientists to inform and refine model parameterizations. As polar regions continue to undergo profound changes marked by increased precipitation, reduced sea and land ice, freshening oceans, and shifting ecosystem dynamics, characterizing present-day primary aerosol populations is vital. Improved understanding will be essential for anticipating future changes in aerosol-radiation and aerosol-cloud interactions and their implications for polar and global climate systems.

Aerosol-cloud↗

Assimilating Scanning Radar Data into High-Resolution Models

This report documents the findings of a project aiming to improve the initial conditions of km-scale simulations of deep convective storms by assimilating cloud-scale weather radar observations. More accurate numerical analyses will increase the effectiveness of research efforts using LES cloud models as tool to better understand land-atmosphere coupling, boundary layer turbulence, and cloud processes, each used for model parameterization development. Work reported on herein includes: i) assessment of radar data sets for use in data assimilation (‘DA’) experiments, ii) quality control of the radar data set and format conversion to one acceptable to DA schemes utilized by the Weather Research and Forecasting (WRF) model, and iii) examination of the sensitivity of the numerical representation of cloud-scale wind and microphysical features to a variety of tunable DA parameters.

54 ENVIRONMENTAL SCIENCES↗

Parametric Modeling and Economic Analysis of a 2MW th 3-Stream sCO 2 Heat Exchanger

Here, this paper presents the design and cost optimization of a novel 2MW th 3-stream sCO2 plate-fin heat exchanger. This heat exchanger design is unique in that it uses reduced metal oxide particle-to-sCO2 heat exchanger for cost-effective energy storage applications. The design uses low velocity, laminar air as the re-oxidizing reactant to transfer the heat of the re-oxidizing reaction to a sCO2 power loop. The design of the heat exchanger is based on a 2-D, 3-fluid plate/fin heat transfer model. The model parameterizes the size, shape, and number of passages of the heat exchanger to calculate the temperature profile, pressure drop, and fluid velocities of all three fluids. Global heat exchanger parameters such as the effectiveness and total heat transferred to the sCO2 are then calculated for overall performance. Due to the value and increased use of sCO2 heat exchangers in power cycles, a cost model of the system based on the unique high temperature/high pressure operating conditions was created using quotes from reference projects and market analysis. These quoted air-to-sCO2 heat exchangers are then processed using multiple weighting factors pertinent to heat exchanger design, including heat exchanger type, maximum temperature, differential pressures, fluids, duty, and more. These factors are then used in an exponential function in order to generate a parameterized cost curve. The design and cost of the heat exchanger are then optimized using the SMPSO genetic algorithm in Python. The optimization objectives for the system are to maximize the overall system effectiveness, including an air recuperator for preheating, and to minimize unit costs. Additional constraints are added to the system for the sCO2 and air pressure drops, air velocity to reduce particle entrainment, and the length and volume of the heat exchanger.

Cost Model↗

Radiation-induced segregation in dilute Fe–Cr: A rate-theory framework for the Cr enrichment–depletion transition at the grain boundary

Radiation-induced segregation (RIS) poses a significant challenge for ferritic Fe–Cr alloys under irradiation, as it can compromise mechanical integrity and increase susceptibility to intergranular corrosion. Yet, the mechanisms governing Cr segregation remain incompletely understood. Here, in this study, we present a physics-based rate-theory model parameterized using self-consistent mean field theory-based Onsager transport coefficients to investigate RIS at the grain boundary in dilute Fe-(0.1 at. %) Cr. Under equal production rates of vacancies and self-interstitial atoms (SIA), and their equal absorption rates by bulk dislocations, the model simulates the experimentally observed transition from Cr enrichment at low temperatures to depletion at higher temperatures. Under these unbiased conditions, systematic investigation reveals that while temperature-dependent transport properties dictate the segregation direction, dose rate, grain size, and dislocation density only influence the magnitude and spatial extent of Cr segregation. However, under more realistic conditions of preferential vacancy production within damage cascade and/or preferential SIA absorption by bulk dislocations, the enrichment-to-depletion transition shifts to lower temperatures. Our findings demonstrate that RIS predictions based solely on transport coefficients are valid only under symmetric point defect flux conditions, and that biases in defect production and absorption must be considered for accurate predictions. This work provides a mechanistic framework for understanding RIS in ferritic alloys and informs alloy design for advanced nuclear systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evidence for Sub-Chandrasekhar Mass Type Ia Supernovae from an Extensive Survey of Radiative Transfer Models

There exists two classes of viable progenitors for normal Type Ia supernovae (SNe Ia): systems in which a white dwarf explodes at the Chandrasekhar mass (${M}_{\mathrm{ch}}$), and systems in which a white dwarf explodes below the Chandrasekhar mass (sub-${M}_{\mathrm{ch}}$). It is not clear which of these channels is dominant; observations and light-curve modeling have provided evidence for both. In this work we use an extensive grid of 4500 time-dependent, multiwavelength radiation transport simulations to show that the sub-${M}_{\mathrm{ch}}$ model can reproduce the entirety of the width–luminosity relation, while the ${M}_{\mathrm{ch}}$ model can only produce the brighter events $(0.8\lt {\rm{\Delta }}{M}_{15}(B)\lt 1.55)$, implying that fast-declining SNe Ia come from sub-${M}_{\mathrm{ch}}$ explosions. We do not assume a particular theoretical paradigm for the progenitor or explosion mechanism, but instead construct parameterized models that vary the mass, kinetic energy, and compositional structure of the ejecta, thereby realizing a broad range of possible outcomes of white dwarf explosions. We provide fitting functions based on our large grid of detailed simulations that map observable properties of SNe Ia, such as peak brightness and light-curve width, to physical parameters such as ${}^{56}\mathrm{Ni}$ and total ejected mass. These can be used to estimate the physical properties of observed SNe Ia.

79 ASTRONOMY AND ASTROPHYSICS↗

HIERARCHICAL PROBABILISTIC INFERENCE OF COSMIC SHEAR

Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing surveys. The first part of our statistical framework requires specification of a likelihood function for the pixel data in an imaging survey given parameterized models for the galaxies in the images. We derive the lensing shear posterior by marginalizing over all intrinsic galaxy properties that contribute to the pixel data (i.e., not limited to galaxy ellipticities) and learn the distributions for the intrinsic galaxy properties via hierarchical inference with a suitably flexible conditional probabilitiy distribution specification. We use importance sampling to separate the modeling of small imaging areas from the global shear inference, thereby rendering our algorithm computationally tractable for large surveys. With simple numerical examples we demonstrate the improvements in accuracy from our importance sampling approach, as well as the significance of the conditional distribution specification for the intrinsic galaxy properties when the data are generated from an unknown number of distinct galaxy populations with different morphological characteristics.

79 ASTRONOMY AND ASTROPHYSICS↗