Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 973 records · Page 54

Radiation breakage of DNA: a model based on random-walk chromatin structure

Monte Carlo computer software, called DNAbreak, has recently been developed to analyze observed non-random clustering of DNA double strand breaks in chromatin after exposure to densely ionizing radiation. The software models coarse-grained configurations of chromatin and radiation tracks, small-scale details being suppressed in order to obtain statistical results for larger scales, up to the size of a whole chromosome. We here give an analytic counterpart of the numerical model, useful for benchmarks, for elucidating the numerical results, for analyzing the assumptions of a more general but less mechanistic "randomly-located-clusters" formalism, and, potentially, for speeding up the calculations. The equations characterize multi-track DNA fragment-size distributions in terms of one-track action; an important step in extrapolating high-dose laboratory results to the much lower doses of main interest in environmental or occupational risk estimation. The approach can utilize the experimental information on DNA fragment-size distributions to draw inferences about large-scale chromatin geometry during cell-cycle interphase.

Chromatin/radiation effects↗

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming↗

Use of satellite data and modeling to asses the influence of stratospheric processes on the troposphere

The research is comprised of the following tasks: use of simple analytical and numerical models of a coupled troposphere-stratosphere system to examine the effects of radiation and ozone on planetary wave dynamics and the tropospheric circulation; use of satellite data obtained from the Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) instrument and Solar Backscattered Ultraviolet (SBUV) experiment, in conjunction with National Meteorological Center (NMC) data, to determine the planetary wave vertical structures, dominant wave spectra, ozone spectra, and time variations in diabatic heating rate; and synthesis of the modeling and observational results to provide a better understanding of the effects that stratospheric processes have on tropospheric dynamics.

Nathan, Terrence R.↗

Assessing the Usefulness of AIRS Radiance Observations in a 4D-Var Assimilation Scheme Using the Penn State/NCAR Mesoscale Model Version 5 (MM5) and a Stand Alone Radiative Transfer Algorithm (SARTA)

Short-term precipitation forecast skill remains poor despite further advances in numerical modeling and data analysis. Utilizing indirect observations, such as the radiance observations obtained from the Atmospheric Infrared Sounder (AIRS), in a 4D-Var assimilation scheme may lead to a rapid improvement in the skill of 12-h QPF. In order to investigate the usefulness of AIRS observations, the adjoint code of a radiative transfer model (SARTA) is developed and linked to MM5 4D-Var system. A statistical comparison of the AIRS radiances with those calculated using the linked MMS-SARTA model is carried out to assess errors in the radiative transfer model, MM5 and AIRS data, and to develop a quality control scheme. Outliers and any bias that may exist in radiative transfer model, MM5 and AIRS data are identified. An adjoint sensitivity study, utilizing the SARTA adjoint model, is then done to identify the relative sensitivities of each spectral channel to such parameters as atmospheric temperature, pressure, and mixing ratio. Another sensitivity study is done using the MM5 alone to determine the sensitivity of the 12-h QPF to model parameters such as atmospheric temperature, pressure, and mixing ratio under various weather regimes (i.e. summertime convection, winter snowfall). Results from these two types of adjoint sensitivity studies are used to determine which AIRS channels contribute most to 12-h QPF under differing weather regimes. This information is useful when conducting 4D-Var experiments using a subset of the 2378 AIRS spectral channels. Finally, 4D-Var experiments are done to analyze the impact of the AIRS observations on the quality of the 12-h QPF results. This project falls under the data assimilation general session topic as specified by the conference organizers.

Carrier, M.↗

Analyzing Electric Field Morphology Through Data-Model Comparisons of the GEM IM/S Assessment Challenge Events

The storm-time inner magnetospheric electric field morphology and dynamics are assessed by comparing numerical modeling results of the plasmasphere and ring current with many in situ and remote sensing data sets. Two magnetic storms are analyzed, April 22,2001 and October 21-23,2001, which are the events selected for the Geospace Environment Modeling (GEM) Inner Magnetosphere/Storms (IM/S) Assessment Challenge (IMSAC). The IMSAC seeks to quantify the accuracy of inner magnetospheric models as well as synthesize our understanding of this region. For each storm, the ring current-atmosphere interaction model (RAM) and the dynamic global core plasma model (DGCPM) were run together with various settings for the large-scale convection electric field and the nightside ionospheric conductance. DGCPM plasmaspheric parameters were compared with IMAGE-EUV plasmapause extractions and LANL-MPA plume locations and velocities. RAM parameters were compared with Dst*, LANL-MPA fluxes and moments, IMAGE-MENA images, and IMAGE-HENA images. Both qualitative and quantitative comparisons were made to determine the electric field morphology that allows the model results to best fit the plasma data at various times during these events. The simulations with self-consistent electric fields were, in general, better than those with prescribed field choices. This indicates that the time-dependent modulation of the inner magnetospheric electric fields by the nightside ionosphere is quite significant for accurate determination of these fields (and their effects). It was determined that a shielded Volland-Stern field description driven by the 3-hour Kp index yields accurate results much of the time, but can be quite inconsistent. The modified Mcllwain field description clearly lagged in overall accuracy compared to the other fields, but matched some data sets (like Dst*) quite well. The rankings between the simulations varied depending on the storm and the individual data sets, indicating that each field description did well for some place, time, and energy range during the events, as well as doing less well in other places, times, and energies. Several unresolved issues regarding the storm-time inner magnetospheric electric field are discussed.

Liemohn, Michael W.↗

Numerical Simulation of a Tornado Generating Supercell

The development of tornadoes from a tornado generating supercell is investigated with a large eddy simulation weather model. Numerical simulations are initialized with a sounding representing the environment of a tornado producing supercell that affected North Carolina and Virginia during the Spring of 2011. The structure of the simulated storm was very similar to that of a classic supercell, and compared favorably to the storm that affected the vicinity of Raleigh, North Carolina. The presence of mid-level moisture was found to be important in determining whether a supercell would generate tornadoes. The simulations generated multiple tornadoes, including cyclonic-anticyclonic pairs. The structure and the evolution of these tornadoes are examined during their lifecycle.

Proctor, Fred H.↗

A numerical simulation of the full two-dimensional electrothermal de-icer pad

The ability to predict the time-temperature history of electrothermal de-icer pads is important in the subsequent design of improved and more efficient versions. These de-icer pads are installed near the surface of aircraft components, for the specific purpose of removing accreted ice. The proposed numerical model can incorporate the full 2-D geometry through a section of a region (i.e., section of an airfoil), that current 1-D numerical codes are unable to do. Thus, the effects of irregular layers, curvature, etc., can now be accounted for in the thermal transients. Each layer in the actual geometry is mapped via a body-fitted coordinate transformation into uniform, rectangular computational grids. The relevant heat transfer equations are transformed and discretized. To model the phase change that might occur in any accreted ice, in an enthalpy formulation the phase change equations are likewise transformed and discretized. The code developed was tested against numerous classical numerical solutions, as well as against experimental de-icing data on a UH1H rotor blade obtained from the NASA Lewis Research Center. The excellent comparisons obtained show that this code can be a useful tool in predicting the performance of current de-icer models, as well as in the designing of future models.

Masiulaniec, Konstanty C.↗

Using Hough harmonics to validate and assess nonlinear shallow-water models

The implementation of a technique for locating programming errors in shallow-water codes, establishing the correctness of the code, and assessing the performance of the numerical model under various flow conditions is described. The right-hand side of the differential equations is modified in such a way that the exact solution of the nonlinear initial-value problem is known, so that the truncation errors of the numerical scheme can be studied in detail. The exact solution is prescribed to be any linear combination of Hough harmonics which propagate in time according to their natural frequencies.

Dee, Dick P.↗

The signature of H/+/ in the equatorial anomaly - An empirical model

Ion composition data from the OGO-6 satellite are used to develop an empirical model of the H(+) signature in the equatorial anomaly. The empirical model is an attempt to represent the generally repeatable details of the satellite data by a parameterized numerical model accounting for such variables as the altitude, local time and longitude of the observations. Although the resulting model is in general agreement with past findings on the variability of the equatorial anomaly phenomenon, there appears to be a complexity of short-term or narrowly localized variations which cannot be adequately understood on the basis of data from a single satellite.

Taylor, H. A., Jr.↗

Evaluation of an aerodynamic-load prediction method on a STOL fighter configuration

A three-dimensional panel method was used to compute the aerodynamic loads on a large-scale, powered wind tunnel model. The size of the computation grid was selected such that turnaround times, on present computers remained satisfactory. With this constraint, the STOL aircraft configuration was modeled and the capabilities and limitations of the prediction method were investigated. After this study, a similar numerical model was established for a forward-swept-wing aircraft configuration which is planned for large-scale, low-speed testing. The resulting aerodynamic-load predictions are discussed and will be utilized, together with the future experimental results, to obtain a broader validation of this method as well as to reduce the necessary size of the matrix of the test parameters.

Katz, J.↗

Evaluation of an aerodynamic-load prediction method on a STOL fighter configuration

A three-dimensional panel method was used to compute the aerodynamic loads on a large-scale, powered wind tunnel model. The size of the computation grid was selected such that turnaround times, on present computers remained satisfactory. With this constraint, the STOL aircraft configuration was modeled and the capabilities and limitations of the prediction method were investigated. After this study, a similar numerical model was established for a forward-swept-wing aircraft configuration which is planned for large-scale, low-speed testing. The resulting aerodynamic-load predictions are discussed and will be utilized, together with the future experimental results, to obtain a broader validation of this method as well as to reduce the necessary size of the matrix of the test parameters.

Katz, J.↗

Investigations of Clouds and Aerosols on Mars, Venus and Titan

This final report was included as part of a new proposal. This new proposal was selected for funding on 9 Nov. 1999. For Titan our interest during the past few years was to explain the observed asymmetry in the albedo. We suggested earlier, from one-dimensional modeling studies, that vertical transport rates were comparable to particle fall speeds. Since heating of the upper atmosphere, which drives dynamical motions, is largely due to the aerosols, a nonlinear interaction between dynamics, radiative heating and particle microphysics is possible. We pursued this interaction in a two-dimensional model. We showed that the observed variations in the albedo between the two hemispheres and over an orbital cycle, could be due to dynamical motions suspending particles so that particle sizes and optical depths vary across the planet. In the Cassini time frame, future studies of this interaction between dynamics, radiation and microphysics may be worthwhile using, the strong modeling base that others, and we have developed. In our recently approved proposal, however, we plan to extend our modeling to hydrocarbon clouds that lie at lower levels. We know very little about such clouds, and numerical models for their properties are non-existent. These clouds may be observed by the Huygen's Probe, and by the Cassini orbiter, so predictions of their properties should help in the analysis of Cassini data. We have also developed a sophisticated model for the lower, condensational, clouds on Venus. In this model we explored the water vapor budget on Venus, and the properties of the clouds such as particle size distribution. During the past few years we have applied our model of the water ice clouds on Mars to new data sets from Pathfinder, and Mars Global Surveyor. We have compared our predictions of cloud properties with those seen by Pathfinder, and found reasonable agreement.

Toon, Owen B.↗

Investigations of Clouds and Aerosols on Mars, Venus and Titan

For Titan our interest during the past few years was to explain the observed asymmetry in the albedo. We suggested earlier, from one-dimensional modeling studies, that vertical transport rates were comparable to particle fall speeds. Since heating of the upper atmosphere, which drives dynamical motions, is largely due to the aerosols, a non- linear interaction between dynamics, radiative heating and particle microphysics is possible. We pursued this interaction in a two-dimensional model. We showed that the observed variations in the albedo between the two hemispheres and over an orbital cycle, could be due to dynamical motions suspending particles so that particle sizes and optical depths vary across the planet. In the Cassini time frame, future studies of this interaction between dynamics, radiation and microphysics may be worthwhile using the strong modeling base that others, and we have developed. In our recently approved proposal, however, we plan to extend our modeling to hydrocarbon clouds that lie at lower levels. We know very little about such clouds, and numerical models for their properties are non-existent. These clouds may be observed by the Huygen's Probe, and by the Cassini orbiter, so predictions of their properties should help in the analysis of Cassini data. We have also developed a sophisticated model for the lower, condensational, clouds on Venus. In this model we explored the water vapor budget on Venus, and the properties of the clouds such as particle size distribution. Most researchers have investigated the upper clouds on Venus, which are essentially a photochemical aerosol with a long lifetime. The lower clouds, however, are similar to stratus clouds on Earth. These clouds have short lifetimes, and are tightly coupled to the dynamics at the base of the Venus cloud deck. We believe that the holes in these clouds seen at near infrared wavelengths by Galileo are related to some interaction between dynamics and cloud physics. One goal of our recently approved proposed work is to better understand this interaction, and to better understand how these clouds may vary if the climate of Venus were to vary. During the past few years we have applied our model of the water ice clouds on Mars to new data sets from Pathfinder, and Mars Global Surveyor. We have compared our predictions of cloud properties with those seen by Pathfinder, and found reasonable agreement. More recently we have improved our model by including a radiative transfer algorithm and a boundary layer transport scheme . The goal of this work is to understand the thermal inversions seen in Mars Global Surveyor data. We find that the inversions are created by clouds radiatively cooling the atmospheric layers in which the clouds form. We have also initiated laboratory work on the physics of carbon dioxide particle formation. In our recently approved work we propose to use these laboratory data, and to extend our modeling to carbon dioxide clouds in the Martian atmosphere.

Toon, Owen B.↗

Managing Analysis Models in the Design Process

Design of large, complex space systems depends on significant model-based support for exploration of the design space. Integrated models predict system performance in mission-relevant terms given design descriptions and multiple physics-based numerical models. Both the design activities and the modeling activities warrant explicit process definitions and active process management to protect the project from excessive risk. Software and systems engineering processes have been formalized and similar formal process activities are under development for design engineering and integrated modeling. JPL is establishing a modeling process to define development and application of such system-level models.

design↗

Tip Vortex and Wake Characteristics of a Counterrotating Open Rotor

One of the primary noise sources for Open Rotor systems is the interaction of the forward rotor tip vortex and blade wake with the aft rotor. NASA has collaborated with General Electric on the testing of a new generation of low noise, counterrotating Open Rotor systems. Three-dimensional particle image velocimetry measurements were acquired in the intra-rotor gap of the Historical Baseline blade set. The velocity measurements are of sufficient resolution to characterize the tip vortex size and trajectory as well as the rotor wake decay and turbulence character. The tip clearance vortex trajectory is compared to results from previously developed models. Forward rotor wake velocity profiles are shown. Results are presented in a form as to assist numerical modeling of Open Rotor system aerodynamics and acoustics.

VanZante, Dale E.↗

Tip Vortex and Wake Characteristics of a Counterrotating Open Rotor

One of the primary noise sources for Open Rotor systems is the interaction of the forward rotor tip vortex and blade wake with the aft rotor. NASA has collaborated with General Electric on the testing of a new generation of low noise, counterrotating Open Rotor systems. Three-dimensional particle image velocimetry measurements were acquired in the intra-rotor gap of the Historical Baseline blade set. The velocity measurements are of sufficient resolution to characterize the tip vortex size and trajectory as well as the rotor wake decay and turbulence character. The tip clearance vortex trajectory is compared to results from previously developed models. Forward rotor wake velocity profiles are shown. Results are presented in a form as to assist numerical modeling of Open Rotor system aerodynamics and acoustics.

Van Zante, Dale E.↗

Hybrid-Vlasov Simulation of Soft X-Ray Emissions at the Earth’S Dayside Magnetospheric Boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause, the magnetosheath and the polar cusps by remote sensing techniques. The Solar wind–Magnetosphere–Ionosphere Link Explorer (SMILE) and Lunar Environment heliospheric X-ray Imager (LEXI) missions aim to obtain soft X-ray images of near-Earth space thanks to their Soft X-ray Imager (SXI) instruments. While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission, the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma. To investigate the possible signatures of ion-kinetic-scale processes in soft X-ray images, we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model. The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field. We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations, providing meridional and equatorial views. We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirror mode wave structures in the magnetosheath and flux transfer events (FTEs) at the magnetopause. Our results suggest that, although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers, mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images. For instance, a local increase by 30% in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12% enhancement in the line-of-sight- and time integrated soft X-ray emissivity originating from this region. Likewise, a proton density increase by 14% in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%. These are likely conservative estimates, given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind. These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale, transient phenomena occurring on the dayside.

magnetosphere↗

Jupiter radiation belt models (July 1974)

Flux profiles which were derived from data returned by Pioneer 10 during Jupiter encounter, form the basis for a new set of numerical models for the energy spectra of electrons and protons in Jupiter's inner magnetosphere.

Divine, N.↗