Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model weights”

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,603 records · Page 89

Vibroacoustic Tailoring of a Rod-Stiffened Composite Fuselage Panel with Multidisciplinary Considerations

An efficient multi-objective design tailoring procedure seeking to improve the vibroacoustic performance of a fuselage panel while maintaining or reducing weight is presented. The structure considered is the pultruded rod stitched efficient unitized structure, a highly integrated composite structure concept designed for a noncylindrical, next-generation flight vehicle fuselage. Modifications to a baseline design are evaluated within a six-parameter design space including spacing, flange width, and web height for both frame and stringer substructure components. The change in sound power radiation attributed to a design change is predicted using finite-element models sized and meshed for analyses in the 500 Hz, 1 kHz, and 2 kHz octave bands. Three design studies are carried out in parallel while considering a diffuse acoustic field excitation and two types of turbulent boundary-layer excitation. Kriging surrogate models are used to reduce the computational costs of resolving the vibroacoustic and weight objective Pareto fronts. The resulting Pareto optimal designs are then evaluated under a static pressurization ultimate load to assess structural strength and stability. Results suggest that choosing alternative configurations within the considered design space can reduce weight and improve vibroacoustic performance without compromising strength and stability of the structure under the static load condition considered, but the tradeoffs are significantly influenced by the spatial characteristics of the assumed excitation field.

Allen, Albert R.↗

Development and Implementation of an Automatic Continuous Online Monitoring and Control Platform for Polymerization Reactions to Sharply Boost Energy and Resource Efficiency in Polymer Manufacturing

The project goal was to create an energy saving paradigm shift in how polymers are manufactured in the 21st century. It used Automatic Continuous Online Monitoring of Polymerization reactions (ACOMP) integrated for the first time with automatic active control to create the innovative ‘ACOMP/Control Interface’, or ‘ACOMP/CI’. ACOMP/CI will begin the transformation from old, inefficient processes into highly evolved, energy and resource efficient ones. The ACOMP platform is broadly applicable to many types of reactions and processes throughout the vast polymer industry. The industry provides materials for sectors such as automotive, aerospace, oil recovery, agriculture, paints, resins, adhesives, pharmaceuticals and therapeutic proteins, optics, electronics, lightweight building materials, and many more. The U.S. chemical industry is one of the last major sectors in which the U.S. has top global stature. It consumes 24% of all U.S. manufacturing energy, produces over $800B of product annually, supports 25% of the U.S. GDP and employs over 6 million people. It is also a major source of GHG emissions. Polymers make up approximately 30% of this sector. It is estimated that annually 60 TBtu of energy could be saved and 3 million tons less of GHG emissions produced by optimizing production in the polyolefin manufacturing sector alone. The project scope included first time design and prototyping of an ACOMP/CI, creation of active reaction controllers, and demonstration of control capabilities on ideal, low concentration polymerization reactions. All these elements of the scope were met, including advances and findings not originally anticipated. Extensions to more complex reactions, beyond the reactor capabilities of the current project ACOMP/CI, such as polyolefins and other high pressure/high temperature reactions, are being proposed in Fall 2017 to CESMII, a DoE based NNMI. The initial proposal was for a three year funded project, but this was reduced to a two year project and budget due to funding constraints. Hence, some of the original plans, such as adaptation of the ACOMP/CI to more relevant industrial processes, such as emulsion and dispersion technologies, could not be carried out. A third year of funding was requested at the end of the project, but DoE did not have resources to grant this. The sub-contractor Fluence Analytics (previously Advanced Polymer Monitoring Technologies, Inc) designed, prototyped, and commissioned a working ACOMP/CI by June 2015. The reaction characteristics to be automatically controlled were i) conversion kinetics, ii) molecular weight, iii) copolymer composition, and iv) simultaneous molecular weight and composition. A two pronged control strategy was used. The Tulane/Fluence group took a basic principles approach that did not rely on kinetic models. The LSU group took a more complex, non-linear model-oriented approach involving complete kinetic descriptions of the reaction system Each of these approaches proved successful in their own way. By April 2016 fully automatic control of conversion and weight average molecular weight, Mw, trajectories was achieved using the Tulane/Fluence (TF) basic principles controller. Similar results were obtained by the LSU non-linear model controller by August 2016. The demonstration system was aqueous free radical polymerization of acrylamide, Am. The control variables were temperature and semi-batch feed to the reactor of Am monomer and initiator. A demonstration of active manual conversion control in an industrial process using high solids in inverse emulsion polymerization of Am was achieved. During Summer 2016 the TF controller was used in conjunction with a chain transfer agent, another control variable, to automatically produce multi-modal molecular weight distributions, MWD, in a single reactor. Industrially, multi-modal MWD are produced by mixing products made in separate reactors, requiring significant extra time, energy, and reactor resources. Recognizing the industrial potential a patent on automatic production of multi-modal polymers was filed, and DoE acknowledged. In Fall 2016 the TF team developed a basic principles controller for copolymer composition and demonstrated it on aqueous free radical copolymerization of comonomers Am and styrene sulfonate, SS. TF then fused the Mw and composition controllers to achieve simultaneous control of both Mw and copolymer composition trajectories. Numerous simultaneous trajectories were demonstrated, including a trimodal composition distribution with constant Mw. Meanwhile, the LSU group developed a Kalman filter to improve the results of their automatic Mw and conversion controller and successful tests were made. During the project the TF team developed a means of computing full MWD during polymer synthesis without need for any chromatographic separation, based on model distributions. This means the polymer product is ‘born characterized’ and this can eliminate post-manufacture analytical laboratory quality control. TF filed a joint patent application on this new approach to chromatography-free determination of MWD with acknowledgment to DoE. The Tulane group obtained a 60MHz NMR during the project and recently completed the first work on separating three comonomers, Am, SS, and Na-acrylate, with a first demonstration of terpolymer composition control with the TF basic principles controller. Widespread dissemination of ACOMP/CI in the polymer manufacturing sector will bolster DoE goals of energy efficiency and reduced GHG emissions: The ability to monitor and actively control polymerization reactions will lead to more efficient use of energy and non-renewable resources, plant and labor time, increase the safety of manufacturing personnel, and will enhance product quality and lead to feasibility of manufacturing of polymers currently too complex for industrial scale production, while leading to less GHG emissions per kilo of product, and allowing for increased U.S. competitiveness in this enormous manufacturing sector. When ACOMP/CI is expanded to the polyolefin industry it is estimated that 60 TeraBTU/year of energy can be saved. Much of this saving is anticipated to come from optimized control of grade changeovers in steady state reactors and maintenance of steady states. Conclusions: ACOMP’s ability to provide continuous realtime data streams of measured polymer and reaction characteristics made it possible, for the first time, to directly and automatically control free radical polymerization reactions. An industrial client of Fluence Analytics has requisitioned the first ACOMP/CI which uses the TF basic principles controller. This sets the stage for FA to add control features to the ACOMP systems it has begun to install on the industrial scale beginning in 2014. Recommendations: This successful project was mainly limited to ideal polymerizations not of an industrial sort. The most energy intensive portion of polymer manufacturing is polyolefins. Adoption of ACOMP/CI to this enormous industrial sector faces the enormous challenges of high temperature, high pressure continuous sampling and high temperature sensor operation to obtain the continuous data needed for direct reaction control. The project team has a strategy for achieving this ambitious goal and will present it in Fall 2017 as a proposal to CESMII/DoE. It is recommended that this upcoming proposal be funded in order to make full use of the achievements of this just ended DoE project as the next step towards making polyolefin ACOMP/CI an energy saving reality. It is projected that ACOMP/CI can have its first polyolefin testbed demonstrations within two years of beginning the proposed project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Impact of Advection Schemes on Tracer Interrelationships in Large-Eddy Simulations of Deep Convection

This study investigates the preservation of tracer interrelationships during advection in large-eddy simulations of an idealized deep convective cloud, which is particularly relevant to chemistry, aerosol, and cloud microphysics models. Employing the Cloud Model 1, advection is represented using third-, fifth-, and seventh-order weighted essentially non-oscillatory schemes. As a simplified analogy for cloud hydrometeors and aerosols, several inert passive tracers following linear and nonlinear relationships are initialized after the cloud reaches ~6-km depth. Numerical mixing in the simulated turbulent convective clouds leads to significant deviations from the initial nonlinear relationships between tracers. In these simulations, a considerable fraction of the grid points where the tracers’ nonlinear relationships are altered from advection are classified as unrealistic (e.g., ~13% for the environmental tracers on average), including errors from range-preserving unmixing and overshooting. Errors in the sum of three tracers are also relatively large, ranging between ~1% and 16% for 5% of the grid points in and near the cloud. The magnitude of unrealistic mixing and errors in the sum of three tracers generally increase with the order of accuracy of the advection scheme. These results are consistent across model grid spacings ranging from 50 to 200 m, and across three different flow realizations for each combination of grid spacing and advection scheme tested. Tests employing a previously proposed scalar normalization procedure show substantially reduced errors in the sum of three tracers with a relatively small negative impact on other tracer relationships. This analysis, therefore, suggests efficacy of the normalization procedure when applied to turbulent three-dimensional cloud simulations.

54 ENVIRONMENTAL SCIENCES↗

The use of high temperature superconductors to levitate lunar telescope

The objective of this paper was to assist in the construction of a lunar telescope mirror model by conducting research on composite materials and other lightweight, rigid materials, and by determining how much weight can be levitated by available superconductors. It is believed that with the construction of four magnets suspended over four bulk superconductors (or vice versa), there should be no problems lifting a model mirror and stabilizing it at different positions. It may be necessary to increase the size and quality of the superconductors and/or magnets in order to achieve this.

Brown, Beth A.↗

Magmatic volatiles and the weathering of Mars

The sources for volatiles on Mars have been the subject of many hypotheses for exogenous influences including late accretion of volatile-enriched material, impact devolatilization to create massive early atmospheres, and even major bombardment by comets. However, the inventory of chemically active volatiles observable at the contemporary surface of Mars is consistent with domination by endogenous, subsequent planetary processes, viz., persistent magmatic outgassing. Volcanism on Mars has been widespread in both space and time. Notwithstanding important specific differences between the mantles of Earth and Mars, the geochemical similarities are such that the suite of gases emitted from Martian volcanic activity should include H2O, CO2, S-containing gases (e.g. H2S and/or SO2), and Cl-containing gases (e.g., Cl2 and/or HCl). H2O and CO2 exist in the atmosphere of Mars. Both are also present as surface condensates. However, spectroscopic observations of the Martian atmosphere clearly show that the S- and Cl-containing gases are severely depleted, with upper limits of less than or equal to 10(exp -7) the abundance of CO2. Likewise, there is no evidence of polar condensates of compounds of these elements as there is for CO2 and H2O. Within the soil, on the other hand, there has been direct measurement of incorporated H2O and abundant compounds containing S and Cl. Barring some as yet implausible geochemical sequestering process, the S/Cl ratio of about 6:1 in Martian soils implies a limit of 5% on the contribution of matter of solarlike composition (e.g., carbonaceous chondrite or cometary material) to these volatiles. Hence, exogenous sources are minor or not yet observed. From analysis of elemental trends in Martian soils, it has been recently shown that a simple two-component model can satisfy the Viking in situ measurements. Component A includes Si and most or all the Al, Ca, Ti, and Fe. Component B, taken as 16 +/- 3% by weight of the total, contains S and most or all the Cl and Mg. These results constrain several models of Martian soil mineralogy but are consistent with a mixture of silicates (such as Fe-rich clays and accessory minerals and soluble salts). The overall element profile is notably like shergottites, with significant incorporation of chemically reactive atmospheric gases from magmatic degassing.

Clark, B. C.↗

MULTIFRACTAL STRUCTURES DETECTED BY VOYAGER 1 AT THE HELIOSPHERIC BOUNDARIES

To better understand the dynamics of turbulent systems, we have proposed a phenomenological model based on a generalized Cantor set with two rescaling and one weight parameters. In this Letter, using recent Voyager 1 magnetic field data, we extend our two-scale multifractal analysis further in the heliosheath beyond the heliospheric termination shock, and even now near the heliopause, when entering the interstellar medium for the first time in human history. We have identified the scaling inertial region for magnetized heliospheric plasma between the termination shock and the heliopause. We also show that the degree of multifractality decreases with the heliocentric distance and is still modulated by the phases of the solar cycle in the entire heliosphere including the heliosheath. Moreover, we observe the change of scaling toward a nonintermittent (nonmultifractal) behavior in the nearby interstellar medium, just beyond the heliopause. We argue that this loss of multifractal behavior could be a signature of the expected crossing of the heliopause by Voyager 2 in the near future. The results obtained demonstrate that our phenomenological multifractal model exhibits some properties of intermittent turbulence in the solar system plasmas, and we hope that it could shed light on universal characteristics of turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Gaussian FLOWERS: Wind-rose-based analytical integration of Gaussian wake model for extremely fast AEP estimation

A major cost in the study of wind farm layout optimization is the repeated evaluation of the annual energy production (AEP). The current approach to estimating AEP requires a large set of flow simulations to be performed that cover each discrete wind speed and direction combination contained within the wind rose, followed by a probability-weighted sum of the power production resulting from each simulation. Even with inexpensive engineering wake models, this numerical integration scheme can lead to high computational costs. In this paper, we derive an analytical formulation for estimating farm AEP across every wind direction, based on a Gaussian wake velocity model, which reduces the number of wind farm simulations to a single function evaluation. As a result, we find that the Gaussian-FLOWERS approach reduces the time for AEP calculations by more than two orders of magnitude with a small trade-off in accuracy when compared to a conventional approach. This massive reduction in computation cost is useful to reduce overall costs in wind farm layout optimization studies.

17 WIND ENERGY↗

Weak-charge form-factor determination at the electron-ion collider

Determining the weak charge form factor, 𝐹 𝑊 ⁡(𝑄 2 ), of nuclei over a continuous range of momentum transfers, 0 ≲ 𝑄 2 ≲ 0.1 GeV 2 , is essential for mapping out the distribution of neutrons in nuclei. The neutron density distribution has significant implications for a broad range of areas, including studies of nuclear structure, neutron stars, and physics beyond the Standard Model. Currently, our knowledge of 𝐹 𝑊 ⁡(𝑄 2 ) comes primarily from fixed target experiments that measure the parity-violating asymmetry in coherent elastic electron-ion scattering. Fixed target experiments, such as CREX and PREX-1,2, have provided high-precision weak charge form factor extractions for the 48 Ca and 208 Pb nuclei, respectively. However, a major limitation of fixed target experiments is that they each provide data only at a single value of 𝑄 2 . With the proposed electron-ion collider (EIC) on the horizon, we explore its potential to impact the determination of the weak charge form factor. While it cannot compete with the precision of fixed target experiments, it can provide data over a wide and continuous range of 𝑄 2 values, and for a wide variety of nuclei. We show that with data corresponding to an integrated luminosity of ℒ ∼ 500/𝐴 fb −1 , where 𝐴 is the nucleus atomic weight, the EIC can significantly impact constraints by lifting degeneracies in theoretical models of the neutron density distribution. Ensuring EIC detector coverage at low 𝑄 2 and large negative pseudorapidities will be essential for such 𝐹 𝑊 ⁡(𝑄 2 ) measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Monte Carlo-based simulation of dynamic jaws tomotherapy

Purpose: Original TomoTherapy systems may involve a trade-off between conformity and treatment speed, the user being limited to three slice widths (1.0, 2.5, and 5.0 cm). This could be overcome by allowing the jaws to define arbitrary fields, including very small slice widths (<1 cm), which are challenging for a beam model. The aim of this work was to incorporate the dynamic jaws feature into a Monte Carlo (MC) model called TomoPen, based on the MC code PENELOPE, previously validated for the original TomoTherapy system. Methods: To keep the general structure of TomoPen and its efficiency, the simulation strategy introduces several techniques: (1) weight modifiers to account for any jaw settings using only the 5 cm phase-space file; (2) a simplified MC based model called FastStatic to compute the modifiers faster than pure MC; (3) actual simulation of dynamic jaws. Weight modifiers computed with both FastStatic and pure MC were compared. Dynamic jaws simulations were compared with the convolution/superposition (C/S) of TomoTherapy in the ''cheese'' phantom for a plan with two targets longitudinally separated by a gap of 3 cm. Optimization was performed in two modes: asymmetric jaws-constant couch speed (''running start stop,'' RSS) and symmetric jaws-variable couch speed (''symmetric running start stop,'' SRSS). Measurements with EDR2 films were also performed for RSS for the formal validation of TomoPen with dynamic jaws. Results: Weight modifiers computed with FastStatic were equivalent to pure MC within statistical uncertainties (0.5% for three standard deviations). Excellent agreement was achieved between TomoPen and C/S for both asymmetric jaw opening/constant couch speed and symmetric jaw opening/variable couch speed, with deviations well within 2%/2 mm. For RSS procedure, agreement between C/S and measurements was within 2%/2 mm for 95% of the points and 3%/3 mm for 98% of the points, where dose is greater than 30% of the prescription dose (gamma analysis). Dose profiles acquired in transverse and longitudinal directions through the center of the phantom were also compared with excellent agreement (2%/2 mm) between all modalities. Conclusions: The combination of weights modifiers and interpolation allowed implementing efficiently dynamic jaws and dynamic couch features into TomoPen at a minimal cost in terms of efficiency (simulation around 8 h on a single CPU).

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Evaluation of CMIP6 GCMs Over the CONUS for Downscaling Studies

Despite the necessity of Global Climate Models (GCMs) sub-selection in downscaling studies, an objective approach for their selection is currently lacking. Building on the previously established concepts in GCMs evaluation frameworks, we develop a weighted averaging technique to remove the redundancy in the evaluation criteria and rank 37 GCMs from the sixth phase of the Coupled Models Intercomparison Project over the contiguous United States. GCMs are rated based on their average performance across 66 evaluation measures in the historical period (1981–2014) after each metric is weighted between zero and one, depending on its uniqueness. The robustness of the outcome is tested by repeating the process with the empirical orthogonal function analysis in which each GCM is ranked based on its sum of distances from the reference in the principal component space. The two methodologies work in contrasting ways to remove the metrics redundancy but eventually develop similar GCMs rankings. A disparity in GCMs' behavior related to their sensitivity to the size of the evaluation suite is observed, highlighting the need for comprehensive multi-variable GCMs evaluation at varying timescales for determining their skillfulness over a region. The sub-selection goal is to use a representative set of skillful models over the region of interest without substantial overlap in their future climate responses and modeling errors in representing historical climate. Additional analyses of GCMs' independence and spread in their future projections provide the necessary information to objectively select GCMs while keeping all aspects of necessity in view.

54 ENVIRONMENTAL SCIENCES↗

Experiment K-310: The effect of space flight on ostenogenesis and dentinogenesis in the mandible of rats. Supplement 1: The effects of space flight on alveolar bone modeling and remodeling in the rat mandible

The histomorphometric study of alveolar bone, a non-weight-bearing bone submitted mainly to the mechanical stimulations of mastication, showed that space flight decreases the remodeling activity but does not induce a negative balance between resorption and formation. The most dramatic effect of space flight has been observed along the periosteal surface, and especially in areas not covered with masticatory muscles, where bone formation almost stopped completely during the flight period. This bone, having been submitted to the same mechanical forces in the flight animals and the controls, leads to the conclusion that factors other than mechanical loading might be involved in the decreased bone formation during flight.

Van, P. T.↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction, aircraft performance model,↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction↗

Linking strain anisotropy and plasticity in copper metallization

The elastic anisotropy of copper leads to significant variation in the x-ray elastic constants (XEC), which link diffraction-based strain measurements to stress. An accurate depiction of the mechanical response in copper thin films requires a determination of an appropriate grain interaction model that lies between Voigt and Reuss limits. It is shown that the associated XEC weighting fraction, x*, between these limits provides a metric by which strain anisotropy can be quantified. Experimental values of x*, as determined by a linear regression scheme of diffraction data collected from multiple reflections, reveal the degree of strain anisotropy and its dependence on plastic deformation induced during in-situ and ex-situ thermal treatments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of CuO{sub 2} Lattice Strain on the Electronic Structure and Properties of High-T{sub c} Cuprate Family

Doping and strain dependences of the electronic structure of the CuO{sub 6}-octahedra layer within LDA+GTB method in the frameworks of six-band p-d model are calculated. Band structure and Fermi surface of the quasiparticle excitations in the effective Hubbard model are characterized by inhomogeneous distribution of the k-dependent spectral weight. Doping results in reconstruction of the band structure, redistribution of the spectral weight over dispersion surface, and reconstruction of Fermi surface from small hole pockets in the underdoped region to large hole contour in the overdoped region through two quantum phase transitions (QPT). Strain increasing leads to displacement of the valence and conductivity bands, bandwidths decreasing, and shift of the concentrations corresponding to QPTs. Strain dependences of the antiferromagnetic exchange coupling and DOS at the Fermi level determining superconducting temperature T{sub c} are obtained. Effective exchange coupling in the equation for T{sub c} demonstrates monotonic strain dependence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Investigating the age distribution of fracture discharge using multiple environmental tracers, Bedrichov Tunnel, Czech Republic

We investigated the transit time distribution (TTD) of discharge collected from fractures in the Bedrichov Tunnel, Czech Republic, using lumped parameter models and multiple environmental tracers. We then utilize time series of δ 18 O, δ 2 H and 3 H along with CFC measurements from individual fractures in the Bedrichov Tunnel of the Czech Republic to investigate the TTD, and the uncertainty in estimated mean travel time in several fracture networks of varying length and discharge. We also compare several TTDs, including the dispersion distribution, the exponential distribution, and a developed TTD which includes the effects of matrix diffusion. The effect of seasonal recharge is explored by comparing several seasonal weighting functions to derive the historical recharge concentration. We identify best fit mean ages for each TTD by minimizing the error-weighted, multi-tracer χ2 residual for each seasonal weighting function. We use this methodology to test the ability of each TTD and seasonal input function to fit the observed tracer concentrations, and the effect of choosing different TTD and seasonal recharge functions on the mean age estimation. We find that the estimated mean transit time is a function of both the assumed TTD and seasonal weighting function. Best fits as measured by the χ2 value were achieved for the dispersion model using the seasonal input function developed here for two of the three modeled sites, while at the third site, equally good fits were achieved with the exponential model and the dispersion model and our seasonal input function. The average mean transit time for all TTDs and seasonal input functions converged to similar values at each location. The sensitivity of the estimated mean transit time to the seasonal weighting function was equal to that of the TTD. These results indicated that understanding seasonality of recharge is at least as important as the uncertainty in the flow path distribution in fracture networks and that unique identification of the TTD and mean transit time is difficult given the uncertainty in the recharge function. But, the mean transit time appears to be relatively robust to the structural model uncertainty. The results presented here should be applicable to other studies using environmental tracers to constrain flow and transport properties in fractured rock systems.

54 ENVIRONMENTAL SCIENCES↗

Bilevel Parameter Optimization for Nonlocal Image Denoising Models

We propose a bilevel optimization approach for the determination of parameters in nonlocal image denoising. We consider both spatial weights in front of the fidelity term, as well as weights within the kernel of the nonlocal operator. In both cases we investigate the differentiability of the solution operator in function spaces and derive a first order optimality system that characterizes local minima. For the numerical solution of the problems, we propose a second-order optimization algorithm in combination with a finite element discretization of the nonlocal denoising models and a computational strategy for the solution of the resulting dense linear systems. Several experiments are run in order to show the suitability of our approach.

97 MATHEMATICS AND COMPUTING↗

Plan, formulate, and discuss a NASTRAN finite element model of the UH-60A helicopter airframe

Under a rotorcraft structural dynamics program sponsored by the NASA Langley Research Center, Sikorsky Aircraft, together with the other major helicopter airframe manufacturers, is engaged in a study to improve the use of finite element analysis to predict the dynamic behavior of helicopter airframes. This program, which was designated DAMVIBS (Design Analysis Methods for VIBrationS), includes activities in the areas of: planning, creating, and documenting finite element models of helicopter airframes; the performance of ground vibration tests; and the correlation of test and analysis. The work performed at Sikorsky Aircraft for planning, creating, and documenting a finite element model of the UH-60A BLACK HAWK helicopter airframe is summarized. A complete description of the components of the helicopter which are to be represented in the model is presented and includes: the structural arrangement, the identification of primary and secondary structure, the components of the drive and power trains, and the attachment of large weight items to the structure. Also presented are the techniques which were used to formulate the structural finite element model for static analysis, for forming the mass and vibration models for dynamic analysis, and the procedures which were used to check out and verify the integrity of the model. Initial predictions for the vibration modes for the helicopter are included.

Dinyovszky, P.↗