Search NASA⌕ Search

SEARCH · Search NASA

Results for “Step height”

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 379 records · Page 21

Determining the Molecular Growth Mechanisms of Protein Crystal faces by Atomic Force Microscopy

A high resolution atomic force microscopy (AFM) study had shown that the molecular packing on the tetragonal lysozyme (110) face corresponded to only one of two possible packing arrangements, suggesting that growth layers on this face were of bimolecular height (Li et al., 1998). Theoretical analyses of the packing had also indicated that growth of this face should proceed by the addition of growth units of at least tetramer size corresponding to the 43 helices in the crystal. In this study an AFM linescan technique was devised to measure the dimensions of individual growth units on protein crystal faces. The growth process of tetragonal lysozyme crystals was slowed down by employing very low supersaturations. As a result images of individual growth events on the (110) face were observed, shown by jump discontinuities in the growth step in the linescan images. The growth unit dimension in the scanned direction was obtained by suitably averaging these images. A large number of scans in two directions on the (110) face were performed and the distribution of lysozyme aggregate sizes were obtained. A variety of growth units, all of which were 43 helical lysozyme aggregates, were shown to participate in the growth process with a 43 tetramer being the minimum observed size. This technique represents a new application for AFM allowing time resolved studies of molecular process to be carried out.

Li, Huayu↗

Scanning Backscatter Lidar Observations for Characterizing 4-D Cloud and Aerosol Fields to Improve Radiative Transfer Parameterizations

Clouds have a powerful influence on atmospheric radiative transfer and hence are crucial to understanding and interpreting the exchange of radiation between the Earth's surface, the atmosphere, and space. Because clouds are highly variable in space, time and physical makeup, it is important to be able to observe them in three dimensions (3-D) with sufficient resolution that the data can be used to generate and validate parameterizations of cloud fields at the resolution scale of global climate models (GCMs). Simulation of photon transport in three dimensionally inhomogeneous cloud fields show that spatial inhomogeneities tend to decrease cloud reflection and absorption and increase direct and diffuse transmission, Therefore it is an important task to characterize cloud spatial structures in three dimensions on the scale of GCM grid elements. In order to validate cloud parameterizations that represent the ensemble, or mean and variance of cloud properties within a GCM grid element, measurements of the parameters must be obtained on a much finer scale so that the statistics on those measurements are truly representative. High spatial sampling resolution is required, on the order of 1 km or less. Since the radiation fields respond almost instantaneously to changes in the cloud field, and clouds changes occur on scales of seconds and less when viewed on scales of approximately 100m, the temporal resolution of cloud properties should be measured and characterized on second time scales. GCM time steps are typically on the order of an hour, but in order to obtain sufficient statistical representations of cloud properties in the parameterizations that are used as model inputs, averaged values of cloud properties should be calculated on time scales on the order of 10-100 s. The Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE) provides exceptional temporal (100 ms) and spatial (30 m) resolution measurements of aerosol and cloud backscatter in three dimensions. HARLIE was used in a ground-based configuration in several recent field campaigns. Principal data products include aerosol backscatter profiles, boundary layer heights, entrainment zone thickness, cloud fraction as a function of altitude and horizontal wind vector profiles based on correlating the motions of clouds and aerosol structures across portions of the scan. Comparisons will be made between various cloud detecting instruments to develop a baseline performance metric.

Schwemmer, Geary K.↗

Wave Energy Power Take-off Validation with a Hydraulicly Actuated Rotary Dynamometer and a Bi-Directional High-Power DC Supply

There are many organizations working towards the commercialization of wave energy converter technologies and are advancing their designs through the technology readiness levels (TRLs). A critical step before the field deployment of prototype wave energy converters is the validation of the subsystems and components that are contained in the wave energy converter through laboratory testing and performance characterization. In 2021, the National Renewable Energy Laboratory (NREL) developed and demonstrated a system for testing power takeoffs (PTO) with a low-speed, high-torque dynamometer and a grid-tied high-power DC power source and sink before field deployment. The hydraulic dynamometer allows for the simulation of PTO actuation from wave motion and is capable of a wide range of wave periods and heights which are represented as various speeds and torques from the dynamometer. The high-power bidirectional power supply allows for hardware in the loop and controller in the loop testing to be conducted on WEC power electronics. This presentation was made to describe the methods used by NREL research staff to test all components and sub-systems in the PTO of a novel wave energy converter before field deployment.

dynamometer↗

Wave Energy Power Take-Off Validation with a Hydraulicly Actuated Rotary Dynamometer and a Bidirectional High-Power DC Supply: Preprint

There are many organizations working towards the commercialization of wave energy converter technologies and are advancing their designs through the technology readiness levels (TRLs). A critical step before the field deployment of prototype wave energy converters is the validation of the subsystems and components that are contained in the wave energy converter through laboratory testing and performance characterization. In 2021, the National Renewable Energy Laboratory (NREL) developed and demonstrated a system for testing power takeoffs (PTO) with a low-speed, high-torque dynamometer and a grid-tied high-power DC power source and sink before field deployment. The hydraulic dynamometer allows for the simulation of PTO actuation from wave motion and is capable of a wide range of wave periods and heights which are represented as various speeds and torques from the dynamometer. The high-power bidirectional power supply allows for hardware in the loop and controller in the loop testing to be conducted on WEC power electronics. This presentation was made to describe the methods used by NREL research staff to test all components and sub-systems in the PTO of a novel wave energy converter before field deployment.

dynamometer↗

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗

Reaction mechanisms of a cyclic ether intermediate: Ethyloxirane

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β -hydroperoxyalkyl radicals ( β -˙QOOH) or from reactions of HOÓ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n -butane oxidation. While rate coefficients for β -˙QOOH → ethyloxirane + ÓH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. Here we provide fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ˙R + O 2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ˙QOOH → cyclic ether + ÓH and ˙R + O 2 → conjugate alkene + HOÓ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ˙QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ˙QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Throughput Field Plant Phenotyping: A Self-Supervised Sequential CNN Method to Segment Overlapping Plants

High-throughput plant phenotyping—the use of imaging and remote sensing to record plant growth dynamics—is becoming more widely used. The first step in this process is typically plant segmentation, which requires a well-labeled training dataset to enable accurate segmentation of overlapping plants. However, preparing such training data is both time and labor intensive. To solve this problem, we propose a plant image processing pipeline using a self-supervised sequential convolutional neural network method for in-field phenotyping systems. This first step uses plant pixels from greenhouse images to segment nonoverlapping in-field plants in an early growth stage and then applies the segmentation results from those early-stage images as training data for the separation of plants at later growth stages. The proposed pipeline is efficient and self-supervising in the sense that no human-labeled data are needed. We then combine this approach with functional principal components analysis to reveal the relationship between the growth dynamics of plants and genotypes. We show that the proposed pipeline can accurately separate the pixels of foreground plants and estimate their heights when foreground and background plants overlap and can thus be used to efficiently assess the impact of treatments and genotypes on plant growth in a field environment by computer vision techniques. This approach should be useful for answering important scientific questions in the area of high-throughput phenotyping.

59 BASIC BIOLOGICAL SCIENCES↗

Determination of surface albedos and aerosol extinction characteristics from satellite imagery

A procedure is developed for calculating atmospheric extinction characteristics (optical thickness, scattering height, single scattering albedo) and surface albedo from radiometric images made at satellite altitudes. The procedure - a fast, computerized method - is suitable for the high-volume processing of satellite imagery data and thus can be used to map temporal and spatial distributions of aerosol parameters. Based on an analytical approximate solution to the equation of radiative transfer in a plane parallel atmosphere, the procedure is primarily applicable to the 0.4-micron to 0.8-micron wavelength range for solar zenith angles from 10 deg to 60 deg, surface albedos between 0.03 and about 0.5, and atmospheric optical thicknesses from 0.2 to 2. With a step-like change in the surface albedo (for example, a seashore or river bank or other similar change of terrain), both surface reflectivity and atmospheric optical thickness can be derived from radiance measurements. In this case, the resultant optical thickness is not based on a known surface reflectivity and is in essence independent of the radiometer calibration.

Kaufman, Y. J.↗

Requirements and Conceptual Design of Off-gas Systems for the Reprocessing of Metallic Fuels

An assessment has been conducted to determine how key regulations regarding volatile radionuclide emissions to the atmosphere may apply to the off-gas streams associated with electrochemical reprocessing. The scope of this assessment was based upon a generic electrochemical reprocessing scheme with a throughput rate of 200 MTIHM/y applied to metallic fuel discharged from a sodium fast reactor (SFR), but the findings are able to be translated to other advanced nuclear scenarios as merited. Air dispersion modeling was performed using the EPA CAP-88 model and evaluated the uncontrolled decontamination factors (DFs) that would be required to achieve regulatory compliance with the dose-based limits set forth by EPA regulation 40 CFR 190.10(a). These DFs were compared to those required by fuel cycle–based limits set forth by EPA regulation 40 CFR 190.10(b). Two theoretical sites with disparate climatological conditions were selected for air dispersion modeling (Idaho and Tennessee). The radionuclides modeled included 3 H, 85 Kr, 129 I, and selected alpha-emitting transuranic isotopes (referred to here as 239 Pu-TRU <1y ). It was found that the fuel cycle-based limits in 40 CFR 190.10(b) are most restrictive for 85 Kr and 239 Pu-TRU <1y , with DFs of 3 and 6.1E+09, respectively. The dose-based limit as derived from 40 CFR 190.10(a) could require mitigation of tritium in some scenarios, with an estimated DF of about 3 for the reference scenarios. The fuel cycle-based limit for 129 I resulted in a DF of about 240 for the reference scenario. The need for iodine mitigation based on dose to the public depended upon the physical form of iodine as either particulate or vapor-phase species. Emission of iodine from the facility as a vapor necessitated DFs of about 2 but emission as a particulate would require DFs >6,000 to meet thyroid dose-based limits. Effects of physical form on needed iodine mitigation are significant, but the understanding of speciation of iodine both during electrochemical reprocessing and after release to the atmosphere is limited. The electrochemical processing unit operations were evaluated to identify potential release points for the volatile radionuclides and to assess the potential for retention of the radionuclides within the process (thus decreasing the need for mitigation). Mitigation strategies for 3 H, 85 Kr, 129 I, and 239 Pu-TRU <1y were identified. In all cases, there are reasonably achievable pathways to regulatory compliance, although in some cases additional R&D is merited to verify the chemical speciation of these isotopes and to develop and demonstrate potential treatment technologies for this application. Whether or not additional off-gas controls (beyond common operations such as HEPA filtration and oxygen and moisture control) are needed for any of these regulated or volatile radionuclides depends on the (a) type of facility (NRC-regulated or DOE), (b) used fuel process rate, (c) used fuel burnup and composition, (d) speciation and retention of volatile radionuclides in the process and in the cell gas cleanup system, (e) site-specific parameters such as location, meteorology, stack height, and site boundaries, and (f) levels of conservatism and safety factors used in assessing compliance to air emissions regulations. Performance of this assessment revealed several areas where information is lacking or additional research is required in order to better determine if or what kinds of off-gas control might be needed. First, and most significantly, the understanding of the chemical speciation and physical form and partitioning of iodine during electrochemical processing operations is lacking and prevents the ability to accurately assess the potential iodine mitigation requirements. Future research in this area should be multifaceted and include thermodynamic modeling of iodine speciation in different process steps, experiments to quantify the kinetics of vapor-phase and melt-phase transitions, bench-scale experiments to determine the potential chemical and physical form of iodine emissions from the electrorefining process, and verification of iodine behavior with experiments utilizing operational facilities. Similarly, an improved understanding of iodine behavior in the environment after release from the facility stack will be required to refine dose estimations, as particulate and vapor-phase emissions result in significantly different doses to the MEI.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Understanding Uncertainty in Microstructure Evolution and Constitutive Properties in Additive Process Modeling

Coupled process–microstructure–property modeling, and understanding the sources of uncertainty and their propagation toward error in part property prediction, are key steps toward full utilization of additive manufacturing (AM) for predictable quality part development. The OpenFOAM model for process conditions, the ExaCA model for as-solidified grain structure, and the ExaConstit model for constitutive mechanical properties are used as part of the ExaAM modeling framework to examine a few of the various sources of uncertainty in the modeling workflow. In addition to “random” uncertainty (due to random number generation in the orientations and locations of grains present), the heterogeneous nucleation density N 0 and the mean substrate grain spacing S 0 are varied to examine their impact of grain area development as a function of build height in the simulated microstructure. While mean grain area after 1 mm of build is found to be sensitive to N 0 and S 0 , particularly at small N 0 and large S 0 (despite some convergence toward similar values), the resulting grain shapes and overall textures develop in a reasonably similar manner. As a result of these similar textures, ExaConstit simulation using ExaCA representative volume elements (RVEs) from various permutations of N 0 , S 0 , and location within the build resulted in similar yield stress, stress–strain curve shape, and stress triaxiality distributions. It is concluded that for this particular material and scan pattern, 15 layers is sufficient for ExaCA texture and ExaConstit predicted properties to become relatively independent of additional layer simulation, provided that reasonable estimates for N 0 and S 0 are used. However, additional layers of ExaCA will need to be run to obtain mean grain areas independent of build height and baseplate structure.

36 MATERIALS SCIENCE↗

Probabilistic Machine Learning Estimation of Ocean Mixed Layer Depth from Dense Satellite and Sparse In-Situ Observations

The ocean mixed layer plays an important role in the coupling between the upper ocean and atmosphere across a wide range of time scales. Estimation of the variability of the ocean mixed layer is therefore important for atmosphere-ocean prediction and analysis. The increasing coverage of in situ Argo profile data allows for an increasingly accurate analysis of the mixed layer depth (MLD) variability associated with deviations from the seasonal climatology. However, sampling rates are not sufficient to fully resolve subseasonal (<90 day) MLD variability. Yet, many multivariate observations-based analyses include implicit modeled subseasonal MLD variability. One analysis method is optimal interpolation of in situ data, but the interior analysis can be improved by leveraging surface data with regression or variational approaches. Here, we demonstrate how machine learning methods and satellite sea surface temperature, salinity, and height facilitate MLD estimation in a pilot study of two regions: the mid-latitude southern Indian and the eastern equatorial Pacific Oceans. We construct multiple machine learning architectures to produce weekly 1/2° gridded MLD anomaly fields (relative to a monthly climatology) with uncertainty estimates. We test multiple traditional and probabilistic machine learning techniques to compare both accuracy and probabilistic calibration. We validate our methodology by applying it to ocean model simulations. We find that incorporating sea surface data through a machine learning model improves the performance of spatiotemporal MLD variability estimation compared to optimal interpolation of Argo observations alone. These preliminary results are a promising first step for the application of machine learning to MLD prediction.

Machine Learning↗

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k‐ ε turbulence model, a complicated three‐dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

Lin, Shu-Chao↗

Mechanisms of dissolution from gibbsite step edges elucidated by ab initio molecular dynamics with enhanced sampling

Surface reactions underpin our collective understanding of mineral dissolution. This impacts a range of predictive geochemical models (for weathering, the fate and transport of metals) as well as industrial utilization of minerals. Yet atomistic details are rarely known due to the complex mineral/fluid interfacial environment. There is significant need to understand the mechanistic details of dissolution reactions and how they depend on surface morphology and solution conditions. This work utilizes surface pit models in conjunction with changes to solution composition that mimic pH to explore surface detachment of aluminate from gibbsite, which is a primary source of Al in soils and within the industrial processing of aluminum. Ab initio molecular dynamics simulations with enhanced sampling has been used to explore the detailed process of the detachment, the results of which indicate two potential pathways that are differentiated based upon the extent of water hydration. Here, the heights of the energy barriers depend upon the local morphology which influence the number of bridges (quasi-)simultaneously broken (1 or 2) or the Al-O coordination of the neighboring aluminum atoms (5 or 6) at the armchair edge. pH effects are significant, with a nearly 50% reduction in barrier height under alkaline conditions that are relevant to geothermal fluids and Al extraction from minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Spatial-Resolution OH and CH2O PLIF Visualization in a Dual-Mode Scramjet Combustor

A high-spatial-resolution planar laser-induced fluorescence (PLIF) imaging system was constructed and used to image a cavity-stabilized, premixed ethylene-air flame. The flame was created within a continuous flow, electrically-heated supersonic combustion facility consisting of a Mach 2 nozzle, an isolator with flush-wall fuel injectors, a combustor with a cavity flameholder of height 9 mm and optical access, and an extender. Tests were conducted at total temperature 1200 K, total pressure 300 kPa, equivalence ratio near 0.4 in the combustor, and Mach number near 0.6 in the combustor. A frequency-doubled Nd:YAG laser pumped a dye laser, which produced light at 283.55 nm. The beam was shaped into a light sheet with full width half-maximum 25 microns, which illuminated a streamwise plane that bisected the cavity. An intensified camera system imaged OH in this plane with a square 6.67 mm field of view and in-plane resolution 39 microns. Images were taken between the backward-facing step and 120 mm downstream of the step. OH structures as small as 110 microns were observed. CH2O was excited using 352.48 nm light; the smallest observed CH2O structures were approximately 200 microns wide. Approximately 15,000 images per species were processed and used to compute composite images.

Geipel, Clayton M.↗

Isomer-Dependent Reaction Mechanisms of Cyclic Ether Intermediates: cis -23-Dimethyloxirane and trans -23-Dimethyloxirane.

Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β-hydroperoxyalkyl radicals (β-˙QOOH) or from reactions of HÒO with alkenes. The cis- and trans-isomers of 2,3-dimethyloxirane are intermediates of n-butane oxidation, and while rate coefficients for β-˙QOOH → 2,3-dimethyloxirane + OH are reported extensively, subsequent reaction mechanisms of the cyclic ethers are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe the consumption of 2,3-dimethyloxirane by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present research examines the isomerdependence of 2,3-dimethyloxirane reaction mechanisms in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred via the detection of products from Cl-initiated oxidation of both cis-2,3-dimethyloxirane and trans-2,3-dimethyloxirane using multiplexed photoionization mass spectrometry (MPIMS). The experiments were conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, the enthalpies of stationary points on the ˙R + O 2 surfaces were computed at the ccCA-PS3 level of theory. In total, 28 barrier heights were computed on the 2,3-dimethyloxiranylperoxy surfaces. Two notable aspects are low-lying pathways that form resonance-stabilized ketohydroperoxide-type radicals caused by ˙QOOH ring-opening when the unpaired electron is localized adjacent to the ether group, and cis-trans isomerization of ˙R and ˙QOOH radicals, via inversion, which enable reaction pathways otherwise restricted by stereochemistry. Several species were identified in the MPIMS experiments from ring opening of 2,3-dimethyloxiranyl radicals. Neither of the two conjugate alkene isomers prototypical of ˙R + O 2 reactions were detected. Products were also identified from decomposition of ketohydroperoxide-type radicals. The present work provides the first analysis of 2,3-dimethyloxirane oxidation chemistry and reveals that consumption pathways are complex and require the expansion of submechanisms in chemical kinetics mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BSEC flux towers: CSAT3B and TRH

The data were collected as part of the BSEC project, during the period from June 2025 to May 2026. Directory "broadway" contains data collected on a multi-level flux tower (US-BWf) in the Broadway East neighborhood (1808 North Patterson Park Ave., Baltimore City, MD 21213; LAT: 39o18'40.31'' N; LONG: 76o35'12.43'' W). At each of the four measurement heights (8.5 m, 11.1 m, 13.4 m, 15.9 m), a Campbell Scientific CSAT3B sonic anemometer was operated at 50 Hz to measure virtual temperature (tc) and three velocity components (u: 270 degrees; v: 180 degrees; w: vertical), and a RM Young temperature sensor (model 41382VC) was operated at 1 Hz inside a compact aspirated radiation shield (model 43502) to measure absolute temperature (T) and relative humidity (RH). Inside directory "broadway", directory "netcdf" contains data collected each day in 5-minute chunks that have been converted to NetCDF format (before quality checking), while "4hr" contains data arranged into 4-hour chunks (also in NetCDF format) that have been through basic quality checking steps (treating data points with nonzero diagnostic codes as missing data; fixing six or fewer consecutive missing data points using linear interpolation). Users are recommended to start with data in directory "4hr", while data in directory "netcdf" can be used for reference purposes.

Baltimore↗

Comparison of Two Alternate Methods for Tracking Toe Clearance

Analyses of toe clearance during the swing phase of locomotion has often been utilized in determining a subject s propensity to trip while either walking or stepping over an obstacle. In the literature, toe clearance has been studied using a marker on the superior aspect of the second toe (rtoe), a marker on the lateral aspect of the fifth metatarsal head (mth5), or a virtual marker positioned at the anterior tip of the toe (vtoe). The purpose of this study was to compute toe clearance and associated parameters using a fifth metatarsal marker and a virtual toe marker, and compare the results with those of the standard toe marker. Subjects walked on a motorized treadmill at five different speeds while performing a visual acuity task at two separate target distances (ten 60-second trials). The minimum vertical height (TCl) was determined for each stride, along with its point of occurence in the gait cycle, and the angles of the foot and ankle at that time. A regression analysis was performed on the vtoe and mth5 results versus rtoe individually. For all TCl parameters, the mth5 marker did not correlate well with rtoe; the vtoe marker showed better agreement. Most importantly, the mth5 marker predicted a later occurence of TCl than rtoe and vtoe - thereby missing the most dangerous point in swing phase for a trip. From this analysis, the vtoe marker proved to be a better analog to rtoe than mth5, especially for determining a subject s propensity to trip.

Miller, Christopher A.↗

An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables

We derive a formulation of the nonhydrostatic equations in spherical geometry with a Lorenz staggered vertical discretization. The combination conserves a discrete energy in exact time integration when coupled with a mimetic horizontal discretization. The formulation is a version of Dubos and Tort (2014, https://doi.org/10.1175/MWR-D-14-00069.1) rewritten in terms of primitive variables. It is valid for terrain following mass or height coordinates and for both Eulerian or vertically Lagrangian discretizations. The discretization relies on an extension to Simmons and Burridge (1981, https://doi.org/10.1175/1520-0493(1981)109<0758:AEAAMC>2.0.CO;2) vertical differencing, which we show obeys a discrete derivative product rule. This product rule allows us to simplify the treatment of the vertical transport terms. Energy conservation is obtained via a term-by-term balance in the kinetic, internal, and potential energy budgets, ensuring an energy-consistent discretization up to time truncation error with no spurious sources of energy. We demonstrate convergence with respect to time truncation error in a spectral element code with a horizontal explicit vertically implicit implicit-explicit time stepping algorithm.

97 MATHEMATICS AND COMPUTING↗