Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic timing analysis”

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 307 records · Page 17

Adaptive methods for nonlinear structural dynamics and crashworthiness analysis

The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Belytschko, Ted↗

Analysis of sea ice dynamics

The ongoing work has established the basis for using multiyear sea ice concentrations from SMMR passive microwave for studies of largescale advection and convergence/divergence of the Arctic sea ice pack. Comparisons were made with numerical model simulations and buoy data showing qualitative agreement on daily to interannual time scales. Analysis of the 7-year SMMR data set shows significant interannual variations in the total area of multiyear ice. The scientific objective is to investigate the dynamics, mass balance, and interannual variability of the Arctic sea ice pack. The research emphasizes the direct application of sea ice parameters derived from passive microwave data (SMMR and SSMI) and collaborative studies using a sea ice dynamics model. The possible causes of observed interannual variations in the multiyear ice area are being examined. The relative effects of variations in the large scale advection and convergence/divergence within the ice pack on a regional and seasonal basis are investigated. The effects of anomolous atmospheric forcings are being examined, including the long-lived effects of synoptic events and monthly variations in the mean geostrophic winds. Estimates to be made will include the amount of new ice production within the ice pack during winter and the amount of ice exported from the pack.

Zwally, J.↗

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗

Time and Frequency-Domain Cross-Verification of SLS 6DOF Trajectory Simulations

The SLS GNC team and its partners have developed several time- and frequency-based simulations for development and analysis of the proposed SLS launch vehicle. The simulations differ in fidelity and some have unique functionality that allows them to perform specific analyses. Some examples of the purposes of the various models are: trajectory simulation, multi-body separation, Monte Carlo, hardware in the loop, loads, and frequency domain stability analyses. While no two simulations are identical, many of the models are essentially six degree-of-freedom (6DOF) representations of the SLS plant dynamics, hardware implementation, and flight software. Thus at a high level all of those models should be in agreement. Comparison of outputs from several SLS trajectory and stability analysis tools are ongoing as part of the program's current verification effort. The purpose of these comparisons is to highlight modeling and analysis differences, verify simulation data sources, identify inconsistencies and minor errors, and ultimately to verify output data as being a good representation of the vehicle and subsystem dynamics. This paper will show selected verification work in both the time and frequency domain from the current design analysis cycle of the SLS for several of the design and analysis simulations. In the time domain, the tools that will be compared are MAVERIC, CLVTOPS, SAVANT, STARS, ARTEMIS, and POST 2. For the frequency domain analysis, the tools to be compared are FRACTAL, SAVANT, and STARS. The paper will include discussion of these tools including their capabilities, configurations, and the uses to which they are put in the SLS program. Determination of the criteria by which the simulations are compared (matching criteria) requires thoughtful consideration, and there are several pitfalls that may occur that can severely punish a simulation if not considered carefully. The paper will discuss these considerations and will present a framework for responding to these issues when they arise. For example, small event timing differences can lead to large differences in mass properties if the criteria are to measure those properties at the same time, or large differences in altitude if the criteria are to measure those properties when the simulation experiences a staging event. Similarly, a tiny difference in phase can lead to large gain margin differences for frequency-domain comparisons of gain margins.

VanZwieten, Tannen↗

Time and Frequency-Domain Cross-Verification of SLS 6DOF Trajectory Simulations

The Space Launch System (SLS) Guidance, Navigation, and Control (GNC) team and its partners have developed several time- and frequency-based simulations for development and analysis of the proposed SLS launch vehicle. The simulations differ in fidelity and some have unique functionality that allows them to perform specific analyses. Some examples of the purposes of the various models are: trajectory simulation, multi-body separation, Monte Carlo, hardware in the loop, loads, and frequency domain stability analyses. While no two simulations are identical, many of the models are essentially six degree-of-freedom (6DOF) representations of the SLS plant dynamics, hardware implementation, and flight software. Thus at a high level all of those models should be in agreement. Comparison of outputs from several SLS trajectory and stability analysis tools are ongoing as part of the program's current verification effort. The purpose of these comparisons is to highlight modeling and analysis differences, verify simulation data sources, identify inconsistencies and minor errors, and ultimately to verify output data as being a good representation of the vehicle and subsystem dynamics. This paper will show selected verification work in both the time and frequency domain from the current design analysis cycle of the SLS for several of the design and analysis simulations. In the time domain, the tools that will be compared are MAVERIC, CLVTOPS, SAVANT, STARS, ARTEMIS, and POST 2. For the frequency domain analysis, the tools to be compared are FRACTAL, SAVANT, and STARS. The paper will include discussion of these tools including their capabilities, configurations, and the uses to which they are put in the SLS program. Determination of the criteria by which the simulations are compared (matching criteria) requires thoughtful consideration, and there are several pitfalls that may occur that can severely punish a simulation if not considered carefully. The paper will discuss these considerations and will present a framework for responding to these issues when they arise. For example, small event timing differences can lead to large differences in mass properties if the criteria are to measure those properties at the same time, or large differences in altitude if the criteria are to measure those properties when the simulation experiences a staging event. Similarly, a tiny difference in phase can lead to large gain margin differences for frequency-domain comparisons of gain margins.

Johnson, Matthew↗

Three-Dimensional Electromagnetic and Structural Analysis of Disruptions in the COMPASS Upgrade Vacuum Vessel

Structural analyses of disruptions and the electromagnetic (EM) analyses to support them typically use symmetry to produce a manageable model size. Nonaxisymmetric halo loads require larger analysis models. Here, in this paper, the results of transient EM analyses of halo strikes during a vertical displacement event are presented. The plasma motions and halo characteristics are prescribed based on analyses performed by the Institute of Plasma Physics of the Czech Republic. The time transient EM solution provides loads to a full three-dimensional transient dynamic analysis of the vacuum vessel. The responses to large lateral halo loads are altered and mitigated by the mode responses of the vessel. Dynamic load factors are computed for specific locations in the vessel and supports.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rising Above the Cloud - Toward High-Rate Delay-Tolerant Networking in Low-Earth Orbit

The High Data Rate Architecture (HiDRA) project is implementing a High-rate Delay Tolerant Networking (HDTN) capability that can support Low Earth Orbit (LEO) applications and environments. The present state of the effort, future work, and other elements of the work to date are described in this paper. This implementation is intended to support applications that run at 1+ Gbps, per the requirements of modern optical and high-frequency RF links. Uniquely, this implementation is also tuned to support relay and data trunking applications, which might require support for large numbers of small bundles per second. The design for this platform is based entirely on commercial-off-the-shelf (COTS) components, and possesses buffering capabilities in the 5 TB range. This document takes results from previous individual tests and integrates them to demonstrate results in the presence of a coherent use-case: consider a network aboard the ISS which intends to utilize an upcoming optical communications capability. For this use-case, orbital analysis software is used to analyze orbital dynamics, from which a list of access times are generated that might take in to account weather, schedule competition, etc. A variant of Contact Graph Routing (CGR) is applied to these windows to determine an optimal schedule. This schedule is then loaded into the HDTN prototype and, in conjunction with various measurement tools, a complete end-to-end analysis of HDTN's performance is conducted. Various bottlenecks (including storage) are identified: these bottlenecks are expected to help us focus our future work on the elements of the system that are most likely to present issues moving forward. Finally, we discuss possible paths for evolution beyond the present rates supported by the system, including (but not limited to) hardware acceleration.

Hylton, Alan↗

A Categorization of Dynamic Analyzers

Program analysis techniques and tools are essential to the development process because of the support they provide in detecting errors and deficiencies at different phases of development. The types of information rendered through analysis includes the following: statistical measurements of code, type checks, dataflow analysis, consistency checks, test data,verification of code, and debugging information. Analyzers can be broken into two major categories: dynamic and static. Static analyzers examine programs with respect to syntax errors and structural properties., This includes gathering statistical information on program content, such as the number of lines of executable code, source lines. and cyclomatic complexity. In addition, static analyzers provide the ability to check for the consistency of programs with respect to variables. Dynamic analyzers in contrast are dependent on input and the execution of a program providing the ability to find errors that cannot be detected through the use of static analysis alone. Dynamic analysis provides information on the behavior of a program rather than on the syntax. Both types of analysis detect errors in a program, but dynamic analyzers accomplish this through run-time behavior. This paper focuses on the following broad classification of dynamic analyzers: 1) Metrics; 2) Models; and 3) Monitors. Metrics are those analyzers that provide measurement. The next category, models, captures those analyzers that present the state of the program to the user at specified points in time. The last category, monitors, checks specified code based on some criteria. The paper discusses each classification and the techniques that are included under them. In addition, the role of each technique in the software life cycle is discussed. Familiarization with the tools that measure, model and monitor programs provides a framework for understanding the program's dynamic behavior from different, perspectives through analysis of the input/output data.

Lujan, Michelle R.↗

Variational Assimilation of Global Microwave Rainfall Retrievals: Physical and Dynamical Impact on GEOS Analyses and Forecasts

Global microwave rainfall retrievals from a 5-satellite constellation, including TMI from TRMM, SSWI from DMSP F13, F14 and F15, and AMSR-E from EOS-AQUA, are assimilated into the NASA Goddard Earth Observing System (GEOS) Data Assimilation System (DAS) using a 1-D variational continuous assimilation (VCA) algorithm. The physical and dynamical impact of rainfall assimilation on GEOS analyses and forecasts is examined at various temporal and spatial scales. This study demonstrates that the 1-D VCA algorithm, which was originally developed and evaluated for rainfall assimilations over tropical oceans, can effectively assimilate satellite microwave rainfall retrievals and improve GEOS analyses over both the Tropics and the extratropics where the atmospheric processes are dominated by different large-scale dynamics and moist physics, and also over the land, where rainfall estimates from passive microwave radiometers are believed to be less accurate. Results show that rainfall assimilation renders the GEOS analysis physically and dynamically more consistent with the observed precipitation at the monthly-mean and 6-hour time scales. Over regions where the model precipitation tends to misbehave in distinctly different rainy regimes, the 1-D VCA algorithm, by compensating for errors in the model s moist time-tendency in a 6-h analysis window, is able to bring the rainfall analysis closer to the observed. The radiation and cloud fields also tend to be in better agreement with independent satellite observations in the rainfall-assimilation m especially over regions where rainfall analyses indicate large improvements. Assimilation experiments with and without rainfall data for a midlatitude frontal system clearly indicates that the GEOS analysis is improved through changes in the thermodynamic and dynamic fields that respond to the rainfall assimilation. The synoptic structures of temperature, moisture, winds, divergence, and vertical motion, as well as vorticity are more realistically captured across the front. Short-term forecasts using initial conditions assimilated with rainfall data also show slight improvements. 1

Lin, Xin↗

Revealing nanoscale dynamics during an epoxy curing reaction with x-ray photon correlation spectroscopy

In this work the evolution of nanoscale properties is measured during the thermally triggered curing of an industrial epoxy adhesive. We use x-ray photon correlation spectroscopy (XPCS) to track the progression of the curing reaction through the local dynamics of filler particles that reflect the formation of a thermoset network. Out-of-equilibrium dynamics are resolved through identification and analysis of the intensity– intensity autocorrelation functions obtained from XPCS. The characteristic time scale and local velocity of the filler is calculated as functions of time and temperature. We find that the dynamics speed up when approaching the curing temperature (T cure ), and decay rapidly once T cure is reached. We compare the results from XPCS to conventional macroscale characterization by differential scanning calorimetry (DSC). The demonstration and implementation of nanoscale characterization of curing reactions by XPCS proves useful for future development and optimization of epoxy thermoset materials and other industrial adhesive systems.

36 MATERIALS SCIENCE↗

Synchronization of Human Autonomic Nervous System Rhythms with Geomagnetic Activity in Human Subjects

A coupling between geomagnetic activity and the human nervous system's function was identified by virtue of continuous monitoring of heart rate variability (HRV) and the time-varying geomagnetic field over a 31-day period in a group of 10 individuals who went about their normal day-to-day lives. A time series correlation analysis identified a response of the group's autonomic nervous systems to various dynamic changes in the solar, cosmic ray, and ambient magnetic field. Correlation coefficients and p values were calculated between the HRV variables and environmental measures during three distinct time periods of environmental activity. There were significant correlations between the group's HRV and solar wind speed, Kp, Ap, solar radio flux, cosmic ray counts, Schumann resonance power, and the total variations in the magnetic field. In addition, the time series data were time synchronized and normalized, after which all circadian rhythms were removed. It was found that the participants' HRV rhythms synchronized across the 31-day period at a period of approximately 2.5 days, even though all participants were in separate locations. Overall, this suggests that daily autonomic nervous system activity not only responds to changes in solar and geomagnetic activity, but is synchronized with the time-varying magnetic fields associated with geomagnetic field-line resonances and Schumann resonances.

Magnetic Fields↗

Space shuttle main engine numerical modeling code modifications and analysis

The user of computational fluid dynamics (CFD) codes must be concerned with the accuracy and efficiency of the codes if they are to be used for timely design and analysis of complicated three-dimensional fluid flow configurations. A brief discussion of how accuracy and efficiency effect the CFD solution process is given. A more detailed discussion of how efficiency can be enhanced by using a few Cray Research Inc. utilities to address vectorization is presented and these utilities are applied to a three-dimensional Navier-Stokes CFD code (INS3D).

Ziebarth, John P.↗

ThO 2 and Th 1– x U x O 2 Nanoscale Materials and Thin Films for Nuclear Science Applications

This study investigates the dynamics and mechanisms of solution combustion synthesis (SCS) for the preparation of nanoscale ThO 2 and Th 1–x U x O 2 materials, utilizing metal nitrates (Th(NO 3 ) 4 and UO 2 (NO 3 ) 2 ) and acetylacetone (C 5 H 8 O 2 ) as reactants dissolved in a 2-methoxyethanol (C 3 H 8 O 2 ) solvent. By combining thermodynamic calculations, dynamic time–temperature profile measurements with differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA), this research reveals how variations in acetylacetone concentration and uranium content influence the structural parameters of the synthesized oxides. The time–temperature measurements show that the heating rate and maximum combustion temperatures are sensitive to acetylacetone concentration. DSC-TGA results indicate shifts in exothermic peak temperatures as the uranium content changes. The complexation between thorium and acetylacetone emerges as a critical factor, impacting combustion parameters and the structural characteristics of the final products. The uniform distribution of Th and U in the Th 1–x U x O 2 solid solution and the formation of nanoscale particles with strained crystallites are considered essential for the low-temperature densification of these materials for nuclear fuel pellet applications. Additionally, high-quality ThO 2 and Th 1–x U x O 2 thin (100–150 nm) films are successfully synthesized via electrospray deposition of combustible solutions followed by a brief period of heat treatment. Furthermore, these films exhibit excellent structural and morphological uniformity, making them ideal candidates for nuclear measurements, irradiation damage studies, and investigations into the physical properties of both pure and mixed oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of non-Gaussian fluctuations in model A

Motivated by the experimental search for the QCD critical point, we perform simulations of a stochastic field theory with purely relaxational dynamics (model A). We verify the expected dynamic scaling of correlation functions. Using a finite size scaling analysis, we obtain the dynamic critical exponent z = 2.026(56). We investigate time dependent correlation functions of higher moments M n (t) of the order parameter M(t) for n = 1, 2, 3, 4. We obtain dynamic scaling with the same critical exponent z for all n, but the relaxation constant depends on n. We also study the relaxation of M n (t) after a quench, where the simulation is initialized in the high temperature phase, and the dynamics is studied at the critical temperature T c . Finally, we find that the evolution does not follow simple scaling with the dynamic exponent z, and that it involves an early time rise followed by late stage relaxation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗

Parallel simulation today

This paper surveys topics that presently define the state of the art in parallel simulation. Included in the tutorial are discussions on new protocols, mathematical performance analysis, time parallelism, hardware support for parallel simulation, load balancing algorithms, and dynamic memory management for optimistic synchronization.

Nicol, David↗

Cometary Evolution: Clues on Physical Properties from Chondritic Interplanetary Dust Particles

The degree of diversity or similarity detected in comets depends primarily on the lifetimes of the individual cometary nuclei at the time of analysis. It is inherent in our understanding of cometary orbital dynamics and the seminal model of comet origins by Oort that cometary evolution is the natural order of events in our Solar System. Thus, predictions of cometary behaviour in terms of bulk physical, mineralogical or chemical parameters should contain an appreciation of temporal variation(s). Previously, Rietmeijer and Mackinnon developed mineralogical bases for the chemical evolution of cometary nuclei primarily with regard to the predominantly silicate fraction of comet nuclei. We suggested that alteration of solids in cometary nuclei should be expected and that indications of likely reactants and products can be derived from judicious comparison with terrestrial diagenetic environments which include hydrocryogenic and low-temperature aqueous alterations. In a further development of this concept, Rietmeijer provides indirect evidence for the formation of sulfides and oxides in comet nuclei. Furthermore, Rietmeijer noted that timescales for hydrocryogenic and low-temperature reactions involving liquid water are probably adequate for relatively mature comets, e.g. P/comet Halley. In this paper, we will address the evolution of comet nuclei physical parameters such as solid particle grain size, porosity and density. In natural environments, chemical evolution (e.g. mineral reactions) is often accompanied by changes in physical properties. These concurrent changes are well-documented in the terrestrial geological literature, especially in studies of sediment diagenesis and we suggest that similar basic principles apply within the upper few meters of active comet nuclei. The database for prediction of comet nuclei physical parameters is, in principle, the same as used for the proposition of chemical evolution. We use detailed mineralogical studies of chondritic interplanetary dust particles (IDPS) as a guide to the likely constitution of mature comets traversing the inner Solar System. While there is, as yet, no direct proof that a specific sub-group or type of chondritic IDP is derived from a specific comet, it is clear that these particles are extraterrestrial in origin and that a certain portion of the interplanetary flux received by the Earth is cometary in origin. Two chondritic porous (CP) MPs, sample numbers W7010A2 and W7029Cl, from the Johnson Space Center Cosmic Dust Collection have been selected for this study of putative cometary physical parameters. This particular type of particle is considered a likely candidate for a cometary origin on the basis of mineralogy, bulk composition and morphology. While many IDPs have been subjected to intensive study over the past decade, we can develop a physical parameter model on only these two CP IDPs because few others have been studied in sufficient detail.

Rietmeijer, Frans J. M.↗

Cometary Evolution: Clues on Physical Properties from Chondritic Interplanetary Dust Particles

The degree of diversity or similarity detected in comets depends primarily on the lifetimes of the individual cometary nuclei at the time of analysis. It is inherent in our understanding of cometary orbital dynamics and the seminal model of comet origins that cometary evolution is the natural order of events in our Solar System. Thus, predictions of cometary behaviour in terms of bulk physical, mineralogical or chemical parameters should contain an appreciation of temporal variation(s). Previously, Rietmeijer and Mackinnon [1987] developed mineralogical bases for the chemical evolution of cometary nuclei primarily with regard to the predominantly silicate fraction of comet nuclei. We suggested that alteration of solids in cometary nuclei should be expected and that indications of likely reactants and products can be derived from judicious comparison with terrestrial diagenetic environments which include hydrocryogenic and low-temperature aqueous alterations. In a further development of this concept, Rietmeijer [1988] provides indirect evidence for the formation of sulfides and oxides in comet nuclei. Furthermore, Rietmeijer [1988] noted that timescales for hydrocryogenic and low-temperature reactions involving liquid water are probably adequate for relatively mature comets, e.g. P/comet Halley. In this paper, we will address the evolution of comet nuclei physical parameters such as solid particle grain size, porosity and density. In natural environments, chemical evolution (e.g. mineral reactions) is often accompanied by changes in physical properties. These concurrent changes are well-documented in the terrestrial geological literature, especially in studies of sediment diagenesis and we suggest that similar basic principles apply within the upper few meters of active comet nuclei. The database for prediction of comet nuclei physical parameters is, in principle, the same as used for the proposition of chemical evolution. We use detailed mineralogical studies of chondritic interplanetary dust particles (IDPS) as a guide to the likely constitution of mature comets traversing the inner Solar System. While there is, as yet, no direct proof that a specific sub-group or type of chondritic IDP is derived from a specific comet, it is clear that these particles are extraterrestrial in origin and that a certain portion of the interplanetary flux received by the Earth is cometary in origin. Two chondritic porous (CP) IDPS, sample numbers W701OA2 and W7029CI, from the Johnson Space Center Cosmic Dust Collection have been selected for this study of putative cometary physical parameters. This particular type of particle is considered a likely candidate for a cometary origin on the basis of mineralogy, bulk composition and morphology. While many IDPs have been subjected to intensive study over the past decade, we can develop a physical parameter model on only these two CP IDPs because few others have been studied in sufficient detail.

Reitmeijer, Frans J. M.↗