Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear estimation”

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

Theoretical Characterization of C{sub 6}N, C{sub 6}N{sup −}, and C{sub 6}N{sup +}

This work emphasizes the stability of C{sub 6}N linear carbon chains and carbon clusters containing three-body rings. C{sub 6}N possesses at least 44 neutral isomers and 38 and 35 isomers with a negative or a positive charge. The lowest-energy structures, which can be candidates for laboratory and astrophysical detection, were studied with RCCSD(T)-F12 and MRCI/CASSCF, specifying properties for various electronic states. Neutral C{sub 6}N displays two prominent equilibrium structures, a nitrogen-terminated linear form (X{sup 2}Π) and a C{sub 2v} form (X{sup 2}B{sub 2}) containing a three-carbon ring. They are separated by 0.21 eV. For the linear one, Renner–Teller and spin–orbit effects are expected. Its equilibrium spin–orbit constant ∣ A{sub so,e} ∣was predicted to be 29.09 cm{sup −1}, the dipole moment was computed to be 0.2059 D, and B{sub e} was computed to be 868.20 MHz, very close to the experimental value (B{sub 0}(exp)-B{sub e}(calc) ∼ 4.9 MHz). The anion presents two stable forms, a linear triplet (X{sup 3}Σ{sup −}) and a C{sub 2v}(X{sup 1}A{sub 1}) isomer, whereas for the cation most stable structures are both linear. For the linear anion and the two linear cations B{sub 0} was estimated to be 870.05 MHz, 875.72 MHz, and 917.36 MHz and μ was estimated to be 4.4200 D, 4.2013D, and 0.7828D. For the neutral C{sub 2V} form, Ae = 45048.74 MHz, Be = 1062.78 MHz, Ce = 1038.30 MHz, and μ = 0.8840 D, whereas for the negative C{sub 2V} form, Ae = 40087.48 MHz, Be = 1067.35 MHz, Ce = 1039.67 MHz, and μ = 6.2854 D. Below the electron affinity of the most stable isomer (EA = −3.42 eV), the linear anion displays three probably existing electronic states. Detectability is discussed in terms of the symmetry and spin multiplicity of the ground electronic states.

79 ASTRONOMY AND ASTROPHYSICS↗

Compact Polarimetry in a Low Frequency Spaceborne Context

Compact polarimetry has been shown to be an interesting alternative mode to full polarimetry when global coverage and revisit time are key issues. It consists on transmitting a single polarization, while receiving on two. Several critical points have been identified, one being the Faraday rotation (FR) correction and the other the calibration. When a low frequency electromagnetic wave travels through the ionosphere, it undergoes a rotation of the polarization plane about the radar line of sight for a linearly polarized wave, and a simple phase shift for a circularly polarized wave. In a low frequency radar, the only possible choice of the transmit polarization is the circular one, in order to guaranty that the scattering element on the ground is illuminated with a constant polarization independently of the ionosphere state. This will allow meaningful time series analysis, interferometry as long as the Faraday rotation effect is corrected for the return path. In full-polarimetric (FP) mode, two techniques allow to estimate the FR: Freeman method using linearly polarized data, and Bickel and Bates theory based on the transformation of the measured scattering matrix to a circular basis. In CP mode, an alternate procedure is presented which relies on the bare surface scattering properties. These bare surfaces are selected by the conformity coefficient, invariant with FR. This coefficient is compared to other published classifications to show its potential in distinguishing three different scattering types: surface, doublebounce and volume. The performances of the bare surfaces selection and FR estimation are evaluated on PALSAR and airborne data. Once the bare surfaces are selected and Faraday angle estimated over them, the correction can be applied over the whole scene. The algorithm is compared with both FP techniques. In the last part of the paper, the calibration of a CP system from the point of view of classical matrix transformation methods in polarimetry is proposed.

radar polarization diversity↗

TH-CD-206-12: Image-Based Motion Estimation for Plaque Visualization in Coronary Computed Tomography Angiography

Purpose: Visualization and quantification of coronary artery calcification and atherosclerotic plaque benefits from coronary artery motion (CAM) artifact elimination. This work applies a rigid linear motion model to a Volume of Interest (VoI) for estimating motion estimation and compensation of image degradation in Coronary Computed Tomography Angiography (CCTA). Methods: In both simulation and testbench experiments, translational CAM was generated by displacement of the imaging object (i.e. simulated coronary artery and explanted human heart) by ∼8 mm, approximating the motion of a main coronary branch. Rotation was assumed to be negligible. A motion degraded region containing a calcification was selected as the VoI. Local residual motion was assumed to be rigid and linear over the acquisition window, simulating motion observed during diastasis. The (negative) magnitude of the image gradient of the reconstructed VoI was chosen as the motion estimation objective and was minimized with Covariance Matrix Adaptation Evolution Strategy (CMAES). Results: Reconstruction incorporated the estimated CAM yielded signification recovery of fine calcification structures as well as reduced motion artifacts within the selected local region. The compensated reconstruction was further evaluated using two image similarity metrics, the structural similarity index (SSIM) and Root Mean Square Error (RMSE). At the calcification site, the compensated data achieved a 3% increase in SSIM and a 91.2% decrease in RMSE in comparison with the uncompensated reconstruction. Conclusion: Results demonstrate the feasibility of our image-based motion estimation method exploiting a local rigid linear model for CAM compensation. The method shows promising preliminary results for the application of such estimation in CCTA. Further work will involve motion estimation of complex motion corrupted patient data acquired from clinical CT scanner.

60 APPLIED LIFE SCIENCES↗

Estimation of abatement potentials and costs of air pollution emissions in China

Understanding the air pollution emission abatement potential and associated control cost is a prerequisite to design cost ef?cient control policies. In this study, a linear programming algorithm model, International Control Cost Estimate Tool, was updated with cost data for applications of 56 types of end-of-pipe technologies and ?ve types of renewable energy in 10 major sectors namely power generation, industry combustion, cement pro-duction, iron and steel production, other industry processes, domestic combustion, transportation, solvent use, livestock rearing, and fertilizer use. The updated model was implemented to estimate the abatement potential and marginal cost of multiple pollutants in China. The total maximum abatement potentials of sulfur dioxide (SO2), nitrogen oxides (NOx), primary particulate matter (PM2.5), non-volatile organic compounds (NMVOCs), and ammonia (NH3) in China were estimated to be 19.2, 20.8, 9.1, 17.2 and 8.6 Mt, respectively, which accounted for 89.7%, 89.9%, 94.6%, 74.0%, and 80.2% of their total emissions in 2014, respectively. The associated control cost of such reductions was estimated as 92.5, 469.7, 75.7, 449.0, and 361.8 billion CNY in SO2, NOx, primary PM2.5, NMVOCs and NH3, respectively. Shandong, Jiangsu, Henan, Zhejiang, and Guangdong provinces exhibited large abatement potentials for all pollutants. Provincial disparity analysis shows that high GDP regions tend to have higher reduction potential and total abatement costs. End-of-pipe technologies tended be a cost-ef?cient way to control pollution in industries processes (i.e., cement plants, iron and steel plants, lime production, building ceramic production, glass and brick production), whereas such technologies were less cost- effective in fossil fuel-related sectors (i.e., power plants, industry combustion, domestic combustion, and transportation) compared with renewable energy. The abatement potentials and marginal abatement cost curves developed in this study can further be used as a crucial component in an integrated model to design optimized cost-ef?cient control policies.

Zhang, Fenfen↗

Model-Based Engine Control Architecture with an Extended Kalman Filter

This paper discusses the design and implementation of an extended Kalman filter (EKF) for model-based engine control (MBEC). Previously proposed MBEC architectures feature an optimal tuner Kalman Filter (OTKF) to produce estimates of both unmeasured engine parameters and estimates for the health of the engine. The success of this approach relies on the accuracy of the linear model and the ability of the optimal tuner to update its tuner estimates based on only a few sensors. Advances in computer processing are making it possible to replace the piece-wise linear model, developed off-line, with an on-board nonlinear model running in real-time. This will reduce the estimation errors associated with the linearization process, and is typically referred to as an extended Kalman filter. The non-linear extended Kalman filter approach is applied to the Commercial Modular Aero-Propulsion System Simulation 40,000 (C-MAPSS40k) and compared to the previously proposed MBEC architecture. The results show that the EKF reduces the estimation error, especially during transient operation.

engine control↗

Input-output characterization of an ultrasonic testing system by digital signal analysis

Ultrasonic test system input-output characteristics were investigated by directly coupling the transmitting and receiving transducers face to face without a test specimen. Some of the fundamentals of digital signal processing were summarized. Input and output signals were digitized by using a digital oscilloscope, and the digitized data were processed in a microcomputer by using digital signal-processing techniques. The continuous-time test system was modeled as a discrete-time, linear, shift-invariant system. In estimating the unit-sample response and frequency response of the discrete-time system, it was necessary to use digital filtering to remove low-amplitude noise, which interfered with deconvolution calculations. A digital bandpass filter constructed with the assistance of a Blackman window and a rectangular time window were used. Approximations of the impulse response and the frequency response of the continuous-time test system were obtained by linearly interpolating the defining points of the unit-sample response and the frequency response of the discrete-time system. The test system behaved as a linear-phase bandpass filter in the frequency range 0.6 to 2.3 MHz. These frequencies were selected in accordance with the criterion that they were 6 dB below the maximum peak of the amplitude of the frequency response. The output of the system to various inputs was predicted and the results were compared with the corresponding measurements on the system.

Williams, J. H., Jr.↗

Model-Based Engine Control Architecture with an Extended Kalman Filter

This paper discusses the design and implementation of an extended Kalman filter (EKF) for model-based engine control (MBEC). Previously proposed MBEC architectures feature an optimal tuner Kalman Filter (OTKF) to produce estimates of both unmeasured engine parameters and estimates for the health of the engine. The success of this approach relies on the accuracy of the linear model and the ability of the optimal tuner to update its tuner estimates based on only a few sensors. Advances in computer processing are making it possible to replace the piece-wise linear model, developed off-line, with an on-board nonlinear model running in real-time. This will reduce the estimation errors associated with the linearization process, and is typically referred to as an extended Kalman filter. The nonlinear extended Kalman filter approach is applied to the Commercial Modular Aero-Propulsion System Simulation 40,000 (C-MAPSS40k) and compared to the previously proposed MBEC architecture. The results show that the EKF reduces the estimation error, especially during transient operation.

engine control↗

wtDigiTwin (Wind Turbine Digital Twin) [SWR-21-15]

This wind turbine digital twin software (wtDigiTwin) provides a digital twin solution for wind turbine applications. The focus of wtDigiTwin is to estimate loads, motions and environmental conditions for an operating wind turbine. The program uses supervisory control and data acquisition (SCADA) measurements as inputs, together with a wind turbine model. The wind industry is currently challenged by the high cost of operation and maintenance. These costs could be mitigated if component failures are predicted, but such predictions are difficult unless the turbines are equipped with expensive measuring devices. The alternative is to use a digital twin such as wtDigiTwin to estimate the necessary signals. wtDigiTwin can perform online prediction of signals that are otherwise not measured, using a limited set of reliable measurements and a physics-based model. The predicted signals can be used in applications that have direct cost benefits: 1) real-time estimation of the fatigue consumption of key components of the wind turbine; 2) root cause analyses and failure detections ; 3) lifetime reassessments ; 4) improvements to follow-on designs. The current version provides examples to estimate wind speed, thrust, torque, tower-top position, and tower loads on an onshore wind turbine using the following measurements tower top acceleration, generator torque, pitch, and rotational speed. The model combines a linear state-space model, a wind speed estimator, and a Kalman filter algorithm that integrates measurements with the state model to perform state estimations. The state space model is obtained either using OpenFAST linearizations, or using the yams package provided with the software.

Branlard, Emmanuel↗

Fermionic Partial Tomography via Classical Shadows

Here we propose a tomographic protocol for estimating any k-body reduced density matrix (k-RDM) of an n-mode fermionic state, a ubiquitous step in near-term quantum algorithms for simulating many-body physics, chemistry, and materials. Our approach extends the framework of classical shadows, a randomized approach to learning a collection of quantum-state properties, to the fermionic setting. Our sampling protocol uses randomized measurement settings generated by a discrete group of fermionic Gaussian unitaries, implementable with linear-depth circuits. We prove that estimating all k-RDM elements to additive precision ϵ requires on the order of ($^{n}_{k}$)k 3/2 log(n)/ϵ 2 repeated state preparations, which is optimal up to the logarithmic factor. Furthermore, numerical calculations show that our protocol offers a substantial improvement in constant overheads for k ≥ 2, as compared to prior deterministic strategies. We also adapt our method to particle-number symmetry, wherein the additional circuit depth may be halved at the cost of roughly 2–5 times more repetitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimation and filter stability of stochastic delay systems

Linear and nonlinear filtering for stochastic delay systems are studied. A representation theorem for conditional moment functionals is obtained, which, in turn, is used to derive stochastic differential equations describing the optimal linear or nonlinear filter. A complete characterization of the optimal filter is given for linear systems with Gaussian noise. Stability of the optimal filter is studied in the case where there are no delays in the observations. Using the duality between linear filtering and control, asymptotic stability of the optimal filter is proved. Finally, the cascade of the optimal filter and the deterministic optimal quadratic control system is shown to be asymptotically stable as well.

Kwong, R. H.↗

On the error in crop acreage estimation using satellite (LANDSAT) data

The problem of crop acreage estimation using satellite data is discussed. Bias and variance of a crop proportion estimate in an area segment obtained from the classification of its multispectral sensor data are derived as functions of the means, variances, and covariance of error rates. The linear discriminant analysis and the class proportion estimation for the two class case are extended to include a third class of measurement units, where these units are mixed on ground. Special attention is given to the investigation of mislabeling in training samples and its effect on crop proportion estimation. It is shown that the bias and variance of the estimate of a specific crop acreage proportion increase as the disparity in mislabeling rates between two classes increases. Some interaction is shown to take place, causing the bias and the variance to decrease at first and then to increase, as the mixed unit class varies in size from 0 to 50 percent of the total area segment.

Chhikara, R.↗

Robust detection, isolation and accommodation for sensor failures

The objective is to extend the recent advances in robust control system design of multivariable systems to sensor failure detection, isolation, and accommodation (DIA), and estimator design. This effort provides analysis tools to quantify the trade-off between performance robustness and DIA sensitivity, which are to be used to achieve higher levels of performance robustness for given levels of DIA sensitivity. An innovations-based DIA scheme is used. Estimators, which depend upon a model of the process and process inputs and outputs, are used to generate these innovations. Thresholds used to determine failure detection are computed based on bounds on modeling errors, noise properties, and the class of failures. The applicability of the newly developed tools are demonstrated on a multivariable aircraft turbojet engine example. A new concept call the threshold selector was developed. It represents a significant and innovative tool for the analysis and synthesis of DiA algorithms. The estimators were made robust by introduction of an internal model and by frequency shaping. The internal mode provides asymptotically unbiased filter estimates.The incorporation of frequency shaping of the Linear Quadratic Gaussian cost functional modifies the estimator design to make it suitable for sensor failure DIA. The results are compared with previous studies which used thresholds that were selcted empirically. Comparison of these two techniques on a nonlinear dynamic engine simulation shows improved performance of the new method compared to previous techniques

Emami-Naeini, A.↗

Variations in solar Lyman alpha irradiance on short time scales

Variations in solar UV irradiance at Lyman alpha are studied on short time scales (from days to months) after removing the long-term changes over the solar cycle. The SME/Lyman alpha irradiance is estimated from various solar indices using linear regression analysis. In order to study the nonlinear effects, Lyman alpha irradiance is modeled with a 5th-degree polynomial as well. It is shown that the full-disk equivalent width of the He line at 1083 nm, which is used as a proxy for the plages and active network, can best reproduce the changes observed in Lyman alpha. Approximately 72 percent of the solar-activity-related changes in Lyman alpha irradiance arise from plages and the network. The network contribution is estimated by the correlation analysis to be about 19 percent. It is shown that significant variability remains in Lyman alpha irradiance, with periods around 300, 27, and 13.5d, which is not explained by the solar activity indices. It is shown that the nonlinear effects cannot account for a significant part of the unexplained variation in Lyman alpha irradiance. Therefore, additional events (e.g., large-scale motions and/or a systematic difference in the area and intensity of the plages and network observed in the lines of Ca-K, He 1083, and Lyman alpha) may explain the discrepancies found between the observed and estimated irradiance values.

Pap, J. M.↗

Determining Crop Coefficients Using Remote Sensing for the Maipo River Valley Basin in Chile

Agriculture is the primary use of water in the Maipo River Valley of Central Chile, accounting for ~ 75% of the total demand. Assessment of irrigation needs for agricultural production has commonly relied on reference crop coefficients (Kc) derived from geographic and climatic conditions that differ from those of Chile. This work focused on calculating site-specific Kc values tailored to crop production in the water-stressed Maipo River Valley over consecutive growing seasons from 2019 to 2022. Two distinct approaches were implemented, each relying on remotely-sensed NASA Earth observation datasets. The first method estimated Kc values based on the linear relationship with the Normalized Difference Vegetation Index (NDVI) obtained from either Terra Moderate Resolution Imaging Spectroradiometer (MODIS) or Landsat 8 Operational Land Imager (OLI) surface reflectance. The second technique leveraged information from the ISS Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) by computing the ratio between actual crop evapotranspiration and potential evapotranspiration. The former approach best captured vegetation signals of annual crops while the latter appeared best suited for perennials. Overall, this study provides a strong basis and novel way to accurately estimate Kc using remote sensing, with the potential for improved irrigation management and reduction in water consumption.

Rishudh Thadur↗

An improved method for direct incident solar radiation calculation from hourly solar insolation data in building energy simulation

Solar radiation considerably influences the energy consumption of buildings and the power production of building integrated photovoltaic (BIPV) systems. Hourly solar insolation (Wh/m 2 ), represented as the amount of solar irradiance collected on the ground during a 1-h period, is the most common solar radiation data available and widely used in weather files applied in building energy modeling programs (BEMPs). Because the solar beam and position vary over time, the use of hourly insolation data as the input might result in errors in the estimation of the direct incident solar radiation on a particular surface. In this study, methods used in BEMPs for direct incident solar radiation calculations are first analyzed, and an improved method adopting a new algorithm for estimating the solar irradiance is proposed. Furthermore, the algorithm assumes that the solar irradiance changes linearly within a 1-h period and can be estimated based on the solar irradiance at the half clock and slope. The collected direct normal solar irradiance data of 2016 from eight solar radiation stations in China were used to demonstrate the proposed method and evaluate its performance by comparing the results with those from three conventional methods used in BEMPs along with the ground truth measurements. In addition, in this study, factors affecting the accuracy of the calculation results are explored. The results of the estimated direct incident solar radiation show that the proposed method achieves the best accuracy, followed by the methods used in DOE-2, EnergyPlus, and DeST. The proposed method guarantees that the hourly direct solar insolation will remain the same and reflects the variation in the direct solar irradiance across a 1-h time frame. The proposed method can be adopted in BEMPs to improve the accuracy of the solar radiation calculation, thereby improving the accuracy of the simulated building performance and the BIPV production.

14 SOLAR ENERGY↗

Exchange field effect in the crystal-field ground state of Ce M Al 4 Si 2

The crystal-field ground-state wave functions of the tetragonal, magnetically ordering Kondo lattice materials CeMAl 4 Si 2 (M = Rh, Ir, and Pt) are determined in this paper with low-temperature linearly polarized soft-x-ray absorption spectroscopy, and estimates for the crystal-field splittings are given from the temperature evolution of the linear dichroism. Values for the dominant exchange field in the magnetically ordered phases can be obtained from fitting the influence of magnetic order on the linear dichroism. The direction of the required exchange field is || c for the antiferromagnetic Rh and Ir compounds, with the corresponding strength of the order of λ ex ≈ 6 meV (65 K). Finally and furthermore, the presence of Kondo screening in the Rh and Ir compound is demonstrated on the basis of the absorption due to f 0 in the initial state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measuring Success: A Refined Methodology for Estimating Long-term Continuous Improvement

Successful resource management systems require current and detailed feedback on operational and corporate-level performance. As corporate accountability concerns intensify, the precision and reliability of these performance metrics have become crucial. Traditional savings estimation methods can be difficult to understand, particularly linear regression, and can provide varying results. This paper reviews common efficiency metrics and highlights underlying mathematical inconsistencies when estimating total and percent savings with current methods. A refined approach to calculating long-term utility savings is proposed that simplifies current methodologies utilizing ratios to define an adjusted baseline, allowing for consistent and fair aggregation of results across multiple scales from resources to corporate performance. A simplified example demonstrates how the proposed methodology improves upon existing methods, especially in intermediate years. This paper’s major contributions are the simplified approach for converting modeled utility usage into estimated savings and the consistent roll-up methodology enabling more comparable, aggregable, and actionable results across scales.

Price, Chris [ORNL] (ORCID:0000000202007906)↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗