Plasmon excitation in metal slab by fast point charge: The role of additional boundary conditions in quantum hydrodynamic model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Ying-Ying; Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario N2L 3G1; An, Sheng-Bai

2014-10-15

We study the wake effect in the induced potential and the stopping power due to plasmon excitation in a metal slab by a point charge moving inside the slab. Nonlocal effects in the response of the electron gas in the metal are described by a quantum hydrodynamic model, where the equation of electronic motion contains both a quantum pressure term and a gradient correction from the Bohm quantum potential, resulting in a fourth-order differential equation for the perturbed electron density. Thus, besides using the condition that the normal component of the electron velocity should vanish at the impenetrable boundary ofmore » the metal, a consistent inclusion of the gradient correction is shown to introduce two possibilities for an additional boundary condition for the perturbed electron density. We show that using two different sets of boundary conditions only gives rise to differences in the wake potential at large distances behind the charged particle. On the other hand, the gradient correction in the quantum hydrodynamic model is seen to cause a reduction in the depth of the potential well closest to the particle, and a reduction of its stopping power. Even for a particle moving in the center of the slab, we observe nonlocal effects in the induced potential and the stopping power due to reduction of the slab thickness, which arise from the gradient correction in the quantum hydrodynamic model.« less

Experimental Determination of Jet Boundary Corrections for Airfoil Tests in Four Open Wind Tunnel Jets of Different Shapes

NASA Technical Reports Server (NTRS)

Knight, Montgomery; Harris, Thomas A

1931-01-01

This experimental investigation was conducted primarily for the purpose of obtaining a method of correcting to free air conditions the results of airfoil force tests in four open wind tunnel jets of different shapes. Tests were also made to determine whether the jet boundaries had any appreciable effect on the pitching moments of a complete airplane model. Satisfactory corrections for the effect of the boundaries of the various jets were obtained for all the airfoils tested, the span of the largest being 0.75 of the jet width. The corrections for angle of attack were, in general, larger than those for drag. The boundaries had no appreciable effect on the pitching moments of either the airfoils or the complete airplane model. Increasing turbulence appeared to increase the minimum drag and maximum lift and to decrease the pitching moment.

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

NASA Technical Reports Server (NTRS)

Parsani, Matteo; Carpenter, Mark H.; Nielsen, Eric J.

2015-01-01

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Statistical Calibration and Validation of a Homogeneous Ventilated Wall-Interference Correction Method for the National Transonic Facility

NASA Technical Reports Server (NTRS)

Walker, Eric L.

2005-01-01

Wind tunnel experiments will continue to be a primary source of validation data for many types of mathematical and computational models in the aerospace industry. The increased emphasis on accuracy of data acquired from these facilities requires understanding of the uncertainty of not only the measurement data but also any correction applied to the data. One of the largest and most critical corrections made to these data is due to wall interference. In an effort to understand the accuracy and suitability of these corrections, a statistical validation process for wall interference correction methods has been developed. This process is based on the use of independent cases which, after correction, are expected to produce the same result. Comparison of these independent cases with respect to the uncertainty in the correction process establishes a domain of applicability based on the capability of the method to provide reasonable corrections with respect to customer accuracy requirements. The statistical validation method was applied to the version of the Transonic Wall Interference Correction System (TWICS) recently implemented in the National Transonic Facility at NASA Langley Research Center. The TWICS code generates corrections for solid and slotted wall interference in the model pitch plane based on boundary pressure measurements. Before validation could be performed on this method, it was necessary to calibrate the ventilated wall boundary condition parameters. Discrimination comparisons are used to determine the most representative of three linear boundary condition models which have historically been used to represent longitudinally slotted test section walls. Of the three linear boundary condition models implemented for ventilated walls, the general slotted wall model was the most representative of the data. The TWICS code using the calibrated general slotted wall model was found to be valid to within the process uncertainty for test section Mach numbers less than or equal to 0.60. The scatter among the mean corrected results of the bodies of revolution validation cases was within one count of drag on a typical transport aircraft configuration for Mach numbers at or below 0.80 and two counts of drag for Mach numbers at or below 0.90.

A stable penalty method for the compressible Navier-Stokes equations. 1: Open boundary conditions

NASA Technical Reports Server (NTRS)

Hesthaven, J. S.; Gottlieb, D.

1994-01-01

The purpose of this paper is to present asymptotically stable open boundary conditions for the numerical approximation of the compressible Navier-Stokes equations in three spatial dimensions. The treatment uses the conservation form of the Navier-Stokes equations and utilizes linearization and localization at the boundaries based on these variables. The proposed boundary conditions are applied through a penalty procedure, thus ensuring correct behavior of the scheme as the Reynolds number tends to infinity. The versatility of this method is demonstrated for the problem of a compressible flow past a circular cylinder.

Identifying the groundwater basin boundaries, using environmental isotopes: a case study

NASA Astrophysics Data System (ADS)

Demiroğlu, Muhterem

2017-06-01

Groundwater, which is renewable under current climatic conditions separately from other natural sources, in fact is a finite resource in terms of quality and fossil groundwater. Researchers have long emphasized the necessity of exploiting, operating, conserving and managing groundwater in an efficient and sustainable manner with an integrated water management approach. The management of groundwater needs reliable information about changes on groundwater quantity and quality. Environmental isotopes are the most important tools to provide this support. No matter which method we use to calculate the groundwater budget and flow equations, we need to determine boundary conditions or the physical boundaries of the domain. The Groundwater divide line or basin boundaries that separate the two adjacent basin recharge areas from each other must be drawn correctly to be successful in defining complex groundwater basin boundary conditions. Environmental isotope data, as well as other methods provide support for determining recharge areas of the aquifers, especially for karst aquifers, residence time and interconnections between aquifer systems. This study demonstrates the use of environmental isotope data to interpret and correct groundwater basin boundaries giving as an example the Yeniçıkrı basin within the main Sakarya basin.

Theory of Irregular Waveguides with Slowly Changing Parameters

DTIC Science & Technology

1979-04-05

different waves are orthogcnal between themselves. The conditions of orthogonality we will record in such a way that it would be correct during any...replace the strain of surface, on which it is correct (3.2), by a small change in this condition on the undisturbed surface. Let us establish/install this...it is possible to show which (6.’ is correct with any sign L. The replacement of strain by boundary condition (6.1) introduces into calculation soca

Assimilation of Wave and Current Data for Prediction of Inlet and River Mouth Dynamics

DTIC Science & Technology

2013-07-01

onto the Delft3D computational grid and the specification of Riemann -type boundary conditions for the boundary-normal velocity and surface elevation...conditions from time- history data from in situ tide gages. The corrections are applied to the surface-elevation contribution to the Riemann boundary...The algorithms described above are all of the strong-constraint variational variety, and make use of adjoint solvers corresponding to the various

Some new results for the one-loop mass correction to the compactified λϕ4 theory

NASA Astrophysics Data System (ADS)

Fucci, Guglielmo; Kirsten, Klaus

2018-03-01

In this work, we consider the one-loop effective action of a self-interacting λϕ4 field propagating in a D dimensional Euclidean space endowed with d ≤ D compact dimensions. The main purpose of this paper is to compute the corrections to the mass of the field due to the presence of the compactified dimensions. Although the results of the one-loop correction to the mass of a λϕ4 field are very well known for compactified toroidal spaces, where the field obeys periodic boundary conditions, similar results do not appear to be readily available for cases in which the scalar field is subject to Dirichlet and Neumann boundary conditions. We apply the results of the one-loop mass correction to the study of the critical temperature in Ginzburg-Landau models.

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Silva, Goncalo, E-mail: goncalo.nuno.silva@gmail.com; Talon, Laurent, E-mail: talon@fast.u-psud.fr; Ginzburg, Irina, E-mail: irina.ginzburg@irstea.fr

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies themore » effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.« less

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

NASA Astrophysics Data System (ADS)

Silva, Goncalo; Talon, Laurent; Ginzburg, Irina

2017-04-01

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes-Brinkman-Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes-Brinkman-Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

CDUCT-LaRC Status - Shear Layer Refraction and Noise Radiation

NASA Technical Reports Server (NTRS)

Nark, Douglas M.; Farassat, F.

2006-01-01

A proposed boundary condition accounting for shear layer effects within the Ffowcs Williams-Hawkings radiation module of the CDUCT-LaRC code is investigated. The development and numerical justification of the boundary condition formulation are reviewed. An initial assessment of the effectiveness of the shear layer correction is conducted through comparison with experimental data. Preliminary results indicate that the correction provides physically meaningful modifications of the baseline predicted directivity patterns. Trends of peak directivity steepening and shifting that appeared in predicted patterns were found to follow similar structures in measured data, particularly at higher radiation angles.

Compressibility Considerations for kappa-omega Turbulence Models in Hypersonic Boundary Layer Applications

NASA Technical Reports Server (NTRS)

Rumsey, C. L.

2009-01-01

The ability of kappa-omega models to predict compressible turbulent skin friction in hypersonic boundary layers is investigated. Although uncorrected two-equation models can agree well with correlations for hot-wall cases, they tend to perform progressively worse - particularly for cold walls - as the Mach number is increased in the hypersonic regime. Simple algebraic models such as Baldwin-Lomax perform better compared to experiments and correlations in these circumstances. Many of the compressibility corrections described in the literature are summarized here. These include corrections that have only a small influence for kappa-omega models, or that apply only in specific circumstances. The most widely-used general corrections were designed for use with jet or mixing-layer free shear flows. A less well-known dilatation-dissipation correction intended for boundary layer flows is also tested, and is shown to agree reasonably well with the Baldwin-Lomax model at cold-wall conditions. It exhibits a less dramatic influence than the free shear type of correction. There is clearly a need for improved understanding and better overall physical modeling for turbulence models applied to hypersonic boundary layer flows.

Assessment of Cracks in Stress Concentration Regions with Localized Plastic Zones

DOE Office of Scientific and Technical Information (OSTI.GOV)

Friedman, E.

1998-11-25

Marty brittle fracture evaluation procedures include plasticity corrections to elastically computed stress intensity factors. These corrections, which are based on the existence of a plastic zone in the vicinity of the crack tip, can overestimate the plasticity effect for a crack embedded in a stress concentration region in which the elastically computed stress exceeds the yield strength of the material in a localized zone. The interactions between the crack, which acts to relieve the high stresses driving the crack, plasticity effects in the stress concentration region, and the nature and source of the loading are examined by formulating explicit flawmore » finite element models for a crack emanating from the root of a notch located in a panel subject to an applied tensile stress. The results of these calculations provide conditions under which a crack-tip plasticity correction based on the Irwin plastic zone size overestimates the plasticity effect. A failure assessment diagram (FAD) curve is used to characterize the effect of plasticity on the crack driving force and to define a less restrictive plasticity correction for cracks at notch roots when load-controlled boundary conditions are imposed. The explicit flaw finite element results also demonstrate that stress intensity factors associated with load-controlled boundary conditions, such as those inherent in the ASME Boiler and Pressure Vessel Code as well as in most handbooks of stress intensity factors, can be much higher than those associated with displacement-controlled conditions, such as those that produce residual or thermal stresses. Under certain conditions, the inclusion of plasticity effects for cracks loaded by displacement-controlled boundary conditions reduces the crack driving force thus justifying the elimination of a plasticity correction for such loadings. The results of this study form the basis for removing unnecessary conservatism from flaw evaluation procedures that utilize plasticity corrections.« less

Numerical analysis of hypersonic turbulent film cooling flows

NASA Technical Reports Server (NTRS)

Chen, Y. S.; Chen, C. P.; Wei, H.

1992-01-01

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

Jet-boundary and Plan-form Corrections for Partial-Span Models with Reflection-Plane, End-Plate, or No End-Plate in a Closed Circular Wind Tunnel

NASA Technical Reports Server (NTRS)

Sivells, James C; Deters, Owen J

1946-01-01

A method is presented for determining the jet-boundary and plan-form corrections necessary for application to test data for a partial-span model with a reflection plane, an end plate, or no end plate in a closed circular wind tunnel. Examples are worked out for a partial-span model with each of the three end conditions in the Langley 19-foot pressure tunnel and the corrections are applied to measured values of lift, drag, pitching-moment, rolling-moment, and yawing-moment coefficients.

Charged hadrons in local finite-volume QED+QCD with C⋆ boundary conditions

NASA Astrophysics Data System (ADS)

Lucini, B.; Patella, A.; Ramos, A.; Tantalo, N.

2016-02-01

In order to calculate QED corrections to hadronic physical quantities by means of lattice simulations, a coherent description of electrically-charged states in finite volume is needed. In the usual periodic setup, Gauss's law and large gauge transformations forbid the propagation of electrically-charged states. A possible solution to this problem, which does not violate the axioms of local quantum field theory, has been proposed by Wiese and Polley, and is based on the use of C⋆ boundary conditions. We present a thorough analysis of the properties and symmetries of QED in isolation and QED coupled to QCD, with C⋆ boundary conditions. In particular we learn that a certain class of electrically-charged states can be constructed in a fully consistent fashion without relying on gauge fixing and without peculiar complications. This class includes single particle states of most stable hadrons. We also calculate finite-volume corrections to the mass of stable charged particles and show that these are much smaller than in non-local formulations of QED.

Impact of bias-corrected reanalysis-derived lateral boundary conditions on WRF simulations

NASA Astrophysics Data System (ADS)

Moalafhi, Ditiro Benson; Sharma, Ashish; Evans, Jason Peter; Mehrotra, Rajeshwar; Rocheta, Eytan

2017-08-01

Lateral and lower boundary conditions derived from a suitable global reanalysis data set form the basis for deriving a dynamically consistent finer resolution downscaled product for climate and hydrological assessment studies. A problem with this, however, is that systematic biases have been noted to be present in the global reanalysis data sets that form these boundaries, biases which can be carried into the downscaled simulations thereby reducing their accuracy or efficacy. In this work, three Weather Research and Forecasting (WRF) model downscaling experiments are undertaken to investigate the impact of bias correcting European Centre for Medium range Weather Forecasting Reanalysis ERA-Interim (ERA-I) atmospheric temperature and relative humidity using Atmospheric Infrared Sounder (AIRS) satellite data. The downscaling is performed over a domain centered over southern Africa between the years 2003 and 2012. The sample mean and the mean as well as standard deviation at each grid cell for each variable are used for bias correction. The resultant WRF simulations of near-surface temperature and precipitation are evaluated seasonally and annually against global gridded observational data sets and compared with ERA-I reanalysis driving field. The study reveals inconsistencies between the impact of the bias correction prior to downscaling and the resultant model simulations after downscaling. Mean and standard deviation bias-corrected WRF simulations are, however, found to be marginally better than mean only bias-corrected WRF simulations and raw ERA-I reanalysis-driven WRF simulations. Performances, however, differ when assessing different attributes in the downscaled field. This raises questions about the efficacy of the correction procedures adopted.

Inferring unknow boundary conditions of the Greenland Ice Sheet by assimilating ICESat-1 and IceBridge altimetry intothe Ice Sheet System Model.

NASA Astrophysics Data System (ADS)

Larour, E. Y.; Khazendar, A.; Seroussi, H. L.; Schlegel, N.; Csatho, B. M.; Schenk, A. F.; Rignot, E. J.; Morlighem, M.

2014-12-01

Altimetry signals from missions such as ICESat-1, CryoSat, EnviSat, as well as altimeters onboard Operation IceBridge provide vital insights into processes such as surface mass balance, mass transport and ice-flow dynamics. Historically however, ice-flow models have been focused on assimilating surface velocities from satellite-based radar observations, to infer properties such as basal friction or the position of the bedrock. Here, we leverage a new methodology based on automatic differentation of the Ice Sheet System Model to assimilate surface altimetry data into a reconstruction of the past decade of ice flow on the North Greenland area. We infer corrections to boundary conditions such as basal friction and surface mass balance, as well as corrections to the ice hardness, to best-match the observed altimetry record. We compare these corrections between glaciers such as Petermann Glacier, 79 North and Zacchariae Isstrom. The altimetry signals exhibit very different patterns between East and West, which translate into very different signatures for the inverted boundary conditions. This study gives us greater insights into what differentiates different basins, both in terms of mass transport and ice-flow dynamics, and what could bethe controlling mechanisms behind the very different evolutions of these basins.

Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices

NASA Technical Reports Server (NTRS)

Nance, Robert P.; Hash, David B.; Hassan, H. A.

1997-01-01

A study is made of the issues surrounding prediction of microchannel flows using the direct simulation Monte Carlo method. This investigation includes the introduction and use of new inflow and outflow boundary conditions suitable for subsonic flows. A series of test simulations for a moderate-size microchannel indicates that a high degree of grid under-resolution in the streamwise direction may be tolerated without loss of accuracy. In addition, the results demonstrate the importance of physically correct boundary conditions, as well as possibilities for reducing the time associated with the transient phase of a simulation. These results imply that simulations of longer ducts may be more feasible than previously envisioned.

Efficient Fluid Dynamic Design Optimization Using Cartesian Grids

NASA Technical Reports Server (NTRS)

Dadone, A.; Grossman, B.; Sellers, Bill (Technical Monitor)

2004-01-01

This report is subdivided in three parts. The first one reviews a new approach to the computation of inviscid flows using Cartesian grid methods. The crux of the method is the curvature-corrected symmetry technique (CCST) developed by the present authors for body-fitted grids. The method introduces ghost cells near the boundaries whose values are developed from an assumed flow-field model in vicinity of the wall consisting of a vortex flow, which satisfies the normal momentum equation and the non-penetration condition. The CCST boundary condition was shown to be substantially more accurate than traditional boundary condition approaches. This improved boundary condition is adapted to a Cartesian mesh formulation, which we call the Ghost Body-Cell Method (GBCM). In this approach, all cell centers exterior to the body are computed with fluxes at the four surrounding cell edges. There is no need for special treatment corresponding to cut cells which complicate other Cartesian mesh methods.

Computation of the shock-wave boundary layer interaction with flow separation

NASA Technical Reports Server (NTRS)

Ardonceau, P.; Alziary, T.; Aymer, D.

1980-01-01

The boundary layer concept is used to describe the flow near the wall. The external flow is approximated by a pressure displacement relationship (tangent wedge in linearized supersonic flow). The boundary layer equations are solved in finite difference form and the question of the presence and unicity of the solution is considered for the direct problem (assumed pressure) or converse problem (assumed displacement thickness, friction ratio). The coupling algorithm presented implicitly processes the downstream boundary condition necessary to correctly define the interacting boundary layer problem. The algorithm uses a Newton linearization technique to provide a fast convergence.

Automated Boundary Conditions for Wind Tunnel Simulations

NASA Technical Reports Server (NTRS)

Carlson, Jan-Renee

2018-01-01

Computational fluid dynamic (CFD) simulations of models tested in wind tunnels require a high level of fidelity and accuracy particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure the proper characterization of both the physical geometry of the wind tunnel and recreating the correct flow conditions inside the wind tunnel. The typical trial-and-error effort used for determining the boundary condition values for a particular tunnel configuration are time and computer resource intensive. This paper describes a method for calculating and updating the back pressure boundary condition in wind tunnel simulations by using a proportional-integral-derivative controller. The controller methodology and equations are discussed, and simulations using the controller to set a tunnel Mach number in the NASA Langley 14- by 22-Foot Subsonic Tunnel are demonstrated.

An empirical method to correct for temperature-dependent variations in the overlap function of CHM15k ceilometers

NASA Astrophysics Data System (ADS)

Hervo, Maxime; Poltera, Yann; Haefele, Alexander

2016-07-01

Imperfections in a lidar's overlap function lead to artefacts in the background, range and overlap-corrected lidar signals. These artefacts can erroneously be interpreted as an aerosol gradient or, in extreme cases, as a cloud base leading to false cloud detection. A correct specification of the overlap function is hence crucial in the use of automatic elastic lidars (ceilometers) for the detection of the planetary boundary layer or of low cloud. In this study, an algorithm is presented to correct such artefacts. It is based on the assumption of a homogeneous boundary layer and a correct specification of the overlap function down to a minimum range, which must be situated within the boundary layer. The strength of the algorithm lies in a sophisticated quality-check scheme which allows the reliable identification of favourable atmospheric conditions. The algorithm was applied to 2 years of data from a CHM15k ceilometer from the company Lufft. Backscatter signals corrected for background, range and overlap were compared using the overlap function provided by the manufacturer and the one corrected with the presented algorithm. Differences between corrected and uncorrected signals reached up to 45 % in the first 300 m above ground. The amplitude of the correction turned out to be temperature dependent and was larger for higher temperatures. A linear model of the correction as a function of the instrument's internal temperature was derived from the experimental data. Case studies and a statistical analysis of the strongest gradient derived from corrected signals reveal that the temperature model is capable of a high-quality correction of overlap artefacts, in particular those due to diurnal variations. The presented correction method has the potential to significantly improve the detection of the boundary layer with gradient-based methods because it removes false candidates and hence simplifies the attribution of the detected gradients to the planetary boundary layer. A particularly significant benefit can be expected for the detection of shallow stable layers typical of night-time situations. The algorithm is completely automatic and does not require any on-site intervention but requires the definition of an adequate instrument-specific configuration. It is therefore suited for use in large ceilometer networks.

Boundary condition at a two-phase interface in the lattice Boltzmann method for the convection-diffusion equation.

PubMed

Yoshida, Hiroaki; Kobayashi, Takayuki; Hayashi, Hidemitsu; Kinjo, Tomoyuki; Washizu, Hitoshi; Fukuzawa, Kenji

2014-07-01

A boundary scheme in the lattice Boltzmann method (LBM) for the convection-diffusion equation, which correctly realizes the internal boundary condition at the interface between two phases with different transport properties, is presented. The difficulty in satisfying the continuity of flux at the interface in a transient analysis, which is inherent in the conventional LBM, is overcome by modifying the collision operator and the streaming process of the LBM. An asymptotic analysis of the scheme is carried out in order to clarify the role played by the adjustable parameters involved in the scheme. As a result, the internal boundary condition is shown to be satisfied with second-order accuracy with respect to the lattice interval, if we assign appropriate values to the adjustable parameters. In addition, two specific problems are numerically analyzed, and comparison with the analytical solutions of the problems numerically validates the proposed scheme.

On the boundary treatment in spectral methods for hyperbolic systems

NASA Technical Reports Server (NTRS)

Canuto, C.; Quarteroni, A.

1986-01-01

Spectral methods were successfully applied to the simulation of slow transients in gas transportation networks. Implicit time advancing techniques are naturally suggested by the nature of the problem. The correct treatment of the boundary conditions are clarified in order to avoid any stability restriction originated by the boundaries. The Beam and Warming and the Lerat schemes are unconditionally linearly stable when used with a Chebyshev pseudospectral method. Engineering accuracy for a gas transportation problem is achieved at Courant numbers up to 100.

On the boundary treatment in spectral methods for hyperbolic systems

NASA Technical Reports Server (NTRS)

Canuto, Claudio; Quarteroni, Alfio

1987-01-01

Spectral methods were successfully applied to the simulation of slow transients in gas transportation networks. Implicit time advancing techniques are naturally suggested by the nature of the problem. The correct treatment of the boundary conditions is clarified in order to avoid any stability restriction originated by the boundaries. The Beam and Warming and the Lerat schemes are unconditionally linearly stable when used with a Chebyshev pseudospectral method. Engineering accuracy for a gas transportation problem is achieved at Courant numbers up to 100.

Comment on: "Corrections to the Mathematical Formulation of a Backwards Lagrangian Particle Dispersion Model" by Gibson and Sailor (2012: Boundary-Layer Meteorology 145, 399-406)

NASA Astrophysics Data System (ADS)

Stöckl, Stefan; Rotach, Mathias W.; Kljun, Natascha

2018-01-01

We discuss the results of Gibson and Sailor (Boundary-Layer Meteorol 145:399-406, 2012) who suggest several corrections to the mathematical formulation of the Lagrangian particle dispersion model of Rotach et al. (Q J R Meteorol Soc 122:367-389, 1996). While most of the suggested corrections had already been implemented in the 1990s, one suggested correction raises a valid point, but results in a violation of the well-mixed criterion. Here we improve their idea and test the impact on model results using a well-mixed test and a comparison with wind-tunnel experimental data. The new approach results in similar dispersion patterns as the original approach, while the approach suggested by Gibson and Sailor leads to erroneously reduced concentrations near the ground in convective and especially forced convective conditions.

Calculation of sidewall boundary-layer parameters from rake measurements for the Langley 0.3-meter transonic cryogenic tunnel

NASA Technical Reports Server (NTRS)

Murthy, A. V.

1987-01-01

Correction of airfoil data for sidewall boundary-layer effects requires a knowledge of the boundary-layer displacement thickness and the shape factor with the tunnel empty. To facilitate calculation of these quantities under various test conditions for the Langley 0.3 m Transonic Cryogenic Tunnel, a computer program was written. This program reads the various tunnel parameters and the boundary-layer rake total head pressure measurements directly from the Engineering Unit tapes to calculate the required sidewall boundary-layer parameters. Details of the method along with the results for a sample case are presented.

Skin Friction and Transition Location Measurement on Supersonic Transport Models

NASA Technical Reports Server (NTRS)

Kennelly, Robert A., Jr.; Goodsell, Aga M.; Olsen, Lawrence E. (Technical Monitor)

2000-01-01

Flow visualization techniques were used to obtain both qualitative and quantitative skin friction and transition location data in wind tunnel tests performed on two supersonic transport models at Mach 2.40. Oil-film interferometry was useful for verifying boundary layer transition, but careful monitoring of model surface temperatures and systematic examination of the effects of tunnel start-up and shutdown transients will be required to achieve high levels of accuracy for skin friction measurements. A more common technique, use of a subliming solid to reveal transition location, was employed to correct drag measurements to a standard condition of all-turbulent flow on the wing. These corrected data were then analyzed to determine the additional correction required to account for the effect of the boundary layer trip devices.

Transport synthetic acceleration with opposing reflecting boundary conditions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zika, M.R.; Adams, M.L.

2000-02-01

The transport synthetic acceleration (TSA) scheme is extended to problems with opposing reflecting boundary conditions. This synthetic method employs a simplified transport operator as its low-order approximation. A procedure is developed that allows the use of the conjugate gradient (CG) method to solve the resulting low-order system of equations. Several well-known transport iteration algorithms are cast in a linear algebraic form to show their equivalence to standard iterative techniques. Source iteration in the presence of opposing reflecting boundary conditions is shown to be equivalent to a (poorly) preconditioned stationary Richardson iteration, with the preconditioner defined by the method of iteratingmore » on the incident fluxes on the reflecting boundaries. The TSA method (and any synthetic method) amounts to a further preconditioning of the Richardson iteration. The presence of opposing reflecting boundary conditions requires special consideration when developing a procedure to realize the CG method for the proposed system of equations. The CG iteration may be applied only to symmetric positive definite matrices; this condition requires the algebraic elimination of the boundary angular corrections from the low-order equations. As a consequence of this elimination, evaluating the action of the resulting matrix on an arbitrary vector involves two transport sweeps and a transmission iteration. Results of applying the acceleration scheme to a simple test problem are presented.« less

Accurate artificial boundary conditions for the semi-discretized linear Schrödinger and heat equations on rectangular domains

NASA Astrophysics Data System (ADS)

Ji, Songsong; Yang, Yibo; Pang, Gang; Antoine, Xavier

2018-01-01

The aim of this paper is to design some accurate artificial boundary conditions for the semi-discretized linear Schrödinger and heat equations in rectangular domains. The Laplace transform in time and discrete Fourier transform in space are applied to get Green's functions of the semi-discretized equations in unbounded domains with single-source. An algorithm is given to compute these Green's functions accurately through some recurrence relations. Furthermore, the finite-difference method is used to discretize the reduced problem with accurate boundary conditions. Numerical simulations are presented to illustrate the accuracy of our method in the case of the linear Schrödinger and heat equations. It is shown that the reflection at the corners is correctly eliminated.

Exact solution for the layered convection of a viscous incompressible fluid at specified temperature gradients and tangential forces on the free boundary

NASA Astrophysics Data System (ADS)

Burmasheva, N. V.; Prosviryakov, E. Yu.

2017-12-01

A new exact analytical solution of a system of thermal convection equations in the Boussinesq approximation describing layered flows in an incompressible viscous fluid is obtained. A fluid flow in an infinite layer is considered. Convection in the fluid is induced by tangential stresses specified on the upper non-deformable boundary. At the fixed lower boundary, the no-slip condition is satisfied. Temperature corrections are given on the both boundaries of the fluid layer. The possibility of physical field stratification is investigated.

Temporal Gain Correction for X-Ray Calorimeter Spectrometers

NASA Technical Reports Server (NTRS)

Porter, F. S.; Chiao, M. P.; Eckart, M. E.; Fujimoto, R.; Ishisaki, Y.; Kelley, R. L.; Kilbourne, C. A.; Leutenegger, M. A.; McCammon, D.; Mitsuda, K.

2016-01-01

Calorimetric X-ray detectors are very sensitive to their environment. The boundary conditions can have a profound effect on the gain including heat sink temperature, the local radiation temperature, bias, and the temperature of the readout electronics. Any variation in the boundary conditions can cause temporal variations in the gain of the detector and compromise both the energy scale and the resolving power of the spectrometer. Most production X-ray calorimeter spectrometers, both on the ground and in space, have some means of tracking the gain as a function of time, often using a calibration spectral line. For small gain changes, a linear stretch correction is often sufficient. However, the detectors are intrinsically non-linear and often the event analysis, i.e., shaping, optimal filters etc., add additional non-linearity. Thus for large gain variations or when the best possible precision is required, a linear stretch correction is not sufficient. Here, we discuss a new correction technique based on non-linear interpolation of the energy-scale functions. Using Astro-HSXS calibration data, we demonstrate that the correction can recover the X-ray energy to better than 1 part in 104 over the entire spectral band to above 12 keV even for large-scale gain variations. This method will be used to correct any temporal drift of the on-orbit per-pixel gain using on-board calibration sources for the SXS instrument on the Astro-H observatory.

Performance evaluation and bias correction of DBS measurements for a 1290-MHz boundary layer profiler.

PubMed

Liu, Zhao; Zheng, Chaorong; Wu, Yue

2018-02-01

Recently, the government installed a boundary layer profiler (BLP), which is operated under the Doppler beam swinging mode, in a coastal area of China, to acquire useful wind field information in the atmospheric boundary layer for several purposes. And under strong wind conditions, the performance of the BLP is evaluated. It is found that, even though the quality controlled BLP data show good agreement with the balloon observations, a systematic bias can always be found for the BLP data. For the low wind velocities, the BLP data tend to overestimate the atmospheric wind. However, with the increment of wind velocity, the BLP data show a tendency of underestimation. In order to remove the effect of poor quality data on bias correction, the probability distribution function of the differences between the two instruments is discussed, and it is found that the t location scale distribution is the most suitable probability model when compared to other probability models. After the outliers with a large discrepancy, which are outside of 95% confidence interval of the t location scale distribution, are discarded, the systematic bias can be successfully corrected using a first-order polynomial correction function. The methodology of bias correction used in the study not only can be referred for the correction of other wind profiling radars, but also can lay a solid basis for further analysis of the wind profiles.

Performance evaluation and bias correction of DBS measurements for a 1290-MHz boundary layer profiler

NASA Astrophysics Data System (ADS)

Liu, Zhao; Zheng, Chaorong; Wu, Yue

2018-02-01

Recently, the government installed a boundary layer profiler (BLP), which is operated under the Doppler beam swinging mode, in a coastal area of China, to acquire useful wind field information in the atmospheric boundary layer for several purposes. And under strong wind conditions, the performance of the BLP is evaluated. It is found that, even though the quality controlled BLP data show good agreement with the balloon observations, a systematic bias can always be found for the BLP data. For the low wind velocities, the BLP data tend to overestimate the atmospheric wind. However, with the increment of wind velocity, the BLP data show a tendency of underestimation. In order to remove the effect of poor quality data on bias correction, the probability distribution function of the differences between the two instruments is discussed, and it is found that the t location scale distribution is the most suitable probability model when compared to other probability models. After the outliers with a large discrepancy, which are outside of 95% confidence interval of the t location scale distribution, are discarded, the systematic bias can be successfully corrected using a first-order polynomial correction function. The methodology of bias correction used in the study not only can be referred for the correction of other wind profiling radars, but also can lay a solid basis for further analysis of the wind profiles.

Directional bottom roughness associated with waves, currents, and ripples

USGS Publications Warehouse

Sherwood, Christopher R.; Rosati, Julie D.; Wang, Ping; Roberts, Tiffany M.

2011-01-01

Roughness lengths are used in wave-current bottom boundary layer models to parameterize drag associated with grain roughness, the effect of saltating grains during sediment transport, and small-scale bottom topography (ripples and biogenic features). We made field measurements of flow parameters and recorded sonar images of ripples at the boundary of a sorted-bedform at ~12-m depth on the inner shelf for a range of wave and current conditions over two months. We compared estimates of apparent bottom roughness inferred from the flow measurements with bottom roughness calculated using ripple geometry and the Madsen (1994) one-dimensional (vertical) wave-current bottom boundary layer model. One result of these comparisons was that the model over predicted roughness of flow from the dormant large ripples when waves were small. We developed a correction to the ripple-roughness model that incorporates an apparent ripple wavelength related to the combined wave-current flow direction. This correction provides a slight improvement for low-wave conditions, but does not address several other differences between observations and the modeled roughness.

Tables for correcting airfoil data obtained in the Langley 0.3-meter transonic cryogenic tunnel for sidewall boundary-layer effects

NASA Technical Reports Server (NTRS)

Jenkins, R. V.; Adcock, J. B.

1986-01-01

Tables for correcting airfoil data taken in the Langley 0.3-meter Transonic Cryogenic Tunnel for the presence of sidewall boundary layer are presented. The corrected Mach number and the correction factor are minutely altered by a 20 percent change in the boundary layer virtual origin distance. The sidewall boundary layer displacement thicknesses measured for perforated sidewall inserts and without boundary layer removal agree with the values calculated for solid sidewalls.

Large Eddy Simulation of a Film Cooling Technique with a Plenum

NASA Astrophysics Data System (ADS)

Dharmarathne, Suranga; Sridhar, Narendran; Araya, Guillermo; Castillo, Luciano; Parameswaran, Sivapathasund

2012-11-01

Factors that affect the film cooling performance have been categorized into three main groups: (i) coolant & mainstream conditions, (ii) hole geometry & configuration, and (iii) airfoil geometry Bogard et al. (2006). The present study focuses on the second group of factors, namely, the modeling of coolant hole and the plenum. It is required to simulate correct physics of the problem to achieve more realistic numerical results. In this regard, modeling of cooling jet hole and the plenum chamber is highly important Iourokina et al. (2006). Substitution of artificial boundary conditions instead of correct plenum design would yield unrealistic results Iourokina et al. (2006). This study attempts to model film cooling technique with a plenum using a Large Eddy Simulation.Incompressible coolant jet ejects to the surface of the plate at an angle of 30° where it meets compressible turbulent boundary layer which simulates the turbine inflow conditions. Dynamic multi-scale approach Araya (2011) is introduced to prescribe turbulent inflow conditions. Simulations are carried out for two different blowing ratios and film cooling effectiveness is calculated for both cases. Results obtained from LES will be compared with experimental results.

Acceleration of incremental-pressure-correction incompressible flow computations using a coarse-grid projection method

NASA Astrophysics Data System (ADS)

Kashefi, Ali; Staples, Anne

2016-11-01

Coarse grid projection (CGP) methodology is a novel multigrid method for systems involving decoupled nonlinear evolution equations and linear elliptic equations. The nonlinear equations are solved on a fine grid and the linear equations are solved on a corresponding coarsened grid. Mapping functions transfer data between the two grids. Here we propose a version of CGP for incompressible flow computations using incremental pressure correction methods, called IFEi-CGP (implicit-time-integration, finite-element, incremental coarse grid projection). Incremental pressure correction schemes solve Poisson's equation for an intermediate variable and not the pressure itself. This fact contributes to IFEi-CGP's efficiency in two ways. First, IFEi-CGP preserves the velocity field accuracy even for a high level of pressure field grid coarsening and thus significant speedup is achieved. Second, because incremental schemes reduce the errors that arise from boundaries with artificial homogenous Neumann conditions, CGP generates undamped flows for simulations with velocity Dirichlet boundary conditions. Comparisons of the data accuracy and CPU times for the incremental-CGP versus non-incremental-CGP computations are presented.

The effects of boundary conditions on the steady-state response of three hypothetical ground-water systems; results and implications of numerical experiments

USGS Publications Warehouse

Franke, O. Lehn; Reilly, Thomas E.

1987-01-01

The most critical and difficult aspect of defining a groundwater system or problem for conceptual analysis or numerical simulation is the selection of boundary conditions . This report demonstrates the effects of different boundary conditions on the steady-state response of otherwise similar ground-water systems to a pumping stress. Three series of numerical experiments illustrate the behavior of three hypothetical groundwater systems that are rectangular sand prisms with the same dimensions but with different combinations of constant-head, specified-head, no-flow, and constant-flux boundary conditions. In the first series of numerical experiments, the heads and flows in all three systems are identical, as are the hydraulic conductivity and system geometry . However, when the systems are subjected to an equal stress by a pumping well in the third series, each differs significantly in its response . The highest heads (smallest drawdowns) and flows occur in the systems most constrained by constant- or specified-head boundaries. These and other observations described herein are important in steady-state calibration, which is an integral part of simulating many ground-water systems. Because the effects of boundary conditions on model response often become evident only when the system is stressed, a close match between the potential distribution in the model and that in the unstressed natural system does not guarantee that the model boundary conditions correctly represent those in the natural system . In conclusion, the boundary conditions that are selected for simulation of a ground-water system are fundamentally important to groundwater systems analysis and warrant continual reevaluation and modification as investigation proceeds and new information and understanding are acquired.

Evaluation of different flamelet tabulation methods for laminar spray combustion

NASA Astrophysics Data System (ADS)

Luo, Yujuan; Wen, Xu; Wang, Haiou; Luo, Kun; Fan, Jianren

2018-05-01

In this work, three different flamelet tabulation methods for spray combustion are evaluated. Major differences among these methods lie in the treatment of the temperature boundary conditions of the flamelet equations. Particularly, in the first tabulation method ("M1"), both the fuel and oxidizer temperature boundary conditions are set to be fixed. In the second tabulation method ("M2"), the fuel temperature boundary condition is varied while the oxidizer temperature boundary condition is fixed. In the third tabulation method ("M3"), both the fuel and oxidizer temperature boundary conditions are varied and set to be equal. The focus of this work is to investigate whether the heat transfer between the droplet phase and gas phase can be represented by the studied tabulation methods through a priori analyses. To this end, spray flames stabilized in a three-dimensional counterflow are first simulated with detailed chemistry. Then, the trajectory variables are calculated from the detailed chemistry solutions. Finally, the tabulated thermo-chemical quantities are compared to the corresponding values from the detailed chemistry solutions. The comparisons show that the gas temperature cannot be predicted by "M1" with only a mixture fraction and reaction progress variable being the trajectory variables. The gas temperature can be correctly predicted by both "M2" and "M3," in which the total enthalpy is introduced as an additional manifold. In "M2," variations of the oxidizer temperature are considered with a temperature modification technique, which is not required in "M3." Interestingly, it is found that the mass fractions of the reactants and major products are not sensitive to the representation of the interphase heat transfer in the flamelet chemtables, and they can be correctly predicted by all tabulation methods. By contrast, the intermediate species CO and H2 in the premixed flame reaction zone are over-predicted by all tabulation methods.

Self-correcting quantum memory with a boundary

NASA Astrophysics Data System (ADS)

Hutter, Adrian; Wootton, James R.; Röthlisberger, Beat; Loss, Daniel

2012-11-01

We study the two-dimensional toric-code Hamiltonian with effective long-range interactions between its anyonic excitations induced by coupling the toric code to external fields. It has been shown that such interactions allow an arbitrary increase in the lifetime of the stored quantum information by making L, the linear size of the memory, larger [Chesi , Phys. Rev. APLRAAN1050-294710.1103/PhysRevA.82.022305 82, 022305 (2010)]. We show that for these systems the choice of boundary conditions (open boundaries as opposed to periodic boundary conditions) is not a mere technicality; the influence of anyons produced at the boundaries becomes in fact dominant for large enough L. This influence can be either beneficial or detrimental. In particular, we study an effective Hamiltonian proposed by Pedrocchi [Phys. Rev. BPRBMDO1098-012110.1103/PhysRevB.83.115415 83, 115415 (2011)] that describes repulsion between anyons and anyon holes. For this system, we find a lifetime of the stored quantum information that grows exponentially in L2 for both periodic and open boundary conditions, although the exponent in the latter case is found to be less favorable. However, L is upper bounded through the breakdown of the perturbative treatment of the underlying Hamiltonian.

On a difficulty in eigenfunction expansion solutions for the start-up of fluid flow

NASA Astrophysics Data System (ADS)

Christov, Ivan C.

2015-11-01

Most mathematics and engineering textbooks describe the process of ``subtracting off'' the steady state of a linear parabolic partial differential equation as a technique for obtaining a boundary-value problem with homogeneous boundary conditions that can be solved by separation of variables (i.e., eigenfunction expansions). While this method produces the correct solution for the start-up of the flow of, e.g., a Newtonian fluid between parallel plates, it can lead to erroneous solutions to the corresponding problem for a class of non-Newtonian fluids. We show that the reason for this is the non-rigorous enforcement of the start-up condition in the textbook approach, which leads to a violation of the principle of causality. Nevertheless, these boundary-value problems can be solved correctly using eigenfunction expansions, and we present the formulation that makes this possible (in essence, an application of Duhamel's principle). The solutions obtained by this new approach are shown to agree identically with those obtained by using the Laplace transform in time only, a technique that enforces the proper start-up condition implicitly (hence, the same error cannot be committed). Supported, in part, by NSF Grant DMS-1104047 and the U.S. DOE (Contract No. DE-AC52-06NA25396) through the LANL/LDRD Program.

A new dynamical downscaling approach with GCM bias corrections and spectral nudging

NASA Astrophysics Data System (ADS)

Xu, Zhongfeng; Yang, Zong-Liang

2015-04-01

To improve confidence in regional projections of future climate, a new dynamical downscaling (NDD) approach with both general circulation model (GCM) bias corrections and spectral nudging is developed and assessed over North America. GCM biases are corrected by adjusting GCM climatological means and variances based on reanalysis data before the GCM output is used to drive a regional climate model (RCM). Spectral nudging is also applied to constrain RCM-based biases. Three sets of RCM experiments are integrated over a 31 year period. In the first set of experiments, the model configurations are identical except that the initial and lateral boundary conditions are derived from either the original GCM output, the bias-corrected GCM output, or the reanalysis data. The second set of experiments is the same as the first set except spectral nudging is applied. The third set of experiments includes two sensitivity runs with both GCM bias corrections and nudging where the nudging strength is progressively reduced. All RCM simulations are assessed against North American Regional Reanalysis. The results show that NDD significantly improves the downscaled mean climate and climate variability relative to other GCM-driven RCM downscaling approach in terms of climatological mean air temperature, geopotential height, wind vectors, and surface air temperature variability. In the NDD approach, spectral nudging introduces the effects of GCM bias corrections throughout the RCM domain rather than just limiting them to the initial and lateral boundary conditions, thereby minimizing climate drifts resulting from both the GCM and RCM biases.

Computational Aeroacoustics: An Overview

NASA Technical Reports Server (NTRS)

Tam, Christopher K. W.

2003-01-01

An overview of recent advances in computational aeroacoustics (CAA) is presented. CAA algorithms must not be dispersive and dissipative. It should propagate waves supported by the Euler equations with the correct group velocities. Computation domains are inevitably finite in size. To avoid the reflection of acoustic and other outgoing waves at the boundaries of the computation domain, it is required that special boundary conditions be imposed at the boundary region. These boundary conditions either absorb all the outgoing waves without reflection or allow the waves to exit smoothly. High-order schemes, invariably, supports spurious short waves. These spurious waves tend to pollute the numerical solution. They must be selectively damped or filtered out. All these issues and relevant computation methods are briefly reviewed. Jet screech tones are known to have caused structural fatigue in military combat aircrafts. Numerical simulation of the jet screech phenomenon is presented as an example of a successful application of CAA.

Vanishing Corrections for the Position in a Linear Model of FKPP Fronts

NASA Astrophysics Data System (ADS)

Berestycki, Julien; Brunet, Éric; Harris, Simon C.; Roberts, Matt

2017-02-01

Take the linearised FKPP equation {partialth = partial2xh + h} with boundary condition h( m( t), t) = 0. Depending on the behaviour of the initial condition h 0( x) = h( x, 0) we obtain the asymptotics—up to a o(1) term r( t)—of the absorbing boundary m( t) such that {ω(x) := lim_{tto∞} h(x + m(t) ,t)} exists and is non-trivial. In particular, as in Bramson's results for the non-linear FKPP equation, we recover the celebrated {-3/2 log t} correction for initial conditions decaying faster than {x^{ν}e^{-x}} for some {ν < -2}. Furthermore, when we are in this regime, the main result of the present work is the identification (to first order) of the r( t) term, which ensures the fastest convergence to {ω(x)}. When h 0( x) decays faster than {x^{ν}e^{-x}} for some {ν < -3}, we show that r( t) must be chosen to be {-3√{π/t}}, which is precisely the term predicted heuristically by Ebert-van Saarloos (Phys. D Nonlin. Phenom. 146(1): 1-99, 2000) in the non-linear case (see also Mueller and Munier Phys Rev E 90(4):042143, 2014, Henderson, Commun Math Sci 14(4):973-985, 2016, Brunet and Derrida Stat Phys 1-20, 2015). When the initial condition decays as {x^{ν}e^{-x}} for some {ν in [-3, -2)}, we show that even though we are still in the regime where Bramson's correction is {-3/2 log t}, the Ebert-van Saarloos correction has to be modified. Similar results were recently obtained by Henderson CommunMath Sci 14(4):973-985, 2016 using an analytical approach and only for compactly supported initial conditions.

Anomalous bulk behavior in the free parafermion Z (N ) spin chain

NASA Astrophysics Data System (ADS)

Alcaraz, Francisco C.; Batchelor, Murray T.

2018-06-01

We demonstrate using direct numerical diagonalization and extrapolation methods that boundary conditions have a profound effect on the bulk properties of a simple Z (N ) model for N ≥3 for which the model Hamiltonian is non-Hermitian. For N =2 the model reduces to the well-known quantum Ising model in a transverse field. For open boundary conditions, the Z (N ) model is known to be solved exactly in terms of free parafermions. Once the ends of the open chain are connected by considering the model on a ring, the bulk properties, including the ground-state energy per site, are seen to differ dramatically with increasing N . Other properties, such as the leading finite-size corrections to the ground-state energy, the mass gap exponent, and the specific-heat exponent, are also seen to be dependent on the boundary conditions. We speculate that this anomalous bulk behavior is a topological effect.

An operational large-scale marine planetary boundary layer model

NASA Technical Reports Server (NTRS)

Brown, R. A.; Liu, W. T.

1982-01-01

A marine planetary boundary layer (PBL) model is presented and compared with data from sea-based experiments. The PBL model comprises two layers, the outer an Ekman-Taylor layer with stratification-dependent secondary flow, and the logarithmic surface layer corrected for stratification and humidity effects and variable surface roughness. Corrections are noted for air much warmer than water in stable conditions and for low wind speeds. The layers are analytically defined along with similarity relations and a resistance law for inclusion in a program. An additional interfacial layer correction is developed and shown to be significant for heat flux calculations. Experimental data from GOASEX were used to predict the windfield in the Gulf of Alaska, and JASIN data was used for windfields SE of Iceland. The JASIN-derived wind field predictions were accurate to within 1 m/sec and 10 deg in a 200 km triangle.

Ego Boundary Deficits and the Negative Therapeutic Inter-Action: A Tale of a Whale, a Whale of a Tale.

PubMed

Peterson, Charles A

2017-06-01

Defined variously and unsatisfactorily as a worsening of the patient's condition following a correct interpretation, the negative therapeutic reaction is typically blamed on the patient: "the operation was a success but the patient died." For most neurotic patients unconscious guilt objects to progress and activates the need to suffer. For most character-disturbed patients envy cannot bear the analyst's cleverness. However, patients with ego boundary problems-even sectors of psychosis-may require a different explanatory mechanism, where a correct interpretation may be experienced as a penetration and an engulfment, threatening the intactness of the self. A short-term, time-limited, psychoanalytic psychotherapy that went off the rails following a correct but ill-timed interpretation is presented as an opportunity to amend analytic theory, here favoring the interactional over the intrapsychic. Herman Melville helps tell the tale.

Wave dynamics in an extended macroscopic traffic flow model with periodic boundaries

NASA Astrophysics Data System (ADS)

Wang, Yu-Qing; Chu, Xing-Jian; Zhou, Chao-Fan; Yan, Bo-Wen; Jia, Bin; Fang, Chen-Hao

2018-06-01

Motivated by the previous traffic flow model considering the real-time traffic state, a modified macroscopic traffic flow model is established. The periodic boundary condition is applied to the car-following model. Besides, the traffic state factor R is defined in order to correct the real traffic conditions in a more reasonable way. It is a key step that we introduce the relaxation time as a density-dependent function and provide corresponding evolvement of traffic flow. Three different typical initial densities, namely the high density, the medium one and the low one, are intensively investigated. It can be found that the hysteresis loop exists in the proposed periodic-boundary system. Furthermore, the linear and nonlinear stability analyses are performed in order to test the robustness of the system.

Simulation of charge breeding of rubidium using Monte Carlo charge breeding code and generalized ECRIS model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhao, L.; Cluggish, B.; Kim, J. S.

2010-02-15

A Monte Carlo charge breeding code (MCBC) is being developed by FAR-TECH, Inc. to model the capture and charge breeding of 1+ ion beam in an electron cyclotron resonance ion source (ECRIS) device. The ECRIS plasma is simulated using the generalized ECRIS model which has two choices of boundary settings, free boundary condition and Bohm condition. The charge state distribution of the extracted beam ions is calculated by solving the steady state ion continuity equations where the profiles of the captured ions are used as source terms. MCBC simulations of the charge breeding of Rb+ showed good agreement with recentmore » charge breeding experiments at Argonne National Laboratory (ANL). MCBC correctly predicted the peak of highly charged ion state outputs under free boundary condition and similar charge state distribution width but a lower peak charge state under the Bohm condition. The comparisons between the simulation results and ANL experimental measurements are presented and discussed.« less

FlexibleSUSY-A spectrum generator generator for supersymmetric models

NASA Astrophysics Data System (ADS)

Athron, Peter; Park, Jae-hyeon; Stöckinger, Dominik; Voigt, Alexander

2015-05-01

We introduce FlexibleSUSY, a Mathematica and C++ package, which generates a fast, precise C++ spectrum generator for any SUSY model specified by the user. The generated code is designed with both speed and modularity in mind, making it easy to adapt and extend with new features. The model is specified by supplying the superpotential, gauge structure and particle content in a SARAH model file; specific boundary conditions e.g. at the GUT, weak or intermediate scales are defined in a separate FlexibleSUSY model file. From these model files, FlexibleSUSY generates C++ code for self-energies, tadpole corrections, renormalization group equations (RGEs) and electroweak symmetry breaking (EWSB) conditions and combines them with numerical routines for solving the RGEs and EWSB conditions simultaneously. The resulting spectrum generator is then able to solve for the spectrum of the model, including loop-corrected pole masses, consistent with user specified boundary conditions. The modular structure of the generated code allows for individual components to be replaced with an alternative if available. FlexibleSUSY has been carefully designed to grow as alternative solvers and calculators are added. Predefined models include the MSSM, NMSSM, E6SSM, USSM, R-symmetric models and models with right-handed neutrinos.

An energy-based equilibrium contact angle boundary condition on jagged surfaces for phase-field methods.

PubMed

Frank, Florian; Liu, Chen; Scanziani, Alessio; Alpak, Faruk O; Riviere, Beatrice

2018-08-01

We consider an energy-based boundary condition to impose an equilibrium wetting angle for the Cahn-Hilliard-Navier-Stokes phase-field model on voxel-set-type computational domains. These domains typically stem from μCT (micro computed tomography) imaging of porous rock and approximate a (on μm scale) smooth domain with a certain resolution. Planar surfaces that are perpendicular to the main axes are naturally approximated by a layer of voxels. However, planar surfaces in any other directions and curved surfaces yield a jagged/topologically rough surface approximation by voxels. For the standard Cahn-Hilliard formulation, where the contact angle between the diffuse interface and the domain boundary (fluid-solid interface/wall) is 90°, jagged surfaces have no impact on the contact angle. However, a prescribed contact angle smaller or larger than 90° on jagged voxel surfaces is amplified. As a remedy, we propose the introduction of surface energy correction factors for each fluid-solid voxel face that counterbalance the difference of the voxel-set surface area with the underlying smooth one. The discretization of the model equations is performed with the discontinuous Galerkin method. However, the presented semi-analytical approach of correcting the surface energy is equally applicable to other direct numerical methods such as finite elements, finite volumes, or finite differences, since the correction factors appear in the strong formulation of the model. Copyright © 2018 Elsevier Inc. All rights reserved.

Wall Interference Study of the NTF Slotted Tunnel Using Bodies of Revolution Wall Signature Data

NASA Technical Reports Server (NTRS)

Iyer, Venkit; Kuhl, David D.; Walker, Eric L.

2004-01-01

This paper is a description of the analysis of blockage corrections for bodies of revolution for the slotted-wall configuration of the National Transonic Facility (NTF) at the NASA Langley Research Center (LaRC). A wall correction method based on the measured wall signature is used. Test data from three different-sized blockage bodies and four wall ventilation settings were analyzed at various Mach numbers and unit Reynolds numbers. The results indicate that with the proper selection of the boundary condition parameters, the wall correction method can predict blockage corrections consistent with the wall measurements for Mach numbers as high as 0.95.

Wall effects in wind tunnels

NASA Technical Reports Server (NTRS)

Chevallier, J. P.; Vaucheret, X.

1986-01-01

A synthesis of current trends in the reduction and computation of wall effects is presented. Some of the points discussed include: (1) for the two-dimensional, transonic tests, various control techniques of boundary conditions are used with adaptive walls offering high precision in determining reference conditions and residual corrections. A reduction in the boundary layer effects of the lateral walls is obtained at T2; (2) for the three-dimensional tests, the methods for the reduction of wall effects are still seldom applied due to a lesser need and to their complexity; (3) the supports holding the model of the probes have to be taken into account in the estimation of perturbatory effects.

Experimental verification of electrostatic boundary conditions in gate-patterned quantum devices

NASA Astrophysics Data System (ADS)

Hou, H.; Chung, Y.; Rughoobur, G.; Hsiao, T. K.; Nasir, A.; Flewitt, A. J.; Griffiths, J. P.; Farrer, I.; Ritchie, D. A.; Ford, C. J. B.

2018-06-01

In a model of a gate-patterned quantum device, it is important to choose the correct electrostatic boundary conditions (BCs) in order to match experiment. In this study, we model gated-patterned devices in doped and undoped GaAs heterostructures for a variety of BCs. The best match is obtained for an unconstrained surface between the gates, with a dielectric region above it and a frozen layer of surface charge, together with a very deep back boundary. Experimentally, we find a  ∼0.2 V offset in pinch-off characteristics of 1D channels in a doped heterostructure before and after etching off a ZnO overlayer, as predicted by the model. Also, we observe a clear quantised current driven by a surface acoustic wave through a lateral induced n-i-n junction in an undoped heterostructure. In the model, the ability to pump electrons in this type of device is highly sensitive to the back BC. Using the improved boundary conditions, it is straightforward to model quantum devices quite accurately using standard software.

Accurate Projection Methods for the Incompressible Navier–Stokes Equations

DOE PAGES

Brown, David L.; Cortez, Ricardo; Minion, Michael L.

2001-04-10

This paper considers the accuracy of projection method approximations to the initial–boundary-value problem for the incompressible Navier–Stokes equations. The issue of how to correctly specify numerical boundary conditions for these methods has been outstanding since the birth of the second-order methodology a decade and a half ago. It has been observed that while the velocity can be reliably computed to second-order accuracy in time and space, the pressure is typically only first-order accurate in the L ∞-norm. Here, we identify the source of this problem in the interplay of the global pressure-update formula with the numerical boundary conditions and presentsmore » an improved projection algorithm which is fully second-order accurate, as demonstrated by a normal mode analysis and numerical experiments. In addition, a numerical method based on a gauge variable formulation of the incompressible Navier–Stokes equations, which provides another option for obtaining fully second-order convergence in both velocity and pressure, is discussed. The connection between the boundary conditions for projection methods and the gauge method is explained in detail.« less

Modeling Plasma Turbulence and Flows in LAPD using BOUT++

NASA Astrophysics Data System (ADS)

Friedman, B.; Carter, T. A.; Schaffner, D.; Popovich, P.; Umansky, M. V.; Dudson, B.

2010-11-01

A Braginskii fluid model of plasma turbulence in the BOUT code has recently been applied to LAPD at UCLA [1]. While these initial simulations with a reduced model and periodic axial boundary conditions have shown good agreement with measurements (e.g. power spectrum, correlation lengths), these simulations have lacked physics essential for modeling self-consistent, quantitatively correct flows. In particular, the model did not contain parallel plasma flow induced by sheath boundary conditions, and the axisymmetric radial electric field was not consistent with experiment. This work addresses these issues by extending the simulation model in the BOUT++ code [2], a more advanced version of BOUT. Specifically, end-plate sheath boundary conditions are added, as well as equations to evolve electron temperature and parallel ion velocity. Finally, various techniques are used to attempt to match the experimental electric potential profile, including fixing an equilibrium profile, fixing the radial boundaries, and adding an angular momentum source. [4pt] [1] Popovich et al., http://arxiv.org/abs/1005.2418 (2010).[0pt] [2] Dudson et al., Computer Physics Communications 180 (2009).

Fractional diffusion on bounded domains

DOE PAGES

Defterli, Ozlem; D'Elia, Marta; Du, Qiang; ...

2015-03-13

We found that the mathematically correct specification of a fractional differential equation on a bounded domain requires specification of appropriate boundary conditions, or their fractional analogue. In this paper we discuss the application of nonlocal diffusion theory to specify well-posed fractional diffusion equations on bounded domains.

Measuring nanometre-scale electric fields in scanning transmission electron microscopy using segmented detectors.

PubMed

Brown, H G; Shibata, N; Sasaki, H; Petersen, T C; Paganin, D M; Morgan, M J; Findlay, S D

2017-11-01

Electric field mapping using segmented detectors in the scanning transmission electron microscope has recently been achieved at the nanometre scale. However, converting these results to quantitative field measurements involves assumptions whose validity is unclear for thick specimens. We consider three approaches to quantitative reconstruction of the projected electric potential using segmented detectors: a segmented detector approximation to differential phase contrast and two variants on ptychographical reconstruction. Limitations to these approaches are also studied, particularly errors arising from detector segment size, inelastic scattering, and non-periodic boundary conditions. A simple calibration experiment is described which corrects the differential phase contrast reconstruction to give reliable quantitative results despite the finite detector segment size and the effects of plasmon scattering in thick specimens. A plasmon scattering correction to the segmented detector ptychography approaches is also given. Avoiding the imposition of periodic boundary conditions on the reconstructed projected electric potential leads to more realistic reconstructions. Copyright © 2017 Elsevier B.V. All rights reserved.

Improvement of structural models using covariance analysis and nonlinear generalized least squares

NASA Technical Reports Server (NTRS)

Glaser, R. J.; Kuo, C. P.; Wada, B. K.

1992-01-01

The next generation of large, flexible space structures will be too light to support their own weight, requiring a system of structural supports for ground testing. The authors have proposed multiple boundary-condition testing (MBCT), using more than one support condition to reduce uncertainties associated with the supports. MBCT would revise the mass and stiffness matrix, analytically qualifying the structure for operation in space. The same procedure is applicable to other common test conditions, such as empty/loaded tanks and subsystem/system level tests. This paper examines three techniques for constructing the covariance matrix required by nonlinear generalized least squares (NGLS) to update structural models based on modal test data. The methods range from a complicated approach used to generate the simulation data (i.e., the correct answer) to a diagonal matrix based on only two constants. The results show that NGLS is very insensitive to assumptions about the covariance matrix, suggesting that a workable NGLS procedure is possible. The examples also indicate that the multiple boundary condition procedure more accurately reduces errors than individual boundary condition tests alone.

Thin wing corrections for phase-change heat-transfer data.

NASA Technical Reports Server (NTRS)

Hunt, J. L.; Pitts, J. I.

1971-01-01

Since no methods are available for determining the magnitude of the errors incurred when the semiinfinite slab assumption is violated, a computer program was developed to calculate the heat-transfer coefficients to both sides of a finite, one-dimensional slab subject to the boundary conditions ascribed to the phase-change coating technique. The results have been correlated in the form of correction factors to the semiinfinite slab solutions in terms of parameters normally used with the technique.

Sound radiation of a railway rail in close proximity to the ground

NASA Astrophysics Data System (ADS)

Zhang, Xianying; Squicciarini, Giacomo; Thompson, David J.

2016-02-01

The sound radiation of a railway in close to proximity to a ground (both rigid and absorptive) is predicted by the boundary element method (BEM) in two dimensions (2D). Results are given in terms of the radiation ratio for both vertical and lateral motion of the rail, when the effects of the acoustic boundary conditions due to the sleepers and ballast are taken into account in the numerical models. Allowance is made for the effect of wave propagation along the rail by applying a correction in the 2D modelling. It is shown that the 2D correction is necessary at low frequency, for both vertical and lateral motion of an unsupported rail, especially in the vicinity of the corresponding critical frequency. However, this correction is not applicable for a supported rail; for vertical motion no correction is needed to the 2D result while for lateral motion the corresponding correction would depend on the pad stiffness. Finally, the corresponding numerical predictions of the sound radiation from a rail are verified by comparison with experimental results obtained using a 1/5 scale rail model in different configurations.

Numerical methods for solving moment equations in kinetic theory of neuronal network dynamics

NASA Astrophysics Data System (ADS)

Rangan, Aaditya V.; Cai, David; Tao, Louis

2007-02-01

Recently developed kinetic theory and related closures for neuronal network dynamics have been demonstrated to be a powerful theoretical framework for investigating coarse-grained dynamical properties of neuronal networks. The moment equations arising from the kinetic theory are a system of (1 + 1)-dimensional nonlinear partial differential equations (PDE) on a bounded domain with nonlinear boundary conditions. The PDEs themselves are self-consistently specified by parameters which are functions of the boundary values of the solution. The moment equations can be stiff in space and time. Numerical methods are presented here for efficiently and accurately solving these moment equations. The essential ingredients in our numerical methods include: (i) the system is discretized in time with an implicit Euler method within a spectral deferred correction framework, therefore, the PDEs of the kinetic theory are reduced to a sequence, in time, of boundary value problems (BVPs) with nonlinear boundary conditions; (ii) a set of auxiliary parameters is introduced to recast the original BVP with nonlinear boundary conditions as BVPs with linear boundary conditions - with additional algebraic constraints on the auxiliary parameters; (iii) a careful combination of two Newton's iterates for the nonlinear BVP with linear boundary condition, interlaced with a Newton's iterate for solving the associated algebraic constraints is constructed to achieve quadratic convergence for obtaining the solutions with self-consistent parameters. It is shown that a simple fixed-point iteration can only achieve a linear convergence for the self-consistent parameters. The practicability and efficiency of our numerical methods for solving the moment equations of the kinetic theory are illustrated with numerical examples. It is further demonstrated that the moment equations derived from the kinetic theory of neuronal network dynamics can very well capture the coarse-grained dynamical properties of integrate-and-fire neuronal networks.

Description of Panel Method Code ANTARES

NASA Technical Reports Server (NTRS)

Ulbrich, Norbert; George, Mike (Technical Monitor)

2000-01-01

Panel method code ANTARES was developed to compute wall interference corrections in a rectangular wind tunnel. The code uses point doublets to represent blockage effects and line doublets to represent lifting effects of a wind tunnel model. Subsonic compressibility effects are modeled by applying the Prandtl-Glauert transformation. The closed wall, open jet, or perforated wall boundary condition may be assigned to a wall panel centroid. The tunnel walls can be represented by using up to 8000 panels. The accuracy of panel method code ANTARES was successfully investigated by comparing solutions for the closed wall and open jet boundary condition with corresponding Method of Images solutions. Fourier transform solutions of a two-dimensional wind tunnel flow field were used to check the application of the perforated wall boundary condition. Studies showed that the accuracy of panel method code ANTARES can be improved by increasing the total number of wall panels in the circumferential direction. It was also shown that the accuracy decreases with increasing free-stream Mach number of the wind tunnel flow field.

CoFFEE: Corrections For Formation Energy and Eigenvalues for charged defect simulations

NASA Astrophysics Data System (ADS)

Naik, Mit H.; Jain, Manish

2018-05-01

Charged point defects in materials are widely studied using Density Functional Theory (DFT) packages with periodic boundary conditions. The formation energy and defect level computed from these simulations need to be corrected to remove the contributions from the spurious long-range interaction between the defect and its periodic images. To this effect, the CoFFEE code implements the Freysoldt-Neugebauer-Van de Walle (FNV) correction scheme. The corrections can be applied to charged defects in a complete range of material shapes and size: bulk, slab (or two-dimensional), wires and nanoribbons. The code is written in Python and features MPI parallelization and optimizations using the Cython package for slow steps.

Finite Element Analysis of Free-Edge Delamination in Laminated Composite Specimens

DTIC Science & Technology

1991-06-18

for the degree of Doctor of Philosophy at the Ohio State University. Revision by H. R. Chu corrected some errors and added further studies on...Galerkin’s approach, in which interlaminar stresses and displacements of each layer satisfying geometrica ’ boundary conditions were represented as -series

Development of a two-equation turbulence model for hypersonic flows. Volume 1; Evaluation of a low Reynolds number correction to the Kappa - epsilon two equation compressible turbulence model

NASA Technical Reports Server (NTRS)

Knight, Doyle D.; Becht, Robert J.

1995-01-01

The objective of the current research is the development of an improved k-epsilon two-equation compressible turbulence model for turbulent boundary layer flows experiencing strong viscous-inviscid interactions. The development of an improved model is important in the design of hypersonic vehicles such as the National Aerospace Plane (NASP) and the High Speed Civil Transport (HSCT). Improvements have been made to the low Reynolds number functions in the eddy viscosity and dissipation of solenoidal dissipation of the k-epsilon turbulence mode. These corrections offer easily applicable modifications that may be utilized for more complex geometries. The low Reynolds number corrections are functions of the turbulent Reynolds number and are therefore independent of the coordinate system. The proposed model offers advantages over some current models which are based upon the physical distance from the wall, that modify the constants of the standard model, or that make more corrections than are necessary to the governing equations. The code has been developed to solve the Favre averaged, boundary layer equations for mass, momentum, energy, turbulence kinetic energy, and dissipation of solenoidal dissipation using Keller's box scheme and the Newton spatial marching method. The code has been validated by removing the turbulent terms and comparing the solution with the Blasius solution, and by comparing the turbulent solution with an existing k-epsilon model code using wall function boundary conditions. Excellent agreement is seen between the computed solution and the Blasius solution, and between the two codes. The model has been tested for both subsonic and supersonic flat-plate turbulent boundary layer flow by comparing the computed skin friction with the Van Driest II theory and the experimental data of Weighardt; by comparing the transformed velocity profile with the data of Weighardt, and the Law of the Wall and the Law of the Wake; and by comparing the computed results of an adverse pressure gradient with the experimental data of Fernando and Smits. Good agreement is obtained with the experimental correlations for all flow conditions.

3D filling simulation of micro- and nanostructures in comparison to iso- and variothermal injection moulding trials

NASA Astrophysics Data System (ADS)

Rytka, C.; Lungershausen, J.; Kristiansen, P. M.; Neyer, A.

2016-06-01

Flow simulations can cut down both costs and time for the development of injection moulded polymer parts with functional surfaces used in life science and optical applications. We simulated the polymer melt flow into 3D micro- and nanostructures with Moldflow and Comsol and compared the results to real iso- and variothermal injection moulding trials below, at and above the transition temperature of the polymer. By adjusting the heat transfer coefficient and the transition temperature in the simulation it was possible to achieve good correlation with experimental findings at different processing conditions (mould temperature, injection velocity) for two polymers, namely polymethylmethacrylate and amorphous polyamide. The macroscopic model can be scaled down in volume and number of elements to save computational time for microstructure simulation and to enable first and foremost the nanostructure simulation, as long as local boundary conditions such as flow front speed are transferred correctly. The heat transfer boundary condition used in Moldflow was further evaluated in Comsol. Results showed that the heat transfer coefficient needs to be increased compared to macroscopic moulding in order to represent interfacial polymer/mould effects correctly. The transition temperature is most important in the packing phase for variothermal injection moulding.

An optimized data fusion method and its application to improve lateral boundary conditions in winter for Pearl River Delta regional PM2.5 modeling, China

NASA Astrophysics Data System (ADS)

Huang, Zhijiong; Hu, Yongtao; Zheng, Junyu; Zhai, Xinxin; Huang, Ran

2018-05-01

Lateral boundary conditions (LBCs) are essential for chemical transport models to simulate regional transport; however they often contain large uncertainties. This study proposes an optimized data fusion approach to reduce the bias of LBCs by fusing gridded model outputs, from which the daughter domain's LBCs are derived, with ground-level measurements. The optimized data fusion approach follows the framework of a previous interpolation-based fusion method but improves it by using a bias kriging method to correct the spatial bias in gridded model outputs. Cross-validation shows that the optimized approach better estimates fused fields in areas with a large number of observations compared to the previous interpolation-based method. The optimized approach was applied to correct LBCs of PM2.5 concentrations for simulations in the Pearl River Delta (PRD) region as a case study. Evaluations show that the LBCs corrected by data fusion improve in-domain PM2.5 simulations in terms of the magnitude and temporal variance. Correlation increases by 0.13-0.18 and fractional bias (FB) decreases by approximately 3%-15%. This study demonstrates the feasibility of applying data fusion to improve regional air quality modeling.

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 5: Propagation of propeller tone noise through a fuselage boundary layer

NASA Technical Reports Server (NTRS)

Magliozzi, B.; Hanson, D. B.

1991-01-01

An analysis of tone noise propagation through a boundary layer and fuselage scattering effects was derived. This analysis is a three dimensional and the complete wave field is solved by matching analytical expressions for the incident and scattered waves in the outer flow to a numerical solution in the boundary layer flow. The outer wave field is constructed analytically from an incident wave appropriate to the source and a scattered wave in the standard Hankel function form. For the incident wave, an existing function - domain propeller noise radiation theory is used. In the boundary layer region, the wave equation is solved by numerical methods. The theoretical analysis is embodied in a computer program which allows the calculation of correction factors for the fuselage scattering and boundary layer refraction effects. The effects are dependent on boundary layer profile, flight speed, and frequency. Corrections can be derived for any point on the fuselage, including those on the opposite side from the source. The theory was verified using limited cases and by comparing calculations with available measurements from JetStar tests of model prop-fans. For the JetStar model scale, the boundary layer refraction effects produce moderate fuselage pressure reinforcements aft of and near the plane of rotation and significant attenuation forward of the plane of rotation at high flight speeds. At lower flight speeds, the calculated boundary layer effects result in moderate amplification over the fuselage area of interest. Apparent amplification forward of the plane of rotation is a result of effective changes in the source directivity due to boundary layer refraction effects. Full scale effects are calculated to be moderate, providing fuselage pressure amplification of about 5 dB at the peak noise location. Evaluation using available noise measurements was made under high-speed, high-altitude flight conditions. Comparisons of calculations made of free field noise, using a current frequency-domain propeller noise prediction method, and fuselage effects using this new procedure show good agreement with fuselage measurements over a wide range of flight speeds and frequencies. Correction factors for the JetStar measurements made on the fuselage are provided in an Appendix.

A Quadtree-gridding LBM with Immersed Boundary for Two-dimension Viscous Flows

NASA Astrophysics Data System (ADS)

Yao, Jieke; Feng, Wenliang; Chen, Bin; Zhou, Wei; Cao, Shikun

2017-07-01

An un-uniform quadtree grids lattice Boltzmann method (LBM) with immersed boundary is presented in this paper. In overlapping for different level grids, temporal and spatial interpolation are necessary to ensure the continuity of physical quantity. In order to take advantage of the equation for temporal and spatial step in the same level grids, equal interval interpolation, which is simple to apply to any refined boundary grids in the LBM, is adopted in temporal and spatial aspects to obtain second-order accuracy. The velocity correction, which can guarantee more preferably no-slip boundary condition than the direct forcing method and the momentum exchange method in the traditional immersed-boundary LBM, is used for solid boundary to make the best of Cartesian grid. In present quadtree-gridding immersed-boundary LBM, large eddy simulation (LES) is adopted to simulate the flows over obstacle in higher Reynolds number (Re). The incompressible viscous flows over circular cylinder are carried out, and a great agreement is obtained.

Including fluid shear viscosity in a structural acoustic finite element model using a scalar fluid representation

PubMed Central

Cheng, Lei; Li, Yizeng; Grosh, Karl

2013-01-01

An approximate boundary condition is developed in this paper to model fluid shear viscosity at boundaries of coupled fluid-structure system. The effect of shear viscosity is approximated by a correction term to the inviscid boundary condition, written in terms of second order in-plane derivatives of pressure. Both thin and thick viscous boundary layer approximations are formulated; the latter subsumes the former. These approximations are used to develop a variational formation, upon which a viscous finite element method (FEM) model is based, requiring only minor modifications to the boundary integral contributions of an existing inviscid FEM model. Since this FEM formulation has only one degree of freedom for pressure, it holds a great computational advantage over the conventional viscous FEM formulation which requires discretization of the full set of linearized Navier-Stokes equations. The results from thick viscous boundary layer approximation are found to be in good agreement with the prediction from a Navier-Stokes model. When applicable, thin viscous boundary layer approximation also gives accurate results with computational simplicity compared to the thick boundary layer formulation. Direct comparison of simulation results using the boundary layer approximations and a full, linearized Navier-Stokes model are made and used to evaluate the accuracy of the approximate technique. Guidelines are given for the parameter ranges over which the accurate application of the thick and thin boundary approximations can be used for a fluid-structure interaction problem. PMID:23729844

Procedures for offline grid nesting in regional ocean models

NASA Astrophysics Data System (ADS)

Mason, Evan; Molemaker, Jeroen; Shchepetkin, Alexander F.; Colas, Francois; McWilliams, James C.; Sangrà, Pablo

One-way offline nesting of a primitive-equation regional ocean numerical model (ROMS) is investigated, with special attention to the boundary forcing file creation process. The model has a modified open boundary condition which minimises false wave reflections, and is optimised to utilise high-frequency boundary updates. The model configuration features a previously computed solution which supplies boundary forcing data to an interior domain with an increased grid resolution. At the open boundaries of the interior grid (the child) the topography is matched to that of the outer grid (the parent), over a narrow transition region. A correction is applied to the normal baroclinic and barotropic velocities at the open boundaries of the child to ensure volume conservation. It is shown that these steps, together with a carefully constructed interpolation of the parent data, lead to a high-quality child solution, with minimal artifacts such as persistent rim currents and wave reflections at the boundaries. Sensitivity experiments provide information about the robustness of the model open boundary condition to perturbations in the surface wind stress forcing field, to the perturbation of the volume conservation enforcement in the boundary forcing, and to perturbation of the vertical density structure in the boundary forcing. This knowledge is important when extending the nesting technique to include external data from alien sources, such as ocean models with physics and/or numerics different from ROMS, or from observed climatologies of temperature, salinity and sea level.

Including fluid shear viscosity in a structural acoustic finite element model using a scalar fluid representation.

PubMed

Cheng, Lei; Li, Yizeng; Grosh, Karl

2013-08-15

An approximate boundary condition is developed in this paper to model fluid shear viscosity at boundaries of coupled fluid-structure system. The effect of shear viscosity is approximated by a correction term to the inviscid boundary condition, written in terms of second order in-plane derivatives of pressure. Both thin and thick viscous boundary layer approximations are formulated; the latter subsumes the former. These approximations are used to develop a variational formation, upon which a viscous finite element method (FEM) model is based, requiring only minor modifications to the boundary integral contributions of an existing inviscid FEM model. Since this FEM formulation has only one degree of freedom for pressure, it holds a great computational advantage over the conventional viscous FEM formulation which requires discretization of the full set of linearized Navier-Stokes equations. The results from thick viscous boundary layer approximation are found to be in good agreement with the prediction from a Navier-Stokes model. When applicable, thin viscous boundary layer approximation also gives accurate results with computational simplicity compared to the thick boundary layer formulation. Direct comparison of simulation results using the boundary layer approximations and a full, linearized Navier-Stokes model are made and used to evaluate the accuracy of the approximate technique. Guidelines are given for the parameter ranges over which the accurate application of the thick and thin boundary approximations can be used for a fluid-structure interaction problem.

Rarefaction and Non-equilibrium Effects in Hypersonic Flows about Leading Edges of Small Bluntness

NASA Astrophysics Data System (ADS)

Ivanov, Mikhail; Khotyanovsky, Dmitry; Kudryavtsev, Alexey; Shershnev, Anton; Bondar, Yevgeniy; Yonemura, Shigeru

2011-05-01

A hypersonic flow about a cylindrically blunted thick plate at a zero angle of attack is numerically studied with the kinetic (DSMC) and continuum (Navier-Stokes equations) approaches. The Navier-Stokes equations with velocity slip and temperature jump boundary conditions correctly predict the flow fields and surface parameters for values of the Knudsen number (based on the radius of leading edge curvature) smaller than 0.1. The results of computations demonstrate significant effects of the entropy layer on the boundary layer characteristics.

Groundwater Flow Model of Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nevada Test Site, Nye County, Nevada, Revision 0

DOE Office of Scientific and Technical Information (OSTI.GOV)

Greg Ruskauff

2006-06-01

The Pahute Mesa groundwater flow model supports the FFACO UGTA corrective action strategy objective of providing an estimate of the vertical and horizontal extent of contaminant migration for each CAU in order to predict contaminant boundaries. A contaminant boundary is the model-predicted perimeter that defines the extent of radionuclide-contaminated groundwater from underground nuclear testing above background conditions exceeding Safe Drinking Water Act (SDWA) standards. The contaminant boundary will be composed of both a perimeter boundary and a lower hydrostratigraphic unit (HSU) boundary. Additional results showing contaminant concentrations and the location of the contaminant boundary at selected times will also bemore » presented. These times may include the verification period, the end of the five-year proof-of-concept period, as well as other times that are of specific interest. The FFACO (1996) requires that the contaminant transport model predict the contaminant boundary at 1,000 years and “at a 95% level of confidence.” The Pahute Mesa Phase I flow model described in this report provides, through the flow fields derived from alternative hydrostratigraphic framework models (HFMs) and recharge models, one part of the data required to compute the contaminant boundary. Other components include the simplified source term model, which incorporates uncertainty and variability in the factors that control radionuclide release from an underground nuclear test (SNJV, 2004a), and the transport model with the concomitant parameter uncertainty as described in Shaw (2003). The uncertainty in all the above model components will be evaluated to produce the final contaminant boundary. This report documents the development of the groundwater flow model for the Central and Western Pahute Mesa CAUs.« less

Direct Validation of the Wall Interference Correction System of the Ames 11-Foot Transonic Wind Tunnel

NASA Technical Reports Server (NTRS)

Ulbrich, Norbert; Boone, Alan R.

2003-01-01

Data from the test of a large semispan model was used to perform a direct validation of a wall interference correction system for a transonic slotted wall wind tunnel. At first, different sets of uncorrected aerodynamic coefficients were generated by physically changing the boundary condition of the test section walls. Then, wall interference corrections were computed and applied to all data points. Finally, an interpolation of the corrected aerodynamic coefficients was performed. This interpolation made sure that the corrected Mach number of a given run would be constant. Overall, the agreement between corresponding interpolated lift, drag, and pitching moment coefficient sets was very good. Buoyancy corrections were also investigated. These studies showed that the accuracy goal of one drag count may only be achieved if reliable estimates of the wall interference induced buoyancy correction are available during a test.

Thomas-Fermi model electron density with correct boundary conditions: Application to atoms and ions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Patil, S.H.

1999-01-01

The author proposes an electron density in atoms and ions, which has the Thomas-Fermi-Dirac form in the intermediate region of r, satisfies the Kato condition for small r, and has the correct asymptotic behavior at large values of r, where r is the distance from the nucleus. He also analyzes the perturbation in the density produced by multipolar fields. He uses these densities in the Poisson equation to deduce average values of r{sup m}, multipolar polarizabilities, and dispersion coefficients of atoms and ions. The predictions are in good agreement with experimental and other theoretical values, generally within about 20%. Hemore » tabulates here the coefficient A in the asymptotic density; radial expectation values (r{sup m}) for m = 2, 4, 6; multipolar polarizabilities {alpha}{sub 1}, {alpha}{sub 2}, {alpha}{sub 3}; expectation values {l_angle}r{sup 0}{r_angle} and {l_angle}r{sup 2}{r_angle} of the asymptotic electron density; and the van der Waals coefficient C{sub 6} for atoms and ions with 2 {le} Z {le} 92. Many of the results, particularly the multipolar polarizabilities and the higher order dispersion coefficients, are the only ones available in the literature. The variation of these properties also provides interesting insight into the shell structure of atoms and ions. Overall, the Thomas-Fermi-Dirac model with the correct boundary conditions provides a good global description of atoms and ions.« less

Aeroheating Measurement of Apollo Shaped Capsule with Boundary Layer Trip in the Free-piston Shock Tunnel HIEST

NASA Technical Reports Server (NTRS)

Hideyuki, TANNO; Tomoyuki, KOMURO; Kazuo, SATO; Katsuhiro, ITOH; Lillard, Randolph P.; Olejniczak, Joseph

2013-01-01

An aeroheating measurement test campaign of an Apollo capsule model with laminar and turbulent boundary layer was performed in the free-piston shock tunnel HIEST at JAXA Kakuda Space Center. A 250mm-diameter 6.4%-scaled Apollo CM capsule model made of SUS-304 stainless steel was applied in this study. To measure heat flux distribution, the model was equipped with 88 miniature co-axial Chromel-Constantan thermocouples on the heat shield surface of the model. In order to promote boundary layer transition, a boundary layer trip insert with 13 "pizza-box" isolated roughness elements, which have 1.27mm square, were placed at 17mm below of the model geometric center. Three boundary layer trip inserts with roughness height of k=0.3mm, 0.6mm and 0.8mm were used to identify the appropriate height to induce transition. Heat flux records with or without roughness elements were obtained for model angles of attack 28º under stagnation enthalpy between H(sub 0)=3.5MJ/kg to 21MJ/kg and stagnation pressure between P(sub 0)=14MPa to 60MPa. Under the condition above, Reynolds number based on the model diameter was varied from 0.2 to 1.3 million. With roughness elements, boundary layer became fully turbulent less than H(sub 0)=9MJ/kg condition. However, boundary layer was still laminar over H(sub 0)=13MJ/kg condition even with the highest roughness elements. An additional experiment was also performed to correct unexpected heat flux augmentation observed over H(sub 0)=9MJ/kg condition.

Predicting Fog in the Nocturnal Boundary Layer

NASA Astrophysics Data System (ADS)

Izett, Jonathan; van de Wiel, Bas; Baas, Peter; van der Linden, Steven; van Hooft, Antoon; Bosveld, Fred

2017-04-01

Fog is a global phenomenon that presents a hazard to navigation and human safety, resulting in significant economic impacts for air and shipping industries as well as causing numerous road traffic accidents. Accurate prediction of fog events, however, remains elusive both in terms of timing and occurrence itself. Statistical methods based on set threshold criteria for key variables such as wind speed have been developed, but high rates of correct prediction of fog events still lead to similarly high "false alarms" when the conditions appear favourable, but no fog forms. Using data from the CESAR meteorological observatory in the Netherlands, we analyze specific cases and perform statistical analyses of event climatology, in order to identify the necessary conditions for correct prediction of fog. We also identify potential "missing ingredients" in current analysis that could help to reduce the number of false alarms. New variables considered include the indicators of boundary layer stability, as well as the presence of aerosols conducive to droplet formation. The poster presents initial findings of new research as well as plans for continued research.

Time domain numerical calculations of unsteady vortical flows about a flat plate airfoil

NASA Technical Reports Server (NTRS)

Hariharan, S. I.; Yu, Ping; Scott, J. R.

1989-01-01

A time domain numerical scheme is developed to solve for the unsteady flow about a flat plate airfoil due to imposed upstream, small amplitude, transverse velocity perturbations. The governing equation for the resulting unsteady potential is a homogeneous, constant coefficient, convective wave equation. Accurate solution of the problem requires the development of approximate boundary conditions which correctly model the physics of the unsteady flow in the far field. A uniformly valid far field boundary condition is developed, and numerical results are presented using this condition. The stability of the scheme is discussed, and the stability restriction for the scheme is established as a function of the Mach number. Finally, comparisons are made with the frequency domain calculation by Scott and Atassi, and the relative strengths and weaknesses of each approach are assessed.

A refined mixed shear flexible finite element for the nonlinear analysis of laminated plates

NASA Technical Reports Server (NTRS)

Putcha, N. S.; Reddy, J. N.

1986-01-01

The present study is concerned with the development of a mixed shear flexible finite element with relaxed continuity for the geometrically linear and nonlinear analysis of laminated anisotropic plates. The formulation of the element is based on a refined higher-order theory. This theory satisfies the zero transverse shear stress boundary conditions on the top and bottom faces of the plate. Shear correction coefficients are not needed. The developed element consists of 11 degrees-of-freedom per node, taking into account three displacements, two rotations, and six moment resultants. An evaluation of the element is conducted with respect to the accuracy obtained in the bending of laminated anistropic rectangular plates with different lamination schemes, loadings, and boundary conditions.

Experimental Verification of the Theory of Wind-Tunnel Boundary Interference

NASA Technical Reports Server (NTRS)

Theodorsen, Theodore; Silverstein, Abe

1935-01-01

The results of an experimental investigation on the boundary-correction factor are presented in this report. The values of the boundary-correction factor from the theory, which at the present time is virtually completed, are given in the report for all conventional types of tunnels. With the isolation of certain disturbing effects, the experimental boundary-correction factor was found to be in satisfactory agreement with the theoretically predicted values, thus verifying the soundness and sufficiency of the theoretical analysis. The establishment of a considerable velocity distortion, in the nature of a unique blocking effect, constitutes a principal result of the investigation.

Influence of the lower boundary in lysimeter observations

NASA Astrophysics Data System (ADS)

Weller, Ulrich; Richter, Katja; Gubis, Jozef; Vogel, Hans-Jörg

2014-05-01

Lysimeters are a valuable tool to study the water household in soils under close to natural conditions. One major drawback is that they are cut off at the lower boundary. This influences strongly the percolation of water. As long as water is leaching down in the soil, it is stagnating at the lower boundary until saturated conditions are reached and the water can percolate through the gravel filter, and under unsaturated conditions there is no flow at all at the lower boundary. In natural soils the water potential at the same depth differs considerably from the regime in a lysimeter. If the depth of the soil or the soil forming substrate is deep enough, the lower boundary is at the potential that allows the percolation of the long term mean of percolation. In other situations, a water table may influence the matric potential in the natural soil, or a less permeable layer may impede free drainage. In all these situations the matric potential at the depth of the lower boundary of the lysimeter will differ substantially in the natural soil. The latest generation of lysimeter therefore has a controlled lower boundary. The matric potential can be actively adjusted to a desired value over a broad range. Most applications connect the suction in the lysimeter to a reference value obtained in the field at the same depth in order to mimic the correct distribution of the soil water. In this presentation we demonstrate the long term influence of the different lower boundary regimes on percolation and evaporation of water based on soil physical models, and we show first field data on the practical implementations with several months of observations.

Dynamic modeling of porous heterogeneous micro/nanobeams

NASA Astrophysics Data System (ADS)

Ebrahimi, Farzad; Jafari, Ali; Reza Barati, Mohammad

2017-12-01

In the present paper, the thermo-mechanical vibration characteristics of a functionally graded (FG) porous microbeam subjected to various types of thermal loadings are investigated based on modified couple stress theory and exact position of neutral axis. The FG micro/nanobeam is modeled via a refined hyperbolic beam theory in which the shear deformation effect is verified without the shear correction factor. A modified power-law distribution which contains porosity volume fraction is used to describe the graded material properties of the FG micro/nanobeam. The temperature field has uniform, linear and nonlinear distributions across the thickness. The governing equations and the related boundary conditions are derived by Hamilton's principle and they are solved applying an analytical solution which satisfies various boundary conditions. A comparison study is performed to verify the present formulation with the known data in the literature and a good agreement is observed. The parametric study covered in this paper includes several parameters, such as thermal loadings, porosity volume fraction, power-law exponents, slenderness ratio, scale parameter and various boundary conditions on natural frequencies of porous FG micro/nanobeams in detail.

A physical-based gas-surface interaction model for rarefied gas flow simulation

NASA Astrophysics Data System (ADS)

Liang, Tengfei; Li, Qi; Ye, Wenjing

2018-01-01

Empirical gas-surface interaction models, such as the Maxwell model and the Cercignani-Lampis model, are widely used as the boundary condition in rarefied gas flow simulations. The accuracy of these models in the prediction of macroscopic behavior of rarefied gas flows is less satisfactory in some cases especially the highly non-equilibrium ones. Molecular dynamics simulation can accurately resolve the gas-surface interaction process at atomic scale, and hence can predict accurate macroscopic behavior. They are however too computationally expensive to be applied in real problems. In this work, a statistical physical-based gas-surface interaction model, which complies with the basic relations of boundary condition, is developed based on the framework of the washboard model. In virtue of its physical basis, this new model is capable of capturing some important relations/trends for which the classic empirical models fail to model correctly. As such, the new model is much more accurate than the classic models, and in the meantime is more efficient than MD simulations. Therefore, it can serve as a more accurate and efficient boundary condition for rarefied gas flow simulations.

Unweighted least squares phase unwrapping by means of multigrid techniques

NASA Astrophysics Data System (ADS)

Pritt, Mark D.

1995-11-01

We present a multigrid algorithm for unweighted least squares phase unwrapping. This algorithm applies Gauss-Seidel relaxation schemes to solve the Poisson equation on smaller, coarser grids and transfers the intermediate results to the finer grids. This approach forms the basis of our multigrid algorithm for weighted least squares phase unwrapping, which is described in a separate paper. The key idea of our multigrid approach is to maintain the partial derivatives of the phase data in separate arrays and to correct these derivatives at the boundaries of the coarser grids. This maintains the boundary conditions necessary for rapid convergence to the correct solution. Although the multigrid algorithm is an iterative algorithm, we demonstrate that it is nearly as fast as the direct Fourier-based method. We also describe how to parallelize the algorithm for execution on a distributed-memory parallel processor computer or a network-cluster of workstations.

How do Stability Corrections Perform in the Stable Boundary Layer Over Snow?

NASA Astrophysics Data System (ADS)

Schlögl, Sebastian; Lehning, Michael; Nishimura, Kouichi; Huwald, Hendrik; Cullen, Nicolas J.; Mott, Rebecca

2017-10-01

We assess sensible heat-flux parametrizations in stable conditions over snow surfaces by testing and developing stability correction functions for two alpine and two polar test sites. Five turbulence datasets are analyzed with respect to, (a) the validity of the Monin-Obukhov similarity theory, (b) the model performance of well-established stability corrections, and (c) the development of new univariate and multivariate stability corrections. Using a wide range of stability corrections reveals an overestimation of the turbulent sensible heat flux for high wind speeds and a generally poor performance of all investigated functions for large temperature differences between snow and the atmosphere above (>10 K). Applying the Monin-Obukhov bulk formulation introduces a mean absolute error in the sensible heat flux of 6 W m^{-2} (compared with heat fluxes calculated directly from eddy covariance). The stability corrections produce an additional error between 1 and 5 W m^{-2}, with the smallest error for published stability corrections found for the Holtslag scheme. We confirm from previous studies that stability corrections need improvements for large temperature differences and wind speeds, where sensible heat fluxes are distinctly overestimated. Under these atmospheric conditions our newly developed stability corrections slightly improve the model performance. However, the differences between stability corrections are typically small when compared to the residual error, which stems from the Monin-Obukhov bulk formulation.

Double-Higgs boson production in the high-energy limit: planar master integrals

NASA Astrophysics Data System (ADS)

Davies, Joshua; Mishima, Go; Steinhauser, Matthias; Wellmann, David

2018-03-01

We consider the virtual corrections to the process gg → HH at NLO in the high energy limit and compute the corresponding planar master integrals in an expansion for small top quark mass. We provide details on the evaluation of the boundary conditions and present analytic results expressed in terms of harmonic polylogarithms.

Estimating surface NO2 and SO2 mixing ratios from fast-response total column observations and potential application to geostationary missions.

PubMed

Knepp, T; Pippin, M; Crawford, J; Chen, G; Szykman, J; Long, R; Cowen, L; Cede, A; Abuhassan, N; Herman, J; Delgado, R; Compton, J; Berkoff, T; Fishman, J; Martins, D; Stauffer, R; Thompson, A M; Weinheimer, A; Knapp, D; Montzka, D; Lenschow, D; Neil, D

Total-column nitrogen dioxide (NO 2 ) data collected by a ground-based sun-tracking spectrometer system (Pandora) and an photolytic-converter-based in-situ instrument collocated at NASA's Langley Research Center in Hampton, Virginia were analyzed to study the relationship between total-column and surface NO 2 measurements. The measurements span more than a year and cover all seasons. Surface mixing ratios are estimated via application of a planetary boundary-layer (PBL) height correction factor. This PBL correction factor effectively corrects for boundary-layer variability throughout the day, and accounts for up to ≈75 % of the variability between the NO 2 data sets. Previous studies have made monthly and seasonal comparisons of column/surface data, which has shown generally good agreement over these long average times. In the current analysis comparisons of column densities averaged over 90 s and 1 h are made. Applicability of this technique to sulfur dioxide (SO 2 ) is briefly explored. The SO 2 correlation is improved by excluding conditions where surface levels are considered background. The analysis is extended to data from the July 2011 DISCOVER-AQ mission over the greater Baltimore, MD area to examine the method's performance in more-polluted urban conditions where NO 2 concentrations are typically much higher.

Atmospheric monitoring in MAGIC and data corrections

NASA Astrophysics Data System (ADS)

Fruck, Christian; Gaug, Markus

2015-03-01

A method for analyzing returns of a custom-made "micro"-LIDAR system, operated alongside the two MAGIC telescopes is presented. This method allows for calculating the transmission through the atmospheric boundary layer as well as thin cloud layers. This is achieved by applying exponential fits to regions of the back-scattering signal that are dominated by Rayleigh scattering. Making this real-time transmission information available for the MAGIC data stream allows to apply atmospheric corrections later on in the analysis. Such corrections allow for extending the effective observation time of MAGIC by including data taken under adverse atmospheric conditions. In the future they will help reducing the systematic uncertainties of energy and flux.

Bounded solutions in a T-shaped waveguide and the spectral properties of the Dirichlet ladder

NASA Astrophysics Data System (ADS)

Nazarov, S. A.

2014-08-01

The Dirichlet problem is considered on the junction of thin quantum waveguides (of thickness h ≪ 1) in the shape of an infinite two-dimensional ladder. Passage to the limit as h → +0 is discussed. It is shown that the asymptotically correct transmission conditions at nodes of the corresponding one-dimensional quantum graph are Dirichlet conditions rather than the conventional Kirchhoff transmission conditions. The result is obtained by analyzing bounded solutions of a problem in the T-shaped waveguide that the boundary layer phenomenon.

High-order Two-way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering

NASA Technical Reports Server (NTRS)

Fibich, Gadi; Tsynkov, Semyon

2000-01-01

When solving linear scattering problems, one typically first solves for the impinging wave in the absence of obstacles. Then, by linear superposition, the original problem is reduced to one that involves only the scattered waves driven by the values of the impinging field at the surface of the obstacles. In addition, when the original domain is unbounded, special artificial boundary conditions (ABCs) that would guarantee the reflectionless propagation of waves have to be set at the outer boundary of the finite computational domain. The situation becomes conceptually different when the propagation equation is nonlinear. In this case the impinging and scattered waves can no longer be separated, and the problem has to be solved in its entirety. In particular, the boundary on which the incoming field values are prescribed, should transmit the given incoming waves in one direction and simultaneously be transparent to all the outgoing waves that travel in the opposite direction. We call this type of boundary conditions two-way ABCs. In the paper, we construct the two-way ABCs for the nonlinear Helmholtz equation that models the laser beam propagation in a medium with nonlinear index of refraction. In this case, the forward propagation is accompanied by backscattering, i.e., generation of waves in the direction opposite to that of the incoming signal. Our two-way ABCs generate no reflection of the backscattered waves and at the same time impose the correct values of the incoming wave. The ABCs are obtained for a fourth-order accurate discretization to the Helmholtz operator; the fourth-order grid convergence is corroborated experimentally by solving linear model problems. We also present solutions in the nonlinear case using the two-way ABC which, unlike the traditional Dirichlet boundary condition, allows for direct calculation of the magnitude of backscattering.

Phase I Flow and Transport Model Document for Corrective Action Unit 97: Yucca Flat/Climax Mine, Nevada National Security Site, Nye County, Nevada, Revision 1 with ROTCs 1 and 2

DOE Office of Scientific and Technical Information (OSTI.GOV)

Andrews, Robert

The Underground Test Area (UGTA) Corrective Action Unit (CAU) 97, Yucca Flat/Climax Mine, in the northeast part of the Nevada National Security Site (NNSS) requires environmental corrective action activities to assess contamination resulting from underground nuclear testing. These activities are necessary to comply with the UGTA corrective action strategy (referred to as the UGTA strategy). The corrective action investigation phase of the UGTA strategy requires the development of groundwater flow and contaminant transport models whose purpose is to identify the lateral and vertical extent of contaminant migration over the next 1,000 years. In particular, the goal is to calculate themore » contaminant boundary, which is defined as a probabilistic model-forecast perimeter and a lower hydrostratigraphic unit (HSU) boundary that delineate the possible extent of radionuclide-contaminated groundwater from underground nuclear testing. Because of structural uncertainty in the contaminant boundary, a range of potential contaminant boundaries was forecast, resulting in an ensemble of contaminant boundaries. The contaminant boundary extent is determined by the volume of groundwater that has at least a 5 percent chance of exceeding the radiological standards of the Safe Drinking Water Act (SDWA) (CFR, 2012).« less

Viscous wing theory development. Volume 1: Analysis, method and results

NASA Technical Reports Server (NTRS)

Chow, R. R.; Melnik, R. E.; Marconi, F.; Steinhoff, J.

1986-01-01

Viscous transonic flows at large Reynolds numbers over 3-D wings were analyzed using a zonal viscid-inviscid interaction approach. A new numerical AFZ scheme was developed in conjunction with the finite volume formulation for the solution of the inviscid full-potential equation. A special far-field asymptotic boundary condition was developed and a second-order artificial viscosity included for an improved inviscid solution methodology. The integral method was used for the laminar/turbulent boundary layer and 3-D viscous wake calculation. The interaction calculation included the coupling conditions of the source flux due to the wing surface boundary layer, the flux jump due to the viscous wake, and the wake curvature effect. A method was also devised incorporating the 2-D trailing edge strong interaction solution for the normal pressure correction near the trailing edge region. A fully automated computer program was developed to perform the proposed method with one scalar version to be used on an IBM-3081 and two vectorized versions on Cray-1 and Cyber-205 computers.

77 FR 49775 - Boundary Technical Correction for the McKenzie Wild and Scenic River “Lower and Upper Portions...

Federal Register 2010, 2011, 2012, 2013, 2014

2012-08-17

... DEPARTMENT OF AGRICULTURE Forest Service Boundary Technical Correction for the McKenzie Wild and... Service, USDA. ACTION: Notice of availability. SUMMARY: In accordance with Section 3(b) of the Wild and... for the boundary of the McKenzie Wild and Scenic River ``Lower Portion'' and ``Upper Portion'' to...

Elastic modelling in tilted transversely isotropic media with convolutional perfectly matched layer boundary conditions

NASA Astrophysics Data System (ADS)

Han, Byeongho; Seol, Soon Jee; Byun, Joongmoo

2012-04-01

To simulate wave propagation in a tilted transversely isotropic (TTI) medium with a tilting symmetry-axis of anisotropy, we develop a 2D elastic forward modelling algorithm. In this algorithm, we use the staggered-grid finite-difference method which has fourth-order accuracy in space and second-order accuracy in time. Since velocity-stress formulations are defined for staggered grids, we include auxiliary grid points in the z-direction to meet the free surface boundary conditions for shear stress. Through comparisons of displacements obtained from our algorithm, not only with analytical solutions but also with finite element solutions, we are able to validate that the free surface conditions operate appropriately and elastic waves propagate correctly. In order to handle the artificial boundary reflections efficiently, we also implement convolutional perfectly matched layer (CPML) absorbing boundaries in our algorithm. The CPML sufficiently attenuates energy at the grazing incidence by modifying the damping profile of the PML boundary. Numerical experiments indicate that the algorithm accurately expresses elastic wave propagation in the TTI medium. At the free surface, the numerical results show good agreement with analytical solutions not only for body waves but also for the Rayleigh wave which has strong amplitude along the surface. In addition, we demonstrate the efficiency of CPML for a homogeneous TI medium and a dipping layered model. Only using 10 grid points to the CPML regions, the artificial reflections are successfully suppressed and the energy of the boundary reflection back into the effective modelling area is significantly decayed.

Thomas-Fermi approximation for a condensate with higher-order interactions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thoegersen, M.; Jensen, A. S.; Zinner, N. T.

We consider the ground state of a harmonically trapped Bose-Einstein condensate within the Gross-Pitaevskii theory including the effective-range corrections for a two-body zero-range potential. The resulting nonlinear Schroedinger equation is solved analytically in the Thomas-Fermi approximation neglecting the kinetic-energy term. We present results for the chemical potential and the condensate profiles, discuss boundary conditions, and compare to the usual Thomas-Fermi approach. We discuss several ways to increase the influence of effective-range corrections in experiment with magnetically tunable interactions. The level of tuning required could be inside experimental reach in the near future.

A new ghost-node method for linking different models and initial investigations of heterogeneity and nonmatching grids

USGS Publications Warehouse

Dickinson, J.E.; James, S.C.; Mehl, S.; Hill, M.C.; Leake, S.A.; Zyvoloski, G.A.; Faunt, C.C.; Eddebbarh, A.-A.

2007-01-01

A flexible, robust method for linking parent (regional-scale) and child (local-scale) grids of locally refined models that use different numerical methods is developed based on a new, iterative ghost-node method. Tests are presented for two-dimensional and three-dimensional pumped systems that are homogeneous or that have simple heterogeneity. The parent and child grids are simulated using the block-centered finite-difference MODFLOW and control-volume finite-element FEHM models, respectively. The models are solved iteratively through head-dependent (child model) and specified-flow (parent model) boundary conditions. Boundary conditions for models with nonmatching grids or zones of different hydraulic conductivity are derived and tested against heads and flows from analytical or globally-refined models. Results indicate that for homogeneous two- and three-dimensional models with matched grids (integer number of child cells per parent cell), the new method is nearly as accurate as the coupling of two MODFLOW models using the shared-node method and, surprisingly, errors are slightly lower for nonmatching grids (noninteger number of child cells per parent cell). For heterogeneous three-dimensional systems, this paper compares two methods for each of the two sets of boundary conditions: external heads at head-dependent boundary conditions for the child model are calculated using bilinear interpolation or a Darcy-weighted interpolation; specified-flow boundary conditions for the parent model are calculated using model-grid or hydrogeologic-unit hydraulic conductivities. Results suggest that significantly more accurate heads and flows are produced when both Darcy-weighted interpolation and hydrogeologic-unit hydraulic conductivities are used, while the other methods produce larger errors at the boundary between the regional and local models. The tests suggest that, if posed correctly, the ghost-node method performs well. Additional testing is needed for highly heterogeneous systems. ?? 2007 Elsevier Ltd. All rights reserved.

Thermal properties of composite materials : effective conductivity tensor and edge effects

NASA Astrophysics Data System (ADS)

Matine, A.; Boyard, N.; Cartraud, P.; Legrain, G.; Jarny, Y.

2012-11-01

The homogenization theory is a powerful approach to determine the effective thermal conductivity tensor of heterogeneous materials such as composites, including thermoset matrix and fibres. Once the effective properties are calculated, they can be used to solve a heat conduction problem on the composite structure at the macroscopic scale. This approach leads to good approximations of both the heat flux and temperature in the interior zone of the structure, however edge effects occur in the vicinity of the domain boundaries. In this paper, following the approach proposed in [10] for elasticity, it is shown how these edge effects can be corrected. Thus an additional asymptotic expansion is introduced, which plays the role of a edge effect term. This expansion tends to zero far from the boundary, and is assumed to decrease exponentially. Moreover, the length of the edge effect region can be determined from the solution of an eigenvalue problem. Numerical examples are considered for a standard multilayered material. The homogenized solutions computed with a finite element software, and corrected with the edge effect terms, are compared to a heterogeneous finite element solution at the microscopic scale. The influences of the thermal contrast and scale factor are illustrated for different kind of boundary conditions.

The two sides of the C-factor.

PubMed

Fok, Alex S L; Aregawi, Wondwosen A

2018-04-01

The aim of this paper is to investigate the effects on shrinkage strain/stress development of the lateral constraints at the bonded surfaces of resin composite specimens used in laboratory measurement. Using three-dimensional (3D) Hooke's law, a recently developed shrinkage stress theory is extended to 3D to include the additional out-of-plane strain/stress induced by the lateral constraints at the bonded surfaces through the Poisson's ratio effect. The model contains a parameter that defines the relative thickness of the boundary layers, adjacent to the bonded surfaces, that are under such multiaxial stresses. The resulting differential equation is solved for the shrinkage stress under different boundary conditions. The accuracy of the model is assessed by comparing the numerical solutions with a wide range of experimental data, which include those from both shrinkage strain and shrinkage stress measurements. There is good agreement between theory and experiments. The model correctly predicts the different instrument-dependent effects that a specimen's configuration factor (C-factor) has on shrinkage stress. That is, for noncompliant stress-measuring instruments, shrinkage stress increases with the C-factor of the cylindrical specimen; while the opposite is true for compliant instruments. The model also provides a correction factor, which is a function of the C-factor, Poisson's ratio and boundary layer thickness of the specimen, for shrinkage strain measured using the bonded-disc method. For the resin composite examined, the boundary layers have a combined thickness that is ∼11.5% of the specimen's diameter. The theory provides a physical and mechanical basis for the C-factor using principles of engineering mechanics. The correction factor it provides allows the linear shrinkage strain of a resin composite to be obtained more accurately from the bonded-disc method. Published by Elsevier Ltd.

Geometric and shading correction for images of printed materials using boundary.

PubMed

Brown, Michael S; Tsoi, Yau-Chat

2006-06-01

A novel technique that uses boundary interpolation to correct geometric distortion and shading artifacts present in images of printed materials is presented. Unlike existing techniques, our algorithm can simultaneously correct a variety of geometric distortions, including skew, fold distortion, binder curl, and combinations of these. In addition, the same interpolation framework can be used to estimate the intrinsic illumination component of the distorted image to correct shading artifacts. We detail our algorithm for geometric and shading correction and demonstrate its usefulness on real-world and synthetic data.

Using a Quasipotential Transformation for Modeling Diffusion Media inPolymer-Electrolyte Fuel Cells

DOE Office of Scientific and Technical Information (OSTI.GOV)

Weber, Adam Z.; Newman, John

2008-08-29

In this paper, a quasipotential approach along with conformal mapping is used to model the diffusion media of a polymer-electrolyte fuel cell. This method provides a series solution that is grid independent and only requires integration along a single boundary to solve the problem. The approach accounts for nonisothermal phenomena, two-phase flow, correct placement of the electronic potential boundary condition, and multilayer media. The method is applied to a cathode diffusion medium to explore the interplay between water and thermal management and performance, the impact of the rib-to-channel ratio, and the existence of diffusion under the rib and flooding phenomena.

Data pre-processing: Stratospheric aerosol perturbing effect on the remote sensing of vegetation: Correction method for the composite NDVI after the Pinatubo eruption

NASA Technical Reports Server (NTRS)

Vermote, E.; Elsaleous, N.; Kaufman, Y. J.; Dutton, E.

1994-01-01

An operational stratospheric correction scheme used after the Mount Pinatubo (Phillipines) eruption (Jun. 1991) is presented. The stratospheric aerosol distribution is assumed to be only variable with latitude. Each 9 days the latitudinal distribution of the optical thickness is computed by inverting radiances observed in the NOAA AVHRR channel 1 (0.63 micrometers) and channel 2 (0.83 micrometers) over the Pacific Ocean. This radiance data set is used to check the validity of model used for inversion by checking consistency of the optical thickness deduced from each channel as well as optical thickness deduced from different scattering angles. Using the optical thickness profile previously computed and radiative transfer code assuming Lambertian boundary condition, each pixel of channel 1 and 2 are corrected prior to computation of NDVI (Normalized Difference Vegetation Index). Comparison between corrected, non corrected, and years prior to Pinatubo eruption (1989 to 1990) NDVI composite, shows the necessity and the accuracy of the operational correction scheme.

The choice of boundary conditions and mesh for scaffolding FEM model on the basis of natural vibrations measurements

NASA Astrophysics Data System (ADS)

Cyniak, Patrycja; Błazik-Borowa, Ewa; Szer, Jacek; Lipecki, Tomasz; Szer, Iwona

2018-01-01

Scaffolding is a specific construction with high susceptibility to low frequency vibrations. The numerical model of scaffolding presented in this paper contains real imperfections received from geodetic measurements of real construction. Boundary conditions were verified on the basis of measured free vibrations. A simulation of a man walking on penultimate working level as a dynamic load variable in time was made for verified model. The paper presents procedure for a choice of selected parameters of the scaffolding FEM model. The main aim of analysis is the best projection of the real construction and correct modeling of worker walking on the scaffolding. Different boundary conditions are considered, because of their impact on construction vibrations. Natural vibrations obtained from FEM calculations are compared with free vibrations measured during in-situ tests. Structure accelerations caused by walking human are then considered in this paper. Methodology of creating numerical models of scaffoldings and analysis of dynamic effects during human walking are starting points for further considerations about dynamic loads acting on such structures and effects of these loads to construction and workers, whose workplaces are situated on the scaffolding.

Stability and natural vibration analysis of laminated plates by using a mixed element based on a refined plate theory

NASA Technical Reports Server (NTRS)

Putcha, N. S.; Reddy, J. N.

1986-01-01

A mixed shear flexible finite element, with relaxed continuity, is developed for the geometrically linear and nonlinear analysis of layered anisotropic plates. The element formulation is based on a refined higher order theory which satisfies the zero transverse shear stress boundary conditions on the top and bottom faces of the plate and requires no shear correction coefficients. The mixed finite element developed herein consists of eleven degrees of freedom per node which include three displacements, two rotations and six moment resultants. The element is evaluated for its accuracy in the analysis of the stability and vibration of anisotropic rectangular plates with different lamination schemes and boundary conditions. The mixed finite element described here for the higher order theory gives very accurate results for buckling loads and natural frequencies.

The Schwinger Variational Method

NASA Technical Reports Server (NTRS)

Huo, Winifred M.

1995-01-01

Variational methods have proven invaluable in theoretical physics and chemistry, both for bound state problems and for the study of collision phenomena. The application of the Schwinger variational (SV) method to e-molecule collisions and molecular photoionization has been reviewed previously. The present chapter discusses the implementation of the SV method as applied to e-molecule collisions. Since this is not a review of cross section data, cross sections are presented only to server as illustrative examples. In the SV method, the correct boundary condition is automatically incorporated through the use of Green's function. Thus SV calculations can employ basis functions with arbitrary boundary conditions. The iterative Schwinger method has been used extensively to study molecular photoionization. For e-molecule collisions, it is used at the static exchange level to study elastic scattering and coupled with the distorted wave approximation to study electronically inelastic scattering.

Folded supersymmetry with a twist

DOE PAGES

Cohen, Timothy; Craig, Nathaniel; Lou, Hou Keong; ...

2016-03-30

Folded supersymmetry (f-SUSY) stabilizes the weak scale against radiative corrections from the top sector via scalar partners whose gauge quantum numbers differ from their Standard Model counterparts. This non-trivial pairing of states can be realized in extra-dimensional theories with appropriate supersymmetry-breaking boundary conditions. We present a class of calculable f-SUSY models that are parametrized by a non-trivial twist in 5D boundary conditions and can accommodate the observed Higgs mass and couplings. Although the distinctive phenomenology associated with the novel folded states should provide strong evidence for this mechanism, the most stringent constraints are currently placed by conventional supersymmetry searches. Asmore » a result, these models remain minimally fine-tuned in light of LHC8 data and provide a range of both standard and exotic signatures accessible at LHC13.« less

Theoretical analysis and experimental study of constraint boundary conditions for acquiring the beacon in satellite-ground laser communications

NASA Astrophysics Data System (ADS)

Yu, Siyuan; Wu, Feng; Wang, Qiang; Tan, Liying; Ma, Jing

2017-11-01

Acquisition and recognition for the beacon is the core technology of establishing the satellite optical link. In order to acquire the beacon correctly, the beacon image should be recognized firstly, excluding the influence of the background light. In this processing, many factors will influence the recognition precision of the beacon. This paper studies the constraint boundary conditions for acquiring the beacon from the perspective of theory and experiment, and as satellite-ground laser communications, an approach for obtaining the adaptive segmentation method is also proposed. Finally, the long distance laser communication experiment (11.16 km) verifies the validity of this method and the tracking error with the method is the least compared with the traditional approaches. The method helps to greatly improve the tracking precision in the satellite-ground laser communications.

Optimal control in adaptive optics modeling of nonlinear systems

NASA Astrophysics Data System (ADS)

Herrmann, J.

The problem of using an adaptive optics system to correct for nonlinear effects like thermal blooming is addressed using a model containing nonlinear lenses through which Gaussian beams are propagated. The best correction of this nonlinear system can be formulated as a deterministic open loop optimal control problem. This treatment gives a limit for the best possible correction. Aspects of adaptive control and servo systems are not included at this stage. An attempt is made to determine that control in the transmitter plane which minimizes the time averaged area or maximizes the fluence in the target plane. The standard minimization procedure leads to a two-point-boundary-value problem, which is ill-conditioned in the case. The optimal control problem was solved using an iterative gradient technique. An instantaneous correction is introduced and compared with the optimal correction. The results of the calculations show that for short times or weak nonlinearities the instantaneous correction is close to the optimal correction, but that for long times and strong nonlinearities a large difference develops between the two types of correction. For these cases the steady state correction becomes better than the instantaneous correction and approaches the optimum correction.

A groundwater data assimilation application study in the Heihe mid-reach

NASA Astrophysics Data System (ADS)

Ragettli, S.; Marti, B. S.; Wolfgang, K.; Li, N.

2017-12-01

The present work focuses on modelling of the groundwater flow in the mid-reach of the endorheic river Heihe in the Zhangye oasis (Gansu province) in arid north-west China. In order to optimise the water resources management in the oasis, reliable forecasts of groundwater level development under different management options and environmental boundary conditions have to be produced. For this means, groundwater flow is modelled with Modflow and coupled to an Ensemble Kalman Filter programmed in Matlab. The model is updated with monthly time steps, featuring perturbed boundary conditions to account for uncertainty in model forcing. Constant biases between model and observations have been corrected prior to updating and compared to model runs without bias correction. Different options for data assimilation (states and/or parameters), updating frequency, and measures against filter inbreeding (damping factor, covariance inflation, spatial localization) have been tested against each other. Results show a high dependency of the Ensemble Kalman filter performance on the selection of observations for data assimilation. For the present regional model, bias correction is necessary for a good filter performance. A combination of spatial localization and covariance inflation is further advisable to reduce filter inbreeding problems. Best performance is achieved if parameter updates are not large, an indication for good prior model calibration. Asynchronous updating of parameter values once every five years (with data of the past five years) and synchronous updating of the groundwater levels is better suited for this groundwater system with not or slow changing parameter values than synchronous updating of both groundwater levels and parameters at every time step applying a damping factor. The filter is not able to correct time lags of signals.

Efficient algorithms for analyzing the singularly perturbed boundary value problems of fractional order

NASA Astrophysics Data System (ADS)

Sayevand, K.; Pichaghchi, K.

2018-04-01

In this paper, we were concerned with the description of the singularly perturbed boundary value problems in the scope of fractional calculus. We should mention that, one of the main methods used to solve these problems in classical calculus is the so-called matched asymptotic expansion method. However we shall note that, this was not achievable via the existing classical definitions of fractional derivative, because they do not obey the chain rule which one of the key elements of the matched asymptotic expansion method. In order to accommodate this method to fractional derivative, we employ a relatively new derivative so-called the local fractional derivative. Using the properties of local fractional derivative, we extend the matched asymptotic expansion method to the scope of fractional calculus and introduce a reliable new algorithm to develop approximate solutions of the singularly perturbed boundary value problems of fractional order. In the new method, the original problem is partitioned into inner and outer solution equations. The reduced equation is solved with suitable boundary conditions which provide the terminal boundary conditions for the boundary layer correction. The inner solution problem is next solved as a solvable boundary value problem. The width of the boundary layer is approximated using appropriate resemblance function. Some theoretical results are established and proved. Some illustrating examples are solved and the results are compared with those of matched asymptotic expansion method and homotopy analysis method to demonstrate the accuracy and efficiency of the method. It can be observed that, the proposed method approximates the exact solution very well not only in the boundary layer, but also away from the layer.

Areal-averaged trace gas emission rates from long-range open-path measurements in stable boundary layer conditions

NASA Astrophysics Data System (ADS)

Schäfer, K.; Grant, R. H.; Emeis, S.; Raabe, A.; von der Heide, C.; Schmid, H. P.

2012-07-01

Measurements of land-surface emission rates of greenhouse and other gases at large spatial scales (10 000 m2) are needed to assess the spatial distribution of emissions. This can be readily done using spatial-integrating micro-meteorological methods like flux-gradient methods which were evaluated for determining land-surface emission rates of trace gases under stable boundary layers. Non-intrusive path-integrating measurements are utilized. Successful application of a flux-gradient method requires confidence in the gradients of trace gas concentration and wind, and in the applicability of boundary-layer turbulence theory; consequently the procedures to qualify measurements that can be used to determine the flux is critical. While there is relatively high confidence in flux measurements made under unstable atmospheres with mean winds greater than 1 m s-1, there is greater uncertainty in flux measurements made under free convective or stable conditions. The study of N2O emissions of flat grassland and NH3 emissions from a cattle lagoon involves quality-assured determinations of fluxes under low wind, stable or night-time atmospheric conditions when the continuous "steady-state" turbulence of the surface boundary layer breaks down and the layer has intermittent turbulence. Results indicate that following the Monin-Obukhov similarity theory (MOST) flux-gradient methods that assume a log-linear profile of the wind speed and concentration gradient incorrectly determine vertical profiles and thus flux in the stable boundary layer. An alternative approach is considered on the basis of turbulent diffusivity, i.e. the measured friction velocity as well as height gradients of horizontal wind speeds and concentrations without MOST correction for stability. It is shown that this is the most accurate of the flux-gradient methods under stable conditions.

Inflow/Outflow Conditions for Unsteady Aerodynamics and Aeroacoustics in Nonuniform Flow

NASA Technical Reports Server (NTRS)

Atassi, Oliver V.; Grady, Joseph E. (Technical Monitor)

2003-01-01

The effect of a nonuniform mean flow on the normal modes; the inflow/outflow nonreflecting boundary conditions; and the sound power are studied. The normal modes in an annular duct are computed using a spectral method in combination with a shooting method. The swirl causes force imbalance which couples the acoustic and vortical modes. The acoustic modes are distinguished from the vortical modes by their large pressure and small vorticity content. The mean swirl also produces a Doppler shift in frequency. This results in more counter-spinning modes cut-on at a given frequency than modes spinning with the swirl. Nonreflecting boundary conditions are formulated using the normal mode solutions. The inflow/outflow boundary conditions are implemented in a linearized Euler scheme and validated by computing the propagation of acoustic and vortical waves in a duct for a variety of swirling mean flows. Numerical results show that the evolution of the vortical disturbances is sensitive to the inflow conditions and the details of the wake excitations. All three components of the wake velocity must be considered to correctly compute the wake evolution and the blade upwash. For high frequencies, the acoustic-vortical mode coupling is weak and a conservation equation for the acoustic energy can be derived. Sound power calculations show significant mean flow swirl effects, but mode interference effects are small.

Effect of aspect ratio on sidewall boundary-layer influence in two-dimensional airfoil testing

NASA Technical Reports Server (NTRS)

Murthy, A. V.

1986-01-01

The effect of sidewall boundary layers in airfoil testing in two-dimensional wind tunnels is investigated. The non-linear crossflow velocity variation induced because of the changes in the sidewall boundary-layer thickness is represented by the flow between a wavy wall and a straight wall. Using this flow model, a correction for the sidewall boundary-layer effects is derived in terms of the undisturbed sidewall boundary-layer properties, the test Mach number and the airfoil aspect ratio. Application of the proposed correction to available experimental data showed good correlation for the shock location and pressure distribution on airfoils.

Wall function treatment for bubbly boundary layers at low void fractions.

PubMed

Soares, Daniel V; Bitencourt, Marcelo C; Loureiro, Juliana B R; Silva Freire, Atila P

2018-01-01

The present work investigates the role of different treatments of the lower boundary condition on the numerical prediction of bubbly flows. Two different wall function formulations are tested against experimental data obtained for bubbly boundary layers: (i) a new analytical solution derived through asymptotic techniques and (ii) the previous formulation of Troshko and Hassan (IJHMT, 44, 871-875, 2001a). A modified k-e model is used to close the averaged Navier-Stokes equations together with the hypothesis that turbulence can be modelled by a linear superposition of bubble and shear induced eddy viscosities. The work shows, in particular, how four corrections must the implemented in the standard single-phase k-e model to account for the effects of bubbles. The numerical implementation of the near wall functions is made through a finite elements code.

Vacuum energy in Einstein-Gauss-Bonnet anti-de Sitter gravity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kofinas, Georgios; Olea, Rodrigo

2006-10-15

A finite action principle for Einstein-Gauss-Bonnet anti-de Sitter gravity is achieved by supplementing the bulk Lagrangian by a suitable boundary term, whose form substantially differs in odd and even dimensions. For even dimensions, this term is given by the boundary contribution in the Euler theorem with a coupling constant fixed, demanding the spacetime to have constant (negative) curvature in the asymptotic region. For odd dimensions, the action is stationary under a boundary condition on the variation of the extrinsic curvature. A well-posed variational principle leads to an appropriate definition of energy and other conserved quantities using the Noether theorem, andmore » to a correct description of black hole thermodynamics. In particular, this procedure assigns a nonzero energy to anti-de Sitter spacetime in all odd dimensions.« less

Error analysis and correction in wavefront reconstruction from the transport-of-intensity equation

PubMed Central

Barbero, Sergio; Thibos, Larry N.

2007-01-01

Wavefront reconstruction from the transport-of-intensity equation (TIE) is a well-posed inverse problem given smooth signals and appropriate boundary conditions. However, in practice experimental errors lead to an ill-condition problem. A quantitative analysis of the effects of experimental errors is presented in simulations and experimental tests. The relative importance of numerical, misalignment, quantization, and photodetection errors are shown. It is proved that reduction of photodetection noise by wavelet filtering significantly improves the accuracy of wavefront reconstruction from simulated and experimental data. PMID:20052302

Viscous pressure correction in the irrotational flow outside Prandtl's boundary layer

NASA Astrophysics Data System (ADS)

Joseph, Daniel; Wang, Jing

2004-11-01

We argue that boundary layers on solid with irrotational motion outside are like a gas bubble because the shear stress vanishes at the edge of the boundary layer but the irrotational shear stress does not. This discrepancy induces a pressure correction and an additional drag which can be advertised as due to the viscous dissipation of the irrotational flow. Typically, this extra correction to the drag would be relatively small. A much more interesting implication of the extra pressure theory arises from the consideration of the effects of viscosity on the normal stress on a solid boundary which are entirely neglected in Prandtl's theory. It is very well known and easily demonstrated that as a consequence of the continuity equation the viscous normal stress must vanish on a rigid solid. It follows that all the greatly important effects of viscosity on the normal stress are buried in the pressure and the leading order effects of viscosity on the normal stress can be obtained from the viscous correction of viscous potential flow.

Non-CMC solutions to the Einstein constraint equations on asymptotically Euclidean manifolds with apparent horizon boundaries

NASA Astrophysics Data System (ADS)

Holst, Michael; Meier, Caleb

2015-01-01

In this article we further develop the solution theory for the Einstein constraint equations on an n-dimensional, asymptotically Euclidean manifold M with interior boundary Σ. Building on recent results for both the asymptotically Euclidean and compact with boundary settings, we show the existence of far-from-CMC and near-CMC solutions to the conformal formulation of the Einstein constraints when nonlinear Robin boundary conditions are imposed on Σ, similar to those analyzed previously by Dain (2004 Class. Quantum Grav. 21 555-73), by Maxwell (2004, 2005 Commun. Math. Phys. 253 561-83), and by Holst and Tsogtgerel (2013 Class. Quantum Grav. 30 205011) as a model of black holes in various CMC settings, and by Holst et al (2013 Non-CMC solutions to the einstein constraint equations with apparent horizon boundaries arXiv:1310.2302v1) in the setting of far-from-CMC solutions on compact manifolds with boundary. These ‘marginally trapped surface’ Robin conditions ensure that the expansion scalars along null geodesics perpendicular to the boundary region Σ are non-positive, which is considered the correct mathematical model for black holes in the context of the Einstein constraint equations. Assuming a suitable form of weak cosmic censorship, the results presented in this article guarantee the existence of initial data that will evolve into a space-time containing an arbitrary number of black holes. A particularly important feature of our results are the minimal restrictions we place on the mean curvature, giving both near- and far-from-CMC results that are new.

An immersed boundary-simplified sphere function-based gas kinetic scheme for simulation of 3D incompressible flows

NASA Astrophysics Data System (ADS)

Yang, L. M.; Shu, C.; Yang, W. M.; Wang, Y.; Wu, J.

2017-08-01

In this work, an immersed boundary-simplified sphere function-based gas kinetic scheme (SGKS) is presented for the simulation of 3D incompressible flows with curved and moving boundaries. At first, the SGKS [Yang et al., "A three-dimensional explicit sphere function-based gas-kinetic flux solver for simulation of inviscid compressible flows," J. Comput. Phys. 295, 322 (2015) and Yang et al., "Development of discrete gas kinetic scheme for simulation of 3D viscous incompressible and compressible flows," J. Comput. Phys. 319, 129 (2016)], which is often applied for the simulation of compressible flows, is simplified to improve the computational efficiency for the simulation of incompressible flows. In the original SGKS, the integral domain along the spherical surface for computing conservative variables and numerical fluxes is usually not symmetric at the cell interface. This leads the expression of numerical fluxes at the cell interface to be relatively complicated. For incompressible flows, the sphere at the cell interface can be approximately considered to be symmetric as shown in this work. Besides that, the energy equation is usually not needed for the simulation of incompressible isothermal flows. With all these simplifications, the simple and explicit formulations for the conservative variables and numerical fluxes at the cell interface can be obtained. Second, to effectively implement the no-slip boundary condition for fluid flow problems with complex geometry as well as moving boundary, the implicit boundary condition-enforced immersed boundary method [Wu and Shu, "Implicit velocity correction-based immersed boundary-lattice Boltzmann method and its applications," J. Comput. Phys. 228, 1963 (2009)] is introduced into the simplified SGKS. That is, the flow field is solved by the simplified SGKS without considering the presence of an immersed body and the no-slip boundary condition is implemented by the immersed boundary method. The accuracy and efficiency of the present scheme are validated by simulating the decaying vortex flow, flow past a stationary and rotating sphere, flow past a stationary torus, and flows over dragonfly flight.

Giant Ferroelectric Polarization in Ultrathin Ferroelectrics via Boundary-Condition Engineering.

PubMed

Xie, Lin; Li, Linze; Heikes, Colin A; Zhang, Yi; Hong, Zijian; Gao, Peng; Nelson, Christopher T; Xue, Fei; Kioupakis, Emmanouil; Chen, Longqing; Schlom, Darrel G; Wang, Peng; Pan, Xiaoqing

2017-08-01

Tailoring and enhancing the functional properties of materials at reduced dimension is critical for continuous advancement of modern electronic devices. Here, the discovery of local surface induced giant spontaneous polarization in ultrathin BiFeO 3 ferroelectric films is reported. Using aberration-corrected scanning transmission electron microscopy, it is found that the spontaneous polarization in a 2 nm-thick ultrathin BiFeO 3 film is abnormally increased up to ≈90-100 µC cm -2 in the out-of-plane direction and a peculiar rumpled nanodomain structure with very large variation in c/a ratios, which is analogous to morphotropic phase boundaries (MPBs), is formed. By a combination of density functional theory and phase-field calculations, it is shown that it is the unique single atomic Bi 2 O 3 - x layer at the surface that leads to the enhanced polarization and appearance of the MPB-like nanodomain structure. This finding clearly demonstrates a novel route to the enhanced functional properties in the material system with reduced dimension via engineering the surface boundary conditions. © 2017 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Impurity distribution and microstructure of Ga-doped ZnO films grown by molecular beam epitaxy

NASA Astrophysics Data System (ADS)

Kvit, A. V.; Yankovich, A. B.; Avrutin, V.; Liu, H.; Izyumskaya, N.; Özgür, Ü.; Morkoç, H.; Voyles, P. M.

2012-12-01

We report microstructural characterization of heavily Ga-doped ZnO (GZO) thin films on GaN and sapphire by aberration-corrected scanning transmission electron microscopy. Growth under oxygen-rich and metal-rich growth conditions leads to changes in the GZO polarity and different extended defects. For GZO layers on sapphire, the primary extended defects are voids, inversion domain boundaries, and low-angle grain boundaries. Ga doping of ZnO grown under metal-rich conditions causes a switch from pure oxygen polarity to mixed oxygen and zinc polarity in small domains. Electron energy loss spectroscopy and energy dispersive spectroscopy spectrum imaging show that Ga is homogeneous, but other residual impurities tend to accumulate at the GZO surface and at extended defects. GZO grown on GaN on c-plane sapphire has Zn polarity and no voids. There are misfit dislocations at the interfaces between GZO and an undoped ZnO buffer layer and at the buffer/GaN interface. Low-angle grain boundaries are the only threading microstructural defects. The potential effects of different extended defects and impurity distributions on free carrier scattering are discussed.

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

NASA Astrophysics Data System (ADS)

Marques, Alexandre Noll; Nave, Jean-Christophe; Rosales, Rodolfo Ruben

2017-04-01

We present a fast and accurate algorithm to solve Poisson problems in complex geometries, using regular Cartesian grids. We consider a variety of configurations, including Poisson problems with interfaces across which the solution is discontinuous (of the type arising in multi-fluid flows). The algorithm is based on a combination of the Correction Function Method (CFM) and Boundary Integral Methods (BIM). Interface and boundary conditions can be treated in a fast and accurate manner using boundary integral equations, and the associated BIM. Unfortunately, BIM can be costly when the solution is needed everywhere in a grid, e.g. fluid flow problems. We use the CFM to circumvent this issue. The solution from the BIM is used to rewrite the problem as a series of Poisson problems in rectangular domains-which requires the BIM solution at interfaces/boundaries only. These Poisson problems involve discontinuities at interfaces, of the type that the CFM can handle. Hence we use the CFM to solve them (to high order of accuracy) with finite differences and a Fast Fourier Transform based fast Poisson solver. We present 2-D examples of the algorithm applied to Poisson problems involving complex geometries, including cases in which the solution is discontinuous. We show that the algorithm produces solutions that converge with either 3rd or 4th order of accuracy, depending on the type of boundary condition and solution discontinuity.

Rotational Diffusion Depends on Box Size in Molecular Dynamics Simulations.

PubMed

Linke, Max; Köfinger, Jürgen; Hummer, Gerhard

2018-06-07

We show that the rotational dynamics of proteins and nucleic acids determined from molecular dynamics simulations under periodic boundary conditions suffer from significant finite-size effects. We remove the box-size dependence of the rotational diffusion coefficients by adding a hydrodynamic correction k B T/6 ηV with k B Boltzmann's constant, T the absolute temperature, η the solvent shear viscosity, and V the box volume. We show that this correction accounts for the finite-size dependence of the rotational diffusion coefficients of horse-heart myoglobin and a B-DNA dodecamer in aqueous solution. The resulting hydrodynamic radii are in excellent agreement with experiment.

Finite-size effects in simulations of electrolyte solutions under periodic boundary conditions

NASA Astrophysics Data System (ADS)

Thompson, Jeffrey; Sanchez, Isaac

The equilibrium properties of charged systems with periodic boundary conditions may exhibit pronounced system-size dependence due to the long range of the Coulomb force. As shown by others, the leading-order finite-size correction to the Coulomb energy of a charged fluid confined to a periodic box of volume V may be derived from sum rules satisfied by the charge-charge correlations in the thermodynamic limit V -> ∞ . In classical systems, the relevant sum rule is the Stillinger-Lovett second-moment (or perfect screening) condition. This constraint implies that for large V, periodicity induces a negative bias of -kB T(2 V) - 1 in the total Coulomb energy density of a homogeneous classical charged fluid of given density and temperature. We present a careful study of the impact of such finite-size effects on the calculation of solute chemical potentials from explicit-solvent molecular simulations of aqueous electrolyte solutions. National Science Foundation Graduate Research Fellowship Program, Grant No. DGE-1610403.

Proteus: a direct forcing method in the simulations of particulate flows

NASA Astrophysics Data System (ADS)

Feng, Zhi-Gang; Michaelides, Efstathios E.

2005-01-01

A new and efficient direct numerical method for the simulation of particulate flows is introduced. The method combines desired elements of the immersed boundary method, the direct forcing method and the lattice Boltzmann method. Adding a forcing term in the momentum equation enforces the no-slip condition on the boundary of a moving particle. By applying the direct forcing scheme, Proteus eliminates the need for the determination of free parameters, such as the stiffness coefficient in the penalty scheme or the two relaxation parameters in the adaptive-forcing scheme. The method presents a significant improvement over the previously introduced immersed-boundary-lattice-Boltzmann method (IB-LBM) where the forcing term was computed using a penalty method and a user-defined parameter. The method allows the enforcement of the rigid body motion of a particle in a more efficient way. Compared to the "bounce-back" scheme used in the conventional LBM, the direct-forcing method provides a smoother computational boundary for particles and is capable of achieving results at higher Reynolds number flows. By using a set of Lagrangian points to track the boundary of a particle, Proteus eliminates any need for the determination of the boundary nodes that are prescribed by the "bounce-back" scheme at every time step. It also makes computations for particles of irregular shapes simpler and more efficient. Proteus has been developed in two- as well as three-dimensions. This new method has been validated by comparing its results with those from experimental measurements for a single sphere settling in an enclosure under gravity. As a demonstration of the efficiency and capabilities of the present method, the settling of a large number (1232) of spherical particles is simulated in a narrow box under two different boundary conditions. It is found that when the no-slip boundary condition is imposed at the front and rear sides of the box the particles motion is significantly hindered. Under the periodic boundary conditions, the particles move faster. The simulations show that the sedimentation characteristics in a box with periodic boundary conditions at the two sides are very close to those found in the sedimentation of two-dimensional circular particles. In the Greek mythology Proteus is a hero, the son of Poseidon. In addition to his ability to change shapes and take different forms at will, Zeus granted him the power to make correct predictions for the future. One cannot expect better attributes from a numerical code.

A Boundary Condition Relaxation Algorithm for Strongly Coupled, Ablating Flows Including Shape Change

NASA Technical Reports Server (NTRS)

Gnoffo, Peter A.; Johnston, Christopher O.

2011-01-01

Implementations of a model for equilibrium, steady-state ablation boundary conditions are tested for the purpose of providing strong coupling with a hypersonic flow solver. The objective is to remove correction factors or film cooling approximations that are usually applied in coupled implementations of the flow solver and the ablation response. Three test cases are considered - the IRV-2, the Galileo probe, and a notional slender, blunted cone launched at 10 km/s from the Earth's surface. A successive substitution is employed and the order of succession is varied as a function of surface temperature to obtain converged solutions. The implementation is tested on a specified trajectory for the IRV-2 to compute shape change under the approximation of steady-state ablation. Issues associated with stability of the shape change algorithm caused by explicit time step limits are also discussed.

Thermodynamic properties of small aggregates of rare-gas atoms

NASA Technical Reports Server (NTRS)

Etters, R. D.; Kaelberer, J.

1975-01-01

The present work reports on the equilibrium thermodynamic properties of small clusters of xenon, krypton, and argon atoms, determined from a biased random-walk Monte Carlo procedure. Cluster sizes ranged from 3 to 13 atoms. Each cluster was found to have an abrupt liquid-gas phase transition at a temperature much less than for the bulk material. An abrupt solid-liquid transition is observed for thirteen- and eleven-particle clusters. For cluster sizes smaller than 11, a gradual transition from solid to liquid occurred over a fairly broad range of temperatures. Distribution of number of bond lengths as a function of bond length was calculated for several systems at various temperatures. The effects of box boundary conditions are discussed. Results show the importance of a correct description of boundary conditions. A surprising result is the slow rate at which system properties approach bulk behavior as cluster size is increased.

On the assumption of vanishing temperature fluctuations at the wall for heat transfer modeling

NASA Technical Reports Server (NTRS)

Sommer, T. P.; So, R. M. C.; Zhang, H. S.

1993-01-01

Boundary conditions for fluctuating wall temperature are required for near-wall heat transfer modeling. However, their correct specifications for arbitrary thermal boundary conditions are not clear. The conventional approach is to assume zero fluctuating wall temperature or zero gradient for the temperature variance at the wall. These are idealized specifications and the latter condition could lead to an ill posed problem for fully-developed pipe and channel flows. In this paper, the validity and extent of the zero fluctuating wall temperature condition for heat transfer calculations is examined. The approach taken is to assume a Taylor expansion in the wall normal coordinate for the fluctuating temperature that is general enough to account for both zero and non-zero value at the wall. Turbulent conductivity is calculated from the temperature variance and its dissipation rate. Heat transfer calculations assuming both zero and non-zero fluctuating wall temperature reveal that the zero fluctuating wall temperature assumption is in general valid. The effects of non-zero fluctuating wall temperature are limited only to a very small region near the wall.

Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: An accurate correction scheme for electrostatic finite-size effects

PubMed Central

Rocklin, Gabriel J.; Mobley, David L.; Dill, Ken A.; Hünenberger, Philippe H.

2013-01-01

The calculation of a protein-ligand binding free energy based on molecular dynamics (MD) simulations generally relies on a thermodynamic cycle in which the ligand is alchemically inserted into the system, both in the solvated protein and free in solution. The corresponding ligand-insertion free energies are typically calculated in nanoscale computational boxes simulated under periodic boundary conditions and considering electrostatic interactions defined by a periodic lattice-sum. This is distinct from the ideal bulk situation of a system of macroscopic size simulated under non-periodic boundary conditions with Coulombic electrostatic interactions. This discrepancy results in finite-size effects, which affect primarily the charging component of the insertion free energy, are dependent on the box size, and can be large when the ligand bears a net charge, especially if the protein is charged as well. This article investigates finite-size effects on calculated charging free energies using as a test case the binding of the ligand 2-amino-5-methylthiazole (net charge +1 e) to a mutant form of yeast cytochrome c peroxidase in water. Considering different charge isoforms of the protein (net charges −5, 0, +3, or +9 e), either in the absence or the presence of neutralizing counter-ions, and sizes of the cubic computational box (edges ranging from 7.42 to 11.02 nm), the potentially large magnitude of finite-size effects on the raw charging free energies (up to 17.1 kJ mol−1) is demonstrated. Two correction schemes are then proposed to eliminate these effects, a numerical and an analytical one. Both schemes are based on a continuum-electrostatics analysis and require performing Poisson-Boltzmann (PB) calculations on the protein-ligand system. While the numerical scheme requires PB calculations under both non-periodic and periodic boundary conditions, the latter at the box size considered in the MD simulations, the analytical scheme only requires three non-periodic PB calculations for a given system, its dependence on the box size being analytical. The latter scheme also provides insight into the physical origin of the finite-size effects. These two schemes also encompass a correction for discrete solvent effects that persists even in the limit of infinite box sizes. Application of either scheme essentially eliminates the size dependence of the corrected charging free energies (maximal deviation of 1.5 kJ mol−1). Because it is simple to apply, the analytical correction scheme offers a general solution to the problem of finite-size effects in free-energy calculations involving charged solutes, as encountered in calculations concerning, e.g., protein-ligand binding, biomolecular association, residue mutation, pKa and redox potential estimation, substrate transformation, solvation, and solvent-solvent partitioning. PMID:24320250

Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: an accurate correction scheme for electrostatic finite-size effects.

PubMed

Rocklin, Gabriel J; Mobley, David L; Dill, Ken A; Hünenberger, Philippe H

2013-11-14

The calculation of a protein-ligand binding free energy based on molecular dynamics (MD) simulations generally relies on a thermodynamic cycle in which the ligand is alchemically inserted into the system, both in the solvated protein and free in solution. The corresponding ligand-insertion free energies are typically calculated in nanoscale computational boxes simulated under periodic boundary conditions and considering electrostatic interactions defined by a periodic lattice-sum. This is distinct from the ideal bulk situation of a system of macroscopic size simulated under non-periodic boundary conditions with Coulombic electrostatic interactions. This discrepancy results in finite-size effects, which affect primarily the charging component of the insertion free energy, are dependent on the box size, and can be large when the ligand bears a net charge, especially if the protein is charged as well. This article investigates finite-size effects on calculated charging free energies using as a test case the binding of the ligand 2-amino-5-methylthiazole (net charge +1 e) to a mutant form of yeast cytochrome c peroxidase in water. Considering different charge isoforms of the protein (net charges -5, 0, +3, or +9 e), either in the absence or the presence of neutralizing counter-ions, and sizes of the cubic computational box (edges ranging from 7.42 to 11.02 nm), the potentially large magnitude of finite-size effects on the raw charging free energies (up to 17.1 kJ mol(-1)) is demonstrated. Two correction schemes are then proposed to eliminate these effects, a numerical and an analytical one. Both schemes are based on a continuum-electrostatics analysis and require performing Poisson-Boltzmann (PB) calculations on the protein-ligand system. While the numerical scheme requires PB calculations under both non-periodic and periodic boundary conditions, the latter at the box size considered in the MD simulations, the analytical scheme only requires three non-periodic PB calculations for a given system, its dependence on the box size being analytical. The latter scheme also provides insight into the physical origin of the finite-size effects. These two schemes also encompass a correction for discrete solvent effects that persists even in the limit of infinite box sizes. Application of either scheme essentially eliminates the size dependence of the corrected charging free energies (maximal deviation of 1.5 kJ mol(-1)). Because it is simple to apply, the analytical correction scheme offers a general solution to the problem of finite-size effects in free-energy calculations involving charged solutes, as encountered in calculations concerning, e.g., protein-ligand binding, biomolecular association, residue mutation, pKa and redox potential estimation, substrate transformation, solvation, and solvent-solvent partitioning.

Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: An accurate correction scheme for electrostatic finite-size effects

NASA Astrophysics Data System (ADS)

Rocklin, Gabriel J.; Mobley, David L.; Dill, Ken A.; Hünenberger, Philippe H.

2013-11-01

The calculation of a protein-ligand binding free energy based on molecular dynamics (MD) simulations generally relies on a thermodynamic cycle in which the ligand is alchemically inserted into the system, both in the solvated protein and free in solution. The corresponding ligand-insertion free energies are typically calculated in nanoscale computational boxes simulated under periodic boundary conditions and considering electrostatic interactions defined by a periodic lattice-sum. This is distinct from the ideal bulk situation of a system of macroscopic size simulated under non-periodic boundary conditions with Coulombic electrostatic interactions. This discrepancy results in finite-size effects, which affect primarily the charging component of the insertion free energy, are dependent on the box size, and can be large when the ligand bears a net charge, especially if the protein is charged as well. This article investigates finite-size effects on calculated charging free energies using as a test case the binding of the ligand 2-amino-5-methylthiazole (net charge +1 e) to a mutant form of yeast cytochrome c peroxidase in water. Considering different charge isoforms of the protein (net charges -5, 0, +3, or +9 e), either in the absence or the presence of neutralizing counter-ions, and sizes of the cubic computational box (edges ranging from 7.42 to 11.02 nm), the potentially large magnitude of finite-size effects on the raw charging free energies (up to 17.1 kJ mol-1) is demonstrated. Two correction schemes are then proposed to eliminate these effects, a numerical and an analytical one. Both schemes are based on a continuum-electrostatics analysis and require performing Poisson-Boltzmann (PB) calculations on the protein-ligand system. While the numerical scheme requires PB calculations under both non-periodic and periodic boundary conditions, the latter at the box size considered in the MD simulations, the analytical scheme only requires three non-periodic PB calculations for a given system, its dependence on the box size being analytical. The latter scheme also provides insight into the physical origin of the finite-size effects. These two schemes also encompass a correction for discrete solvent effects that persists even in the limit of infinite box sizes. Application of either scheme essentially eliminates the size dependence of the corrected charging free energies (maximal deviation of 1.5 kJ mol-1). Because it is simple to apply, the analytical correction scheme offers a general solution to the problem of finite-size effects in free-energy calculations involving charged solutes, as encountered in calculations concerning, e.g., protein-ligand binding, biomolecular association, residue mutation, pKa and redox potential estimation, substrate transformation, solvation, and solvent-solvent partitioning.

Thermocapillary Bubble Migration: Thermal Boundary Layers for Large Marangoni Numbers

NASA Technical Reports Server (NTRS)

Balasubramaniam, R.; Subramanian, R. S.

1996-01-01

The migration of an isolated gas bubble in an immiscible liquid possessing a temperature gradient is analyzed in the absence of gravity. The driving force for the bubble motion is the shear stress at the interface which is a consequence of the temperature dependence of the surface tension. The analysis is performed under conditions for which the Marangoni number is large, i.e. energy is transferred predominantly by convection. Velocity fields in the limit of both small and large Reynolds numbers are used. The thermal problem is treated by standard boundary layer theory. The outer temperature field is obtained in the vicinity of the bubble. A similarity solution is obtained for the inner temperature field. For both small and large Reynolds numbers, the asymptotic values of the scaled migration velocity of the bubble in the limit of large Marangoni numbers are calculated. The results show that the migration velocity has the same scaling for both low and large Reynolds numbers, but with a different coefficient. Higher order thermal boundary layers are analyzed for the large Reynolds number flow field and the higher order corrections to the migration velocity are obtained. Results are also presented for the momentum boundary layer and the thermal wake behind the bubble, for large Reynolds number conditions.

Wall touching kink mode calculations with the M3D code

NASA Astrophysics Data System (ADS)

Breslau, J. A.

2014-10-01

In recent years there have been a number of results published concerning the transient vessel currents and forces occurring during a tokamak VDE, as predicted by simulations with the nonlinear MHD code M3D. The nature of the simulations is such that these currents and forces occur at the boundary of the computational domain, making the proper choice of boundary conditions critical to the reliability of the results. The M3D boundary condition includes the prescription that the normal component of the velocity vanish at the wall. It has been argued that this prescription invalidates the calculations because it would seem to rule out the possibility of advection of plasma surface currents into the wall. This claim has been tested by applying M3D to an idealized case - a kink-unstable plasma column - in order to abstract the essential physics from the complications involved in the attempt to model real devices. While comparison of the results is complicated by effects arising from the higher dimensionality and complexity of M3D, we have verified that M3D is capable of reproducing both the correct saturation behavior of the free boundary kink and the ``Hiro'' currents arising when the kink interacts with a conducting tile surface interior to the ideal wall.

A correction function method for the wave equation with interface jump conditions

NASA Astrophysics Data System (ADS)

Abraham, David S.; Marques, Alexandre Noll; Nave, Jean-Christophe

2018-01-01

In this paper a novel method to solve the constant coefficient wave equation, subject to interface jump conditions, is presented. In general, such problems pose issues for standard finite difference solvers, as the inherent discontinuity in the solution results in erroneous derivative information wherever the stencils straddle the given interface. Here, however, the recently proposed Correction Function Method (CFM) is used, in which correction terms are computed from the interface conditions, and added to affected nodes to compensate for the discontinuity. In contrast to existing methods, these corrections are not simply defined at affected nodes, but rather generalized to a continuous function within a small region surrounding the interface. As a result, the correction function may be defined in terms of its own governing partial differential equation (PDE) which may be solved, in principle, to arbitrary order of accuracy. The resulting scheme is not only arbitrarily high order, but also robust, having already seen application to Poisson problems and the heat equation. By extending the CFM to this new class of PDEs, the treatment of wave interface discontinuities in homogeneous media becomes possible. This allows, for example, for the straightforward treatment of infinitesimal source terms and sharp boundaries, free of staircasing errors. Additionally, new modifications to the CFM are derived, allowing compatibility with explicit multi-step methods, such as Runge-Kutta (RK4), without a reduction in accuracy. These results are then verified through numerous numerical experiments in one and two spatial dimensions.

On the Flow of a Compressible Fluid by the Hodograph Method. II - Fundamental Set of Particular Flow Solutions of the Chaplygin Differential Equation

NASA Technical Reports Server (NTRS)

Garrick, I. E.; Kaplan, Carl

1944-01-01

The differential equation of Chaplygin's jet problem is utilized to give a systematic development of particular solutions of the hodograph flow equations, which extends the treatment of Chaplygin into the supersonic range and completes the set of particular solutions. The particular solutions serve to place on a reasonable basis the use of velocity correction formulas for the comparison of incompressible and compressible flows. It is shown that the geometric-mean type of velocity correction formula introduced in part I has significance as an over-all type of approximation in the subsonic range. A brief review of general conditions limiting the potential flow of an adiabatic compressible fluid is given and application is made to the particular solutions, yielding conditions for the existence of singular loci in the supersonic range. The combining of particular solutions in accordance with prescribed boundary flow conditions is not treated in the present paper.

Application of the Shiono and Knight Method in asymmetric compound channels with different side slopes of the internal wall

NASA Astrophysics Data System (ADS)

Alawadi, Wisam; Al-Rekabi, Wisam S.; Al-Aboodi, Ali H.

2018-03-01

The Shiono and Knight Method (SKM) is widely used to predict the lateral distribution of depth-averaged velocity and boundary shear stress for flows in compound channels. Three calibrating coefficients need to be estimated for applying the SKM, namely eddy viscosity coefficient ( λ), friction factor ( f) and secondary flow coefficient ( k). There are several tested methods which can satisfactorily be used to estimate λ, f. However, the calibration of secondary flow coefficients k to account for secondary flow effects correctly is still problematic. In this paper, the calibration of secondary flow coefficients is established by employing two approaches to estimate correct values of k for simulating asymmetric compound channel with different side slopes of the internal wall. The first approach is based on Abril and Knight (2004) who suggest fixed values for main channel and floodplain regions. In the second approach, the equations developed by Devi and Khatua (2017) that relate the variation of the secondary flow coefficients with the relative depth ( β) and width ratio ( α) are used. The results indicate that the calibration method developed by Devi and Khatua (2017) is a better choice for calibrating the secondary flow coefficients than using the first approach which assumes a fixed value of k for different flow depths. The results also indicate that the boundary condition based on the shear force continuity can successfully be used for simulating rectangular compound channels, while the continuity of depth-averaged velocity and its gradient is accepted boundary condition in simulations of trapezoidal compound channels. However, the SKM performance for predicting the boundary shear stress over the shear layer region may not be improved by only imposing the suitable calibrated values of secondary flow coefficients. This is because difficulties of modelling the complex interaction that develops between the flows in the main channel and on the floodplain in this region.

Boundary condition determined wave functions for the ground states of one- and two-electron homonuclear molecules

NASA Astrophysics Data System (ADS)

Patil, S. H.; Tang, K. T.; Toennies, J. P.

1999-10-01

Simple analytical wave functions satisfying appropriate boundary conditions are constructed for the ground states of one-and two-electron homonuclear molecules. Both the asymptotic condition when one electron is far away and the cusp condition when the electron coalesces with a nucleus are satisfied by the proposed wave function. For H2+, the resulting wave function is almost identical to the Guillemin-Zener wave function which is known to give very good energies. For the two electron systems H2 and He2++, the additional electron-electron cusp condition is rigorously accounted for by a simple analytic correlation function which has the correct behavior not only for r12→0 and r12→∞ but also for R→0 and R→∞, where r12 is the interelectronic distance and R, the internuclear distance. Energies obtained from these simple wave functions agree within 2×10-3 a.u. with the results of the most sophisticated variational calculations for all R and for all systems studied. This demonstrates that rather simple physical considerations can be used to derive very accurate wave functions for simple molecules thereby avoiding laborious numerical variational calculations.

Two-Dimensional Thermal Boundary Layer Corrections for Convective Heat Flux Gauges

NASA Technical Reports Server (NTRS)

Kandula, Max; Haddad, George

2007-01-01

This work presents a CFD (Computational Fluid Dynamics) study of two-dimensional thermal boundary layer correction factors for convective heat flux gauges mounted in flat plate subjected to a surface temperature discontinuity with variable properties taken into account. A two-equation k - omega turbulence model is considered. Results are obtained for a wide range of Mach numbers (1 to 5), gauge radius ratio, and wall temperature discontinuity. Comparisons are made for correction factors with constant properties and variable properties. It is shown that the variable-property effects on the heat flux correction factors become significant

Beta value coupled wave theory for nonslanted reflection gratings.

PubMed

Neipp, Cristian; Francés, Jorge; Gallego, Sergi; Bleda, Sergio; Martínez, Francisco Javier; Pascual, Inmaculada; Beléndez, Augusto

2014-01-01

We present a modified coupled wave theory to describe the properties of nonslanted reflection volume diffraction gratings. The method is based on the beta value coupled wave theory, which will be corrected by using appropriate boundary conditions. The use of this correction allows predicting the efficiency of the reflected order for nonslanted reflection gratings embedded in two media with different refractive indices. The results obtained by using this method will be compared to those obtained using a matrix method, which gives exact solutions in terms of Mathieu functions, and also to Kogelnik's coupled wave theory. As will be demonstrated, the technique presented in this paper means a significant improvement over Kogelnik's coupled wave theory.

Beta Value Coupled Wave Theory for Nonslanted Reflection Gratings

PubMed Central

Neipp, Cristian; Francés, Jorge; Gallego, Sergi; Bleda, Sergio; Martínez, Francisco Javier; Pascual, Inmaculada; Beléndez, Augusto

2014-01-01

We present a modified coupled wave theory to describe the properties of nonslanted reflection volume diffraction gratings. The method is based on the beta value coupled wave theory, which will be corrected by using appropriate boundary conditions. The use of this correction allows predicting the efficiency of the reflected order for nonslanted reflection gratings embedded in two media with different refractive indices. The results obtained by using this method will be compared to those obtained using a matrix method, which gives exact solutions in terms of Mathieu functions, and also to Kogelnik's coupled wave theory. As will be demonstrated, the technique presented in this paper means a significant improvement over Kogelnik's coupled wave theory. PMID:24723811

Numerical simulation and comparison of symmetrical/supercritical airfoils for the near tip region of a helicopter in forward flight

NASA Technical Reports Server (NTRS)

Badavi, F. F.

1989-01-01

Aerodynamic loads on a multi-bladed helicopter rotor in forward flight at transonic tip conditions are calculated. The unsteady, three-dimensional, time-accurate compressible Reynolds-averaged thin layer Navier-Stokes equations are solved in a rotating coordinate system on a body-conformed, curvilinear grid of C-H topology. Detailed boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 supercritical airfoils are made under identical forward flight conditions. The rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blade. This correction is obtained by computing the local induced downwash velocity with a free wake analysis program. The calculations are performed on the Numerical Aerodynamic Simulation Cray 2 and the VPS32 (a derivative of a Cyber 205 at the Langley Research Center) for a model helicopter rotor in forward flight.

Wind-tunnel interference with particular reference to off-center positions of the wing and to the downwash at the tail

NASA Technical Reports Server (NTRS)

Silverstein, Abe; White, James A

1937-01-01

The theory of wind tunnel boundary influence on the downwash from a wing has been extended to provide more complete corrections for application to airplane test data. The first section of the report gives the corrections of the lifting line for wing positions above or below the tunnel center line; the second section shows the manner in which the induced boundary influence changes with distance aft of the lifting line. Values of the boundary corrections are given for off-center positions of the wing in circular, square, 2:1 rectangular, and 2:1 elliptical tunnels. Aft of the wing the corrections are presented for only the square and the 2:1 rectangular tunnels, but it is believed that these may be applied to jets of circular and 2:1 elliptical cross sections. In all cases results are included for both open and closed tunnels.

Entrainment of bed material by Earth-surface mass flows: review and reformulation of depth-integrated theory

USGS Publications Warehouse

Iverson, Richard M.; Chaojun Ouyang,

2015-01-01

Earth-surface mass flows such as debris flows, rock avalanches, and dam-break floods can grow greatly in size and destructive potential by entraining bed material they encounter. Increasing use of depth-integrated mass- and momentum-conservation equations to model these erosive flows motivates a review of the underlying theory. Our review indicates that many existing models apply depth-integrated conservation principles incorrectly, leading to spurious inferences about the role of mass and momentum exchanges at flow-bed boundaries. Model discrepancies can be rectified by analyzing conservation of mass and momentum in a two-layer system consisting of a moving upper layer and static lower layer. Our analysis shows that erosion or deposition rates at the interface between layers must in general satisfy three jump conditions. These conditions impose constraints on valid erosion formulas, and they help determine the correct forms of depth-integrated conservation equations. Two of the three jump conditions are closely analogous to Rankine-Hugoniot conditions that describe the behavior of shocks in compressible gasses, and the third jump condition describes shear traction discontinuities that necessarily exist across eroding boundaries. Grain-fluid mixtures commonly behave as compressible materials as they undergo entrainment, because changes in bulk density occur as the mixtures mobilize and merge with an overriding flow. If no bulk density change occurs, then only the shear-traction jump condition applies. Even for this special case, however, accurate formulation of depth-integrated momentum equations requires a clear distinction between boundary shear tractions that exist in the presence or absence of bed erosion.

Phonological Codes Are Assembled before Word Fixation: Evidence from Boundary Paradigm in Sentence Reading

ERIC Educational Resources Information Center

Miellet, Sebastien; Sparrow, Laurent

2004-01-01

This experiment employed the boundary paradigm during sentence reading to explore the nature of early phonological coding in reading. Fixation durations were shorter when the parafoveal preview was the correct word than when it was a spelling control pseudoword. In contrast, there was no significant difference between correct word and…

Three-dimensional Monte Carlo model of pulsed-laser treatment of cutaneous vascular lesions

NASA Astrophysics Data System (ADS)

Milanič, Matija; Majaron, Boris

2011-12-01

We present a three-dimensional Monte Carlo model of optical transport in skin with a novel approach to treatment of side boundaries of the volume of interest. This represents an effective way to overcome the inherent limitations of ``escape'' and ``mirror'' boundary conditions and enables high-resolution modeling of skin inclusions with complex geometries and arbitrary irradiation patterns. The optical model correctly reproduces measured values of diffuse reflectance for normal skin. When coupled with a sophisticated model of thermal transport and tissue coagulation kinetics, it also reproduces realistic values of radiant exposure thresholds for epidermal injury and for photocoagulation of port wine stain blood vessels in various skin phototypes, with or without application of cryogen spray cooling.

Finite-size effects for anisotropic 2D Ising model with various boundary conditions

NASA Astrophysics Data System (ADS)

Izmailian, N. Sh

2012-12-01

We analyze the exact partition function of the anisotropic Ising model on finite M × N rectangular lattices under four different boundary conditions (periodic-periodic (pp), periodic-antiperiodic (pa), antiperiodic-periodic (ap) and antiperiodic-antiperiodic (aa)) obtained by Kaufman (1949 Phys. Rev. 76 1232), Wu and Hu (2002 J. Phys. A: Math. Gen. 35 5189) and Kastening (2002 Phys. Rev. E 66 057103)). We express the partition functions in terms of the partition functions Zα, β(J, k) with (α, β) = (0, 0), (1/2, 0), (0, 1/2) and (1/2, 1/2), J is an interaction coupling and k is an anisotropy parameter. Based on such expressions, we then extend the algorithm of Ivashkevich et al (2002 J. Phys. A: Math. Gen. 35 5543) to derive the exact asymptotic expansion of the logarithm of the partition function for all boundary conditions mentioned above. Our result is f = fbulk + ∑∞p = 0fp(ρ, k)S-p - 1, where f is the free energy of the system, fbulk is the free energy of the bulk, S = MN is the area of the lattice and ρ = M/N is the aspect ratio. All coefficients in this expansion are expressed through analytical functions. We have introduced the effective aspect ratio ρeff = ρ/sinh 2Jc and show that for pp and aa boundary conditions all finite size correction terms are invariant under the transformation ρeff → 1/ρeff. This article is part of ‘Lattice models and integrability’, a special issue of Journal of Physics A: Mathematical and Theoretical in honour of F Y Wu's 80th birthday.

The Boundary Layers in Fluids with Little Friction

NASA Technical Reports Server (NTRS)

Blasius, H.

1950-01-01

The vortices forming in flowing water behind solid bodies are not represented correctly by the solution of the potential theory nor by Helmholtz's jets. Potential theory is unable to satisfy the condition that the water adheres at the wetted bodies, and its solutions of the fundamental hydrodynamic equations are at variance with the observation that the flow separates from the body at a certain point and sends forth a highly turbulent boundary layer into the free flow. Helmholtz's theory attempts to imitate the latter effect in such a way that it joins two potential flows, jet and still water, nonanalytical along a stream curve. The admissibility of this method is based on the fact that, at zero pressure, which is to prevail at the cited stream curve, the connection of the fluid, and with it the effect of adjacent parts on each other, is canceled. In reality, however, the pressure at these boundaries is definitely not zero, but can even be varied arbitrarily. Besides, Helmholtz's theory with its potential flows does not satisfy the condition of adherence nor explain the origin of the vortices, for in all of these problems, the friction must be taken into account on principle, according to the vortex theorem.

New and Topologically Massive Gravity, from the Outside In

NASA Astrophysics Data System (ADS)

Cunliff, Colin

This thesis examines the asymptotically anti-de Sitter solutions of higher-derivative gravity in 2+1 dimensions, using a Fefferman-Graham-like approach that expands solutions from the boundary (at infinity) into the interior. First, solutions of topologically massive gravity (TMG) are analyzed for values of the mass parameter in the range mu ≥ 1. The traditional Fefferman-Graham expansion fails to capture the dynamics of TMG, and new terms in the asymptotic expansion are needed to include the massive graviton modes. The linearized modes of Carlip, Deser, Waldron and Wise map onto the non-Einstein solutions for all μ, with nonlinear corrections appearing at higher order in the expansion. A similar result is found for new massive gravity (NMG), where the asymptotic behavior of massive gravitons is found to depend on the coupling parameter m2. Additionally, new boundary conditions are discovered for a range of values -1 < 2m2 l2 < 1 at which non-Einstein modes decay more slowly than the rate required for Brown-Henneaux boundary conditions. The holographically renormalized stress tensor is computed for these modes, and the relevant counterterms are identified up to unphysical ambiguities.

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 1: Theory document (version 3.0)

NASA Technical Reports Server (NTRS)

Epton, Michael A.; Magnus, Alfred E.

1990-01-01

An outline of the derivation of the differential equation governing linear subsonic and supersonic potential flow is given. The use of Green's Theorem to obtain an integral equation over the boundary surface is discussed. The engineering techniques incorporated in the Panel Aerodynamics (PAN AIR) program (a discretization method which solves the integral equation for arbitrary first order boundary conditions) are then discussed in detail. Items discussed include the construction of the compressibility transformation, splining techniques, imposition of the boundary conditions, influence coefficient computation (including the concept of the finite part of an integral), computation of pressure coefficients, and computation of forces and moments. Principal revisions to version 3.0 are the following: (1) appendices H and K more fully describe the Aerodynamic Influence Coefficient (AIC) construction; (2) appendix L now provides a complete description of the AIC solution process; (3) appendix P is new and discusses the theory for the new FDP module (which calculates streamlines and offbody points); and (4) numerous small corrections and revisions reflecting the MAG module rewrite.

Automatic software correction of residual aberrations in reconstructed HRTEM exit waves of crystalline samples

DOE PAGES

Ophus, Colin; Rasool, Haider I.; Linck, Martin; ...

2016-11-30

We develop an automatic and objective method to measure and correct residual aberrations in atomic-resolution HRTEM complex exit waves for crystalline samples aligned along a low-index zone axis. Our method uses the approximate rotational point symmetry of a column of atoms or single atom to iteratively calculate a best-fit numerical phase plate for this symmetry condition, and does not require information about the sample thickness or precise structure. We apply our method to two experimental focal series reconstructions, imaging a β-Si 3N 4 wedge with O and N doping, and a single-layer graphene grain boundary. We use peak and latticemore » fitting to evaluate the precision of the corrected exit waves. We also apply our method to the exit wave of a Si wedge retrieved by off-axis electron holography. In all cases, the software correction of the residual aberration function improves the accuracy of the measured exit waves.« less

Automatic software correction of residual aberrations in reconstructed HRTEM exit waves of crystalline samples

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ophus, Colin; Rasool, Haider I.; Linck, Martin

We develop an automatic and objective method to measure and correct residual aberrations in atomic-resolution HRTEM complex exit waves for crystalline samples aligned along a low-index zone axis. Our method uses the approximate rotational point symmetry of a column of atoms or single atom to iteratively calculate a best-fit numerical phase plate for this symmetry condition, and does not require information about the sample thickness or precise structure. We apply our method to two experimental focal series reconstructions, imaging a β-Si 3N 4 wedge with O and N doping, and a single-layer graphene grain boundary. We use peak and latticemore » fitting to evaluate the precision of the corrected exit waves. We also apply our method to the exit wave of a Si wedge retrieved by off-axis electron holography. In all cases, the software correction of the residual aberration function improves the accuracy of the measured exit waves.« less

Contact Stress Analysis of Spiral Bevel Gears Using Finite Element Analysis

NASA Technical Reports Server (NTRS)

Bibel, G. D.; Kumar, A; Reddy, S.; Handschuh, R.

1995-01-01

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Voting based object boundary reconstruction

NASA Astrophysics Data System (ADS)

Tian, Qi; Zhang, Like; Ma, Jingsheng

2005-07-01

A voting-based object boundary reconstruction approach is proposed in this paper. Morphological technique was adopted in many applications for video object extraction to reconstruct the missing pixels. However, when the missing areas become large, the morphological processing cannot bring us good results. Recently, Tensor voting has attracted people"s attention, and it can be used for boundary estimation on curves or irregular trajectories. However, the complexity of saliency tensor creation limits its applications in real-time systems. An alternative approach based on tensor voting is introduced in this paper. Rather than creating saliency tensors, we use a "2-pass" method for orientation estimation. For the first pass, Sobel d*etector is applied on a coarse boundary image to get the gradient map. In the second pass, each pixel puts decreasing weights based on its gradient information, and the direction with maximum weights sum is selected as the correct orientation of the pixel. After the orientation map is obtained, pixels begin linking edges or intersections along their direction. The approach is applied to various video surveillance clips under different conditions, and the experimental results demonstrate significant improvement on the final extracted objects accuracy.

Dilatation-dissipation corrections for advanced turbulence models

NASA Technical Reports Server (NTRS)

Wilcox, David C.

1992-01-01

This paper analyzes dilatation-dissipation based compressibility corrections for advanced turbulence models. Numerical computations verify that the dilatation-dissipation corrections devised by Sarkar and Zeman greatly improve both the k-omega and k-epsilon model predicted effect of Mach number on spreading rate. However, computations with the k-gamma model also show that the Sarkar/Zeman terms cause an undesired reduction in skin friction for the compressible flat-plate boundary layer. A perturbation solution for the compressible wall layer shows that the Sarkar and Zeman terms reduce the effective von Karman constant in the law of the wall. This is the source of the inaccurate k-gamma model skin-friction predictions for the flat-plate boundary layer. The perturbation solution also shows that the k-epsilon model has an inherent flaw for compressible boundary layers that is not compensated for by the dilatation-dissipation corrections. A compressibility modification for k-gamma and k-epsilon models is proposed that is similar to those of Sarkar and Zeman. The new compressibility term permits accurate predictions for the compressible mixing layer, flat-plate boundary layer, and a shock separated flow with the same values for all closure coefficients.

Ship Trim Optimization: Assessment of Influence of Trim on Resistance of MOERI Container Ship

PubMed Central

Duan, Wenyang

2014-01-01

Environmental issues and rising fuel prices necessitate better energy efficiency in all sectors. Shipping industry is a stakeholder in environmental issues. Shipping industry is responsible for approximately 3% of global CO2 emissions, 14-15% of global NOX emissions, and 16% of global SOX emissions. Ship trim optimization has gained enormous momentum in recent years being an effective operational measure for better energy efficiency to reduce emissions. Ship trim optimization analysis has traditionally been done through tow-tank testing for a specific hullform. Computational techniques are increasingly popular in ship hydrodynamics applications. The purpose of this study is to present MOERI container ship (KCS) hull trim optimization by employing computational methods. KCS hull total resistances and trim and sinkage computed values, in even keel condition, are compared with experimental values and found in reasonable agreement. The agreement validates that mesh, boundary conditions, and solution techniques are correct. The same mesh, boundary conditions, and solution techniques are used to obtain resistance values in different trim conditions at Fn = 0.2274. Based on attained results, optimum trim is suggested. This research serves as foundation for employing computational techniques for ship trim optimization. PMID:24578649

Study of boundary-layer transition using transonic cone Preston tube data

NASA Technical Reports Server (NTRS)

Reed, T. D.; Abu-Mostafa, A.

1982-01-01

Laminar layer Preston tube data on a sharp nose, ten degree cone obtained in the Ames 11 ft TWT and in flight tests are analyzed. During analyses of the laminar-boundary layer data, errors were discovered in both the wind tunnel and the flight data. A correction procedure for errors in the flight data is recommended which forces the flight data to exhibit some of the orderly characteristics of the wind tunnel data. From corrected wind tunnel data, a correlation is developed between Preston tube pressures and the corresponding values of theoretical laminar skin friction. Because of the uncertainty in correcting the flight data, a correlation for the unmodified data is developed, and, in addition, three other correlations are developed based on different correction procedures. Each of these correlations are used in conjunction with the wind tunnel correlation to define effective freestream unit Reynolds numbers for the 11 ft TWT over a Mach number range of 0.30 to 0.95. The maximum effective Reynolds numbers are approximately 6.5% higher than the normal values. These maximum values occur between freestream Mach numbers of 0.60 and 0.80. Smaller values are found outside this Mach number range. These results indicate wind tunnel noise affects the average laminar skin friction much less than it affects boundary layer transition. Data on the onset, extent, and end of boundary layer transition are summarized. Application of a procedure for studying the relative effects of varying nose radius on a ten degree cone at supercritical speeds indicates that increasing nose radius promotes boundary layer transition and separation of laminar boundary layers.

Transport dissipative particle dynamics model for mesoscopic advection- diffusion-reaction problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhen, Li; Yazdani, Alireza; Tartakovsky, Alexandre M.

2015-07-07

We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. tDPD is an extension of the classic DPD framework with extra variables for describing the evolution of concentration fields. The transport of concentration is modeled by a Fickian flux and a random flux between particles, and an analytical formula is proposed to relate the mesoscopic concentration friction to the effective diffusion coefficient. To validate the present tDPD model and the boundary conditions, we perform three tDPDmore » simulations of one-dimensional diffusion with different boundary conditions, and the results show excellent agreement with the theoretical solutions. We also performed two-dimensional simulations of ADR systems and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, we present an application of the tDPD model to the dynamic process of blood coagulation involving 25 reacting species in order to demonstrate the potential of tDPD in simulating biological dynamics at the mesoscale. We find that the tDPD solution of this comprehensive 25-species coagulation model is only twice as computationally expensive as the DPD simulation of the hydrodynamics only, which is a significant advantage over available continuum solvers.« less

Greenland Regional and Ice Sheet-wide Geometry Sensitivity to Boundary and Initial conditions

NASA Astrophysics Data System (ADS)

Logan, L. C.; Narayanan, S. H. K.; Greve, R.; Heimbach, P.

2017-12-01

Ice sheet and glacier model outputs require inputs from uncertainly known initial and boundary conditions, and other parameters. Conservation and constitutive equations formalize the relationship between model inputs and outputs, and the sensitivity of model-derived quantities of interest (e.g., ice sheet volume above floatation) to model variables can be obtained via the adjoint model of an ice sheet. We show how one particular ice sheet model, SICOPOLIS (SImulation COde for POLythermal Ice Sheets), depends on these inputs through comprehensive adjoint-based sensitivity analyses. SICOPOLIS discretizes the shallow-ice and shallow-shelf approximations for ice flow, and is well-suited for paleo-studies of Greenland and Antarctica, among other computational domains. The adjoint model of SICOPOLIS was developed via algorithmic differentiation, facilitated by the source transformation tool OpenAD (developed at Argonne National Lab). While model sensitivity to various inputs can be computed by costly methods involving input perturbation simulations, the time-dependent adjoint model of SICOPOLIS delivers model sensitivities to initial and boundary conditions throughout time at lower cost. Here, we explore both the sensitivities of the Greenland Ice Sheet's entire and regional volumes to: initial ice thickness, precipitation, basal sliding, and geothermal flux over the Holocene epoch. Sensitivity studies such as described here are now accessible to the modeling community, based on the latest version of SICOPOLIS that has been adapted for OpenAD to generate correct and efficient adjoint code.

Seismic performance of arch dams on non-homogeneous and discontinuous foundations (a case study: Karun 4 Dam)

NASA Astrophysics Data System (ADS)

Ferdousi, A.

2017-06-01

The present study set out to investigate the nonlinear seismic response of the dam-reservoir-rock foundation system, taking into consideration the effects of change in the material properties of discontinuous foundation. To this end, it is important to provide the proper modeling of truncated boundary conditions at the far-end of rock foundation and reservoir fluid domain and to correctly apply the in situ stresses for rock foundation. The nonlinear seismic response of an arch dam mainly depends on the opening and sliding of the dam body's contraction joints and foundation discontinuities, failure of the jointed rock and concrete materials, etc. In this paper, a time domain dynamic analysis of the 3D dam-reservoir-foundation interaction problem was performed by developing a nonlinear Finite Element program. The results of the analysis of Karun-4 Dam revealed the essential role of modeling discontinuities and boundary conditions of rock foundation under seismic excitation.

Modelling the baroclinic circulation with tidal components in the Adriatic Sea

NASA Astrophysics Data System (ADS)

Guarnieri, A.; Pinardi, N.; Oddo, P.; Bortoluzzi, G.; Ravaioli, M.

2012-04-01

The impact of tides in the circulation of the Adriatic sea has been investigated by means of a nested baroclinic numerical ocean model. Tides have been introduced using a modified Flather boundary condition at the open side of the domain. The results show that tidal amplitudes and phases are reproduced correctly by the baroclinic model and the tidal harmonic constants errors are comparable with those resulting from the most consolidated barotropic models. Numerical experiments were conducted to estimate and assess the impact of (i) the modified Flather lateral boundary condition, (ii) the tides on temperature, salinity and stratification structures in the basin, and (iii) the tides on mixing and circulation in general. Tides induce a different momentum advective component in the basin which in turn produces a different distribution of water masses in the basin. Tides impact on mixing and stratification in the Po river region (north-western Adriatic) and induce fluctuations of salinity and temperature on semidiurnal frequencies in all seasons for the first and only winter for the second.

Crystal Growth Simulations To Establish Physically Relevant Kinetic Parameters from the Empirical Kolmogorov-Johnson-Mehl-Avrami Model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dill, Eric D.; Folmer, Jacob C.W.; Martin, James D.

A series of simulations was performed to enable interpretation of the material and physical significance of the parameters defined in the Kolmogorov, Johnson and Mehl, and Avrami (KJMA) rate expression commonly used to describe phase boundary controlled reactions of condensed matter. The parameters k, n, and t 0 are shown to be highly correlated, which if unaccounted for seriously challenge mechanistic interpretation. It is demonstrated that rate measurements exhibit an intrinsic uncertainty without precise knowledge of the location and orientation of nucleation with respect to the free volume into which it grows. More significantly, it is demonstrated that the KJMAmore » rate constant k is highly dependent on sample size. However, under the simulated conditions of slow nucleation relative to crystal growth, sample volume and sample anisotropy correction affords a means to eliminate the experimental condition dependence of the KJMA rate constant, k, producing the material-specific parameter, the velocity of the phase boundary, v pb.« less

Analysis of a high speed civil transport configuration at subsonic flow conditions using a Navier-Stokes solver

NASA Technical Reports Server (NTRS)

Lessard, Victor R.

1993-01-01

Computations of three dimensional vortical flows over a generic High Speed Civil Transport (HSCT) configuration with an aspect ratio of 3.04 are performed using a thin-layer Navier-Stokes solver. The HSCT cruise configuration is modeled without leading or trailing edge flap deflections and without engine nacelles. The flow conditions, which correspond to tests done in the NASA Langley 8-Foot Transonic Pressure Tunnel (TPT), are a subsonic Mach number of 0.3 and Reynolds number of 4.4 million for a range-of-attack (-.23 deg to 17.78 deg). The effects of the farfield boundary location with respect to the body are investigated. The boundary layer is assumed turbulent and simulated using an algebraic turbulence model. The key features of the vortices and their interactions are captured. Grid distribution in the vortex regions is critical for predicting the correct induced lift. Computed forces and surface pressures compare reasonably well with the experimental TPT data.

Pulsatile flow in ventricular catheters for hydrocephalus

NASA Astrophysics Data System (ADS)

Giménez, Á.; Galarza, M.; Thomale, U.; Schuhmann, M. U.; Valero, J.; Amigó, J. M.

2017-05-01

The obstruction of ventricular catheters (VCs) is a major problem in the standard treatment of hydrocephalus, the flow pattern of the cerebrospinal fluid (CSF) being one important factor thereof. As a first approach to this problem, some of the authors studied previously the CSF flow through VCs under time-independent boundary conditions by means of computational fluid dynamics in three-dimensional models. This allowed us to derive a few basic principles which led to designs with improved flow patterns regarding the obstruction problem. However, the flow of the CSF has actually a pulsatile nature because of the heart beating and blood flow. To address this fact, here we extend our previous computational study to models with oscillatory boundary conditions. The new results will be compared with the results for constant flows and discussed. It turns out that the corrections due to the pulsatility of the CSF are quantitatively small, which reinforces our previous findings and conclusions. This article is part of the themed issue `Mathematical methods in medicine: neuroscience, cardiology and pathology'.

The unbiasedness of a generalized mirage boundary correction method for Monte Carlo integration estimators of volume

Treesearch

Thomas B. Lynch; Jeffrey H. Gove

2014-01-01

The typical "double counting" application of the mirage method of boundary correction cannot be applied to sampling systems such as critical height sampling (CHS) that are based on a Monte Carlo sample of a tree (or debris) attribute because the critical height (or other random attribute) sampled from a mirage point is generally not equal to the critical...

A perspective on thirty years of the Webb, Pearman and Leuning density corrections

Treesearch

Xuhui Lee; William J. Massman

2011-01-01

The density correction theory of Webb et al. (1980, Q J Roy Meteorol Soc 106: 85-100, hereafter WPL) is a principle underpinning the experimental investigation of surface fluxes of energy and masses in the atmospheric boundary layer. It has a long-lasting influence in boundary-layer meteorology and micrometeorology, and the year 2010 marks the 30th anniversary of the...

Towards real-time flood forecasting in hydraulics: merits of in situ discharge and water level data assimilation for the modeling of the Marne catchment in France

NASA Astrophysics Data System (ADS)

Ricci, S. M.; Habert, J.; Le Pape, E.; Piacentini, A.; Jonville, G.; Thual, O.; Zaoui, F.

2011-12-01

The present study describes the assimilation of river flow and water level observations and the resulting improvement in flood forecasting. The Kalman Filter algorithm was built on top of the one-dimensional hydraulic model, MASCARET, [1] which describes the Saint-Venant equations. The assimilation algorithm folds in two steps: the first one was based on the assumption that the upstream flow can be adjusted using a three-parameter correction; the second one consisted of directly correcting the hydraulic state. This procedure was previously applied on the Adour Maritime Catchment using water level observations [2]. On average, it was shown that the data assimilation procedure enables an improvement of 80% in the simulated water level over the reanalysis period, 60 % in the forecast water level at a one-hour lead time, and 25% at a twelve-hour lead time. The procedure was then applied on the Marne Catchment, which includes karstic tributaries, located East of the Paris basin, characterized by long flooding periods and strong sensitivity to local precipitations. The objective was to geographically extend and improve the existing model used by the flood forecasting service located in Chalons-en-Champagne. A hydrological study first enabled the specification of boundary conditions (upstream flow or lateral inflow), then the hydraulic model was calibrated using in situ discharge data (adjustment of Strickler coefficients or cross sectional geometry). The assimilation of water level data enabled the reduction of the uncertainty in the hydrological boundary conditions and led to significant improvement of the simulated water level in re-analysis and forecast modes. Still, because of errors in the Strickler coefficients or cross section geometry, the improvement of the simulated water level sometimes resulted in a degradation of discharge values. This problem was overcome by controlling the correction of the hydrological boundary conditions by directly assimilating discharge observations rather than water level observations. As this approach leads to a satisfying simulation of flood events in the Marne catchment in re-analysis and forecast mode, ongoing work aims at controlling Strickler coefficients through data assimilation procedures in order to simultaneously improve the water level and discharge state. [1] N. Goutal, F. Maurel: A finite volume solver for 1D shallow water equations applied to an actual river, Int. J. Numer. Meth. Fluids, 38(2), 1--19, 2002. [2] S. Ricci, A. Piacentini, O. Thual, E. Le Pape, G. Jonville, 2011: Correction of upstream flow and hydraulic state with data assimilation on the context of flood forecasting. Submitted to Hydrol. Earth Syst. Sci, In review.

A dissipative particle dynamics method for arbitrarily complex geometries

NASA Astrophysics Data System (ADS)

Li, Zhen; Bian, Xin; Tang, Yu-Hang; Karniadakis, George Em

2018-02-01

Dissipative particle dynamics (DPD) is an effective Lagrangian method for modeling complex fluids in the mesoscale regime but so far it has been limited to relatively simple geometries. Here, we formulate a local detection method for DPD involving arbitrarily shaped geometric three-dimensional domains. By introducing an indicator variable of boundary volume fraction (BVF) for each fluid particle, the boundary of arbitrary-shape objects is detected on-the-fly for the moving fluid particles using only the local particle configuration. Therefore, this approach eliminates the need of an analytical description of the boundary and geometry of objects in DPD simulations and makes it possible to load the geometry of a system directly from experimental images or computer-aided designs/drawings. More specifically, the BVF of a fluid particle is defined by the weighted summation over its neighboring particles within a cutoff distance. Wall penetration is inferred from the value of the BVF and prevented by a predictor-corrector algorithm. The no-slip boundary condition is achieved by employing effective dissipative coefficients for liquid-solid interactions. Quantitative evaluations of the new method are performed for the plane Poiseuille flow, the plane Couette flow and the Wannier flow in a cylindrical domain and compared with their corresponding analytical solutions and (high-order) spectral element solution of the Navier-Stokes equations. We verify that the proposed method yields correct no-slip boundary conditions for velocity and generates negligible fluctuations of density and temperature in the vicinity of the wall surface. Moreover, we construct a very complex 3D geometry - the "Brown Pacman" microfluidic device - to explicitly demonstrate how to construct a DPD system with complex geometry directly from loading a graphical image. Subsequently, we simulate the flow of a surfactant solution through this complex microfluidic device using the new method. Its effectiveness is demonstrated by examining the rich dynamics of surfactant micelles, which are flowing around multiple small cylinders and stenotic regions in the microfluidic device without wall penetration. In addition to stationary arbitrary-shape objects, the new method is particularly useful for problems involving moving and deformable boundaries, because it only uses local information of neighboring particles and satisfies the desired boundary conditions on-the-fly.

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

NASA Technical Reports Server (NTRS)

Bibel, G. D.; Kumar, A.; Reddy, S.; Handschuh, R.

1993-01-01

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

A discontinuous Poisson-Boltzmann equation with interfacial jump: homogenisation and residual error estimate.

PubMed

Fellner, Klemens; Kovtunenko, Victor A

2016-01-01

A nonlinear Poisson-Boltzmann equation with inhomogeneous Robin type boundary conditions at the interface between two materials is investigated. The model describes the electrostatic potential generated by a vector of ion concentrations in a periodic multiphase medium with dilute solid particles. The key issue stems from interfacial jumps, which necessitate discontinuous solutions to the problem. Based on variational techniques, we derive the homogenisation of the discontinuous problem and establish a rigorous residual error estimate up to the first-order correction.

A Validated All-Pressure Fluid Drop Model and Lewis Number Effects for a Binary Mixture

NASA Technical Reports Server (NTRS)

Harstad, K.; Bellan, J.

1999-01-01

The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Conditions and Linear Stability Analysis at the Transition to Synchronization of Three Coupled Phase Oscillators in a Ring

NASA Astrophysics Data System (ADS)

El-Nashar, Hassan F.

2017-06-01

We consider a system of three nonidentical coupled phase oscillators in a ring topology. We explore the conditions that must be satisfied in order to obtain the phases at the transition to a synchrony state. These conditions lead to the correct mathematical expressions of phases that aid to find a simple analytic formula for critical coupling when the oscillators transit to a synchronization state having a common frequency value. The finding of a simple expression for the critical coupling allows us to perform a linear stability analysis at the transition to the synchronization stage. The obtained analytic forms of the eigenvalues show that the three coupled phase oscillators with periodic boundary conditions transit to a synchrony state when a saddle-node bifurcation occurs.

Exact solution and precise asymptotics of a Fisher-KPP type front

NASA Astrophysics Data System (ADS)

Berestycki, Julien; Brunet, Éric; Derrida, Bernard

2018-01-01

The present work concerns a version of the Fisher-KPP equation where the nonlinear term is replaced by a saturation mechanism, yielding a free boundary problem with mixed conditions. Following an idea proposed in Brunet and Derrida (2015 J. Stat. Phys. 161 801), we show that the Laplace transform of the initial condition is directly related to some functional of the front position μt . We then obtain precise asymptotics of the front position by means of singularity analysis. In particular, we recover the so-called Ebert and van Saarloos correction (Ebert and van Saarloos 2000 Physica D 146 1), we obtain an additional term of order log t /t in this expansion, and we give precise conditions on the initial condition for those terms to be present.

Interferogram conditioning for improved Fourier analysis and application to X-ray phase imaging by grating interferometry.

PubMed

Montaux-Lambert, Antoine; Mercère, Pascal; Primot, Jérôme

2015-11-02

An interferogram conditioning procedure, for subsequent phase retrieval by Fourier demodulation, is presented here as a fast iterative approach aiming at fulfilling the classical boundary conditions imposed by Fourier transform techniques. Interference fringe patterns with typical edge discontinuities were simulated in order to reveal the edge artifacts that classically appear in traditional Fourier analysis, and were consecutively used to demonstrate the correction efficiency of the proposed conditioning technique. Optimization of the algorithm parameters is also presented and discussed. Finally, the procedure was applied to grating-based interferometric measurements performed in the hard X-ray regime. The proposed algorithm enables nearly edge-artifact-free retrieval of the phase derivatives. A similar enhancement of the retrieved absorption and fringe visibility images is also achieved.

A class of renormalised meshless Laplacians for boundary value problems

NASA Astrophysics Data System (ADS)

Basic, Josip; Degiuli, Nastia; Ban, Dario

2018-02-01

A meshless approach to approximating spatial derivatives on scattered point arrangements is presented in this paper. Three various derivations of approximate discrete Laplace operator formulations are produced using the Taylor series expansion and renormalised least-squares correction of the first spatial derivatives. Numerical analyses are performed for the introduced Laplacian formulations, and their convergence rate and computational efficiency are examined. The tests are conducted on regular and highly irregular scattered point arrangements. The results are compared to those obtained by the smoothed particle hydrodynamics method and the finite differences method on a regular grid. Finally, the strong form of various Poisson and diffusion equations with Dirichlet or Robin boundary conditions are solved in two and three dimensions by making use of the introduced operators in order to examine their stability and accuracy for boundary value problems. The introduced Laplacian operators perform well for highly irregular point distribution and offer adequate accuracy for mesh and mesh-free numerical methods that require frequent movement of the grid or point cloud.

Building a Lego wall: Sequential action selection.

PubMed

Arnold, Amy; Wing, Alan M; Rotshtein, Pia

2017-05-01

The present study draws together two distinct lines of enquiry into the selection and control of sequential action: motor sequence production and action selection in everyday tasks. Participants were asked to build 2 different Lego walls. The walls were designed to have hierarchical structures with shared and dissociated colors and spatial components. Participants built 1 wall at a time, under low and high load cognitive states. Selection times for correctly completed trials were measured using 3-dimensional motion tracking. The paradigm enabled precise measurement of the timing of actions, while using real objects to create an end product. The experiment demonstrated that action selection was slowed at decision boundary points, relative to boundaries where no between-wall decision was required. Decision points also affected selection time prior to the actual selection window. Dual-task conditions increased selection errors. Errors mostly occurred at boundaries between chunks and especially when these required decisions. The data support hierarchical control of sequenced behavior. (PsycINFO Database Record (c) 2017 APA, all rights reserved).

Exploratory Calibration of Adjustable-Protrusion Surface-Obstacle (APSO) Skin Friction Vector Gage

NASA Technical Reports Server (NTRS)

Hakkinen, Raimo J.; Neubauer, Jeremy S.; Hamory, Philip J.; Bui, Trong T.; Noffz, Gregory K.; Young, Ron (Technical Monitor)

2003-01-01

The design of an adjustable-protrusion surface-obstacle (APSO) skin friction vector gage is presented. Results from exploratory calibrations conducted in laminar and turbulent boundary layers at the Washington University Low-Speed Wind Tunnel and for turbulent boundary layers at speeds up to Mach 2 on the ceiling of the NASA Glenn Research Center 8- X 6-ft Supersonic Wind Tunnel are also discussed. The adjustable-height gage was designed to yield both the magnitude and direction of the surface shear stress vector and to measure the local static pressure distribution. Results from the NASA test show good correlation for subsonic and low supersonic conditions covering several orders of magnitude in terms of the adopted similarity variables. Recommendations for future work in this area consist of identifying the physical parameters responsible for the disagreement between the university and NASA data sets, developing a compressibility correction specific to the APSO geometry, and examining the effect that static pressure distribution and skewed boundary layers have on the results from the APSO.

Corrected simulations for one-dimensional diffusion processes with naturally occurring boundaries.

PubMed

Shafiey, Hassan; Gan, Xinjun; Waxman, David

2017-11-01

To simulate a diffusion process, a usual approach is to discretize the time in the associated stochastic differential equation. This is the approach used in the Euler method. In the present work we consider a one-dimensional diffusion process where the terms occurring, within the stochastic differential equation, prevent the process entering a region. The outcome is a naturally occurring boundary (which may be absorbing or reflecting). A complication occurs in a simulation of this situation. The term involving a random variable, within the discretized stochastic differential equation, may take a trajectory across the boundary into a "forbidden region." The naive way of dealing with this problem, which we refer to as the "standard" approach, is simply to reset the trajectory to the boundary, based on the argument that crossing the boundary actually signifies achieving the boundary. In this work we show, within the framework of the Euler method, that such resetting introduces a spurious force into the original diffusion process. This force may have a significant influence on trajectories that come close to a boundary. We propose a corrected numerical scheme, for simulating one-dimensional diffusion processes with naturally occurring boundaries. This involves correcting the standard approach, so that an exact property of the diffusion process is precisely respected. As a consequence, the proposed scheme does not introduce a spurious force into the dynamics. We present numerical test cases, based on exactly soluble one-dimensional problems with one or two boundaries, which suggest that, for a given value of the discrete time step, the proposed scheme leads to substantially more accurate results than the standard approach. Alternatively, the standard approach needs considerably more computation time to obtain a comparable level of accuracy to the proposed scheme, because the standard approach requires a significantly smaller time step.

Corrected simulations for one-dimensional diffusion processes with naturally occurring boundaries

NASA Astrophysics Data System (ADS)

Shafiey, Hassan; Gan, Xinjun; Waxman, David

2017-11-01

To simulate a diffusion process, a usual approach is to discretize the time in the associated stochastic differential equation. This is the approach used in the Euler method. In the present work we consider a one-dimensional diffusion process where the terms occurring, within the stochastic differential equation, prevent the process entering a region. The outcome is a naturally occurring boundary (which may be absorbing or reflecting). A complication occurs in a simulation of this situation. The term involving a random variable, within the discretized stochastic differential equation, may take a trajectory across the boundary into a "forbidden region." The naive way of dealing with this problem, which we refer to as the "standard" approach, is simply to reset the trajectory to the boundary, based on the argument that crossing the boundary actually signifies achieving the boundary. In this work we show, within the framework of the Euler method, that such resetting introduces a spurious force into the original diffusion process. This force may have a significant influence on trajectories that come close to a boundary. We propose a corrected numerical scheme, for simulating one-dimensional diffusion processes with naturally occurring boundaries. This involves correcting the standard approach, so that an exact property of the diffusion process is precisely respected. As a consequence, the proposed scheme does not introduce a spurious force into the dynamics. We present numerical test cases, based on exactly soluble one-dimensional problems with one or two boundaries, which suggest that, for a given value of the discrete time step, the proposed scheme leads to substantially more accurate results than the standard approach. Alternatively, the standard approach needs considerably more computation time to obtain a comparable level of accuracy to the proposed scheme, because the standard approach requires a significantly smaller time step.

Validation of the k- ω turbulence model for the thermal boundary layer profile of effusive cooled walls

NASA Astrophysics Data System (ADS)

Hink, R.

2015-09-01

The choice of materials for rocket chamber walls is limited by its thermal resistance. The thermal loads can be reduced substantially by the blowing out of gases through a porous surface. The k- ω-based turbulence models for computational fluid dynamic simulations are designed for smooth, non-permeable walls and have to be adjusted to account for the influence of injected fluids. Wilcox proposed therefore an extension for the k- ω turbulence model for the correct prediction of turbulent boundary layer velocity profiles. In this study, this extension is validated against experimental thermal boundary layer data from the Thermosciences Division of the Department of Mechanical Engineering from the Stanford University. All simulations are performed with a finite volume-based in-house code of the German Aerospace Center. Several simulations with different blowing settings were conducted and discussed in comparison to the results of the original model and in comparison to an additional roughness implementation. This study has permitted to understand that velocity profile corrections are necessary in contrast to additional roughness corrections to predict the correct thermal boundary layer profile of effusive cooled walls. Finally, this approach is applied to a two-dimensional simulation of an effusive cooled rocket chamber wall.

Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to Robin boundary condition

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pan, Wenxiao; Bao, Jie; Tartakovsky, Alexandre M.

2014-02-15

Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to variousmore » forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.« less

Generalization of the Hartree-Fock approach to collision processes

NASA Astrophysics Data System (ADS)

Hahn, Yukap

1997-06-01

The conventional Hartree and Hartree-Fock approaches for bound states are generalized to treat atomic collision processes. All the single-particle orbitals, for both bound and scattering states, are determined simultaneously by requiring full self-consistency. This generalization is achieved by introducing two Ansäauttze: (a) the weak asymptotic boundary condition, which maintains the correct scattering energy and target orbitals with correct number of nodes, and (b) square integrable amputated scattering functions to generate self-consistent field (SCF) potentials for the target orbitals. The exact initial target and final-state asymptotic wave functions are not required and thus need not be specified a priori, as they are determined simultaneously by the SCF iterations. To check the asymptotic behavior of the solution, the theory is applied to elastic electron-hydrogen scattering at low energies. The solution is found to be stable and the weak asymptotic condition is sufficient to produce the correct scattering amplitudes. The SCF potential for the target orbital shows the strong penetration by the projectile electron during the collision, but the exchange term tends to restore the original form. Potential applicabilities of this extension are discussed, including the treatment of ionization and shake-off processes.

Preferred color correction for digital LCD TVs

NASA Astrophysics Data System (ADS)

Kim, Kyoung Tae; Kim, Choon-Woo; Ahn, Ji-Young; Kang, Dong-Woo; Shin, Hyun-Ho

2009-01-01

Instead of colorimetirc color reproduction, preferred color correction is applied for digital TVs to improve subjective image quality. First step of the preferred color correction is to survey the preferred color coordinates of memory colors. This can be achieved by the off-line human visual tests. Next step is to extract pixels of memory colors representing skin, grass and sky. For the detected pixels, colors are shifted towards the desired coordinates identified in advance. This correction process may result in undesirable contours on the boundaries between the corrected and un-corrected areas. For digital TV applications, the process of extraction and correction should be applied in every frame of the moving images. This paper presents a preferred color correction method in LCH color space. Values of chroma and hue are corrected independently. Undesirable contours on the boundaries of correction are minimized. The proposed method change the coordinates of memory color pixels towards the target color coordinates. Amount of correction is determined based on the averaged coordinate of the extracted pixels. The proposed method maintains the relative color difference within memory color areas. Performance of the proposed method is evaluated using the paired comparison. Results of experiments indicate that the proposed method can reproduce perceptually pleasing images to viewers.

Turbulent Transfer Between Street Canyons and the Overlying Atmospheric Boundary Layer

NASA Astrophysics Data System (ADS)

Salizzoni, Pietro; Marro, Massimo; Soulhac, Lionel; Grosjean, Nathalie; Perkins, Richard J.

2011-12-01

The turbulent exchange of momentum between a two-dimensional cavity and the overlying boundary layer has been studied experimentally, using hot-wire anemometry and particle image velocimetry (PIV). Conditions within the boundary layer were varied by changing the width of the canyons upstream of the test canyon, whilst maintaining the square geometry of the test canyon. The results show that turbulent transfer is due to the coupling between the instabilities generated in the shear layer above the canyons and the turbulent structures in the oncoming boundary layer. As a result, there is no single, unique velocity scale that correctly characterizes all the processes involved in the turbulent exchange of momentum across the boundary layer. Similarly, there is no single velocity scale that can characterize the different properties of the turbulent flow within the canyon, which depends strongly on the way in which turbulence from the outer flow is entrained into the cavity and carried round by the mean flow. The results from this study will be useful in developing simple parametrizations for momentum exchange in the urban canopy, in situations where the street geometry consists principally of relatively long, uniform streets arranged in grid-like patterns; they are unlikely to be applicable to sparse geometries composed of isolated three-dimensional obstacles.

Capabilities of wind tunnels with two-adaptive walls to minimize boundary interference in 3-D model testing

NASA Technical Reports Server (NTRS)

Rebstock, Rainer; Lee, Edwin E., Jr.

1989-01-01

An initial wind tunnel test was made to validate a new wall adaptation method for 3-D models in test sections with two adaptive walls. First part of the adaptation strategy is an on-line assessment of wall interference at the model position. The wall induced blockage was very small at all test conditions. Lift interference occurred at higher angles of attack with the walls set aerodynamically straight. The adaptation of the top and bottom tunnel walls is aimed at achieving a correctable flow condition. The blockage was virtually zero throughout the wing planform after the wall adjustment. The lift curve measured with the walls adapted agreed very well with interference free data for Mach 0.7, regardless of the vertical position of the wing in the test section. The 2-D wall adaptation can significantly improve the correctability of 3-D model data. Nevertheless, residual spanwise variations of wall interference are inevitable.

Inflow/Outflow Boundary Conditions with Application to FUN3D

NASA Technical Reports Server (NTRS)

Carlson, Jan-Renee

2011-01-01

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.

Fast words boundaries localization in text fields for low quality document images

NASA Astrophysics Data System (ADS)

Ilin, Dmitry; Novikov, Dmitriy; Polevoy, Dmitry; Nikolaev, Dmitry

2018-04-01

The paper examines the problem of word boundaries precise localization in document text zones. Document processing on a mobile device consists of document localization, perspective correction, localization of individual fields, finding words in separate zones, segmentation and recognition. While capturing an image with a mobile digital camera under uncontrolled capturing conditions, digital noise, perspective distortions or glares may occur. Further document processing gets complicated because of its specifics: layout elements, complex background, static text, document security elements, variety of text fonts. However, the problem of word boundaries localization has to be solved at runtime on mobile CPU with limited computing capabilities under specified restrictions. At the moment, there are several groups of methods optimized for different conditions. Methods for the scanned printed text are quick but limited only for images of high quality. Methods for text in the wild have an excessively high computational complexity, thus, are hardly suitable for running on mobile devices as part of the mobile document recognition system. The method presented in this paper solves a more specialized problem than the task of finding text on natural images. It uses local features, a sliding window and a lightweight neural network in order to achieve an optimal algorithm speed-precision ratio. The duration of the algorithm is 12 ms per field running on an ARM processor of a mobile device. The error rate for boundaries localization on a test sample of 8000 fields is 0.3

Low-to-moderate Reynolds number swirling flow in an annular channel with a rotating end wall.

PubMed

Davoust, Laurent; Achard, Jean-Luc; Drazek, Laurent

2015-02-01

This paper presents a new method for solving analytically the axisymmetric swirling flow generated in a finite annular channel from a rotating end wall, with no-slip boundary conditions along stationary side walls and a slip condition along the free surface opposite the rotating floor. In this case, the end-driven swirling flow can be described from the coupling between an azimuthal shear flow and a two-dimensional meridional flow driven by the centrifugal force along the rotating floor. A regular asymptotic expansion based on a small but finite Reynolds number is used to calculate centrifugation-induced first-order correction to the azimuthal Stokes flow obtained as the solution at leading order. For solving the first-order problem, the use of an integral boundary condition for the vorticity is found to be a convenient way to attribute boundary conditions in excess for the stream function to the vorticity. The annular geometry is characterized by both vertical and horizontal aspect ratios, whose respective influences on flow patterns are investigated. The vertical aspect ratio is found to involve nontrivial changes in flow patterns essentially due to the role of corner eddies located on the left and right sides of the rotating floor. The present analytical method can be ultimately extended to cylindrical geometries, irrespective of the surface opposite the rotating floor: a wall or a free surface. It can also serve as an analytical tool for monitoring confined rotating flows in applications related to surface viscosimetry or crystal growth from the melt.

Advances in Numerical Boundary Conditions for Computational Aeroacoustics

NASA Technical Reports Server (NTRS)

Tam, Christopher K. W.

1997-01-01

Advances in Computational Aeroacoustics (CAA) depend critically on the availability of accurate, nondispersive, least dissipative computation algorithm as well as high quality numerical boundary treatments. This paper focuses on the recent developments of numerical boundary conditions. In a typical CAA problem, one often encounters two types of boundaries. Because a finite computation domain is used, there are external boundaries. On the external boundaries, boundary conditions simulating the solution outside the computation domain are to be imposed. Inside the computation domain, there may be internal boundaries. On these internal boundaries, boundary conditions simulating the presence of an object or surface with specific acoustic characteristics are to be applied. Numerical boundary conditions, both external or internal, developed for simple model problems are reviewed and examined. Numerical boundary conditions for real aeroacoustic problems are also discussed through specific examples. The paper concludes with a description of some much needed research in numerical boundary conditions for CAA.

A fast multipole method combined with a reaction field for long-range electrostatics in molecular dynamics simulations: The effects of truncation on the properties of water

NASA Astrophysics Data System (ADS)

Mathias, Gerald; Egwolf, Bernhard; Nonella, Marco; Tavan, Paul

2003-06-01

We present a combination of the structure adapted multipole method with a reaction field (RF) correction for the efficient evaluation of electrostatic interactions in molecular dynamics simulations under periodic boundary conditions. The algorithm switches from an explicit electrostatics evaluation to a continuum description at the maximal distance that is consistent with the minimum image convention, and, thus, avoids the use of a periodic electrostatic potential. A physically motivated switching function enables charge clusters interacting with a given charge to smoothly move into the solvent continuum by passing through the spherical dielectric boundary surrounding this charge. This transition is complete as soon as the cluster has reached the so-called truncation radius Rc. The algorithm is used to examine the dependence of thermodynamic properties and correlation functions on Rc in the three point transferable intermolecular potential water model. Our test simulations on pure liquid water used either the RF correction or a straight cutoff and values of Rc ranging from 14 Å to 40 Å. In the RF setting, the thermodynamic properties and the correlation functions show convergence for Rc increasing towards 40 Å. In the straight cutoff case no such convergence is found. Here, in particular, the dipole-dipole correlation functions become completely artificial. The RF description of the long-range electrostatics is verified by comparison with the results of a particle-mesh Ewald simulation at identical conditions.

Phase separation in the six-vertex model with a variety of boundary conditions

NASA Astrophysics Data System (ADS)

Lyberg, I.; Korepin, V.; Ribeiro, G. A. P.; Viti, J.

2018-05-01

We present numerical results for the six-vertex model with a variety of boundary conditions. Adapting an algorithm for domain wall boundary conditions, proposed in the work of Allison and Reshetikhin [Ann. Inst. Fourier 55(6), 1847-1869 (2005)], we examine some modifications of these boundary conditions. To be precise, we discuss partial domain wall boundary conditions, reflecting ends, and half turn boundary conditions (domain wall boundary conditions with half turn symmetry). Dedicated to the memory of Ludwig Faddeev

External Peer Review Team Report Underground Testing Area Subproject for Frenchman Flat, Revision 1

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sam Marutzky

2010-09-01

An external peer review was conducted to review the groundwater models used in the corrective action investigation stage of the Underground Test Area (UGTA) subproject to forecast zones of potential contamination in 1,000 years for the Frenchman Flat area. The goal of the external peer review was to provide technical evaluation of the studies and to assist in assessing the readiness of the UGTA subproject to progress to monitoring activities for further model evaluation. The external peer review team consisted of six independent technical experts with expertise in geology, hydrogeology,'''groundwater modeling, and radiochemistry. The peer review team was tasked withmore » addressing the following questions: 1. Are the modeling approaches, assumptions, and model results for Frenchman Flat consistent with the use of modeling studies as a decision tool for resolution of environmental and regulatory requirements? 2. Do the modeling results adequately account for uncertainty in models of flow and transport in the Frenchman Flat hydrological setting? a. Are the models of sufficient scale/resolution to adequately predict contaminant transport in the Frenchman Flat setting? b. Have all key processes been included in the model? c. Are the methods used to forecast contaminant boundaries from the transport modeling studies reasonable and appropriate? d. Are the assessments of uncertainty technically sound and consistent with state-of-the-art approaches currently used in the hydrological sciences? 3. Are the datasets and modeling results adequate for a transition to Corrective Action Unit monitoring studies—the next stage in the UGTA strategy for Frenchman Flat? The peer review team is of the opinion that, with some limitations, the modeling approaches, assumptions, and model results are consistent with the use of modeling studies for resolution of environmental and regulatory requirements. The peer review team further finds that the modeling studies have accounted for uncertainty in models of flow and transport in the Frenchman Flat except for a few deficiencies described in the report. Finally, the peer review team concludes that the UGTA subproject has explored a wide range of variations in assumptions, methods, and data, and should proceed to the next stage with an emphasis on monitoring studies. The corrective action strategy, as described in the Federal Facility Agreement and Consent Order, states that the groundwater flow and transport models for each corrective action unit will consider, at a minimum, the following: • Alternative hydrostratigraphic framework models of the modeling domain. • Uncertainty in the radiological and hydrological source terms. • Alternative models of recharge. • Alternative boundary conditions and groundwater flows. • Multiple permissive sets of calibrated flow models. • Probabilistic simulations of transport using plausible sets of alternative framework and recharge models, and boundary and groundwater flows from calibrated flow models. • Ensembles of forecasts of contaminant boundaries. • Sensitivity and uncertainty analyses of model outputs. The peer review team finds that these minimum requirements have been met. While the groundwater modeling and uncertainty analyses have been quite detailed, the peer review team has identified several modeling-related issues that should be addressed in the next phase of the corrective action activities: • Evaluating and using water-level gradients from the pilot wells at the Area 5 Radioactive Waste Management Site in model calibration. • Re-evaluating the use of geochemical age-dating data to constrain model calibrations. • Developing water budgets for the alluvial and upper volcanic aquifer systems in Frenchman Flat. • Considering modeling approaches in which calculated groundwater flow directions near the water table are not predetermined by model boundary conditions and areas of recharge, all of which are very uncertain. • Evaluating local-scale variations in hydraulic conductivity on the calculated contaminant boundaries. • Evaluating the effects of non-steady-state flow conditions on calculated contaminant boundaries, including the effects of long-term declines in water levels, climatic change, and disruption of groundwater system by potential earthquake faulting along either of the two major controlling fault zones in the flow system (the Cane Spring and Rock Valley faults). • Considering the use of less-complex modeling approaches. • Evaluating the large change in water levels in the vicinity of the Frenchman Flat playa and developing a conceptual model to explain these water-level changes. • Developing a long-term groundwater level monitoring program for Frenchman Flat with regular monitoring of water levels at key monitoring wells. Despite these reservations, the peer review team strongly believes that the UGTA subproject should proceed to the next stage.« less

Regional Model Nesting Within GFS Daily Forecasts Over West Africa

NASA Technical Reports Server (NTRS)

Druyan, Leonard M.; Fulakeza, Matthew; Lonergan, Patrick; Worrell, Ruben

2010-01-01

The study uses the RM3, the regional climate model at the Center for Climate Systems Research of Columbia University and the NASA/Goddard Institute for Space Studies (CCSR/GISS). The paper evaluates 30 48-hour RM3 weather forecasts over West Africa during September 2006 made on a 0.5 grid nested within 1 Global Forecast System (GFS) global forecasts. September 2006 was the Special Observing Period #3 of the African Monsoon Multidisciplinary Analysis (AMMA). Archived GFS initial conditions and lateral boundary conditions for the simulations from the US National Weather Service, National Oceanographic and Atmospheric Administration were interpolated four times daily. Results for precipitation forecasts are validated against Tropical Rainfall Measurement Mission (TRMM) satellite estimates and data from the Famine Early Warning System (FEWS), which includes rain gauge measurements, and forecasts of circulation are compared to reanalysis 2. Performance statistics for the precipitation forecasts include bias, root-mean-square errors and spatial correlation coefficients. The nested regional model forecasts are compared to GFS forecasts to gauge whether nesting provides additional realistic information. They are also compared to RM3 simulations driven by reanalysis 2, representing high potential skill forecasts, to gauge the sensitivity of results to lateral boundary conditions. Nested RM3/GFS forecasts generate excessive moisture advection toward West Africa, which in turn causes prodigious amounts of model precipitation. This problem is corrected by empirical adjustments in the preparation of lateral boundary conditions and initial conditions. The resulting modified simulations improve on the GFS precipitation forecasts, achieving time-space correlations with TRMM of 0.77 on the first day and 0.63 on the second day. One realtime RM3/GFS precipitation forecast made at and posted by the African Centre of Meteorological Application for Development (ACMAD) in Niamey, Niger is shown.

A new real-time guidance strategy for aerodynamic ascent flight

NASA Astrophysics Data System (ADS)

Yamamoto, Takayuki; Kawaguchi, Jun'ichiro

2007-12-01

Reusable launch vehicles are conceived to constitute the future space transportation system. If these vehicles use air-breathing propulsion and lift taking-off horizontally, the optimal steering for these vehicles exhibits completely different behavior from that in conventional rockets flight. In this paper, the new guidance strategy is proposed. This method derives from the optimality condition as for steering and an analysis concludes that the steering function takes the form comprised of Linear and Logarithmic terms, which include only four parameters. The parameter optimization of this method shows the acquired terminal horizontal velocity is almost same with that obtained by the direct numerical optimization. This supports the parameterized Liner Logarithmic steering law. And here is shown that there exists a simple linear relation between the terminal states and the parameters to be corrected. The relation easily makes the parameters determined to satisfy the terminal boundary conditions in real-time. The paper presents the guidance results for the practical application cases. The results show the guidance is well performed and satisfies the terminal boundary conditions specified. The strategy built and presented here does guarantee the robust solution in real-time excluding any optimization process, and it is found quite practical.

A new class of N=2 topological amplitudes

NASA Astrophysics Data System (ADS)

Antoniadis, I.; Hohenegger, S.; Narain, K. S.; Sokatchev, E.

2009-12-01

We describe a new class of N=2 topological amplitudes that compute a particular class of BPS terms in the low energy effective supergravity action. Specifically they compute the coupling F(( where F, λ and ϕ are gauge field strengths, gaugino and holomorphic vector multiplet scalars. The novel feature of these terms is that they depend both on the vector and hypermultiplet moduli. The BPS nature of these terms implies that they satisfy a holomorphicity condition with respect to vector moduli and a harmonicity condition as well as a second order differential equation with respect to hypermultiplet moduli. We study these conditions explicitly in heterotic string theory and show that they are indeed satisfied up to anomalous boundary terms in the world-sheet moduli space. We also analyze the boundary terms in the holomorphicity and harmonicity equations at a generic point in the vector and hyper moduli space. In particular we show that the obstruction to the holomorphicity arises from the one loop threshold correction to the gauge couplings and we argue that this is due to the contribution of non-holomorphic couplings to the connected graphs via elimination of the auxiliary fields.

A continuum treatment of sliding in Eulerian simulations of solid-solid and solid-fluid interfaces

NASA Astrophysics Data System (ADS)

Subramaniam, Akshay; Ghaisas, Niranjan; Lele, Sanjiva

2017-11-01

A novel treatment of sliding is developed for use in an Eulerian framework for simulating elastic-plastic deformations of solids coupled with fluids. In this method, embedded interfacial boundary conditions for perfect sliding are imposed by enforcing the interface normal to be a principal direction of the Cauchy stress and appropriate consistency conditions ensure correct transmission and reflection of waves at the interface. This sliding treatment may be used either to simulate a solid-solid sliding interface or to incorporate an internal slip boundary condition at a solid-fluid interface. Sliding laws like the Coulomb friction law can also be incorporated with relative ease into this framework. Simulations of sliding interfaces are conducted using a 10th order compact finite difference scheme and a Localized Artificial Diffusivity (LAD) scheme for shock and interface capturing. 1D and 2D simulations are used to assess the accuracy of the sliding treatment. The Richmyer-Meshkov instability between copper and aluminum is simulated with this sliding treatment as a demonstration test case. Support for this work was provided through Grant B612155 from the Lawrence Livermore National Laboratory, US Department of Energy.

Scalable improvement of SPME multipolar electrostatics in anisotropic polarizable molecular mechanics using a general short-range penetration correction up to quadrupoles.

PubMed

Narth, Christophe; Lagardère, Louis; Polack, Étienne; Gresh, Nohad; Wang, Qiantao; Bell, David R; Rackers, Joshua A; Ponder, Jay W; Ren, Pengyu Y; Piquemal, Jean-Philip

2016-02-15

We propose a general coupling of the Smooth Particle Mesh Ewald SPME approach for distributed multipoles to a short-range charge penetration correction modifying the charge-charge, charge-dipole and charge-quadrupole energies. Such an approach significantly improves electrostatics when compared to ab initio values and has been calibrated on Symmetry-Adapted Perturbation Theory reference data. Various neutral molecular dimers have been tested and results on the complexes of mono- and divalent cations with a water ligand are also provided. Transferability of the correction is adressed in the context of the implementation of the AMOEBA and SIBFA polarizable force fields in the TINKER-HP software. As the choices of the multipolar distribution are discussed, conclusions are drawn for the future penetration-corrected polarizable force fields highlighting the mandatory need of non-spurious procedures for the obtention of well balanced and physically meaningful distributed moments. Finally, scalability and parallelism of the short-range corrected SPME approach are addressed, demonstrating that the damping function is computationally affordable and accurate for molecular dynamics simulations of complex bio- or bioinorganic systems in periodic boundary conditions. Copyright © 2016 Wiley Periodicals, Inc.

Analytic Corrections to CFD Heating Predictions Accounting for Changes in Surface Catalysis. Part II

NASA Technical Reports Server (NTRS)

Gnoffo, Peter A.; Inger, George R.

1996-01-01

A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (but time intensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations is described. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. It allows a designer at a work-station to address two questions, given a single CFD solution. (1) How much does the heating change for a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change at the interface of two different TPS materials with an abrupt change in catalytic efficiency? The answer to the second question is particularly important, because abrupt changes from low to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. Capabilities of this approach for application to Reusable Launch Vehicle (RLV) design are demonstrated. If the definition of surface catalysis is uncertain early in the design process, results show that fully catalytic wall boundary conditions provide the best baseline for CFD design points.

Numerical investigation of nonlinear fluid-structure interaction dynamic behaviors under a general Immersed Boundary-Lattice Boltzmann-Finite Element method

NASA Astrophysics Data System (ADS)

Gong, Chun-Lin; Fang, Zhe; Chen, Gang

A numerical approach based on the immersed boundary (IB), lattice Boltzmann and nonlinear finite element method (FEM) is proposed to simulate hydrodynamic interactions of very flexible objects. In the present simulation framework, the motion of fluid is obtained by solving the discrete lattice Boltzmann equations on Eulerian grid, the behaviors of flexible objects are calculated through nonlinear dynamic finite element method, and the interactive forces between them are implicitly obtained using velocity correction IB method which satisfies the no-slip conditions well at the boundary points. The efficiency and accuracy of the proposed Immersed Boundary-Lattice Boltzmann-Finite Element method is first validated by a fluid-structure interaction (F-SI) benchmark case, in which a flexible filament flaps behind a cylinder in channel flow, then the nonlinear vibration mechanism of the cylinder-filament system is investigated by altering the Reynolds number of flow and the material properties of filament. The interactions between two tandem and side-by-side identical objects in a uniform flow are also investigated, and the in-phase and out-of-phase flapping behaviors are captured by the proposed method.

Prototypical Concepts and Misconceptions of Plate Tectonic Boundaries

NASA Astrophysics Data System (ADS)

Sibley, D. F.; Patino, L. C.

2003-12-01

Students of geology encounter many prototypical/exemplar concepts* that include representative, but not necessarily defining, features and characteristics. This study of students' prototypical representations of plate tectonic boundaries indicates that their representations are rich sources of information about their misconceptions about plate tectonics. After lectures in plate tectonics and mountain building, 353 students in a general education geology class were asked to draw a continent-continent convergent boundary. For this study, a correct answer is defined as having the major features in correct proportions as depicted in the plate boundary diagrams on the USGS web. Fifty-two percent of the drawings were either incorrect or incomplete such that they could not be interpreted. Only 48% were readily interpretable, and of these 22% drew the boundary correctly, showing a thickening of crust where two continents collide. Thirty-three percent drew the boundary showing concave slabs of continental crust as one might imagine two pieces of firm rubber pushed together on a rigid surface and 45% depicted mountains as one might imagine inverted ice cream cones on a rigid plank. Twenty-one senior class geology majors and graduate students were given the same assignment. Forty-eight percent rendered a correct drawing, whereas 38% drew the same ice cream cone on a plank type picture that 45% of the general education students drew. In a second class of 12 geology majors, only 1 student drew a cross section of a continent-ocean boundary similar to standard representation. Four of 12 drew mountains on the top of continental crust over a subduction zone but did not draw a compensating mass within the crust or lithosphere. Prototypical drawings provide more information about students' concepts than do most multiple-choice questions. For example, sixty-two percent of theses students who drew mountains similar to foam rubber pads pushed together on a desk or ice cream cones on a plank correctly answered a multiple-choice question that would appear to indicate a better understanding than the drawings reveal. Furthermore, 12 interviewed students made statements that could be interpreted to indicate that they understood the concept of mountain building at plate tectonic boundaries better than their drawings suggest. Incoherence of multiple-choice responses, verbal statements and drawings may be common in novice learners. If cognitive scientists are correct in their model of multiple types of mental representations for the same term, then the fact that novices may hold inconsistent representations is not surprising. The fact that students at various academic levels draw very similar prototypes that are incorrect is evidence that students have distinct and persistent prototype misconceptions. * Cognitive scientists define a prototypical/exemplar concept as a mental representation of the best examples or central tendencies of a term.

Toward a virtual building laboratory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Klems, J.H.; Finlayson, E.U.; Olsen, T.H.

1999-03-01

In order to achieve in a timely manner the large energy and dollar savings technically possible through improvements in building energy efficiency, it will be necessary to solve the problem of design failure risk. The most economical method of doing this would be to learn to calculate building performance with sufficient detail, accuracy and reliability to avoid design failure. Existing building simulation models (BSM) are a large step in this direction, but are still not capable of this level of modeling. Developments in computational fluid dynamics (CFD) techniques now allow one to construct a road map from present BSM's tomore » a complete building physical model. The most useful first step is a building interior model (BIM) that would allow prediction of local conditions affecting occupant health and comfort. To provide reliable prediction a BIM must incorporate the correct physical boundary conditions on a building interior. Doing so raises a number of specific technical problems and research questions. The solution of these within a context useful for building research and design is not likely to result from other research on CFD, which is directed toward the solution of different types of problems. A six-step plan for incorporating the correct boundary conditions within the context of the model problem of a large atrium has been outlined. A promising strategy for constructing a BIM is the overset grid technique for representing a building space in a CFD calculation. This technique promises to adapt well to building design and allows a step-by-step approach. A state-of-the-art CFD computer code using this technique has been adapted to the problem and can form the departure point for this research.« less

Generation of Accurate Lateral Boundary Conditions for a Surface-Water Groundwater Interaction Model

NASA Astrophysics Data System (ADS)

Khambhammettu, P.; Tsou, M.; Panday, S. M.; Kool, J.; Wei, X.

2010-12-01

The 106 mile long Peace River in Florida flows south from Lakeland to Charlotte Harbor and has a drainage basin of approximately 2,350 square miles. A long-term decline in stream flows and groundwater potentiometric levels has been observed in the region. Long-term trends in rainfall, along with effects of land use changes on runoff, surface-water storage, recharge and evapotranspiration patterns, and increased groundwater and surface-water withdrawals have contributed to this decline. The South West Florida Water Management District (SWFWMD) has funded the development of the Peace River Integrated Model (PRIM) to assess the effects of land use, water use, and climatic changes on stream flows and to evaluate the effectiveness of various management alternatives for restoring stream flows. The PRIM was developed using MODHMS, a fully integrated surface-water groundwater flow and transport simulator developed by HydroGeoLogic, Inc. The development of the lateral boundary conditions (groundwater inflow and outflow) for the PRIM in both historical and predictive contexts is discussed in this presentation. Monthly-varying specified heads were used to define the lateral boundary conditions for the PRIM. These head values were derived from the coarser Southern District Groundwater Model (SDM). However, there were discrepancies between the simulated SDM heads and measured heads: the likely causes being spatial (use of a coarser grid) and temporal (monthly average pumping rates and recharge rates) approximations in the regional SDM. Finer re-calibration of the SDM was not feasible, therefore, an innovative approach was adopted to remove the discrepancies. In this approach, point discrepancies/residuals between the observed and simulated heads were kriged with an appropriate variogram to generate a residual surface. This surface was then added to the simulated head surface of the SDM to generate a corrected head surface. This approach preserves the trends associated with groundwater pumping / recharge in the SDM and adds the kriged residual surface as variations back to the trend. The variations could be from the scale effects of grid resolution and from the temporal averaging of stresses (pumping, recharge, etc.,). The validity of the approach is demonstrated by visual and statistical comparison of the observed and simulated heads before and after correction. For predictive simulations, an Artificial Neural Network was trained to predict heads at monitoring wells based on precipitation and pumping. These predicted head values could then be used as surrogate observations for correcting the results of the regional SDM. In summary, an appropriate approach to link a regional groundwater model to a detailed surface-water groundwater interaction model is demonstrated with an example.

Transport dissipative particle dynamics model for mesoscopic advection-diffusion-reaction problems

PubMed Central

Yazdani, Alireza; Tartakovsky, Alexandre; Karniadakis, George Em

2015-01-01

We present a transport dissipative particle dynamics (tDPD) model for simulating mesoscopic problems involving advection-diffusion-reaction (ADR) processes, along with a methodology for implementation of the correct Dirichlet and Neumann boundary conditions in tDPD simulations. tDPD is an extension of the classic dissipative particle dynamics (DPD) framework with extra variables for describing the evolution of concentration fields. The transport of concentration is modeled by a Fickian flux and a random flux between tDPD particles, and the advection is implicitly considered by the movements of these Lagrangian particles. An analytical formula is proposed to relate the tDPD parameters to the effective diffusion coefficient. To validate the present tDPD model and the boundary conditions, we perform three tDPD simulations of one-dimensional diffusion with different boundary conditions, and the results show excellent agreement with the theoretical solutions. We also performed two-dimensional simulations of ADR systems and the tDPD simulations agree well with the results obtained by the spectral element method. Finally, we present an application of the tDPD model to the dynamic process of blood coagulation involving 25 reacting species in order to demonstrate the potential of tDPD in simulating biological dynamics at the mesoscale. We find that the tDPD solution of this comprehensive 25-species coagulation model is only twice as computationally expensive as the conventional DPD simulation of the hydrodynamics only, which is a significant advantage over available continuum solvers. PMID:26156459

Correlation-based Transition Modeling for External Aerodynamic Flows

NASA Astrophysics Data System (ADS)

Medida, Shivaji

Conventional turbulence models calibrated for fully turbulent boundary layers often over-predict drag and heat transfer on aerodynamic surfaces with partially laminar boundary layers. A robust correlation-based model is developed for use in Reynolds-Averaged Navier-Stokes simulations to predict laminar-to-turbulent transition onset of boundary layers on external aerodynamic surfaces. The new model is derived from an existing transition model for the two-equation k-omega Shear Stress Transport (SST) turbulence model, and is coupled with the one-equation Spalart-Allmaras (SA) turbulence model. The transition model solves two transport equations for intermittency and transition momentum thickness Reynolds number. Experimental correlations and local mean flow quantities are used in the model to account for effects of freestream turbulence level and pressure gradients on transition onset location. Transition onset is triggered by activating intermittency production using a vorticity Reynolds number criterion. In the new model, production and destruction terms of the intermittency equation are modified to improve consistency in the fully turbulent boundary layer post-transition onset, as well as ensure insensitivity to freestream eddy viscosity value specified in the SA model. In the original model, intermittency was used to control production and destruction of turbulent kinetic energy. Whereas, in the new model, only the production of eddy viscosity in SA model is controlled, and the destruction term is not altered. Unlike the original model, the new model does not use an additional correction to intermittency for separation-induced transition. Accuracy of drag predictions are improved significantly with the use of the transition model for several two-dimensional single- and multi-element airfoil cases over a wide range of Reynolds numbers. The new model is able to predict the formation of stable and long laminar separation bubbles on low-Reynolds number airfoils that is not captured with conventional turbulence models. The validated transition model is successfully applied to rotating blade configurations in axial flow conditions to study the effects of transitional boundary layers on rotor thrust and torque. In helicopter rotors, inclusion of transition effects increased thrust prediction by 2% and decreased torque by as much as 8% at lower collective angles, due to reduced airfoil profile drag. In wind turbine rotors, transition model predicted a 7%--70% increase in generated shaft torque at lower wind speeds, due to lower viscous drag. This has important implications for CFD analysis of small wind turbines operating at low values of rated power. Transition onset locations along upper and lower surfaces of rotor blades are analyzed in detail. A new crossflow transition onset criterion is developed to account for crossflow instability effects in three-dimensional boundary layers. Preliminary results for swept wing and rotating blade flows demonstrate the need to account for crossflow transition in three-dimensional simulations of wings, rotating blades, and airframes. Inclusion of crossflow effects resulted in accelerated transition in the presence of favorable pressure gradients and yawed flow. Finally, a new correction to the wall damping function in the Spalart-Allmaras turbulence model is proposed to improve sensitivity of the model to strong adverse pressure gradients (APG). The correction reduces turbulence production in the boundary layer when the ratio of magnitudes of local turbulent stress to the wall shear stress exceeds a threshold value, therefore enabling earlier separation of boundary layer. Improved prediction of static and dynamic stall on two-dimensional airfoils is demonstrated with the APG correction.

Time-dependent boundary conditions for hyperbolic systems. II

NASA Technical Reports Server (NTRS)

Thompson, Kevin W.

1990-01-01

A general boundary condition formalism is developed for all types of boundary conditions to which hyperbolic systems are subject; the formalism makes possible a 'cookbook' approach to boundary conditions, by means of which novel boundary 'recipes' may be derived and previously devised ones may be consulted as required. Numerous useful conditions are derived for such CFD problems as subsonic and supersonic inflows and outflows, nonreflecting boundaries, force-free boundaries, constant pressure boundaries, and constant mass flux. Attention is given to the computation and integration of time derivatives.

Time-dependent boundary conditions for hyperbolic systems. II

NASA Astrophysics Data System (ADS)

Thompson, Kevin W.

1990-08-01

A general boundary condition formalism is developed for all types of boundary conditions to which hyperbolic systems are subject; the formalism makes possible a 'cookbook' approach to boundary conditions, by means of which novel boundary 'recipes' may be derived and previously devised ones may be consulted as required. Numerous useful conditions are derived for such CFD problems as subsonic and supersonic inflows and outflows, nonreflecting boundaries, force-free boundaries, constant pressure boundaries, and constant mass flux. Attention is given to the computation and integration of time derivatives.

Tectonic predictions with mantle convection models

NASA Astrophysics Data System (ADS)

Coltice, Nicolas; Shephard, Grace E.

2018-04-01

Over the past 15 yr, numerical models of convection in Earth's mantle have made a leap forward: they can now produce self-consistent plate-like behaviour at the surface together with deep mantle circulation. These digital tools provide a new window into the intimate connections between plate tectonics and mantle dynamics, and can therefore be used for tectonic predictions, in principle. This contribution explores this assumption. First, initial conditions at 30, 20, 10 and 0 Ma are generated by driving a convective flow with imposed plate velocities at the surface. We then compute instantaneous mantle flows in response to the guessed temperature fields without imposing any boundary conditions. Plate boundaries self-consistently emerge at correct locations with respect to reconstructions, except for small plates close to subduction zones. As already observed for other types of instantaneous flow calculations, the structure of the top boundary layer and upper-mantle slab is the dominant character that leads to accurate predictions of surface velocities. Perturbations of the rheological parameters have little impact on the resulting surface velocities. We then compute fully dynamic model evolution from 30 and 10 to 0 Ma, without imposing plate boundaries or plate velocities. Contrary to instantaneous calculations, errors in kinematic predictions are substantial, although the plate layout and kinematics in several areas remain consistent with the expectations for the Earth. For these calculations, varying the rheological parameters makes a difference for plate boundary evolution. Also, identified errors in initial conditions contribute to first-order kinematic errors. This experiment shows that the tectonic predictions of dynamic models over 10 My are highly sensitive to uncertainties of rheological parameters and initial temperature field in comparison to instantaneous flow calculations. Indeed, the initial conditions and the rheological parameters can be good enough for an accurate prediction of instantaneous flow, but not for a prediction after 10 My of evolution. Therefore, inverse methods (sequential or data assimilation methods) using short-term fully dynamic evolution that predict surface kinematics are promising tools for a better understanding of the state of the Earth's mantle.

Correction of Excessive Precipitation Over Steep and High Mountains in a General Circulation Model

NASA Technical Reports Server (NTRS)

Chao, Winston C.

2012-01-01

Excessive precipitation over steep and high mountains (EPSM) is a well-known problem in GCMs and meso-scale models. This problem impairs simulation and data assimilation products. Among the possible causes investigated in this study, we found that the most important one, by far, is a missing upward transport of heat out of the boundary layer due to the vertical circulations forced by the daytime upslope winds, which are forced by the heated boundary layer on subgrid-scale slopes. These upslope winds are associated with large subgrid-scale topographic variation, which is found over steep and high mountains. Without such subgridscale heat ventilation, the resolvable-scale upslope flow in the boundary layer generated by surface sensible heat flux along the mountain slopes is excessive. Such an excessive resolvablescale upslope flow combined with the high moisture content in the boundary layer results in excessive moisture transport toward mountaintops, which in turn gives rise to EPSM. Other possible causes of EPSM that we have investigated include 1) a poorly-designed horizontal moisture flux in the terrain-following coordinates, 2) the condition for cumulus convection being too easily satisfied at mountaintops, 3) the presence of conditional instability of the computational kind, and 4) the absence of blocked flow drag. These are all minor or inconsequential. We have parameterized the ventilation effects of the subgrid-scale heated-slope-induced vertical circulation (SHVC) by removing heat from the boundary layer and depositing it in layers higher up when the topographic variance exceeds a critical value. Test results using NASA/Goddard's GEOS-S GCM have shown that this largely solved the EPSM problem.

A closed-form analytical model for predicting 3D boundary layer displacement thickness for the validation of viscous flow solvers

NASA Astrophysics Data System (ADS)

Kumar, V. R. Sanal; Sankar, Vigneshwaran; Chandrasekaran, Nichith; Saravanan, Vignesh; Natarajan, Vishnu; Padmanabhan, Sathyan; Sukumaran, Ajith; Mani, Sivabalan; Rameshkumar, Tharikaa; Nagaraju Doddi, Hema Sai; Vysaprasad, Krithika; Sharan, Sharad; Murugesh, Pavithra; Shankar, S. Ganesh; Nejaamtheen, Mohammed Niyasdeen; Baskaran, Roshan Vignesh; Rahman Mohamed Rafic, Sulthan Ariff; Harisrinivasan, Ukeshkumar; Srinivasan, Vivek

2018-02-01

A closed-form analytical model is developed for estimating the 3D boundary-layer-displacement thickness of an internal flow system at the Sanal flow choking condition for adiabatic flows obeying the physics of compressible viscous fluids. At this unique condition the boundary-layer blockage induced fluid-throat choking and the adiabatic wall-friction persuaded flow choking occur at a single sonic-fluid-throat location. The beauty and novelty of this model is that without missing the flow physics we could predict the exact boundary-layer blockage of both 2D and 3D cases at the sonic-fluid-throat from the known values of the inlet Mach number, the adiabatic index of the gas and the inlet port diameter of the internal flow system. We found that the 3D blockage factor is 47.33 % lower than the 2D blockage factor with air as the working fluid. We concluded that the exact prediction of the boundary-layer-displacement thickness at the sonic-fluid-throat provides a means to correctly pinpoint the causes of errors of the viscous flow solvers. The methodology presented herein with state-of-the-art will play pivotal roles in future physical and biological sciences for a credible verification, calibration and validation of various viscous flow solvers for high-fidelity 2D/3D numerical simulations of real-world flows. Furthermore, our closed-form analytical model will be useful for the solid and hybrid rocket designers for the grain-port-geometry optimization of new generation single-stage-to-orbit dual-thrust-motors with the highest promising propellant loading density within the given envelope without manifestation of the Sanal flow choking leading to possible shock waves causing catastrophic failures.

General Boundary Conditions for a Majorana Single-Particle in a Box in (1 + 1) Dimensions

NASA Astrophysics Data System (ADS)

De Vincenzo, Salvatore; Sánchez, Carlet

2018-05-01

We consider the problem of a Majorana single-particle in a box in (1 + 1) dimensions. We show that the most general set of boundary conditions for the equation that models this particle is composed of two families of boundary conditions, each one with a real parameter. Within this set, we only have four confining boundary conditions—but infinite not confining boundary conditions. Our results are also valid when we include a Lorentz scalar potential in this equation. No other Lorentz potential can be added. We also show that the four confining boundary conditions for the Majorana particle are precisely the four boundary conditions that mathematically can arise from the general linear boundary condition used in the MIT bag model. Certainly, the four boundary conditions for the Majorana particle are also subject to the Majorana condition.

The effects of aging on the speed-accuracy compromise: Boundary optimality in the diffusion model.

PubMed

Starns, Jeffrey J; Ratcliff, Roger

2010-06-01

We evaluated age-related differences in the optimality of decision boundary settings in a diffusion model analysis. In the model, the width of the decision boundary represents the amount of evidence that must accumulate in favor of a response alternative before a decision is made. Wide boundaries lead to slow but accurate responding, and narrow boundaries lead to fast but inaccurate responding. There is a single value of boundary separation that produces the most correct answers in a given period of time, and we refer to this value as the reward rate optimal boundary (RROB). We consistently found across a variety of decision tasks that older adults used boundaries that were much wider than the RROB value. Young adults used boundaries that were closer to the RROB value, although age differences in optimality were smaller with instructions emphasizing speed than with instructions emphasizing accuracy. Young adults adjusted their boundary settings to more closely approach the RROB value when they were provided with accuracy feedback and extensive practice. Older participants showed no evidence of making boundary adjustments in response to feedback or task practice, and they consistently used boundary separation values that produced accuracy levels that were near asymptote. Our results suggest that young adults attempt to balance speed and accuracy to achieve the most correct answers per unit time, whereas older adultts attempt to minimize errors even if they must respond quite slowly to do so. (c) 2010 APA, all rights reserved

A Reconsideration of the Extension Strain Criterion for Fracture and Failure of Rock

NASA Astrophysics Data System (ADS)

Wesseloo, J.; Stacey, T. R.

2016-12-01

The complex behaviours of rocks and rock masses have presented paradoxes to the rock engineer, including the fracturing of seemingly strong rock under low stress conditions, which often occurs near excavation boundaries. The extension strain criterion was presented as a fracture initiation criterion under these conditions (Stacey in Int J Rock Mech Min Sci 18:469-474, 1981). This criterion has been used successfully by some and criticised by others. In this paper, we review the literature on the extension strain criterion and present a case for the correct interpretation of the criterion and the conditions suitable for its use. We argue that the extension strain criterion can also be used to provide an indication of damage level under conditions of relatively low confining stress. We also present an augmentation of the criterion, the ultimate extension strain, which is applicable under extensional loading conditions when σ 2 is similar in magnitude to σ 1.

A simplified fourwall interference assessment procedure for airfoil data obtained in the Langley 0.3-meter transonic cryogenic tunnel

NASA Technical Reports Server (NTRS)

Murthy, A. V.

1987-01-01

A simplified fourwall interference assessment method has been described, and a computer program developed to facilitate correction of the airfoil data obtained in the Langley 0.3-m Transonic Cryogenic Tunnel (TCT). The procedure adopted is to first apply a blockage correction due to sidewall boundary-layer effects by various methods. The sidewall boundary-layer corrected data are then used to calculate the top and bottom wall interference effects by the method of Capallier, Chevallier and Bouinol, using the measured wall pressure distribution and the model force coefficients. The interference corrections obtained by the present method have been compared with other methods and found to give good agreement for the experimental data obtained in the TCT with slotted top and bottom walls.

Mesh refinement in a two-dimensional large eddy simulation of a forced shear layer

NASA Technical Reports Server (NTRS)

Claus, R. W.; Huang, P. G.; Macinnes, J. M.

1989-01-01

A series of large eddy simulations are made of a forced shear layer and compared with experimental data. Several mesh densities were examined to separate the effect of numerical inaccuracy from modeling deficiencies. The turbulence model that was used to represent small scale, 3-D motions correctly predicted some gross features of the flow field, but appears to be structurally incorrect. The main effect of mesh refinement was to act as a filter on the scale of vortices that developed from the inflow boundary conditions.

Damage localization and quantification of composite stratified beam Structures using residual force method

NASA Astrophysics Data System (ADS)

Behtani, A.; Bouazzouni, A.; Khatir, S.; Tiachacht, S.; Zhou, Y.-L.; Abdel Wahab, M.

2017-05-01

In this paper, the problem of using measured modal parameters to detect and locate damage in beam composite stratified structures with four layers of graphite/epoxy [0°/902°/0°] is investigated. A technique based on the residual force method is applied to composite stratified structure with different boundary conditions, the results of damage detection for several damage cases demonstrate that using residual force method as damage index, the damage location can be identified correctly and the damage extents can be estimated as well.

Computational studies of an impulsively started viscous flow

NASA Technical Reports Server (NTRS)

Davis, Sanford S.

1988-01-01

Progress in validating incompressible Navier-Stokes codes is described using a predictor/corrector scheme. The flow field under study is the impulsive start of a circular cylinder and the unsteady evolution of the separation bubble. In the current code, a uniform asymptotic expansion is used as an initial condition in order to correctly capture the initial growth of the vortex sheet. Volocity fields at selected instants of time are decomposed into vectorial representations of Navier-Stokes equations which are then used to analyze dominant contributions in the boundary-layer region.

Influence of boundary conditions on the hydrodynamic forces of an oscillating sphere

NASA Astrophysics Data System (ADS)

Mirauda, Domenica; Negri, Marco; Martinelli, Luca; Malavasi, Stefano

2018-06-01

The design of submerged structures in sea currents presents certain problems that are not only connected to the shape of the obstacle but also to the number of acting forces as well as the correct modelling of the structures dynamic response. Currently, the common approach is that of integrated numerical modelling, which considers the contribution of both current and structure. The reliability of such an approach is better verified with experimental tests performed on models of simple geometry. On the basis of these considerations, the present work analyses the hydrodynamic forces acting on a sphere, which is characterised by a low mass ratio and damping. The sphere is immersed in a free surface flow and can oscillate along the streamwise and transverse flow direction. It is located at three different positions inside the current: close to the channel bottom, near the free surface and in the middle, and equally distant from both the bottom and free surface. The obtained results for different boundaries and flow kinematic conditions show a relevant influence of the free surface on the hydrodynamic forces along both the streamwise and transverse flow directions.

A General Computational Approach for Magnetohydrodynamic Flows Using the CFX Code: Buoyant Flow Through a Vertical Square Channel

DOE Office of Scientific and Technical Information (OSTI.GOV)

Di Piazza, Ivan; Buehler, Leo

2000-09-15

The buoyancy-driven magnetoconvection in the cross section of an infinitely long vertical square duct is investigated numerically using the CFX code package. The implementation of a magnetohydrodynamic (MHD) problem in CFX is discussed, with particular reference to the Lorentz forces and the electric potential boundary conditions for arbitrary electrical conductivity of the walls. The method proposed is general and applies to arbitrary geometries with an arbitrary orientation of the magnetic field. Results for fully developed flow under various thermal boundary conditions are compared with asymptotic analytical solutions. The comparison shows that the asymptotic analysis is confirmed for highly conducting wallsmore » as high velocity jets occur at the side walls. For weakly conducting walls, the side layers become more conducting than the side walls, and strong electric currents flow within these layers parallel to the magnetic field. As a consequence, the velocity jets are suppressed, and the core solution is only corrected by the viscous forces near the wall. The implementation of MHD in CFX is achieved.« less

A finite element model development for simulation of the impact of slab thickness, joints, and membranes on indoor radon concentration.

PubMed

Muñoz, E; Frutos, B; Olaya, M; Sánchez, J

2017-10-01

The focus of this study is broadly to define the physics involved in radon generation and transport through the soil and other materials using different parameter-estimation tools from the literature. The effect of moisture in the soil and radon transport via water in the pore space was accounted for with the application of a porosity correction coefficient. A 2D finite element model is created, which reproduces the diffusion and advection mechanisms resulting from specified boundary conditions. A comparison between the model and several analytical and numerical solutions obtained from the literature and field studies validates the model. Finally, the results demonstrate that the model can predict radon entry through different building boundary conditions, such as concrete slabs with or without joints, variable slab thicknesses and diffusion coefficients, and the use of several radon barrier membranes. Cracks in the concrete or the radon barrier membrane have been studied to understand how indoor concentration is affected by these issues. Copyright © 2017 Elsevier Ltd. All rights reserved.

Impact of tides in a baroclinic circulation model of the Adriatic Sea

NASA Astrophysics Data System (ADS)

Guarnieri, A.; Pinardi, N.; Oddo, P.; Bortoluzzi, G.; Ravaioli, M.

2013-01-01