Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled models”

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 883 records · Page 49

A 7.5-Year Dataset of SSM/I-Derived Surface Turbulent Fluxes Over Global Oceans

The global air-sea turbulent fluxes are needed for driving ocean models and validating coupled ocean-atmosphere global models. A method was developed to retrieve surface air humidity from the radiances measured by the Special Sensor Microwave/Imager (SSM/I) Using both SSM/I-retrieved surface wind and air humidity, they computed daily turbulent fluxes over global oceans with a stability-dependent bulk scheme. Based on this method, we have produced Version 1 of Goddard Satellite-Based Surface Turbulent Fluxes (GSSTF) dataset from the SSM/I data and other data. It provides daily- and monthly-mean surface turbulent fluxes and some relevant parameters over global oceans for individual F8, F10, and F11 satellites covering the period July 1987-December 1994. It also provides 1988-94 annual- and monthly-mean climatologies of the same variables, using only F8 and F1 1 satellite data. It has a spatial resolution of 2.0 degrees x 2.5 degrees lat-long and is archived at the NASA/GSFC DAAC. The purpose of this paper is to present an updated assessment of the GSSTF 1.0 dataset.

Chou, Shu-Hsien↗

Fast-Running Aeroelastic Code Based on Unsteady Linearized Aerodynamic Solver Developed

The NASA Glenn Research Center has been developing aeroelastic analyses for turbomachines for use by NASA and industry. An aeroelastic analysis consists of a structural dynamic model, an unsteady aerodynamic model, and a procedure to couple the two models. The structural models are well developed. Hence, most of the development for the aeroelastic analysis of turbomachines has involved adapting and using unsteady aerodynamic models. Two methods are used in developing unsteady aerodynamic analysis procedures for the flutter and forced response of turbomachines: (1) the time domain method and (2) the frequency domain method. Codes based on time domain methods require considerable computational time and, hence, cannot be used during the design process. Frequency domain methods eliminate the time dependence by assuming harmonic motion and, hence, require less computational time. Early frequency domain analyses methods neglected the important physics of steady loading on the analyses for simplicity. A fast-running unsteady aerodynamic code, LINFLUX, which includes steady loading and is based on the frequency domain method, has been modified for flutter and response calculations. LINFLUX, solves unsteady linearized Euler equations for calculating the unsteady aerodynamic forces on the blades, starting from a steady nonlinear aerodynamic solution. First, we obtained a steady aerodynamic solution for a given flow condition using the nonlinear unsteady aerodynamic code TURBO. A blade vibration analysis was done to determine the frequencies and mode shapes of the vibrating blades, and an interface code was used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor was used to interpolate the mode shapes from the structural dynamic mesh onto the computational dynamics mesh. Then, we used LINFLUX to calculate the unsteady aerodynamic forces for a given mode, frequency, and phase angle. A postprocessor read these unsteady pressures and calculated the generalized aerodynamic forces, eigenvalues, and response amplitudes. The eigenvalues determine the flutter frequency and damping. As a test case, the flutter of a helical fan was calculated with LINFLUX and compared with calculations from TURBO-AE, a nonlinear time domain code, and from ASTROP2, a code based on linear unsteady aerodynamics.

Reddy, T. S. R.↗

Recovery of the Antarctic Ozone Hole

The Antarctic ozone hole develops each year and culminates by early Spring. Antarctic ozone values have been monitored since 1979 using satellite observations from the TOMS and OMI instruments. The severity of the hole has been assessed using the minimum total ozone value from the October monthly mean (depth of the hole), the average size during the September-October period, and the ozone mass deficit. Ozone is mainly destroyed by halogen catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this presentation, we show the relationships of halogens and temperature to both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. We use two methods to estimate ozone hole recovery. First, we use projections of halogen levels combined with age-of-air estimates in a parametric model. Second, we use a coupled chemistry climate model to assess recovery. We find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. Furthermore, full recovery to 1980 levels will not occur until approximately 2068. We will also show some error estimates of these dates and the impact of climate change on the recovery.

Newman, Paul A.↗

Towards Improved Forecasts of Atmospheric and Oceanic Circulations over the Complex Terrain of the Eastern Mediterranean

Forecasting atmospheric and oceanic circulations accurately over the Eastern Mediterranean has proved to be an exceptional challenge. The existence of fine-scale topographic variability (land/sea coverage) and seasonal dynamics variations can create strong spatial gradients in temperature, wind and other state variables, which numerical models may have difficulty capturing. The Hellenic Center for Marine Research (HCMR) is one of the main operational centers for wave forecasting in the eastern Mediterranean. Currently, HCMR's operational numerical weather/ocean prediction model is based on the coupled Eta/Princeton Ocean Model (POM). Since 1999, HCMR has also operated the POSEIDON floating buoys as a means of state-of-the-art, real-time observations of several oceanic and surface atmospheric variables. This study attempts a first assessment at improving both atmospheric and oceanic prediction by initializing a regional Numerical Weather Prediction (NWP) model with high-resolution sea surface temperatures (SST) from remotely sensed platforms in order to capture the small-scale characteristics.

Chronis, Themis↗

Optimal Initial Perturbations for Ensemble Prediction of the Madden-Julian Oscillation during Boreal Winter

An initialization strategy, tailored to the prediction of the Madden-Julian oscillation (MJO), is evaluated using the Goddard Earth Observing System Model, version 5 (GEOS-5), coupled general circulation model (CGCM). The approach is based on the empirical singular vectors (ESVs) of a reduced-space statistically determined linear approximation of the full nonlinear CGCM. The initial ESV, extracted using 10 years (1990-99) of boreal winter hindcast data, has zonal wind anomalies over the western Indian Ocean, while the final ESV (at a forecast lead time of 10 days) reflects a propagation of the zonal wind anomalies to the east over the Maritime Continent an evolution that is characteristic of the MJO. A new set of ensemble hindcasts are produced for the boreal winter season from 1990 to 1999 in which the leading ESV provides the initial perturbations. The results are compared with those from a set of control hindcasts generated using random perturbations. It is shown that the ESV-based predictions have a systematically higher bivariate correlation skill in predicting the MJO compared to those using the random perturbations. Furthermore, the improvement in the skill depends on the phase of the MJO. The ESV is particularly effective in increasing the forecast skill during those phases of the MJO in which the control has low skill (with correlations increasing by as much as 0.2 at 20 25-day lead times), as well as during those times in which the MJO is weak.

Ham, Yoo-Geun↗

A Simple Model of Global Aerosol Indirect Effects

Most estimates of the global mean indirect effect of anthropogenic aerosol on the Earth's energy balance are from simulations by global models of the aerosol lifecycle coupled with global models of clouds and the hydrologic cycle. Extremely simple models have been developed for integrated assessment models, but lack the flexibility to distinguish between primary and secondary sources of aerosol. Here a simple but more physically based model expresses the aerosol indirect effect (AIE) using analytic representations of cloud and aerosol distributions and processes. Although the simple model is able to produce estimates of AIEs that are comparable to those from some global aerosol models using the same global mean aerosol properties, the estimates by the simple model are sensitive to preindustrial cloud condensation nuclei concentration, preindustrial accumulation mode radius, width of the accumulation mode, size of primary particles, cloud thickness, primary and secondary anthropogenic emissions, the fraction of the secondary anthropogenic emissions that accumulates on the coarse mode, the fraction of the secondary mass that forms new particles, and the sensitivity of liquid water path to droplet number concentration. Estimates of present-day AIEs as low as 5 W/sq m and as high as 0.3 W/sq m are obtained for plausible sets of parameter values. Estimates are surprisingly linear in emissions. The estimates depend on parameter values in ways that are consistent with results from detailed global aerosol-climate simulation models, which adds to understanding of the dependence on AIE uncertainty on uncertainty in parameter values.

Clouds↗

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale repeating unit cell model of a single spherulite containing four disparate length scales was developed to predict the thermoelastic behavior of semicrystalline thermoplastic materials for composite aerospace applications. The continuum level scales were fully coupled and modeled using the generalized method of cells and the high fidelity generalized method of cells micromechanics theories. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents in the multiscale continuum models. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, showing good agreement with the available experimental data from the open literature. Moreover, it is shown that predicted properties are fairly insensitive to the fidelity of the micromechanics model used at the highest continuum scale or the assumed shape of the spherulite.

thermoplastics↗

Simulation of Seasonal Snow Microwave TB Using Coupled Multi-Layered Snow Evolution and Microwave Emission Models

The accurate quantification of SWE has important societal benefits, including improving domestic and agricultural water planning, flood forecasting and electric power generation. However, passive-microwave SWE algorithms suffer from variations in TB due to snow metamorphism, difficult to distinguish from those due to SWE variations. Coupled snow evolution-emission models are able to predict snow metamorphism, allowing us to account for emissivity changes. They can also be used to identify weaknesses in the snow evolution model. Moreover, thoroughly evaluating coupled models is a contribution toward the assimilation of TB, which leads to a significant increase in the accuracy of SWE estimates.

microwave↗

LES of a Jet Excited by the Localized Arc Filament Plasma Actuators

The fluid dynamics of a high-speed jet are governed by the instability waves that form in the free-shear boundary layer of the jet. Jet excitation manipulates the growth and saturation of particular instability waves to control the unsteady flow structures that characterize the energy cascade in the jet.The results may include jet noise mitigation or a reduction in the infrared signature of the jet. The Localized Arc Filament Plasma Actuators (LAFPA) have demonstrated the ability to excite a high-speed jets in laboratory experiments. Extending and optimizing this excitation technology, however, is a complex process that will require many tests and trials. Computational simulations can play an important role in understanding and optimizing this actuator technology for real-world applications. Previous research has focused on developing a suitable actuator model and coupling it with the appropriate computational fluid dynamics (CFD) methods using two-dimensional spatial flow approximations. This work is now extended to three-dimensions (3-D) in space. The actuator model is adapted to a series of discrete actuators and a 3-D LES simulation of an excited jet is run. The results are used to study the fluid dynamics near the actuator and in the jet plume.

Brown, Clifford A.↗

Simulating Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi

A microstructure-based FEM model that couples crystal plasticity, crystallographic descriptions of the B2-B19' martensitic phase transformation, and anisotropic elasticity is used to simulate thermal cycling and isothermal deformation in polycrystalline NiTi (49.9at% Ni). The model inputs include anisotropic elastic properties, polycrystalline texture, DSC data, and a subset of isothermal deformation and load-biased thermal cycling data. A key experimental trend is captured.namely, the transformation strain during thermal cycling is predicted to reach a peak with increasing bias stress, due to the onset of plasticity at larger bias stress. Plasticity induces internal stress that affects both thermal cycling and isothermal deformation responses. Affected thermal cycling features include hysteretic width, two-way shape memory effect, and evolution of texture with increasing bias stress. Affected isothermal deformation features include increased hardening during loading and retained martensite after unloading. These trends are not captured by microstructural models that lack plasticity, nor are they all captured in a robust manner by phenomenological approaches. Despite this advance in microstructural modeling, quantitative differences exist, such as underprediction of open loop strain during thermal cycling.

Manchiraju, Sivom↗

Calibration of Two-dimensional Floodplain Modeling in the Atchafalaya River Basin Using SAR Interferometry

Two-dimensional (2D) satellite imagery has been increasingly employed to improve prediction of floodplain inundation models. However, most focus has been on validation of inundation extent, with little attention on the 2D spatial variations of water elevation and slope. The availability of high resolution Interferometric Synthetic Aperture Radar (InSAR) imagery offers unprecedented opportunity for quantitative validation of surface water heights and slopes derived from 2D hydrodynamic models. In this study, the LISFLOOD-ACC hydrodynamic model is applied to the central Atchafalaya River Basin, Louisiana, during high flows typical of spring floods in the Mississippi Delta region, for the purpose of demonstrating the utility of InSAR in coupled 1D/2D model calibration. Two calibration schemes focusing on Manning s roughness are compared. First, the model is calibrated in terms of water elevations at a single in situ gage during a 62 day simulation period from 1 April 2008 to 1 June 2008. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. The best-fit models show that the mean absolute errors are 3.8 cm for a single in situ gage calibration and 5.7 cm/46 days for InSAR water level calibration. The optimum values of Manning's roughness coefficients are 0.024/0.10 for the channel/floodplain, respectively, using a single in situ gage, and 0.028/0.10 for channel/floodplain the using SAR. Based on the calibrated water elevation changes, daily storage changes within the size of approx 230 sq km of the model area are also calculated to be of the order of 107 cubic m/day during high water of the modeled period. This study demonstrates the feasibility of SAR interferometry to support 2D hydrodynamic model calibration and as a tool for improved understanding of complex floodplain hydrodynamics

Hahn Chul Jung↗

Linear Parameter-Varying Control Compatible with Excitation for Real-Time Modeling

The Modeling and Control for Agile Aircraft Development project (MCAAD) seeks to reduce the time spent on wind tunnel testing and computational fluid dynamic modeling when developing control laws for new commercial and military aircraft. MCAAD uses advanced dynamic modeling procedures coupled with real-time control design which can reduce the time and cost spent on traditional aircraft development methods. This is done in-flight by using a real-time modeling procedure where programmed test inputs (PTI) deflect the aircraft control surfaces, and the responses are measured in order to create a dynamic model. This paper presents an 𝐻∞-based linear parameter-varying (LPV) control law designed on a nonlinear rigid-body dynamic model which commands control surface deflections at frequencies above the PTI frequency band. Thus the controller will not suppress the response due to the PTI. The LPV controller is compared with a previously designed Nonlinear Dynamic Inversion (NDI) controller and the coherence between the commanded control surface deflections and the PTI is plotted for both controllers. The LPV controller is shown to distort the PTI less than the NDI controller.

Luke J. Miller↗

Development of a Reactor Model for Chemical Conversion of Lunar Regolith

Lunar regolith will be used for a variety of purposes such as oxygen and propellant production and manufacture of various materials. The design and development of chemical conversion reactors for processing lunar regolith will require an understanding of the coupling among the chemical, mass and energy transport processes occurring at the length and time scales of the overall reactor with those occurring at the corresponding scales of the regolith particles. To this end, a coupled transport model is developed using, as an example, the reduction of ilmenite-containing regolith by a continuous flow of hydrogen in a flow-through reactor. The ilmenite conversion occurs on the surface and within the regolith particles. As the ilmenite reduction proceeds, the hydrogen in the reactor is consumed, and this, in turn, affects the conversion rate of the ilmenite in the particles. Several important quantities are identified as a result of the analysis. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the time for hydrogen to diffuse into the pores of the regolith particles and the chemical reaction time. The paper investigates the relationships between these quantities and their impact on the regolith conversion. Application of the model to various chemical reactor types, such as fluidized-bed, packed-bed, and rotary-bed configurations, are discussed.

Hedge, uday↗

Development of a Reactor Model for Chemical Conversion of Lunar Regolith

Lunar regolith will be used for a variety of purposes such as oxygen and propellant production and manufacture of various materials. The design and development of chemical conversion reactors for processing lunar regolith will require an understanding of the coupling among the chemical, mass and energy transport processes occurring at the length and time scales of the overall reactor with those occurring at the corresponding scales of the regolith particles. To this end, a coupled transport model is developed using, as an example, the reduction of ilmenite-containing regolith by a continuous flow of hydrogen in a flow-through reactor. The ilmenite conversion occurs on the surface and within the regolith particles. As the ilmenite reduction proceeds, the hydrogen in the reactor is consumed, and this, in turn, affects the conversion rate of the ilmenite in the particles. Several important quantities are identified as a result of the analysis. Reactor scale parameters include the void fraction (i.e., the fraction of the reactor volume not occupied by the regolith particles) and the residence time of hydrogen in the reactor. Particle scale quantities include the time for hydrogen to diffuse into the pores of the regolith particles and the chemical reaction time. The paper investigates the relationships between these quantities and their impact on the regolith conversion. Application of the model to various chemical reactor types, such as fluidized-bed, packed-bed, and rotary-bed configurations, are discussed.

Hegde, U.↗

Coupled Attitude and Orbit Dynamics and Control in Formation Flying Systems

Formation flying systems can range from global constellations offering extended service coverage to clusters of highly coordinated vehicles that perform distributed sensing. Recently, the use of groups of micro-satellites in the areas of near Earth explorations, deep space explorations, and military applications has received considerable attention by researchers and practitioners. To date, most proposed control strategies are based on linear models (e.g., Hill-Clohessy-Wiltshire equations) or nonlinear models that are restricted to circular reference orbits. Also, all models in the literature are uncoupled between relative position and relative attitude. In this paper, a generalized dynamic model is proposed. The reference orbit is not restricted to the circular case. In this formulation, the leader or follower satellite can be in either a circular or an elliptic orbit. In addition to maintaining a specified relative position, the satellites are also required to maintain specified relative attitudes. Thus the model presented couples vehicle attitude and orbit requirements. Orbit perturbations are also included. In particular, the J(sub 2) effects are accounted in the model. Finally, a sliding mode controller is developed and used to control the relative attitude of the formation and the simulation results are presented.

Yunjun Xu↗

Propagation studies using a theoretical ionosphere model

The mid-latitude ionospheric and neutral atmospheric models are coupled with an advanced three dimensional ray tracing program to see what success would be obtained in predicting the wave propagation conditions and to study to what extent the use of theoretical ionospheric models is practical. The Penn State MK 1 ionospheric model, the Mitra-Rowe D region model, and the Groves' neutral atmospheric model are used throughout this work to represent the real electron densities and collision frequencies. The Faraday rotation and differential Doppler velocities from satellites, the propagation modes for long distance high frequency propagation, the group delays for each mode, the ionospheric absorption, and the spatial loss are all predicted.

Lee, M.↗

Modeling the Formation and Eruption of Coronal Structures by Linking Data-Driven Magnetofrictional and MHD Simulations for AR 12673

Context .The data-driven and time-dependent modeling of coronal magnetic fields is crucial for understanding solar eruptions. These efforts are complicated by the challenges of finding a balance between physical realism and computing efficiency. One possible technique is to couple two modeling approaches. Aims . Our aim here is to showcase our progress in using time-dependent magnetofrictional model (TMFM) results as input to dynamical magnetohydrodynamic (MHD) simulations. However, due to the different evolution processes in these two models, using TMFM snapshots in an MHD simulation is nontrivial. We address these issues, both physically and numerically, discuss the incompatibility of the TMFM output to serve as the initial condition in MHD simulations, and show our methods of mitigating this. The evolution of the flux systems and the cause of the eruption are investigated. Methods . TMFM is a prevalent approach that has proven to be a very useful tool in the study of the formation of unstable structures in the solar corona. In particular, it is capable of incorporating observational data as initial and boundary conditions and requires shorter computational time compared to MHD simulations. To leverage the efficiency of data-driven TMFM and also to simulate eruptive events in the MHD framework, one can apply TMFM up to a certain time before the expected eruption(s) and then proceed with the simulation in the full or ideal MHD regime in order to more accurately capture the eruption process. Results . We show the results of a benchmark test case with a linked TMFM and MHD simulation to study the evolution of NOAA active region 12673. A rise of a twisted flux bundle through the MHD simulation domain is observed, but we find that the rate of the rise and the altitude reached depends on the time of the TMFM snapshot that was used to initialize the MHD simulation and the helicity injected into the system. The analysis suggested that torus instability and slip-running reconnection could play an important role in the eruption. Conclusions . The results show that the linkage of TMFM and zero-β MHD models can be successfully used to model the eruptive coronal magnetic fields.

F. Daei↗

Generation and Calibration of Linear Models of Aircraft with Highly Coupled Aeroelastic and Flight Dynamics

The lightweight structures and unconventional configurations being considered for the next generation of aircraft mean that any effort to predict or control the flight dynamics is impacted by the structural dynamics. One of the most severe forms of coupling between aeroelasticity and flight dynamics is an instability called body freedom flutter. The existing tools often assume a relatively weak effect of structural dynamics on the flight dynamics, and are therefore incapable of modeling strong interactions like body freedom flutter. A method of combining different sources of data traditionally used for aeroelasticity and flight dynamics is described by reconciling many of the differences between these models. By building upon past modeling efforts, a level of familiarity in the approach is achieved. Generally the differences from the traditional approaches are subtle but significant. The traditional frequency domain flutter model in a modal coordinate system is converted to a form consistent with a time domain flight dynamics model. The time domain rational function approximation about a non-inertial coordinate system and the unique constraints for the conversion between the inertial and non-inertial coordinate systems are discussed. A consistent transformation of the states of aeroelastic models to flight dynamics models is derived, which enables the integration of data from higher fidelity computational fluid dynamics models or wind-tunnel testing. The present method of integrating multidisciplinary data was used to create models that compare well with X-56A flight-test data, including conditions past the flutter speed.

Ouellette, Jeffrey↗