SEARCH · Engineering Papers
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.
Efforts to produce a Monte Carlo computer program for the analysis of crystal growth are briefly discussed. A literature survey was conducted of articles relating to the subject. A list of references reviewed is presented.
This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.
This paper develops a deterministic model of frequency distributions for energy imparted (total energy deposition) in small volumes similar to DNA molecules from high-energy ions of interest for space radiation protection and cancer therapy. Frequency distributions for energy imparted are useful for considering radiation quality and for modeling biological damage produced by ionizing radiation. For high-energy ions, secondary electron (delta-ray) tracks originating from a primary ion track make dominant contributions to energy deposition events in small volumes. Our method uses the distribution of electrons produced about an ion's path and incorporates results from Monte Carlo simulation of electron tracks to predict frequency distributions for ions, including their dependence on radial distance. The contribution from primary ion events is treated using an impact parameter formalism of spatially restricted linear energy transfer (LET) and energy-transfer straggling. We validate our model by comparing it directly to results from Monte Carlo simulations for proton and alpha-particle tracks. We show for the first time frequency distributions of energy imparted in DNA structures by several high-energy ions such as cosmic-ray iron ions. Our comparison with results from Monte Carlo simulations at low energies indicates the accuracy of the method.
Quantifying the uncertainty in model parameters for the purpose of damage prognosis can be accomplished utilizing Bayesian inference and damage diagnosis data from sources such as non-destructive evaluation or structural health monitoring. The number of samples required to solve the Bayesian inverse problem through common sampling techniques (e.g., Markov chain Monte Carlo) renders high-fidelity finite element-based damage growth models unusable due to prohibitive computation times. However, these types of models are often the only option when attempting to model complex damage growth in real-world structures. Here, a recently developed high-fidelity crack growth model is used which, when compared to finite element-based modeling, has demonstrated reductions in computation times of three orders of magnitude through the use of surrogate models and machine learning. The model is flexible in that only the expensive computation of the crack driving forces is replaced by the surrogate models, leaving the remaining parameters accessible for uncertainty quantification. A probabilistic prognosis framework incorporating this model is developed and demonstrated for non-planar crack growth in a modified, edge-notched, aluminum tensile specimen. Predictions of remaining useful life are made over time for five updates of the damage diagnosis data, and prognostic metrics are utilized to evaluate the performance of the prognostic framework. Challenges specific to the probabilistic prognosis of non-planar fatigue crack growth are highlighted and discussed in the context of the experimental results.
The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.
The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.
Aperture-synthesis and single-dish (sub-) millimeter molecular-line and continuum observations reveal in great detail the envelope structure of deeply embedded young stellar objects (SMM 1 = FIRS 1, SMM 2, SMM 3, SMM 4) in the densely star-forming Serpens Molecular Cloud. SMM 1, 3, and 4 show partially resolved (>2" = 800 AU) continuum emission in the beam of the Owens Valley Millimeter Array at lambda = 3.4-1.4 mm. The continuum visibilities accurately constrain the density structure in the envelopes, which can be described by a radial power law with slope -2.0 +/- 0.5 on scales of 300 to 8000 AU. Inferred envelope masses within a radius of 8000 AU are 8.7, 3.0, and 5.3 Msolar for SMM 1, 3, and 4, respectively. A point source with 20%-30% of the total flux at 1.1 mm is required to fit the observations on long baselines, corresponding to warm envelope material within approximately 100 AU or a circumstellar disk. No continuum emission is detected interferometrically toward SMM 2, corresponding to an upper limit of 0.2 Msolar assuming Td = 24 K. The lack of any compact dust emission suggests that the SMM 2 core does not contain a central protostar. Aperture-synthesis observations of the 13CO, C18O, HCO+, H13CO+, HCN, H13CN, N2H+ 1-0, SiO 2-1, and SO 2(2)-1(1) transitions reveal compact emission toward SMM 1, 3, and 4. SMM 2 shows only a number of clumps scattered throughout the primary field of view, supporting the conclusion that this core does not contain a central star. The compact molecular emission around SMM 1, 3, and 4 traces 5"-10" (2000-4000 AU) diameter cores that correspond to the densest regions of the envelopes, as well as material directly associated with the molecular outflow. Especially prominent are the optically thick HCN and HCO+ lines that show up brightly along the walls of the outflow cavities. SO and SiO trace shocked material, where their abundances may be enhanced by 1-2 orders of magnitude over dark-cloud values. A total of 31 molecular transitions have been observed with the James Clerk Maxwell and Caltech Submillimeter telescopes in the 230, 345, 490, and 690 GHz atmospheric windows toward all four sources, containing, among others, lines of CO, HCO+, HCN, H2CO, SiO, SO, and their isotopomers. These lines show 20-30 km s-1 wide line wings, deep and narrow (1-2 km s-1) self-absorption, and 2-3 km s-1 FWHM line cores. The presence of highly excited lines like 12CO 4-3 and 6-5, 13CO 6-5, and several H2CO transitions indicates the presence of material with temperatures > or approximately 100 K. Monte Carlo calculations of the molecular excitation and line transfer show that the envelope model derived from the dust emission can successfully reproduce the observed line intensities. The depletion of CO in the cold gas is modest compared to values inferred in objects like NGC 1333 IRAS 4, suggesting that the phase of large depletions through the entire envelope is short lived and may be influenced by the local star formation density. Emission in high-excitation lines of CO and H2CO requires the presence of a small amount of approximately 100 K material, comprising less than 1% of the total envelope mass and probably associated with the outflow or the innermost region of the envelope. The derived molecular abundances in the warm (Tkin > 20 K) envelope are similar to those found toward other class 0 YSOs like IRAS 16293-2422, though some species appear enhanced toward SMM 1. Taken together, the presented observations and analysis provide the first comprehensive view of the physical and chemical structure of the envelopes of deeply embedded young stellar objects in a clustered environment on scales between 1000 and 10,000 AU.
Logical flow and guidelines are provided for the construction of a low thrust orbit determination computer program. The program, tentatively called FRACAS (filter response analysis for continuously accelerating spacecraft), is capable of generating a reference low thrust trajectory, performing a linear covariance analysis of guidance and navigation processes, and analyzing trajectory nonlinearities in Monte Carlo fashion. The choice of trajectory, guidance and navigation models has been made after extensive literature surveys and investigation of previous software. A key part of program design relied upon experience gained in developing and using Martin Marietta Aerospace programs: TOPSEP (Targeting/Optimization for Solar Electric Propulsion), GODSEP (Guidance and Orbit Determination for SEP) and SIMSEP (Simulation of SEP).
An optical model is developed for imaging optical radars such as the SAR on Seasat and the Shuttle Imaging Radar (SIR-B) by analyzing the Doppler shift of individual speckles in the image. The signal received at the spacecraft is treated in terms of a Fresnel-Kirchhoff integration over all backscattered radiation within a Huygen aperture at the earth. Account is taken of the movement of the spacecraft along the orbital path between emission and reception. The individual points are described by integration of the point source amplitude with a Green's function scattering kernel. Doppler data at each point furnishes the coordinates for visual representations. A Rayleigh-Poisson model of the surface scattering characteristics is used with Monte Carlo methods to generate simulations of Doppler radar speckle that compare well with Seasat SAR data SIR-B data.
The lattice-regularized U(1) Higgs model in an external electromagnetic field is studied by Monte Carlo techniques. In the Coulomb phase, magnetic flux can flow through uniformly. The Higgs phase splits into a region where magnetic flux can penetrate only in the form of vortices and a region where the magnetic flux is completely expelled, the relativistic analog of the Meissner effect in superconductivity. Evidence is presented for symmetry restoration in strong external fields.
The objective of the research is to use the low-frequency wave spectrum measured by the Plasma Wave Instrument on the DE-1 spacecraft to include the wave-particle interaction (WPI) in the different polar wind models that are available at Utah State University. A Monte Carlo technique was used to simulate the ion diffusion in the velocity space due to scattering by the waves. This enabled us to study the effect of WPI on the magnitude as well as the composition of the outflow of the ionospheric ions. In particular, in the first year the 1-D semi-kinetic code developed by Barakat and Schunk (1983) was modified in order to include the effect of the WPI.
The project BEST has been proposed by several laboratories of the French Space Agency (CNES) as part of the GEWEX program. Its objective is the study of the energy budget of the tropical system, and aims at measuring different components of the water cycle by coupling several observing systems. As part of this, Service d'Aeronomie (SA) proposed to study a Differential Absorption Lidar (DIAL) system to measure tropical water vapor profiles and to determine humidity convergence fields, when coupled with the Doppler wind lidar. We have performed simulation analyses for water vapor measurement from the selected orbit of 400 km. For these simulations, we have adopted an atmospheric model developed at SA. We have optimized, by a Monte Carlo statistical method, the system parameters and evaluated bias and standard deviation of the retrieved water vapor mixing ratio.
The detailed flow structure of a normal shock wave for a gas mixture is investigated using the direct-simulation Monte Carlo method. A variable diameter hard-sphere (VDHS) model is employed to investigate the effect of different viscosity temperature exponents (VTE) for each species in a gas mixture. Special attention is paid to the irregular behavior in the density profiles which was previously observed in a helium-xenon experiment. It is shown that the VTE can have substantial effects in the prediction of the structure of shock waves. The variable hard-sphere model of Bird shows good agreement, but with some limitations, with the experimental data if a common VTE is chosen properly for each case. The VDHS model shows better agreement with the experimental data without adjusting the VTE. The irregular behavior of the light-gas component in shock waves of disparate mass mixtures is observed not only in the density profile, but also in the parallel temperature profile. The strength of the shock wave, the type of molecular interactions, and the mole fraction of heavy species have substantial effects on the existence and structure of the irregularities.
A summary of three-dimensional simulations on the hypersonic rarefied flows in an effort to understand the highly nonequilibrium flows about space vehicles entering the Earth's atmosphere for a realistic estimation of the aerothermal loads is presented. Calculations are performed using the direct simulation Monte Carlo method with a five-species reacting gas model, which accounts for rotational and vibrational internal energies. Results are obtained for the external flows about various bodies in the transitional flow regime. For the cases considered, convective heating, flowfield structure and overall aerodynamic coefficients are presented and comparisons are made with the available experimental data. The agreement between the calculated and measured results are very good.
Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration-dissociation recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.
Rarefied flows about hypersonic vehicles entering the upper atmosphere or through nozzles expanding into a near vacuum may only be simulated accurately with a direct simulation Monte Carlo (DSMC) method. Under this grant, researchers enhanced the models employed in the DSMC method and performed simulations in support of existing NASA projects or missions. DSMC models were developed and validated for simulating rotational, vibrational, and chemical relaxation in high-temperature flows, including effects of quantized anharmonic oscillators and temperature-dependent relaxation rates. State-of-the-art advancements were made in simulating coupled vibration - dissociation - recombination for post-shock flows. Models were also developed to compute vehicle surface temperatures directly in the code rather than requiring isothermal estimates. These codes were instrumental in simulating aerobraking of NASA's Magellan spacecraft during orbital maneuvers to assess heat transfer and aerodynamic properties of the delicate satellite. NASA also depended upon simulations of entry of the Galileo probe into the atmosphere of Jupiter to provide drag and flow field information essential for accurate interpretation of an onboard experiment. Finally, the codes have been used extensively to simulate expanding nozzle flows in low-power thrusters in support of propulsion activities at NASA-Lewis. Detailed comparisons between continuum calculations and DSMC results helped to quantify the limitations of continuum CFD codes in rarefied applications.
The accurate pixel-by-pixel retrieval of cloud optical properties from space is influenced by radiative smoothing due to high order photon scattering and radiative roughening due to low order scattering events. Both are caused by cloud heterogeneity and the three-dimensional nature of radiative transfer and can be studied with the aid of computer simulations. We use Monte Carlo simulations on variable 1-D and 2-D model cloud fields to seek for dependencies of smoothing and roughening phenomena on single scattering albedo, solar zenith angle, and cloud characteristics. The results are discussed in the context of high resolution satellite (such as Landsat) retrieval applications. The current work extends the investigation on the inverse NIPA (Non-local Independent Pixel Approximation) as a tool for removing smoothing and improving retrievals of cloud optical depth. This is accomplished by: (1) Delineating the limits of NIPA applicability; (2) Exploring NIPA parameter dependences on cloud macrostructural features, such as mean cloud optical depth and geometrical thickness, degree of extinction and cloud top height variability. We also compare parameter values from empirical and theoretical considerations; (3) Examining the differences between applying NIPA on radiation quantities vs direct application on optical properties; (4) Studying the radiation budget importance of the NIPA corrections as a function of scale. Finally, we discuss fundamental adjustments that need to be considered for successful radiance inversion at non-conservative wavelengths and oblique Sun angles. These adjustments are necessary to remove roughening signatures which become more prominent with increasing absorption and solar zenith angle.
A multi-body flight simulation for the Phoenix Mars Lander has been developed that includes high fidelity six degree-of-freedom rigid-body models for the parachute and lander system. The simulation provides attitude and rate history predictions of all bodies throughout the flight, as well as loads on each of the connecting lines. In so doing, a realistic behavior of the descending parachute/lander system dynamics can be simulated that allows assessment of the Phoenix descent performance and identification of potential sensitivities for landing. This simulation provides a complete end-to-end capability of modeling the entire entry, descent, and landing sequence for the mission. Time histories of the parachute and lander aerodynamic angles are presented. The response of the lander system to various wind models and wind shears is shown to be acceptable. Monte Carlo simulation results are also presented.