Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,117 records · Page 62

A covariance methodology for theoretically generated thermal neutron inelastic scattering data - 14502

The The S(α,β) double-differential thermal neutron scattering law tabulated in ENDF File 7 is, by convention, produced theoretically through fundamental scattering physics models for material-specific interatomic structures. Differential and integral inelastic scattering cross sections are calculated by numerically integrating over S(α,β). This information is supplied to nuclear system simulation codes. For systems with a significant thermal neutron flux, calculated system response parameters can be sensitive to this thermal cross section data. Formats and procedures are currently established for representing uncertainty information, or covariances, for various types of ENDF reaction data. However, no published ENDF evaluations contain covariance data for S(α,β) or associated differential and integral inelastic scattering cross sections. Furthermore, no accepted methodology exists for quantifying or representing covariances for this thermal scattering data. For most solids, the influence of a material's interatomic structure and dynamics on inelastic scattering cross sections can be adequately described through the phonon density of states (DOS). The phonon DOS is a probability density function of the population of available vibrational modes in energy, and it is typically the fundamental input for calculating S(α,β). Probable variation in the spectral shape of the phonon DOS may be established which characterizes uncertainties in the physics models and methodology employed in its production. Through Monte Carlo sampling of perturbations from the reference phonon spectrum, an S(α,β) covariance matrix may be generated. With appropriate sensitivity information, this can be propagated to generate covariance data for differential and integral inelastic scattering cross sections. In this work, hexagonal graphite is used as an example material for the purpose of demonstrating procedures for analyzing, calculating and representing uncertainty information for thermal neutron inelastic scattering data. (authors)

42 ENGINEERING↗

Coupled Inertial Navigation and Flush Air Data Sensing Algorithm for Atmosphere Estimation

This paper describes an algorithm for atmospheric state estimation that is based on a coupling between inertial navigation and flush air data sensing pressure measurements. In this approach, the full navigation state is used in the atmospheric estimation algorithm along with the pressure measurements and a model of the surface pressure distribution to directly estimate atmospheric winds and density using a nonlinear weighted least-squares algorithm. The approach uses a high fidelity model of atmosphere stored in table-look-up form, along with simplified models of that are propagated along the trajectory within the algorithm to provide prior estimates and covariances to aid the air data state solution. Thus, the method is essentially a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, descent, and landing from August 2012. Reasonable estimates of the atmosphere and winds are produced by the algorithm. The observability of winds along the trajectory are examined using an index based on the discrete-time observability Gramian and the pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding information content to the system. The algorithm is then applied to the design of the pressure measurement system for the Mars 2020 mission. The pressure port layout is optimized to maximize the observability of atmospheric states along the trajectory. Linear covariance analysis is performed to assess estimator performance for a given pressure measurement uncertainty. The results indicate that the new tightly-coupled estimator can produce enhanced estimates of atmospheric states when compared with existing algorithms.

Karlgaard, Christopher D.↗

Fabrication and characterization of Bismuth-Cerium composite iron garnet epitaxial films for magneto optical applications

The Bi{sub x}Ce{sub 3-x}Fe{sub 5}O{sub 12} (x = 0.8) epitaxial films of high quality were grown by means of pulsed laser deposition on paramagnetic substrates of Gadolinium Gallium Garnet. We study the modifications of substitutions in the parent garnet Y{sub 3}Fe{sub 5}O{sub 12} that produces a higher magneto-optical response at communication wavelengths. These films displayed a strong in plane textures which are treated in argon as well as reduced atmosphere conditions. The elemental constituents of these films were confirmed by energy dispersive-X ray analysis, elastic recoil detection system, Rutherford backscattering spectroscopy, and X-ray photoelectron spectroscopy measurements. The transmittance spectra were measured and found these films exhibit good transmittance values. The transmittance-spectra were fitted with the theoretical model and the optical constants such as refractive index and absorption edge were evaluated. The highest (negative) Faraday rotation was found for these films treated in the environment of Ar + H{sub 2}. A density matrix theory has been developed for the Faraday rotation and a good agreement between the theory and experiment is found. These epitaxial garnet films can be used in a wide range of frequencies from visible to infrared spectra making them ideal for many magneto optical applications. Therefore, these films may overcome many issues in fabricating all optical isolators which is the viable solution for integrated photonics.

36 MATERIALS SCIENCE↗

First-principles thermodynamics and defect kinetics guidelines for engineering a tailored RRAM device

Resistive Random Access Memories are among the most promising candidates for the next generation of non-volatile memory. Transition metal oxides such as HfOx and TaOx attracted a lot of attention due to their CMOS compatibility. Furthermore, these materials do not require the inclusion of extrinsic conducting defects since their operation is based on intrinsic ones (oxygen vacancies). Using Density Functional Theory, we evaluated the thermodynamics of the defects formation and the kinetics of diffusion of the conducting species active in transition metal oxide RRAM materials. The gained insights based on the thermodynamics in the Top Electrode, Insulating Matrix and Bottom Electrode and at the interfaces are used to design a proper defect reservoir, which is needed for a low-energy reliable switching device. The defect reservoir has also a direct impact on the retention of the Low Resistance State due to the resulting thermodynamic driving forces. The kinetics of the diffusing conducting defects in the Insulating Matrix determine the switching dynamics and resistance retention. The interface at the Bottom Electrode has a significant impact on the low-current operation and long endurance of the memory cell. Our first-principles findings are confirmed by experimental measurements on fabricated RRAM devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A neutronic feasibility study of the AP1000 design loaded with fully ceramic micro-encapsulated fuel

A neutronic feasibility study is performed to evaluate the utilization of fully ceramic microencapsulated (FCM) fuel in the AP1000 reactor design. The widely used Monte Carlo code MCNP is employed to perform the full core analysis at the beginning of cycle (BOC). Both the original AP1000 design and the modified design with the replacement of uranium dioxide fuel pellets with FCM fuel compacts are modeled and simulated for comparison. To retain the original excess reactivity, ranges of fuel particle packing fraction and fuel enrichment in the FCM fuel design are first determined. Within the determined ranges, the reactor control mechanism employed by the original design is directly used in the modified design and the utilization feasibility is evaluated. The worth of control of each type of fuel burnable absorber (discrete/integral fuel burnable absorbers and soluble boron in primary coolant) is calculated for each design and significant differences between the two designs are observed. Those differences are interpreted by the fundamental difference of the fuel form used in each design. Due to the usage of silicon carbide as the matrix material and the fuel particles fuel form in FCM fuel design, neutron slowing down capability is increased in the new design, leading to a much higher thermal spectrum than the original design. This results in different reactivity and fission power density distributions in each design. We conclude that a direct replacement of fuel pellets by the FCM fuel in the AP1000 cannot retain the original optimum reactor core performance. Necessary modifications of the core design should be done and the original control mechanism needs to be re-designed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Dry-pressed lithium nickel cobalt manganese oxide (NCM) cathodes enabled by holey graphene host

There is an increasing demand for electric vehicle batteries with higher energy densities that can be manufactured with low cost utilizing environmentally friendly processes. Lithium nickel cobalt manganese oxides (NCM), such as LiNi(0.5)Co(0.2)Mr(0.3)O(2) (NCM523), are commonly used in commercial lithium ion batteries for electric vehicle applications. These cathodes are typically fabricated using environmentally hazardous, energy intensive, and time-consuming solvent-based slurry casting processes. Polymeric binders, also required for cathode preparation, add weight and can undergo parasitic reactions during battery cycling. Here we report a facile, solvent-free, binder-free, and rapid cold-press electrode fabrication process to prepare high mass loading NCM cathodes within minutes. This process is enabled by the use of holey graphene, a unique lightweight, compressible carbon nanomaterial that serves as both a binder and an electrically conductive matrix. NCM cathodes fabricated by this method are stable and exhibit good active material utilization. A significant advantage for the dry-press process is that the electrode mass loadings can be readily increased simply by adding more composite powder to the compression molding die, resulting in improved areal capacities that are attractive for applications requiring high energy per unit footprint. (C) 2020 Elsevier Ltd. All rights reserved.

25 ENERGY STORAGE↗

Development of an in-situ multi-component reinforced Al-based metal matrix composite by direct metal laser sintering technique — Optimization of process parameters

In the present investigation, an in-situ multi-component reinforced aluminum based metal matrix composite was fabricated by the combination of self-propagating high-temperature synthesis and direct metal laser sintering process. The different mixtures of Al, TiO{sub 2} and B{sub 4}C powders were used to initiate and maintain the self-propagating high-temperature synthesis by laser during the sintering process. It was found from the X-ray diffraction analysis and scanning electron microscopy that the reinforcements like Al{sub 2}O{sub 3}, TiC, and TiB{sub 2} were formed in the composite. The scanning electron microscopy revealed the distribution of the reinforcement phases in the composite and phase identities. The variable parameters such as powder layer thickness, laser power, scanning speed, hatching distance and composition of the powder mixture were optimized for higher density, lower porosity and higher microhardness using Taguchi method. Experimental investigation shows that the density of the specimen mainly depends upon the hatching distance, composition and layer thickness. On the other hand, hatching distance, layer thickness and laser power are the significant parameters which influence the porosity. The composition, laser power and layer thickness are the key influencing parameters for microhardness. - Highlights: • The reinforcements such as Al{sub 2}O{sub 3}, TiC, and TiB{sub 2} were produced in Al-MMC through SHS. • The density is mainly influenced by the material composition and hatching distance. • Hatching distance is the major influencing parameter on porosity. • The material composition is the significant parameter to enhance the microhardness. • The SEM micrographs reveal the distribution of TiC, TiB{sub 2} and Al{sub 2}O{sub 3} in the composite.

36 MATERIALS SCIENCE↗

Knockout beyond the dripline

In this contribution, we will describe neutron and proton removal from {sup 9}C and {sup 7}Be which are two particularly interesting nuclei entering the nucleo-synthesis pp-chain [1, 2]. Neutron and proton removal reactions have been used in the past twenty years to probe the single-particle structure of exotic nuclei. The core parallel-momentum distribution can give information on the angular momentum and spin of the nucleon initial state while the total removal cross section is sensitive to the asymptotic part of the initial wave function and also to the reaction mechanism. Because knockout is a peripheral reaction from which the Asymptotic Normalization Constant (ANC) of the single-particle wave function can be extracted, it has been used as an indirect method to obtain the rate of reactions like {sup 8}B(p,γ){sup 9}C or {sup 7}Be(p,γ){sup 8}B. Nucleon removal has recently been applied by the HiRA collaboration [3] to situations in which the remaining “core” is beyond the drip line, such as {sup 8}C and {sup 6}Be, unbound by one or more protons, and whose excitation-energy spectrum can be obtained by the invariant-mass method. By gating on the ground-state peak, “core” parallel-momentum distributions and total knockout cross sections have been obtained similar to previous studies with well-bound “cores”. In addition for each projectile, knock out to final bound states has also been obtained in several cases. We will report on the theoretical description and comparison to this experimental data for a few cases for which advances in the accuracy of the transfer-to-the continuum model [4, 5] have been made [6]. These include the use, when available, of “ab-initio” overlaps for the initial state [7] and in particular their ANC values [8]. Also, the construction of a nucleus-target folding potential for the treatment of the core-target S-matrix [9] using for the cores “ab-initio” densities [10] and state-of-the-art n−{sup 9}Be optical potentials [11]. Preliminary results and open problems will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A peridynamic failure analysis of fiber-reinforced composite laminates using finite element discontinuous Galerkin approximations

Here, this paper presents a discontinuous Galerkin weak form for bond-based peridynamic models to predict the damage of fiber-reinforced composite laminates. To represent the anisotropy of a laminate in a peridynamic model, a lamina is simplified as a transversely isotropic medium under a plane stress condition. The laminated structure is modeled by stacking the surface mesh layers along the thickness direction according to the laminate sequence. To avoid a mesh dependence on either the fiber orientation or the discretization, the spherical harmonic expansion theory is employed to construct a function for the micro-elastic modulus in terms of the bond-fiber angle. The laminate material is decomposed into an isotropic matrix material part and a transversely isotropic material part. The bond stiffness can be evaluated using the engineering material constants, based on the equivalence between the elastic energy density in the peridynamic theory and the elastic energy density in the classic continuum mechanics theory. Benchmark tests are conducted to verify the proposed model. Numerical results illustrate that the convergence of simulations with different horizon sizes and meshes can be achieved. In terms of damage analysis, the proposed model can capture the dynamic process of the complex coupling of the inner-layer and delamination damage modes.

42 ENGINEERING↗

The Influence of Modification on Crystal Lattice Stability of Austenite in Stainless Steel

Using the methods of electron diffraction microscopy and X-ray diffraction analysis, the influence of alloying of the austenitic steel, Grade 110H13, with chromium and vanadium, as well as high-melting, ultrafine-grained TiO{sub 2}, ZrO{sub 2} powders and Na{sub 3}AlF{sub 6} cryolite on its structural-phase state and microstructure is investigated. It is shown that the matrix of non-modified steel is completely austenitic and consists of an iron-based solid solution and the interstitial (C, N, O and other) and substitutional (Cr, V and other) atoms simultaneously. Alloying with chromium and vanadium changes neither its phase composition nor defect structure, while alloy modification results in qualitatively new structural features: γ → ε-transformation, high-intensity microtwinning, defect structure changes, and a sharp increase in the scalar dislocation density. The features of the deformation-induced microtwinning and ε-martensite plates identified in the modified steel promote revealing additional microtwin systems in the matrix γ-phase, which result in structural changes making it possible to classify it as a γ′-phase. It is found out that an introduction of modifying additions gives rise to the following sequence of structural-phase transformations: γ→γ′→(γ′ +ε). The experimental data obtained demonstrate that as a result of modification the crystal lattice transits into a low-stability state. This transition is accompanied by marked structural-phase changes consisting in the formation of several microtwin systems and γ → ε-transformation. These structural-phase changes in the modified steel are due to the crystal-lattice transition into the low-stability state, followed by new structural-phase alterations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Morphology and Flow Behavior of Cellulose Nanofibers Dispersed in Glycols

Understanding the morphology and flow behavior of cellulose nanofibers (CNFs) dispersed in organic solvents can improve the process of fabricating new cellulose-based nanocomposites. In this study, jute-based 2,2,6,6-tetramethylpiperidinyl-1-oxyl (TEMPO)-oxidized CNFs with two different charge densities (0.64 and 1.03 mmol/g) were dispersed in ethylene glycol (EG) and propylene glycol (PG) using the solvent exchange method. The morphology and dimensions of CNFs in dry and suspension states were characterized using transmission electron microscopy, atomic force microscopy, and small-angle X-ray scattering techniques. The results showed that the cross-sectional dimensions remained the same in different solvents. Rheological measurements showed that CNF suspensions in water or glycol (EG and PG) behaved similar to typical polymer solutions with a solvent-independent overlap concentration corresponding to the crowding factor of about 14. Moreover, a thixotropic behavior was found in the concentrated CNF/glycol systems as determined in typical CNF aqueous suspensions. The fact that TEMPO-oxidized CNFs can be well dispersed in organic solvents opens up new possibilities to improve the CNF–polymer matrix blending, where the use of a viscous solvent can delay the transition to turbulence in processing and improve the control of fiber orientation because of a slower Brownian diffusive motion.

36 MATERIALS SCIENCE↗

Production of an impermeable composite of irradiated graphite and glass by hot isostatic pressing as a long term leach resistant waste form

Around 250,000 tons of irradiated graphite (i-graphite) exists worldwide and can be considered as a current waste or future waste stream. The largest national i-graphite inventory is located in UK (∼ 100,000 tons) with significant quantities also in Russia and France [5]. Most of the i-graphite remains in the cores of shutdown nuclear reactors including the MAGNOX type in UK and the UNGG in France. Whilst there are still operational power reactors with graphite cores, such as the Russian RBMKs and the AGRs in UK, all of them will reach their end of life during the next two decades. The most common reference waste management option of i-graphite is a wet or dry retrieval of the graphite blocks from the reactor core and the grouting of these blocks in a container without further conditioning. This produces large waste package volumes because the encapsulation capacity of the grout is limited and large cavities in the graphite blocks could reduce the packing densities. Packing densities from 0.5 to 1 tons per cubic meter have been assumed for grouting solutions. Furthermore the grout is permeable. This could over time allow the penetration of aqueous phases into the waste block and a potential dissolution and release of radionuclides. As a result particularly highly soluble radionuclides may not be retained by the grout. Vitrification could present an alternative, however a similar waste package volume increase may be expected since the encapsulation capacity of glass is potentially similar to or worse than that of grout. FNAG has developed a process for the production of a graphite-glass composite material called Impermeable Graphite Matrix (IGM) [3]. This process is also applicable to irradiated graphite which allows the manufacturing of an impermeable material without volume increase. Crushed i-graphite is mixed with 20 vol.% of glass and then pressed under vacuum at an elevated temperature in an axial hot vacuum press (HVP). The obtained product has zero or negligible porosity and a water impermeable structure. Structural analysis shows that the glass in the composite has replaced the pores in the graphite structure. The typical pore volume of a graphite material is in the range of 20 vol.%. Therefore no volume increase will occur in comparison with the former graphite material. This IGM material will allow the encapsulation of graphite with package densities larger than 1.5 ton per cubic meter. Therefore a huge volume saving can be achieved by such an alternative encapsulation method. Disposal performance is also enhanced since little or no leaching of radionuclides is observed due to the impermeability of the material NNL and FNAG have proved that IGM can be produced by hot isostatic pressing (HIP) which has several advantages for radioactive materials over the HVP process. - The sealed HIP container avoids the release of any radionuclides. - The outside of the waste package is not contaminated. - The HIP process time is shorter than the HVP process time. The isostatic press avoids anisotropic density distributions. - Simple filling of the HIP container has advantages over the filling of an axial die. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vibratory tactile display for textures

We have developed a tactile display that produces vibratory stimulus to a fingertip in contact with a vibrating tactor matrix. The display depicts tactile surface textures while the user is exploring a virtual object surface. A piezoelectric actuator drives the individual tactor in accordance with both the finger movement and the surface texture being traced. Spatiotemporal display control schemes were examined for presenting the fundamental surface texture elements. The temporal duration of vibratory stimulus was experimentally optimized to simulate the adaptation process of cutaneous sensation. The selected duration time for presenting a single line edge agreed with the time threshold of tactile sensation. Then spatial stimulus disposition schemes were discussed for representation of other edge shapes. As an alternative means not relying on amplitude control, a method of augmented duration at the edge was investigated. Spatial resolution of the display was measured for the lines presented both in perpendicular and parallel to a finger axis. Discrimination of texture density was also measured on random dot textures.

Ikei, Yasushi↗

FAST GENERATION OF ENSEMBLES OF COSMOLOGICAL N-BODY SIMULATIONS VIA MODE RESAMPLING

We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the large-scale (Gaussian-distributed) Fourier modes in a periodic N-body simulation box in a manner that properly accounts for the nonlinear mode coupling between large and small scales. We find that our method for adding new large-scale mode realizations recovers the nonlinear power spectrum to sub-percent accuracy on scales larger than about half the Nyquist frequency of the simulation box. Using 20 N-body simulations, we obtain a power spectrum covariance matrix estimate that matches the estimator from Takahashi et al. (from 5000 simulations) with <20% errors in all matrix elements. Comparing the rates of convergence, we determine that our algorithm requires {approx}8 times fewer simulations to achieve a given error tolerance in estimates of the power spectrum covariance matrix. The degree of success of our algorithm indicates that we understand the main physical processes that give rise to the correlations in the matter power spectrum. Namely, the large-scale Fourier modes modulate both the degree of structure growth through the variation in the effective local matter density and also the spatial frequency of small-scale perturbations through large-scale displacements. We expect our algorithm to be useful for noise modeling when constraining cosmological parameters from weak lensing (cosmic shear) and galaxy surveys, rescaling summary statistics of N-body simulations for new cosmological parameter values, and any applications where the influence of Fourier modes larger than the simulation size must be accounted for.

79 ASTRONOMY AND ASTROPHYSICS↗

Fermionic Partial Tomography via Classical Shadows

Here we propose a tomographic protocol for estimating any k-body reduced density matrix (k-RDM) of an n-mode fermionic state, a ubiquitous step in near-term quantum algorithms for simulating many-body physics, chemistry, and materials. Our approach extends the framework of classical shadows, a randomized approach to learning a collection of quantum-state properties, to the fermionic setting. Our sampling protocol uses randomized measurement settings generated by a discrete group of fermionic Gaussian unitaries, implementable with linear-depth circuits. We prove that estimating all k-RDM elements to additive precision ϵ requires on the order of ($^{n}_{k}$)k 3/2 log(n)/ϵ 2 repeated state preparations, which is optimal up to the logarithmic factor. Furthermore, numerical calculations show that our protocol offers a substantial improvement in constant overheads for k ≥ 2, as compared to prior deterministic strategies. We also adapt our method to particle-number symmetry, wherein the additional circuit depth may be halved at the cost of roughly 2–5 times more repetitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using Redundancy To Reduce Errors in Magnetometer Readings

A method of reducing errors in noisy magnetic-field measurements involves exploitation of redundancy in the readings of multiple magnetometers in a cluster. By "redundancy"is meant that the readings are not entirely independent of each other because the relationships among the magnetic-field components that one seeks to measure are governed by the fundamental laws of electromagnetism as expressed by Maxwell's equations. Assuming that the magnetometers are located outside a magnetic material, that the magnetic field is steady or quasi-steady, and that there are no electric currents flowing in or near the magnetometers, the applicable Maxwell 's equations are delta x B = 0 and delta(raised dot) B = 0, where B is the magnetic-flux-density vector. By suitable algebraic manipulation, these equations can be shown to impose three independent constraints on the values of the components of B at the various magnetometer positions. In general, the problem of reducing the errors in noisy measurements is one of finding a set of corrected values that minimize an error function. In the present method, the error function is formulated as (1) the sum of squares of the differences between the corrected and noisy measurement values plus (2) a sum of three terms, each comprising the product of a Lagrange multiplier and one of the three constraints. The partial derivatives of the error function with respect to the corrected magnetic-field component values and the Lagrange multipliers are set equal to zero, leading to a set of equations that can be put into matrix.vector form. The matrix can be inverted to solve for a vector that comprises the corrected magnetic-field component values and the Lagrange multipliers.

Kulikov, Igor↗

Simulating seawater intrusion in a complex coastal karst aquifer using an improved variable-density flow and solute transport–conduit flow process model

VDFST-CFP (variable-density flow and solute transport–conduit flow process) is a density-dependent discrete-continuum numerical model for simulating seawater intrusion in a dual-permeability coastal karst aquifer. A previous study (Xu and Hu 2017) simulates variable-density flow only in a single conduit, and studies the parameter sensitivities only in the horizontal case (2D domain as horizontal section) by the VDFST-CFP model. This paper focuses on the density-dependent vertical case (2D domain as vertical section) with two major improvements: 1) when implementing double-conduit networks in the domain, simulated intruded plumes in the porous medium are extended in the double-conduit scenario, compared to the single-conduit system; 2) by quantifying micro-textures on the conduit wall by the Goudar-Sonnad equation and considering macro-structures as local head loss. Sensitivity analysis shows that medium hydraulic conductivity, conduit diameter and effective porosity are important parameters for simulating seawater intrusion in the discrete-continuum system. On the other hand, rougher micro-structures and additional macro-structure components on the conduit wall would reduce the distance of seawater intrusion to the conduit system, but, rarely affect salinity distribution in the matrix. Compared to the equivalent mean roughness height, the new method (with more detailed description of structure) simulates seawater intrusion slightly landward in the conduit system. Here, the macro-structure measured by local head loss is more reasonable but needs further study on conduit flow. Xu and Hu (2017) Development of a discrete-continuum VDFST-CFP numerical model for simulating seawater intrusion to a coastal karst aquifer with a conduit system.

58 GEOSCIENCES↗

Effect of porosity, tissue density, and mechanical properties on radial sound speed in human cortical bone

Purpose: The purpose of this study was to investigate the effect of simultaneous changes in cortical porosity, tissue mineral density, and elastic properties on radial speed of sound (SOS) in cortical bone. The authors applied quantitative pulse-echo (PE) ultrasound techniques that hold much potential especially for screening of osteoporosis at primary healthcare facilities. Currently, most PE measurements of cortical thickness, a well-known indicator of fracture risk, use a predefined estimate for SOS in bone to calculate thickness. Due to variation of cortical bone porosity, the use of a constant SOS value propagates to an unknown error in cortical thickness assessment by PE ultrasound. Methods: The authors conducted 2.25 and 5.00 MHz focused PE ultrasound time of flight measurements on femoral diaphyses of 18 cadavers in vitro. Cortical porosities of the samples were determined using microcomputed tomography and related to SOS in the samples. Additionally, the effect of cortical bone porosity and mechanical properties of the calcified matrix on SOS was investigated using numerical finite difference time domain simulations. Results: Both experimental measurements and simulations demonstrated significant negative correlation between radial SOS and cortical porosity (R{sup 2} ≥ 0.493, p < 0.01 and R{sup 2} ≥ 0.989, p < 0.01, respectively). When a constant SOS was assumed for cortical bone, the error due to variation of cortical bone porosity (4.9%–16.4%) was about 6% in the cortical thickness assessment in vitro. Conclusions: Use of a predefined, constant value for radial SOS in cortical bone, i.e., neglecting the effect of measured variation in cortical porosity, propagated to an error of 6% in cortical thickness. This error can be critical as characteristic cortical thinning of 1.10% ± 1.06% per yr decreases bending strength of the distal radius and results in increased fragility in postmenopausal women. Provided that the cortical porosity can be estimated in vivo, the relationship between radial SOS and cortical porosity can be utilized and a porosity based radial SOS estimate could be implemented to determine cortical thickness. This would constitute a step toward individualized quantitative ultrasound diagnostics of osteoporosis.

60 APPLIED LIFE SCIENCES↗