Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Design Extrapolation”

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 451 records · Page 25

Metamodeling Time-Temperature Creep Parameters

Abstract There exist many time-temperature parameter (TTP) models for creep rupture prediction of components including the Larson–Miller (LM), Manson–Haferd (MH), Manson–Brown (MB), Orr–Sherby–Dorn (OSD), Manson–Succop (MS), Graham–Walles (GW), Chitty–Duval (CD), Goldhoff–Sherby (GS) models. It remains a challenge to determine which model is “best”, capable of accurate interpolation and physically realistic extrapolation of creep rupture data for a given material. In this study, metamodeling is applied to create a unified TTP metamodel that combines and regresses into twelve TTP models (eight existing and four newly derived). An analysis of the mathematical problems that exist in TTP models is provided. A matlab code is written that can: (1) calibrate the material constants of any of the twelve TTP models (using the metamodel); (2) determine the most suitable stress-parameter function; (3) and report the normalized mean square error (NMSE) of rupture predictions for a given material database. Using the metamodel, and code, a design engineer can make an intelligent selection of the “best” TTP model for creep resistant design. This process is demonstrated using four isotherms of alloy P91 creep rupture data. To assess the influence of material, further validation is performed on alloys Hastelloy X, 304SS, and 316SS. It is determined that the “best” model is dependent on material type and the quality and quantity of available data.

Engineering↗

Constellation Probabilistic Risk Assessment (PRA): Design Consideration for the Crew Exploration Vehicle

Managed by NASA's Office of Safety and Mission Assurance, a pilot probabilistic risk analysis (PRA) of the NASA Crew Exploration Vehicle (CEV) was performed in early 2006. The PRA methods used follow the general guidance provided in the NASA PRA Procedures Guide for NASA Managers and Practitioners'. Phased-mission based event trees and fault trees are used to model a lunar sortie mission of the CEV - involving the following phases: launch of a cargo vessel and a crew vessel; rendezvous of these two vessels in low Earth orbit; transit to th$: moon; lunar surface activities; ascension &om the lunar surface; and return to Earth. The analysis is based upon assumptions, preliminary system diagrams, and failure data that may involve large uncertainties or may lack formal validation. Furthermore, some of the data used were based upon expert judgment or extrapolated from similar components~systemsT. his paper includes a discussion of the system-level models and provides an overview of the analysis results used to identify insights into CEV risk drivers, and trade and sensitivity studies. Lastly, the PRA model was used to determine changes in risk as the system configurations or key parameters are modified.

Prassinos, Peter G.↗

Effects of land-based wind turbine upsizing on community sound levels and power and energy density

Multiple technological, social, and market factors of wind power are evolving rapidly. Most notably, significant wind turbine scaling is occurring and is forecasted to continue. While the larger turbines expected to be deployed in the future are more powerful and efficient, they are also expected to operate at higher sound levels and require larger setbacks than those installed in the last decade. These sometimes-competing deployment trends and impacts cannot be understood via simple extrapolations of past trends. This study analyzes the effect of these future larger turbines on wind turbine micro-siting, project-level power and energy density, and community noise impacts. Due to their taller heights, larger rotors, and higher sound power levels, future wind turbines will require larger setbacks from homes and greater inter-turbine spacing, resulting in fewer turbines deployed for a given land area. This research finds these changes more than offset the effect of the higher turbine sound emissions, significantly decreasing the average sound levels that wind plant hosting communities experience. Yet, simultaneously, plant layouts using future turbine designs also result in projects with higher installed capacities and annual energy output for a given land area. These increases will likely lead to increased tax benefits and local income in the community. The deployment of fewer turbines on a smaller number of parcels could have implications on siting flexibility and landowner payments.

17 WIND ENERGY↗

Charged-particle spectra from μ – capture on Al

Published data on the emission of charged particles following nuclear muon capture are extremely limited. In addition to its interest as a probe of the nuclear response, these data are important for the design of some current searches for lepton flavor violation. This work presents momentum spectra of protons and deuterons following μ- capture in aluminum. It is the first measurement of a muon capture process performed with a tracking spectrometer. A precision of better than 10% over the momentum range of 100–190MeV/c for protons is obtained; for deuterons of 145–250MeV/c the precision is better than 20%. The observed partial yield of protons with emission momenta above 80MeV/c (kinetic energy 3.4MeV) is 0.0322 ± 0.0007(stat) ± 0.0022(syst) per capture, and for deuterons above 130MeV/c (4.5MeV) it is 0.0122 ± 0.0009(stat) ± 0.0006(syst). Extrapolating to total yields gives 0.045 ± 0.001(stat) ± 0.003(syst) ± 0.001(extrapolation) per capture for protons and 0.018 ± 0.001(stat) ± 0.001(syst)± 0.002(extrapolation) for deuterons, which are the most precise measurements of these quantities to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrated Application of Active Controls (IAAC) technology to an advanced subsonic transport project. ACT/Control/Guidance System study, volume 1

The active control technology (ACT) control/guidance system task of the integrated application of active controls (IAAC) technology project within the NASA energy efficient transport program was documented. The air traffic environment of navigation and air traffic control systems and procedures were extrapolated. An approach to listing flight functions which will be performed by systems and crew of an ACT configured airplane of the 1990s, and a determination of function criticalities to safety of flight, are the basis of candidate integrated ACT/Control/Guidance System architecture. The system mechanizes five active control functions: pitch augmented stability, angle of attack limiting, lateral/directional augmented stability, gust load alleviation, and maneuver load control. The scope and requirements of a program for simulating the integrated ACT avionics and flight deck system, with pilot in the loop, are defined, system and crew interface elements are simulated, and mechanization is recommended. Relationships between system design and crew roles and procedures are evaluated.

Source record↗

Performance of Several Method-of-Characteristics Exhaust Nozzles

Nozzle performance data were obtained with three "method-of-characteristics" nozzles and a 150 conical nozzle at pressure ratios up to 130. Each basic configuration was cut off and tested at expansion ratios of 25, 20, 15, and 10. Unheated dry air was used at nozzle inlet pressures up to 22,000 pounds per square foot absolute. Nozzle thrust data were extrapolated to infinite pressure ratio (zero discharge pressure). As much as 1-percent increase in thrust with no increase in nozzle surface area (weight), can be obtained by using a method-of-characteristics, nozzle instead of a 15 conical nozzle when operating with a nozzle expansion ratio of 25 and nozzle pressure ratios from 200 to infinity. Conversely, for the same thrust, reductions in nozzle divergent surface area in the order of 25 percent are possible. The thrust performance of the method-of-characteristics nozzle was not as good as that of the 150 conical nozzle when operating at pressure ratios considerably below design (below 100 for the expansion ratio 25 nozzles). Theoretical and measured nozzle momentum coefficients agreed within about 0.6 percent. This is the order of accuracy of both the measured and theoretical values.

Farley, John M.↗

BioSentinel: Mission Development of a Radiation Biosensor to Gauge DNA Damage and Repair Beyond Low Earth Orbit on a 6U Nanosatellite.

We are designing and developing a "6U" (10 x 22 x 34 cm; 14 kg) nanosatellite as a secondary payload to fly aboard NASA's Space Launch System (SLS) Exploration Mission (EM) 1, scheduled for launch in late 2017. For the first time in over forty years, direct experimental data from biological studies beyond low Earth orbit (LEO) will be obtained during BioSentinel's 12- to 18- month mission. BioSentinel will measure the damage and repair of DNA in a biological organism and allow us to compare that to information from onboard physical radiation sensors. In order to understand the relative contributions of the space environment's two dominant biological perturbations, reduced gravity and ionizing radiation, results from deep space will be directly compared to data obtained in LEO (on ISS) and on Earth. These data points will be available for validation of existing biological radiation damage and repair models, and for extrapolation to humans, to assist in mitigating risks during future long-term exploration missions beyond LEO. The BioSentinel Payload occupies 4U of the spacecraft and will utilize the monocellular eukaryotic organism Saccharomyces cerevisiae (yeast) to report DNA double-strand-break (DSB) events that result from ambient space radiation. DSB repair exhibits striking conservation of repair proteins from yeast to humans. Yeast was selected because of 1) its similarity to cells in higher organisms, 2) the well-established history of strains engineered to measure DSB repair, 3) its spaceflight heritage, and 4) the wealth of available ground and flight reference data. The S. cerevisiae flight strain will include engineered genetic defects to prevent growth and division until a radiation-induced DSB activates the yeast's DNA repair mechanisms. The triggered culture growth and metabolic activity directly indicate a DSB and its successful repair. The yeast will be carried in the dry state within the 1-atm P/L container in 18 separate fluidics cards with each card having 16 independent culture microwells, with integral microchannels and filters to supply nutrients and reagents, confine the yeast to the wells, and enable optical measurement. The measurement subsystem will monitor each subgroup of culture wells continuously for several weeks, optically tracking DSBtriggered cell growth and metabolism. BioSentinel will also include physical radiation sensors based on the TimePix sensor, as implemented by JSC's RadWorks group, which record individual radiation events including estimates of their linear-energytransfer (LET) values. Radiation-dose and LET data will be compared directly to the rate of DSB-and-repair events measured by the S. cerevisiae biosentinels. The spacecraft bus will operate in a deep space environment with functions that include command and data handling, communications, power generation (via deployable solar panels) and storage, and attitude determination-and-control system with micropropulsion. Development of the BioSentinel spacecraft will mature and prove multiple nanosatellite advances in order to function well beyond LEO: Communications from distances of ≥ 500,000 km; Autonomous attitude control, momentum management, and safe mode of nanosatellites in deep space; Shielding-, hardening-, design-, and software-derived radiation tolerance for electronics; Reliable functionality for 12 - 18 months of key subsystems for biofluidics, memory, communications, power, etc.; Close integration of living biological radiation event monitors with miniature physical radiation spectrometers; Biological measurement of solar particle events beyond Earth orbit In addition to providing the first biological results from beyond LEO in over 4 decades, BioSentinel will provide an adaptable small-satellite instrument platform to perform a range of human-exploration-relevant measurements that characterize the biological consequences of multiple outer space environments. BioSentinel is being developed under NASA's Advanced Exploration Systems program.

DNA damage↗

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗

Ultraviolet interstellar linear polarization. 2: The wavelength dependence

We present new ultraviolet (UV) polarimetry of the well-studied interstellar line of sight toward HD 204827 obtained with the Faint Object Spectrograph on Hubble Space Telescope. HD 204827 is of great interest because the dust along this line of sight has extremely low values of both lambda(sub max) and R(sub V). Its far-UV extinction is very large, reflecting its small R(sub V) value. In addition, we reexamine the entire sample of 14 interstellar lines of sight for which there are now UV polarization data. We find that the previously suggested relationship between lambda(sub max) and the wavelength dependence of the polarization in the UV is strongly supported by the data for this larger sample including HD 204827. Seven stars with lambda(sub max) greater than or equal to 0.54 micron agree well with an extrapolation of the Serkowski relation into the UV while seven stars with lambda(sub max) less than or equal to 0.53 micron show polarization in excess of the Serkowski extrapolation (super-Serkowski). However, the division of the observed lines of sight into Serkowski and super-Serkowski categories is artificial. In fact, the amount of polarization in the UV is correlated with a single parameter, lambda(sub max). This may indicate that there is a mean interstellar polarization law analogous to the mean interstellar extinction law of Cardelli, Calyton, & Mathis which is based on R(sub V). The data are consistent with a linear relationship between 1/lambda(sub max) and rho(UV)/rho(sub max) but more data are needed to define the functional form. We suggest that the Serkowski and super-Serkowski designations be replaced by high and low lambda(sub max) which are more physically descriptive. At the same time, we note that all seven super-Serkowski (low lambda(sub max)) stars lie in a relatively small region of the sky between l(sup II) = 90 deg - 150 deg and b = -5 deg - 15 deg. These stars all lie in or behind a spur of the local Orion spiral arm. Similarly, most of the Serkowski (high lambda(sub max) stars lie in or near the Scorpio-Centaurus OB Association. So lines of sight covering larger areas of the sky are needed to test the universality of the lambda(sub max)/UV polarization relationship. The recent discovery of warm dust near HD 197770 suggests the possibility that a mechanism other than the traditional alignment to the Galactic magnetic field may be invoked to explain its 2175 A polarization bump.

Clayton, Geoffrey C.↗

Industry Facing PV Degradation Prediction Tool and Database to Enable a 50-Year Life Module

The Photovoltaic (PV) industry constantly aims for lower costs, higher-efficiency cells, and improved module designs. These trends lead to using new materials, designs, and manufacturing processes, resulting in a continually changing technological landscape. These changes can potentially introduce new, unknown degradation mechanisms and failure modes that are difficult to diagnose, analyze, test, and model. This introduces uncertainty into the expected lifetime of PV modules of 25 years. Furthermore, research efforts aim for up to 50 years of service life while keeping performance degradation at a minimum for decades of outdoor weathering - putting additional pressure on improving the accuracy of long-term durability and reliability assessments. Here there is a need to organize the existing degradation data into an accessible format and to provide industry relevant tools for extrapolation from laboratory to field conditions. Because the core of this type of analysis involves calculations that are complicated but ubiquitous for many degradation processes, an enhanced predictive modeling framework will facilitate the analysis to help researchers keep up with the rapid pace of technological changes. In this work, we present an online tool that can be used to search for and analyze degradation information and extrapolate PV module performance and durability to field exposure. A graphical user interface will aid in the understanding of the results. The prediction tool will be built modular and published as open source, enabling users to expand on the existing framework. We use an integration pipeline approach that allows us to leverage weather data from the National Solar Radiation Database to perform geospatial degradation analysis in the US and worldwide. Our repository will contain various degradation models and material parameters suitable for the reliability and durability assessment of materials and components deployed outdoors. We hope to become a repository that can be used for weathering and degradation analysis for various applications beyond the PV industry.

degradation↗

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗

Metrology of Large Parts

As discussed in the first chapter of this book, there are many different methods to measure a part using optical technology. Chapter 2 discussed the use of machine vision to measure macroscopic features such as length and position, which was extended to the use of interferometry as a linear measurement tool in chapter 3, and laser or other trackers to find the relation of key points on large parts in chapter 4. This chapter looks at measuring large parts to optical tolerances in the sub-micron range using interferometry, ranging, and optical tools discussed in the previous chapters. The purpose of this chapter is not to discuss specific metrology tools (such as interferometers or gauges), but to describe a systems engineering approach to testing large parts. Issues such as material warpage and temperature drifts that may be insignificant when measuring a part to micron levels under a microscope, as will be discussed in later chapters, can prove to be very important when making the same measurement over a larger part. In this chapter, we will define a set of guiding principles for successfully overcoming these challenges and illustrate the application of these principles with real world examples. While these examples are drawn from specific large optical testing applications, they inform the problems associated with testing any large part to optical tolerances. Manufacturing today relies on micrometer level part performance. Fields such as energy and transportation are demanding higher tolerances to provide increased efficiencies and fuel savings. By looking at how the optics industry approaches sub-micrometer metrology, one can gain a better understanding of the metrology challenges for any larger part specified to micrometer tolerances. Testing large parts, whether optical components or precision structures, to optical tolerances is just like testing small parts, only harder. Identical with what one does for small parts, a metrologist tests large parts and optics in particular to quantify their mechanical properties (such as dimensions, mass, etc); their optical prescription or design (i.e. radius of curvature, conic constant, vertex location, size); and their full part shape. And, just as with small parts, a metrologist accomplishes these tests using distance measuring instruments such as tape measures, inside micrometers, coordinate measuring machines, distance measuring interferometers; angle measuring instruments such as theodolites, autocollimators; and surface measuring instruments including interferometers, stylus profilers, interference microscopes, photogrammetric cameras, or other tools. However, while the methodology may be similar, it is more difficult to test a large object for the simple reason that most metrologists do not have the necessary intuition. The skills used to test small parts or optics in a laboratory do not extrapolate to testing large parts in an industrial setting any more than a backyard gardener might successfully operate a farm. But first, what is a large part? A simple definition might be the part's size or diameter. For optics and diffuse surface parts alike, the driving constraint is ability to illuminate the part's surface. For reflective convex mirrors, large is typically anything greater than 1 meter. But, for refractive optics, flats or convex mirrors, large is typically greater than 0.5 meter. While a size definition is simple, it may be less than universal. A more nuanced definition might be that a large part is any component which cannot be easily tested in a standard laboratory environment, on a standard vibration isolated table using standard laboratory infrastructure. A micro-switch or a precision lens might be easily measured to nanometer levels under a microscope in a lab, but a power turbine spline or a larger telescope mirror will not fit under that microscope and may not even fit on the table.

Stahl, H. Philip↗

Design Study for the Asteroid Redirect Vehicle (ARV) Composite Primary Bulkhead

A design study was undertaken of a carbon fiber primary bulkhead for a large solar electric propulsion (SEP) spacecraft. The bulkhead design, supporting up to 16 t of xenon propellant, progressed from one consisting of many simple parts with many complex joints, to one consisting of a few complex parts with a few simple joints. The unique capabilities of composites led to a topology that transitioned loads from bending to in-plane tension and shear, with low part count. This significantly improved bulkhead manufacturability, cost, and mass. The stiffness-driven structure utilized high-modulus M55J fiber unidirectional prepregs. A full-scale engineering demonstration unit (EDU) of the concept was used to demonstrate manufacturability of the concept. Actual labor data was obtained, which could be extrapolated to a full bulkhead. The effort demonstrated the practicality of using high-modulus fiber (HMF) composites for unique shape topologies that minimize mass and cost. The lessons are applicable to primary and secondary aerospace structures that are stiffness driven.

Cressman, Thomas O.↗

The ESA Hera Mission: Detailed Investigation of the NASA Dart Impact Outcome and Characterization of the Binary Asteroid Didymos

The Hera mission is in development in the Space Safety Program of the European Space Agency (ESA) in collaboration with JAXA for launch in October 2024. It will rendezvous with the binary asteroid (65803) Didymos in early 2027 and over the following 6 months will investigate the properties of Didymos and its small moon Dirmorphos, including the effects of the successful impact of the NASA DART probe. DART was launched on 24 November 2021 at 06:21 UTC and successfully impacted at about 6.1 km/s on the 160 meter-size Dimorphos on 26 September 2022 at 23:14 UTC. Two weeks before impact, it deployed the Light Italian Cubesat for Imaging of Asteroids (LICIACube) that provided images during the few minutes following the impact. A campaign of observations from Earth and space provided distant images of the event, which were used to determine a reduction of 33 min in the orbital period of Dimorphos around Didymos. The Hera mission will investigate the outcome of the DART impact in detail. In particular, it will determine the effects of the impact on Dimorphos’ surface, e.g., whether it produced a crater, and if so its properties, and/or whether it led to global reshaping of Dimorphos. It will accurately measure the mass of Dimorphos, which will allow quantification of the momentum enhancement factor, which is a crucial parameter to evaluate the efficiency of the kinetic impactor technique. Hera will also determine in detail the physical, thermal and compositional properties of the asteroid, including first time measurements of the internal properties, which have a critical influence on the impact outcome. With Hera and DART together, we will have a fully documented deflection test and impact experiment that will allow us to verify our numerical impact models, to be able to extrapolate this deflection test to other asteroids and to improve our understanding of the complex geological properties and response of small asteroids. Hera consists of a mother spacecraft and two cubesats. It will perform the first rendezvous with a binary asteroid, arriving at Didymos about 4 years after DART. Hera will not only contribute significantly to the planetary defense effort by documenting entirely the outcome of the DART impact, it will also greatly advance our understanding of processes in the very low gravity regime of a small asteroidv (Dimorphos) as well as our knowledge of binary asteroids and asteroid geophysics. The mission development is progressing nominally, with the Critical Design Review having been passed in Nov. 2022. The various working groups of the Hera Science Team are working intensively to support this development and be best prepared for the surprises that Hera will offer during its visit to Didymos in 2027.

Asteroid deflection↗

Plenum response to simulated disturbances of the model and fan inlet guide vanes in a transonic tunnel

In order to aid in the design of the National Transonic Facility (NTF) control system, test section/plenum response studies were carried out in a 0.186 scale model of the NTF high speed duct. Two types of disturbances, those induced by the model and those induced by the compressor inlet guide vanes were simulated. Some observations with regard to the test section/plenum response tests are summarized as follows. A resonance frequency for the test section/plenum area of the tunnel of approximately 50 Hz was observed for Mach numbers from 0.40 to 0.90. However, since the plenum is 3.1 times (based on volume) too large for the scaled size of the test section, care must be taken in extrapolating these data to NTF conditions. The plenum pressure data indicate the existence of pressure gradients in the plenum. The test results indicate that the difference between test section static pressure and plenum pressure is dependent on test section flow conditions. Plenum response to inlet guide vane type disturbances appears to be slower than plenum response to test section disturbances.

Gloss, B. B.↗

Evaluation of dual-weighted residual and machine learning error estimation for projection-based reduced-order models of steady partial differential equations

Projection-based reduced-order models (pROMs) show great promise as a means to accelerate many-query applications such as forward error propagation, solving inverse problems, and design optimization. In order to deploy pROMs in the context of high-consequence decision making, accurate error estimates are required to determine the region(s) of applicability in the parameter space. The following paper considers the dual-weighted residual (DWR) error estimate for pROMs and compares it to another promising pROM error estimate, machine learned error models (MLEM). Here, we show how DWR can be applied to ROMs and then evaluate DWR on two partial differential equations (PDEs): a two-dimensional linear convection–reaction–diffusion equation, and a three-dimensional static hyper-elastic beam. It is shown that DWR is able to estimate errors for pROMs extrapolating outside of their training set while MLEM is best suited for pROMs used to interpolate within the pROM training set.

42 ENGINEERING↗

Perfectly Matched Layers and Characteristic Boundaries in Lattice Boltzmann: Accuracy vs Cost

Artificial boundary conditions (BCs) play a ubiquitous role in numerical simulations of transport phenomena in several diverse fields, such as fluid dynamics, electromagnetism, acoustics, geophysics, and many more. They are essential for accurately capturing the behavior of physical systems whenever the simulation domain is truncated for computational efficiency purposes. Ideally, an artificial BC would allow relevant information to enter or leave the computational domain without introducing artifacts or unphysical effects. Boundary conditions designed to control spurious wave reflections are referred to as nonreflective boundary conditions (NRBCs). Another approach is given by the perfectly matched layers (PMLs), in which the computational domain is extended with multiple dampening layers, where outgoing waves are absorbed exponentially in time. Here, in this work, the definition of PML is revised in the context of the lattice Boltzmann method. The impact of adopting different types of BCs at the edge of the dampening zone is evaluated and compared, in terms of both accuracy and computational costs. It is shown that for sufficiently large buffer zones, PMLs allow stable and accurate simulations even when using a simple zeroth-order extrapolation BC. Moreover, employing PMLs in combination with NRBCs potentially offers significant gains in accuracy at a modest computational overhead, provided the parameters of the BC are properly tuned to match the properties of the underlying fluid flow.

97 MATHEMATICS AND COMPUTING↗

Depth Dose Distribution Study within a Phantom Torso after Irradiation with a Simulated Solar Particle Event at NSRL

The adequate knowledge of the radiation environment and the doses incurred during a space mission is essential for estimating an astronaut's health risk. The space radiation environment is complex and variable, and exposures inside the spacecraft and the astronaut's body are compounded by the interactions of the primary particles with the atoms of the structural materials and with the body itself Astronauts' radiation exposures are measured by means of personal dosimetry, but there remains substantial uncertainty associated with the computational extrapolation of skin dose to organ dose, which can lead to over- or underestimation of the health risk. Comparisons of models to data showed that the astronaut's Effective dose (E) can be predicted to within about a +10% accuracy using space radiation transport models for galactic cosmic rays (GCR) and trapped radiation behind shielding. However for solar particle event (SPE) with steep energy spectra and for extra-vehicular activities on the surface of the moon where only tissue shielding is present, transport models predict that there are large differences in model assumptions in projecting organ doses. Therefore experimental verification of SPE induced organ doses may be crucial for the design of lunar missions. In the research experiment "Depth dose distribution study within a phantom torso" at the NASA Space Radiation Laboratory (NSRL) at BNL, Brookhaven, USA the large 1972 SPE spectrum was simulated using seven different proton energies from 50 up to 450 MeV. A phantom torso constructed of natural bones and realistic distributions of human tissue equivalent materials, which is comparable to the torso of the MATROSHKA phantom currently on the ISS, was equipped with a comprehensive set of thermoluminescence detectors and human cells. The detectors are applied to assess the depth dose distribution and radiation transport codes (e.g. GEANT4) are used to assess the radiation field and interactions of the radiation field with the phantom torso. Lymphocyte cells are strategically embedded at selected locations at the skin and internal organs and are processed after irradiation to assess the effects of shielding on the yield of chromosome damage. The initial focus of the present experiment is to correlate biological results with physical dosimetry measurements in the phantom torso. Further on, the results of the passive dosimetry within the anthropomorphic phantoms represent the best tool to generate reliable data to benchmark computational radiation transport models in a radiation field of interest. The presentation will give first results of the physical dose distribution, the comparison with GEANT4 computer simulations based on a Voxel model of the phantom, and a comparison with the data from the chromosome aberration study.

Berger, Thomas↗