Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation modeling”

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 1,315 records · Page 73

Role of Troposphere-Convection-Land Coupling in the Southwestern Amazon Precipitation Bias of the Community Earth System Model version 1 (CESM1)

The Amazon forest plays an important role in the tropical circulation and global water and biogeochemical cycles. However, the common Amazon precipitation biases in global climate models undermine our ability to understand future changes in the Amazon and their global consequences. The Community Earth System Model exhibits a persistent wet season dry bias over the southwestern Amazon. The bias is also characterized by inaccurate diurnal timing and intensity of rainfall events. Analysis of a diverse set of simulations suggests that the bias originates from local processes. A known land model bias in the partitioning of surface energy fluxes is found not to contribute significantly to the rainfall bias. Key model processes are identified to be the low-level divergent circulation controlling horizontal moisture flow and the sensitivity of the modeled deep convection to the lower tropospheric moisture. The lower tropospheric moisture divergence dictates the spatiotemporal distributions of convective precipitation through parameterized entrainment mixing, which generally improves the rainfall diurnal cycle. However, it also limits convective precipitation over the areas with frequent and large-scale divergence in the lower troposphere. The southwestern Amazon is such an area, particularly in the late afternoon, although uncertainty exists in the low-level flow in the simulation and reanalyses. One way to sustain rainfall in the presence of low-level divergence is a convective trigger based on both the moisture and temperature changes by large-scale processes, which is shown to help in reducing the dry bias.

Sakaguchi, Koichi↗

Numerical Investigation of Plasma Detachment in Magnetic Nozzle Experiments

At present there exists no generally accepted theoretical model that provides a consistent physical explanation of plasma detachment from an externally-imposed magnetic nozzle. To make progress towards that end, simulation of plasma flow in the magnetic nozzle of an arcjet experiment is performed using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The aim is to compare the computational results with various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver (see Fig. I), and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. This geometry is modeled because there is a substantial amount of experimental data that can be compared to the computational results, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters. Further studies will focus on modeling literature experiments both for the purpose of additional code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗

A basis for solid modeling of gear teeth with application in design and manufacture

A new approach to modeling gear tooth surfaces is discussed. A computer graphics solid modeling procedure is used to simulate the tooth fabrication process. This procedure is based on the principles of differential geometry that pertain to envelopes of curves and surfaces. The procedure is illustrated with the modeling of spur, helical, bevel, spiral bevel, and hypoid gear teeth. Applications in design and manufacturing are discussed. Extensions to nonstandard tooth forms, to cams, and to rolling element bearings are proposed.

Huston, Ronald L.↗

A Basis for Solid Modeling of Gear Teeth with Application in Design and Manufacture

This paper discusses a new approach to modeling gear tooth surfaces. A computer graphics solid modeling procedure is used to simulate the tooth fabrication processes. This procedure is based on the principles of differential geometry that pertain to envelopes of curves and surfaces. The procedure is illustrated with the modeling of spur, helical, bevel, spiral bevel and hypoid gear teeth. Applications in design and manufacturing arc discussed. Extensions to nonstandard tooth forms, to cams, and to rolling element hearings are proposed.

Huston, Ronald L.↗

The Global Perspective on the Evolution of Solids in a Protoplanetary Disk

It is currently thought that planets around solar-type stars form by the accumulation of solid matter entrained in a gaseous, turbulent protoplanetary disk. We have developed a model designed to simulate the part of this process that starts from small particles suspended in the gaseous disk at the end of the formation stage, and ends up with most of the solid material aggregated into 1-10-km planetesimals. The major novelty of our approach is its emphasis on the global, comprehensive treatment of the problem, as our model simultaneously keeps track of the evolution of gas and solid particles due to gas-solid coupling, coagulation, sedimentation, and evaporation/condensation. The result of our calculations is the radial distribution of solid material circumnavigating a star in the form of a planetesimal swarm. Such a distribution should well approximate the radial apportionment of condensed components of the planets spread over the radial extent of the mature planetary system. Therefore we view our calculations as an attempt to predict the large-scale architecture of planetary systems and to assess their potential diversity. In particular, we have found that some initial conditions lead to all solids being lost to the star, but we can also identify initial conditions leading to a radial distribution of solid material quite reminiscent of what is found in our solar system.

Stepinski, T. F.↗

Matching ecohydrological processes and scales of banded vegetation patterns in semiarid catchments

While the claim that water-carbon interactions result in spatially coherent vegetation patterning is rarely disputed in many arid and semiarid regions, the significance of the detailed water pathways and other high frequency variability remain an open question. How the short temporal scale meteorological fluctuations form the long-term spatial variability of available soil water in complex terrains due to the various hydrological, land surface, and vegetation dynamic feedbacks frames the scope of the work here. Knowledge of the detailed mechanistic feedbacks among soil, plants, and the atmosphere will lead to advances in our understanding of plant water availability in arid and semiarid ecosystems and will provide insights for future model development concerning vegetation pattern formation. In this paper, quantitative estimates of water fluxes and vegetation productivity are provided for a semiarid ecosystem with established vegetation bands on hillslopes using numerical simulations. A state-of-the-science process based ecohydrological model is used, which resolves hydrological and plant physiological processes at the relevant space and time scales, for relatively small periods (e.g., decades) of mature ecosystems (i.e., spatially static vegetation distribution). To unfold the mechanisms that shape the spatial distribution of soil moisture, plant productivity and the relevant surface/subsurface and atmospheric water fluxes, idealized hillslope numerical experiments are constructed, where the effects of soil type, slope steepness, and overland flow accumulation area are quantified. Those mechanisms are also simulated in the presence of complex topography features on landscapes. Finally, the main results are (a) short temporal scale meteorological variability and accurate representation of the scales at which each ecohydrological process operates are crucial for the estimation of the spatial variability of soil water availability to the plant root zone; (b) water fluxes such as evapotranspiration, infiltration, runoff-run-on, and subsurface soil water movement have a dynamic short temporal scale behavior that determines the long-term spatial organization of plant soil water availability in ecosystems with established vegetation patterns; and (c) hypotheses concerning the hydrological responses that can lead to vegetation pattern formation have to accommodate realistic and physically based representations of the fast dynamics of key ecohydrological fluxes.

54 ENVIRONMENTAL SCIENCES↗

Time-of-Flight Estimation for Nondestructive Evaluation

Ultrasonic testing (UT) for nondestructive evaluation (NDE) is a critical entity necessary to resolve both the quality and precision questions of complex parts evolving from the classical approaches to the innovative additive manufacturing (AM) process. This modality provides the essential quantitative information for acceptance and potential flaw detectionof apart under investigation. A primary ingredient in UT besides the required precision robotic hardware for theacquisition of high-qualitymeasurement data is the underlying signal processing. It is here that much of the system performance capability resides. In this report,we discuss the basic steps in UT signal processing along with current capabilities that must be achieved in order to satisfy the criticaldemands of the various LLNL programs.In itsbasic foundations, signal processing is essentially the “extraction of critical information (signals) from noisy, uncertain data.” Clearly, in a perfect world the best signal processing is none---just make a reliable, noise-free, uncertainty-free, measurement. Unfortunately, even the best of systems is still straddled within the confines and limitations of physical instrumentation. With this in mind, we discuss the development of signal processing techniques toextract the desired UT information from noisy measurement data. We start with a discussion of a simple homogeneous representation of a “part” under investigation and its insonification by an ultrasonicexcitation from a physics-based perspective. Once developed, we briefly discussthe various choices of excitation signals and their tradeoffs.Next,we discuss a simulation approach employing simple models from the signal processing perspective and then move on to the development of the basic signal processing approach to ultrasonic signal processing all based on estimating the pulse arrival estimation. Starting with simple peak detection techniques, progressing to the processing workhorse so-called “matched-filter” and finally to the more sophisticated “model-based matched filter.”The underlying pre-and post-processing of noisy measurement data along with the application of these techniques are subsequently demonstrated on experimental ultrasonic data.Finally, we present some results for processing “weld” data when the ultrasonic excitation signal is not directly available.

42 ENGINEERING↗

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗

A stochastic model for particle impingements on orbiting spacecraft

A general methodology for simulating particle impingements on orbiting spacecraft is developed. Major steps in the modeling process are presented as (1) modeling objective, (2) construction of the spacecraft geometrical model, (3) simulation of the particles in the space environment, (4) particle impact and subsequent events of interest, and (5) results of the simulation. A simulation of the expected meteoroid impingements on the Hubble Space Telescope and the resulting angular momentum transfers which can cause telescope pointing disturbances is given to illustrate these methods.

Howell, L. W., Jr.↗

A stochastic model for particle impingements on orbiting spacecraft

A general methodology for simulating particle impingements on orbiting spacecraft is developed. Major steps in the modeling process are presented as (1) modeling objective, (2) construction of the spacecraft geometrical model, (3) simulation of the particles in the space environment, (4) particle impact and subsequent events of interest, and (5) results of the simulation. A simulation of the expected meteoroid impingements on the Hubble Space Telescope and the resulting angular momentum transfers which can cause telescope pointing disturbances is given to illustrate these methods.

Howell, Leonard W., Jr.↗

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗

A dynamic parameterization of sulfuric acid–dimethylamine nucleation and its application in three-dimensional modeling

Sulfuric acid (SA) is a governing gaseous precursor for atmospheric new particle formation (NPF), a major source of global ultrafine particles, in environments studied around the world. In polluted urban atmospheres with high condensation sinks (CSs), the formation of stable SA–amine clusters, such as SA–dimethylamine (DMA) clusters, usually initializes intense NPF events. Coagulation scavenging and cluster evaporation are dominant sink processes of SA–amine clusters in urban atmospheres, yet these loss processes are not quantitatively included in the present parameterizations of SA–amine nucleation. We herein report a parameterization of SA–DMA nucleation, based on cluster dynamic simulations and quantum chemistry calculations, with certain simplifications to greatly reduce the computational costs. Compared with previous SA–DMA nucleation parameterizations, this new parameterization was able to reproduce the dependences of particle formation rates on temperature and CSs. We then incorporated it in a three-dimensional (3-D) chemical transport model to simulate the evolution of the particle number size distributions. Simulation results showed good consistency with the observations in the occurrence of NPF events and particle number size distributions in wintertime Beijing and represented a significant improvement compared to that using a parameterization without coagulation scavenging. Quantitative analysis shows that SA–DMA nucleation contributes significantly to nucleation rates and aerosol population during the 3-D simulations in Beijing (>99 % and >60 %, respectively). These results broaden the understanding of NPF in urban atmospheres and stress the necessity of including the effects of coagulation scavenging and cluster stability in simulating SA–DMA nucleation in 3-D simulations. Representing these processes is thus likely to improve model performance in particle source apportionment and quantification of aerosol effects on air quality, human health, and climate.

Li, Yuyang↗

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗

Simulation of multi-photon emission isotopes using time-resolved SimSET multiple photon history generator

Abstract-Multiple-photon emitters, such as In-111 or Se-75, have enormous potential in the field of nuclear medicine imaging. For example, Se-75 can be used to investigate the bile acid malabsorption and measure the bile acid pool loss. The simulation system for emission tomography (SimSET) is a well-known Monte Carlo simulation (MCS) code in nuclear medicine for its high computational efficiency. However, current SimSET cannot simulate these isotopes due to the lack of modeling of complex decay scheme and the time-dependent decay process. To extend the versatility of SimSET for simulation of those multi-photon emission isotopes, a time-resolved multiple photon history generator based on SimSET codes is developed in present study. For developing the time-resolved SimSET (trSimSET) with radionuclide decay process, the new MCS model introduce new features, including decay time information and photon time-of-flight information, into this new code. The half-life of energy states were tabulated from the Evaluated Nuclear Structure Data File (ENSDF) database. The MCS results indicate that the overall percent difference is less than 8.5% for all simulation trials as compared to GATE. To sum up, we demonstrated that time-resolved SimSET multiple photon history generator can have comparable accuracy with GATE and keeping better computational efficiency. The new MCS code is very useful to study the multi-photon imaging of novel isotopes that needs the simulation of lifetime and the time-of-fight measurements. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Multibody Based Digital Astronaut Dynamics Simulation

BACKGROUND: This study provides the Software, Robotics, & Simulation Division at the NASA Johnson Space Center with a verification tool for multibody dynamics simulation requiring human motion. The motivation stems from current studies of several Vibration Isolation & Stabilization(VIS)system designs that attenuate the moments and forces which would be transmitted to a spacecraft during an exercise. A multibody dynamics model for a proposed VIS was available previously[1], therefore modeling of the VIS was not needed for this work. The interest here is in creating the multibody dynamics model of an astronaut in motion which may be utilized independently or while attached to a mechanism. An existing simulation [2] that utilizes OpenSim [3,4] and an in-house multibody dynamics package (MBDyn) [5] is used in order to verify the astronaut model. The main advantage this model will have over the existing simulation is that everything will be processed in one tool. METHODS AND RESULTS: Creating the simulation required; estimation of Body Segment Inertial Parameters (BSIP),a multibody model of the human-VIS system, joint acceleration profiles, and input files for MBDyn, which is used for this analysis. The scaling factors provided by Dumas et al. [6] are utilized in estimating the BSIP. Anthropometric data are used for estimating these parameters, the Anthropometric Survey of US Army Personnel (ANSUR II) [7] was the source. The astronaut model consists of 15 bodies, 14 joints and 32 degrees of freedom, with the dynamics topology generated using the center of mass locations and anthropometric data. MBDyn has an option for prescribed joint motion (PJM), which requires joint acceleration data as input. The joint angle data is first obtained from a motion capture system and then processed through code that has been created to generate approximate joint acceleration profiles. The topology tree for the astronaut model begins at the right foot up to the pelvis where there is one branch for going down the left leg and another for the torso. The torso branch leads to branches for the arms and a leaf body for the head/neck segment. For attachment to the VIS, the heel of the right foot is connected to the VIS platform through a fixed joint, resembling a foot restraint. The left foot does not attach to the platform in order to prevent a system with a closed loop. Topology and symmetry of the astronaut model were verified through kinematic analysis. Further verification of the forces and moments transmitted to the VIS were verified against the existing simulation. There was a satisfactory level of agreement when testing a simple motion, for example, rocking back and forth. Full exercise motions are to be tested soon. The main outcome has been a novel application of MBDyn for biomechanics modeling that is now available for dynamic simulations involving human motion. The estimation of BSIP was another useful result of this study, requiring only 15 inputs for generating mass properties of a theoretical astronaut model. Expansion on this work is possible by going through an alternative route in obtaining the joint motion data. Instead of high-tech and often expensive motion capture systems, an individual may watch videos with high focus and at a slow motion for each individual segment in order to determine the initial and final time and angle for that specific degree of freedom. Synthetic trajectories may also be created if there is no video reference available.

F N Matari↗

Electromagnetic emission from long-lived binary neutron star merger remnants. I. Formulation of the problem

Binary neutron star (BNS) mergers are the leading model to explain the phenomenology of short gamma-ray bursts (SGRBs). Recent observations of long-lasting X-ray afterglows of SGRBs challenge standard paradigms and indicate that in a large fraction of events a long-lived neutron star (NS) may be formed rather than a black hole. Understanding the mechanisms underlying these afterglows is necessary in order to address the open questions concerning the nature of SGRB central engines. However, recent theoretical progress has been hampered by the fact that the timescales of interest for the afterglow emission are inaccessible to numerical relativity simulations. Here we present a detailed model to bridge the gap between numerical simulations of the merger process and the relevant timescales for the afterglows, assuming that the merger results in a long-lived NS. This model is formulated in terms of a set of coupled differential equations that follow the evolution of the post-merger system and predict its electromagnetic (EM) emission in a self-consistent way, starting from initial data that can be extracted from BNS merger simulations. The model presented here also allows us to search for suitable EM counterparts for multimessenger astronomy, which is expected to become reality within the next few years thanks to ground-based GW detectors such as advanced LIGO and Virgo. This paper discusses the formulation and implementation of the model. In a companion paper, we employ this model to predict the EM emission from ∼10{sup −2} to ∼10{sup 7} s after a BNS merger and discuss the implications in the context of SGRBs and multimessenger astronomy.

79 ASTRONOMY AND ASTROPHYSICS↗