SEARCH · Engineering Papers

Results for “Microscale 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1991

High altitude aircraft direction using real-time scientific analysis of telemetered data

Microscale and some mesoscale scientific targets are often too small and fast moving in a dynamic atmosphere. In this environment the ER-2 may spend only a small time sampling the phenomena of interest. The ER-2 scientists currently rely on post-flight analysis to determine the success of the mission and to provide input for planning the next flight. Real time analysis, combined with realistic modeling, can provide important clues to significantly enhance the probability of success in the investigation of specific scientific targets.

Wegener, Steven↗

NASA NTRS2017

Global Statistics of Microphysical Properties of Cloud-Top Ice Crystals

Ice crystals in clouds are highly complex. Their sizes, macroscale shape (i.e., habit), mesoscale shape (i.e., aspect ratio of components) and microscale shape (i.e., surface roughness) determine optical properties and affect physical properties such as fall speeds, growth rates and aggregation efficiency. Our current understanding on the formation and evolution of ice crystals under various conditions can be considered poor. Commonly, ice crystal size and shape are related to ambient temperature and humidity, but global observational statistics on the variation of ice crystal size and particularly shape have not been available. Here we show results of a project aiming to infer ice crystal size, shape and scattering properties from a combination of MODIS measurements and POLDER-PARASOL multi-angle polarimetry. The shape retrieval procedure infers the mean aspect ratios of components of ice crystals and the mean microscale surface roughness levels, which are quantifiable parameters that mostly affect the scattering properties, in contrast to a habit. We present global statistics on the variation of ice effective radius, component aspect ratio, microscale surface roughness and scattering asymmetry parameter as a function of cloud top temperature, latitude, location, cloud type, season, etc. Generally, with increasing height, sizes decrease, roughness increases, asymmetry parameters decrease and aspect ratios increase towards unity. Some systematic differences are observed for clouds warmer and colder than the homogeneous freezing level. Uncertainties in the retrievals will be discussed. These statistics can be used as observational targets for modeling efforts and to better constrain other satellite remote sensing applications and their uncertainties.

Van Diedenhoven, Bastiaan↗

DOE OSTI2026

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Paper or Plastic? Multiscale Material Handling Properties of Two Model Municipal Solid Waste Streams

Purpose: Municipal Solid Waste (MSW) is a potentially valuable sustainable feedstock for fuel and chemical production due to its carbon-rich content and low cost. This study aims to assess the material handling properties of paperand plastic-rich MSW feedstocks to mitigate equipment failure and processing downtime. Methods: The material handling properties of crumbled MSW feedstocks were measured using apowder rheometer with mass flow hopper calculations to assess handling performance. Inverse gas chromatography was use to measure the surface energy differences between feedstocks. Electron microscopy and Raman spectroscopy was used to evaluate microscale features that may contribute to material handling differences. Results: Plastic and paper rich feedstocks crumbled to a nominal 2 mm particle size were observed to have similar flow and handling characteristics with reasonable hopper outlets. 2 mm plastic rich crumbles, with their higher bulk density, exhibited superior flow performance. By contrast, 4 mm material required significantly larger hopper outlets, indicating poor flowability. Paper rich and 4 mm plastic rich samples displayed broad particle size distributions, which contributed to particle interlocking, jamming, and other flow issues. Electron microscopy revealed that plastic rich samples were significantly smoother, enhancing their flowability compared to the rougher, paper rich materials. Conclusions: This study establishes critical material handling baselines for processing MSW as a viable feedstock for fuel and chemical production. The findings highlight the importance of optimizing particle size and feedstock composition to improve flowability and handling performance.

09 BIOMASS FUELS↗

DOE OSTI2025

Forecasting Battery Electrode Performance via Electrochemical Fluorescence Microscopy and Machine-Learning

Predicting lithium-ion battery performance is hindered by microscale electrode heterogeneities invisible to conventional diagnostics. Here, we combine electrochemical fluorescence microscopy (EFM), which maps electronic connectivity by visualizing an electrofluorophore reaction distribution, with a multitask ElasticNet regression to forecast discharge capacity from spatial heterogeneity. Analyzing 196 images from six pilot-scale LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathodes with varying carbon loadings, we extract 62 descriptors that capture morphology and texture. A compact five-feature model predicts capacity across eight discharge rates, achieving a per-target R 2 of up to 0.63 and an overall R 2 of 0.92, with a mean absolute percentage error of less than 2%. This performance rivals impedance-based approaches while avoiding their reliance on postformation data and incomplete electronic network information. Our facile and rapid, image-driven method may enable electrode quality control upstream of costly cell assembly to offer a transformative tool for data-driven battery research and manufacturing.

battery electrodes↗

DOE OSTI2026

Electrospinning of ultra-high molecular weight polymers into aligned nanofibers and their application in chemical separations

Ultra-high molecular weight (UHMW, >10 6 g mol −1 ) vinyl polymers are increasingly accessible through photoiniferter polymerization, yet their behavior in electrospinning has received comparatively little attention. Here, we explore the electrospinning of UHMW poly(methyl acrylate), poly(methyl methacrylate), poly(N,N-dimethylacrylamide), and poly(styrene-co-pentafluorostyrene) synthesized by photoiniferter polymerization, with number-average molecular weights (M n ) greater than 1000 kDa and dispersities below 1.4. UHMW polymers form stable Taylor cones and discrete fibers at solution concentrations as low as 0.5 wt% due to their large number of entanglements, a threshold unachievable with low molecular weight (M n = 20–54 kDa) polymers of the same chemistry. Using poly(methyl methacrylate) and poly(styrene-co-pentafluorostyrene) as model systems, we demonstrate that fiber diameter can be tuned from 350 nm to 20 μm through systematic variation of polymer concentration (0.01 to 4 wt%) and flow rate (0.5 to 10 mL h −1 ). Additionally, we find that collector rotation rate governs fiber alignment, with highly aligned mats obtained with drum rotation speeds of 1000 RPM. Reversing the polarity of the electric field during electrospinning reorients the surface chemistry of these fiber mats without altering their microscale morphology, shifting the water contact angle by 7° (89° to 96°) and providing a two-fold binding capacity enhancement in cationic dye uptake. These results show that the high entanglement densities of UHMW vinyl polymers shift the accessible electrospinning window to dilute solutions and that electric-field polarity is a useful lever to control the surface chemistry in the resulting fibers.

Marquez, Joshua D. [University of Florida, Gainesv↗

NASA NTRS2018

Composite Overwrapped Pressure Vessel (COPV) Life Test

Recent attention has been focused on how NASA and its commercial partners have been placing less emphasis on testing and have become more dependent on analytical methods when evaluating design margins for fracture-critical components. Of high concern is the possible misuse of such analytic models for thin-walled metallic liners for composite overwrapped pressure vessels (COPVs). The NASA Engineering and Safety Center (NESC) has initiated an assessment to understand the limitations of linear-elastic fracture mechanics (LEFM) computational methods used to predict fatigue crack growth rate behavior of small cracks in thin metal liners for COPVs. It has been observed that fabrication of some thin metallic liners results in a wide variation in microstructure morphology, which results in varying microscale crack growth mechanisms. The ultimate goal is to develop and demonstrate a test-based methodology validating the safe-life requirements for COPVs with thin, elastically responding metal liners where LEFM methods are not appropriate.

Russell, Richard↗

DOE OSTI2024

Data Generation Workflow for Meso/Microscale Coupled Offshore Wind Farm Simulations

A robust, simple-to-use workflow is developed in this study which allows mesoscale information from the NOW-23 database to be easily incorporated into microscale wind farm simulations. This process will enable many different wind farm configurations to be simulated under realistic inflow conditions spanning a variety of atmospheric phenomena.

AMR-Wind↗

NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗

NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

processing↗

NASA NTRS2009

Multiscale Failure Analysis of Laminated Composite Panels Subjected to Blast Loading Using FEAMAC/Explicit

This preliminary report demonstrates the capabilities of the recently developed software implementation that links the Generalized Method of Cells to explicit finite element analysis by extending a previous development which tied the generalized method of cells to implicit finite elements. The multiscale framework, which uses explicit finite elements at the global-scale and the generalized method of cells at the microscale is detailed. This implementation is suitable for both dynamic mechanics problems and static problems exhibiting drastic and sudden changes in material properties, which often encounter convergence issues with commercial implicit solvers. Progressive failure analysis of stiffened and un-stiffened fiber-reinforced laminates subjected to normal blast pressure loads was performed and is used to demonstrate the capabilities of this framework. The focus of this report is to document the development of the software implementation; thus, no comparison between the results of the models and experimental data is drawn. However, the validity of the results are assessed qualitatively through the observation of failure paths, stress contours, and the distribution of system energies.

Pineda, Evan J.↗

NASA NTRS2026

Proactive Wildfire Management: A Remote Sensing and Multimodal CNN-MLP Architecture for Ignition Risk Forecasting

As the frequency and intensity of wildfires increase, with fire seasons now starting earlier and ending later than they have over the past decades, current monitoring systems, such as lookout towers and satellites, are hindered by cloud cover, low-resolution imagery, and static data gaps that fail to track vegetation moisture levels fast enough to catch rapid pre-ignition changes. This report proposes a Machine Learning-enabled Wildfire Ignition Prediction framework that combines satellite monitoring with dynamic and high-resolution remote sensing from Unmanned Aerial Vehicle (UAV) swarms. The method would use multispectral and thermal data from the Landsat program to create a baseline for vegetation health, calculating a two-band Enhanced Vegetation Index (EVI2) and the moisture content of the vegetation. These inputs will later be fused with microscale UAV weather data, including thermal hotspots found through thick canopies, hyperspectral chemical signatures of pre-visual combustion, and local weather streams. The multispectral satellite, multispectral Light Detection and Ranging (LiDAR), and thermal data would then be processed through a Convolutional Neural Network (CNN), alongside a Multilayer Perceptron (MLP) for the micro-weather telemetry. The outputs of these networks would be fused into a single feature representation and passed through a final prediction network to generate real-time ignition risk scores and hotspot alerts. Model performance would be assessed using standard classification metrics, including a Receiver Operating Characteristic - Area Under the Curve (ROC AUC) and F1 score. This system would allow first responders to identify high-risk zones and intervene before ignition occurs, improving emergency response time compared to current approaches.

machine learning↗

DOE OSTI2024

Increasing aggregate size reduces single-cell organic carbon incorporation by hydrogel-embedded wetland microbes

Abstract Microbial degradation of organic carbon in sediments is impacted by the availability of oxygen and substrates for growth. To better understand how particle size and redox zonation impact microbial organic carbon incorporation, techniques that maintain spatial information are necessary to quantify elemental cycling at the microscale. In this study, we produced hydrogel microspheres of various diameters (100, 250, and 500 μm) and inoculated them with an aerobic heterotrophic bacterium isolated from a freshwater wetland (Flavobacterium sp.), and in a second experiment with a microbial community from an urban lacustrine wetland. The hydrogel-embedded microbial populations were incubated with 13C-labeled substrates to quantify organic carbon incorporation into biomass via nanoSIMS. Additionally, luminescent nanosensors enabled spatially explicit measurements of oxygen concentrations inside the microspheres. The experimental data were then incorporated into a reactive-transport model to project long-term steady-state conditions. Smaller (100 μm) particles exhibited the highest microbial cell-specific growth per volume, but also showed higher absolute activity near the surface compared to the larger particles (250 and 500 μm). The experimental results and computational models demonstrate that organic carbon availability was not high enough to allow steep oxygen gradients and as a result, all particle sizes remained well-oxygenated. Our study provides a foundational framework for future studies investigating spatially dependent microbial activity in aggregates using isotopically labeled substrates to quantify growth.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

NASA NTRS1990

Pressure-balanced structures between 1 AU and 24 AU and their implications for solar wind electrons and interstellar pickup ions

'Pressure-balanced structures' (PBSs) in the heliosphere are microscale phenomena with a characteristic length along the radial dimension of the order of 0.05 AU across which the sum of the thermal pressures and the magnetic pressure is constant. PBSs have been identified in the Voyager data from 1 AU to 24 AU. If the density of pickup ions were proportional to the proton density across a PBS, and if the electron pressure were negligible, then from measurements of the magnetic field and the proton density and temperature across pressure balanced structures a pickup ion density is derived which is significantly smaller than the density predicted by a model based on the same assumption. This suggests that there are processes which 'smear out' the pressure of the pickup ions so that it is nearly constant on the scale of a PBS. If the pressure of the pickup ions is uniform across a PBS, then one can determine the electron temperature from the pressure balance condition.

Burlaga, L. F.↗

NASA NTRS2024

Global Sensitivity Analysis of Simulated Remote Sensing Polarimetric Observations Over Snow

This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors. The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum. The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance. The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration. The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.

remote sensing↗

NASA NTRS2003

DNS and LES of a Shear-Free Mixing Layer

The purpose of this work is twofold. First, given the computational resources available today, it is possible to reach, using DNS, higher Reynolds numbers than in Briggs et al.. In the present study, the microscale Reynolds numbers reached in the low- and high-energy homogeneous regions are, respectively, 32 and 69. The results reported earlier can thus be complemented and their robustness in the presence of increased turbulence studied. The second aim of this work is to perform a detailed and documented LES of the shear-free mixing layer. In that respect, the creation of a DNS database at higher Reynolds number is necessary in order to make meaningful LES assessments. From the point of view of LES, the shear-free mixing-layer is interesting since it allows one to test how traditional LES models perform in the presence of an inhomogeneity without having to deal with difficult numerical issues. Indeed, as argued in Briggs et al., it is possible to use a spectral code to study the shear-free mixing layer and one can thus focus on the accuracy of the modelling while avoiding contamination of the results by commutation errors etc. This paper is organized as follows. First we detail the initialization procedure used in the simulation. Since the flow is not statistically stationary, this initialization procedure has a fairly strong influence on the evolution. Although we will focus here on the shear-free mixing layer, the method proposed in the present work can easily be used for other flows with one inhomogeneous direction. The next section of the article is devoted to the description of the DNS. All the relevant parameters are listed and comparison with the Veeravalli & Warhaft experiment is performed. The section on the LES of the shear-free mixing layer follows. A detailed comparison between the filtered DNS data and the LES predictions is presented. It is shown that simple eddy viscosity models perform very well for the present test case, most probably because the flow seems to be almost isotropic in the small-scale range that is not resolved by the LES.

Knaepen, B.↗

NASA NTRS1980

A magnetohydrodynamic model for corotating interplanetary structures

A model is developed which treats interplanetary structures as a small-scale incompressible Alfvenic fluctuation superposed on a large-scale background flow, and uses the complete magnetohydrodynamic equations and the full energy equation to describe the interplanetary structures. For steady-state corotating structures the model assumes the background flow to be field aligned, and the nonlinear system is treated as being composed of two subsets: one dealing with the expansion of the solar wind in stream tubes, the other with the interaction between neighboring stream tubes. The expansion subset has a singularity at Parker's critical radius, and the interaction subset has a singularity at the Alfvenic point. Concerning the region where the flow velocity is super-Alfvenic, the characteristic equation is derived, and a general method of numerical solution for the three-dimensional problems based on the characteristics theory is discussed. In addition, a viscous WKB solution of the microscale system is presented.

Whang, Y. C.↗