Numerical Optimization Using Computer Experiments
NASA Technical Reports Server (NTRS)
Trosset, Michael W.; Torczon, Virginia
1997-01-01
Engineering design optimization often gives rise to problems in which expensive objective functions are minimized by derivative-free methods. We propose a method for solving such problems that synthesizes ideas from the numerical optimization and computer experiment literatures. Our approach relies on kriging known function values to construct a sequence of surrogate models of the objective function that are used to guide a grid search for a minimizer. Results from numerical experiments on a standard test problem are presented.
Comparison of integrated numerical experiments with accelerator and FEL experiments
Thode, L.E.; Carlsten, B.E.; Chan, K.C.D.; Cooper, R.K.; Elliott, J.C.; Gitomer, S.J.; Goldstein, J.C.; Jones, M.E.; McVey, B.D.; Schmitt, M.J.; Takeda, H.; Tokar, R.L.; Wang, T.S.; Young, L.M.
1991-01-01
Even at the conceptual level the strong coupling between the laser subsystem elements, such as the accelerator, wiggler, optics, and control, greatly complicates the understanding and design of an FEL. Given the requirements for a high-performance FEL, the coupling between the laser subsystems must be included in the design approach. To address the subsystem coupling the concept of an integrated numerical experiment (INEX) has been implemented. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostic. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. A complete INEX model has been applied to the 10{mu}m high-extraction-efficiency experiment at Los Alamos and the 0.6-{mu}m Burst Mode experiment at Boeing Aerospace. In addition, various subsets of the INEX model have been compared with a number of other experiments. Overall, the agreement between INEX and the experiments is very good. With the INEX approach, it now appears possible to design high-performance FELS for numerous applications. The first full-scale test of the INEX approach is the Los Alamos HIBAF experiment. The INEX concept, implementation, and validation with experiments are discussed. 28 refs., 13 figs., 1 tab.
Numerical experiments in homogeneous turbulence
NASA Technical Reports Server (NTRS)
Rogallo, R. S.
1981-01-01
The direct simulation methods developed by Orszag and Patternson (1972) for isotropic turbulence were extended to homogeneous turbulence in an incompressible fluid subjected to uniform deformation or rotation. The results of simulations for irrotational strain (plane and axisymmetric), shear, rotation, and relaxation toward isotropy following axisymmetric strain are compared with linear theory and experimental data. Emphasis is placed on the shear flow because of its importance and because of the availability of accurate and detailed experimental data. The computed results are used to assess the accuracy of two popular models used in the closure of the Reynolds-stress equations. Data from a variety of the computed fields and the details of the numerical methods used in the simulation are also presented.
Numerical experiments on the theta pinch
NASA Technical Reports Server (NTRS)
Volosevich, P. P.; Zukakishyili, G. G.
1979-01-01
Numerical calculation of theta pinch problems are presented. Physical processes in theta pinch systems are considered in a one dimensional, two temperature magnetohydrodynamic, approximation with allowance for end losses by longitudinal heat conductivity. The numerical calculations are compared with results of earlier experiments.
French and English Together: An "Additive" Experience
ERIC Educational Resources Information Center
Wiltshire, Jessica; Harbon, Lesley
2010-01-01
This paper examines the nature of the "additive" experience of a bilingual French-English curriculum at Killarney Heights Public School in New South Wales. Predictably, the well-supported "additive" nature of the languages program model elicited positive reactions regarding educational success. The paper also explores issues for administration,…
Numerical Simulations of Falling Sphere Viscometry Experiments.
NASA Astrophysics Data System (ADS)
O Dwyer, L.; Kellogg, L. H.; Lesher, C. E.
2007-12-01
The falling sphere technique based on Stokes' law is widely used to determine the viscosities of geologically relevant melts at high pressures. Stokes' law is valid when a rigid sphere falls slowly and steadily through a stationary and infinite Newtonian medium of uniform properties. High-pressure falling sphere experiments however, usually involve dropping a dense, refractory sphere through a liquid contained by a cylindrical capsule of finite size. The sphere velocity is influenced by the walls (Faxen correction) and ends of the capsule, and possible convective motion of the fluid. Efforts are made to minimize thermal gradients in laboratory experiments, but small temperature differences within the capsule can lead to convection complicating interpretation. We utilize GALE (Moresi et al., 2003;), a finite element particle-in-cell code, to examine these factors in numerical models of conditions similar to those of high-pressure experiments. Our modeling considers a three- dimensional box or cylinder containing a cluster of particles that represent the dense sphere in laboratory experiments surrounded by low viscosity particles representing the melt. GALE includes buoyancy forces, heat flow, and viscosity variations so our model can be used to assess the effects of the capsule's walls and ends, and the consequences of thermal gradients on the sphere's velocity and trajectory. Comparisons between our numerical simulations and real-time falling sphere experiments involving lower viscosity molten komatiite are made to assess the validity of Stokes' law with the standard Faxen correction included, and formulations considering end effects. The modeling also permits an evaluation of the uncertainties in recovering accurate liquid viscosities from Stokes' law when a dense sphere falls through a convecting low viscosity melt. It also allows us to assess acceleration to a terminal velocity that can provide constraints on melt viscosity in experiments in which the terminal
Submarine sand volcanos: experiments and numerical modelling
NASA Astrophysics Data System (ADS)
Philippe, P.; Ngoma, J.; Delenne, J.
2012-12-01
Fluid overpressure at the bottom of a soil layer may generate fracturation in preferential paths for a cohesive material. But the case of sandy soils is rather different: a significant internal flow is allowed within the material and can potentially induce hydro-mechanical instabilities whose most common example is fluidization. Many works have been devoted to fluidization but very few have the issue of initiation and development of a fluidized zone inside a granular bed, prior entire fluidization of the medium. In this contribution, we report experimental results and numerical simulations on a model system of immersed sand volcanos generated by a localized upward spring of liquid, injected at constant flow-rate at the bottom of a granular layer. Such a localized state of fluidization is relevant for some industrial processes (spouted bed, maintenance of navigable waterways,…) and for several geological issues (kimberlite volcano conduits, fluid venting, oil recovery in sandy soil, More precisely, what is presented here is a comparison between experiments, carried out by direct visualization throughout the medium, and numerical simulations, based on DEM modelling of the grains coupled to resolution of NS equations in the liquid phase (LBM). There is a very good agreement between the experimental phenomenology and the simulation results. When the flow-rate is increased, three regimes are successively observed: static bed, fluidized cavity that does not extend to the top of the layer, and finally fluidization over the entire height of layer that creates a fluidized chimney. A very strong hysteretic effect is present here with an extended range of stability for fluidized cavities when flow-rate is decreased back. This can be interpreted in terms force chains and arches. The influences of grain diameter, layer height and injection width are studied and interpreted using a model previously developed by Zoueshtiagh [1]. Finally, growing rate of the fluidized zone and
Seafloor weathering buffering climate: numerical experiments
NASA Astrophysics Data System (ADS)
Farahat, N. X.; Archer, D. E.; Abbot, D. S.
2013-12-01
Continental silicate weathering is widely held to consume atmospheric CO2 at a rate controlled in part by temperature, resulting in a climate-weathering feedback [Walker et al., 1981]. It has been suggested that weathering of oceanic crust of warm mid-ocean ridge flanks also has a CO2 uptake rate that is controlled by climate [Sleep and Zahnle, 2001; Brady and Gislason, 1997]. Although this effect might not be significant on present-day Earth [Caldeira, 1995], seafloor weathering may be more pronounced during snowball states [Le Hir et al., 2008], during the Archean when seafloor spreading rates were faster [Sleep and Zahnle, 2001], and on waterworld planets [Abbot et al., 2012]. Previous studies of seafloor weathering have made significant contributions using qualitative, generally one-box, models, and the logical next step is to extend this work using a spatially resolved model. For example, experiments demonstrate that seafloor weathering reactions are temperature dependent, but it is not clear whether the deep ocean temperature affects the temperature at which the reactions occur, or if instead this temperature is set only by geothermal processes. Our goal is to develop a 2-D numerical model that can simulate hydrothermal circulation and resulting alteration of oceanic basalts, and can therefore address such questions. A model of diffusive and convective heat transfer in fluid-saturated porous media simulates hydrothermal circulation through porous oceanic basalt. Unsteady natural convection is solved for using a Darcy model of porous media flow that has been extensively benchmarked. Background hydrothermal circulation is coupled to mineral reaction kinetics of basaltic alteration and hydrothermal mineral precipitation. In order to quantify seafloor weathering as a climate-weathering feedback process, this model focuses on hydrothermal reactions that influence carbon uptake as well as ocean alkalinity: silicate rock dissolution, calcium and magnesium leaching
Processing biobased polymers using plasticizers: Numerical simulations versus experiments
NASA Astrophysics Data System (ADS)
Desplentere, Frederik; Cardon, Ludwig; Six, Wim; Erkoç, Mustafa
2016-03-01
In polymer processing, the use of biobased products shows lots of possibilities. Considering biobased materials, biodegradability is in most cases the most important issue. Next to this, bio based materials aimed at durable applications, are gaining interest. Within this research, the influence of plasticizers on the processing of the bio based material is investigated. This work is done for an extrusion grade of PLA, Natureworks PLA 2003D. Extrusion through a slit die equipped with pressure sensors is used to compare the experimental pressure values to numerical simulation results. Additional experimental data (temperature and pressure data along the extrusion screw and die are recorded) is generated on a dr. Collin Lab extruder producing a 25mm diameter tube. All these experimental data is used to indicate the appropriate functioning of the numerical simulation tool Virtual Extrusion Laboratory 6.7 for the simulation of both the industrial available extrusion grade PLA and the compound in which 15% of plasticizer is added. Adding the applied plasticizer, resulted in a 40% lower pressure drop over the extrusion die. The combination of different experiments allowed to fit the numerical simulation results closely to the experimental values. Based on this experience, it is shown that numerical simulations also can be used for modified bio based materials if appropriate material and process data are taken into account.
Spray combustion experiments and numerical predictions
NASA Technical Reports Server (NTRS)
Mularz, Edward J.; Bulzan, Daniel L.; Chen, Kuo-Huey
1993-01-01
The next generation of commercial aircraft will include turbofan engines with performance significantly better than those in the current fleet. Control of particulate and gaseous emissions will also be an integral part of the engine design criteria. These performance and emission requirements present a technical challenge for the combustor: control of the fuel and air mixing and control of the local stoichiometry will have to be maintained much more rigorously than with combustors in current production. A better understanding of the flow physics of liquid fuel spray combustion is necessary. This paper describes recent experiments on spray combustion where detailed measurements of the spray characteristics were made, including local drop-size distributions and velocities. Also, an advanced combustor CFD code has been under development and predictions from this code are compared with experimental results. Studies such as these will provide information to the advanced combustor designer on fuel spray quality and mixing effectiveness. Validation of new fast, robust, and efficient CFD codes will also enable the combustor designer to use them as additional design tools for optimization of combustor concepts for the next generation of aircraft engines.
Numerical experiments on quantum chaotic billiards
NASA Astrophysics Data System (ADS)
de Menezes, D. D.; Jar e Silva, M.; de Aguiar, F. M.
2007-06-01
A recently proposed numerical technique for generation of high-quality unstructured meshes is combined with a finite-element method to solve the Helmholtz equation that describes the quantum mechanics of a particle confined in two-dimensional cavities. Different shapes are treated on equal footing, including Sinai, stadium, annular, threefold symmetric, mushroom, cardioid, triangle, and coupled billiards. The results are shown to be in excellent agreement with available measurements in flat microwave resonator counterparts with nonintegrable geometries.
Numerical experiments on quantum chaotic billiards.
de Menezes, D D; Jar e Silva, M; de Aguiar, F M
2007-06-01
A recently proposed numerical technique for generation of high-quality unstructured meshes is combined with a finite-element method to solve the Helmholtz equation that describes the quantum mechanics of a particle confined in two-dimensional cavities. Different shapes are treated on equal footing, including Sinai, stadium, annular, threefold symmetric, mushroom, cardioid, triangle, and coupled billiards. The results are shown to be in excellent agreement with available measurements in flat microwave resonator counterparts with nonintegrable geometries. PMID:17614670
Numerical experiments on unstructured PIC stability.
Day, David Minot
2011-04-01
Particle-In-Cell (PIC) is a method for plasmas simulation. Particles are pushed with Verlet time integration. Fields are modeled using finite differences on a tensor product mesh (cells). The Unstructured PIC methods studied here use instead finite element discretizations on unstructured (simplicial) meshes. PIC is constrained by stability limits (upper bounds) on mesh and time step sizes. Numerical evidence (2D) and analysis will be presented showing that similar bounds constrain unstructured PIC.
Numerical experiments with flows of elongated granules
NASA Technical Reports Server (NTRS)
Elrod, Harold G.; Brewe, David E.
1992-01-01
Theory and numerical results are given for a program simulating two dimensional granular flow (1) between two infinite, counter-moving, parallel, roughened walls, and (2) for an infinitely wide slider. Each granule is simulated by a central repulsive force field ratcheted with force restitution factor to introduce dissipation. Transmission of angular momentum between particles occurs via Coulomb friction. The effect of granular hardness is explored. Gaps from 7 to 28 particle diameters are investigated, with solid fractions ranging from 0.2 to 0.9. Among features observed are: slip flow at boundaries, coagulation at high densities, and gross fluctuation in surface stress. A videotape has been prepared to demonstrate the foregoing effects.
Numerical Investigation of Plasma Detachment in Magnetic Nozzle Experiments
NASA Technical Reports Server (NTRS)
Sankaran, Kamesh; Polzin, Kurt A.
2008-01-01
At present there exists no generally accepted theoretical model that provides a consistent physical explanation of plasma detachment from an externally-imposed magnetic nozzle. To make progress towards that end, simulation of plasma flow in the magnetic nozzle of an arcjet experiment is performed using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The aim is to compare the computational results with various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver (see Fig. I), and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. This geometry is modeled because there is a substantial amount of experimental data that can be compared to the computational results, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters. Further studies will focus on modeling literature experiments both for the purpose of additional code validation and to extract physical insight regarding the mechanisms driving detachment.
Surface nanobubbles: Theory, numerics and experiments
NASA Astrophysics Data System (ADS)
Weijs, Joost H.
2013-11-01
When a solid is brought into contact with water, surface nanobubbles can be formed at the solid-liquid interface. These nanobubbles are small; their height is of order 10nm and their lateral sizes vary from 10-100 nm. Initially, the only proof of the existence of surface nanobubbles was delivered by atomic force microscopy. Later, additional techniques such as infrared attenuated total reflectance have confirmed the existence of gaseous domains on the solid-liquid interface. Before this overwhelming evidence, the existence of surface nanobubbles was controversial because they possess some unusual properties. For example, nanobubbles are surprisingly robust against dissolution by diffusion and Laplace pressure: Instead of the expected lifetime of about a microsecond, nanobubbles are found to survive for several hours and in some cases even several days. Additionally, surface nanobubbles are flatter than predicted by Young's law and are able to resist strong tensile stresses (~-6 MPa), rather than serving as a nucleation site for a macroscopic bubble. A deep understanding of surface nanobubbles is crucial for practical applications (e.g. drag reduction in microfluidic devices) but nanobubbles also pose fundamental questions on the validity of continuum models at the nanoscale. In this talk, we will discuss these open questions in detail by considering theoretical efforts and molecular dynamics simulations. Theoretically, we study the consequences of a pinned contact line. We find that the pinned contact line can explain the long lifetimes and many other nanobubble properties. From molecular dynamics results, we clarify the influence of the gas species on the contact angle. Finally, we will discuss some very recent experimental and theoretical work on the effects of an acoustic field on nanobubbles. We provide experimental data combined with a theoretical analysis and find that the acoustic driving can cause the nanobubbles to grow by rectified diffusion.
Numerical modeling of Deep Impact experiment
NASA Astrophysics Data System (ADS)
Sultanov, V. G.; Kim, V. V.; Lomonosov, I. V.; Shutov, A. V.; Fortov, V. E.
2007-06-01
The Deep Impact active space experiment has been done [1,2] to study a hypervelocity collision of a metal impactor with the comet 9P/Temple 1. The modeling of impact on solid or porous ice made it possible to conclude: the form and size of crater depends strongly on the density of comet material; the copper impactor does not melt and remains in the solid state; the temperature of ejecta varies from 5000 K for solid ice to 15000 K for porous ice. The impact on moist water- saturated sand demonstrated different results. In this case, the copper impactor practically does not penetrate the comet surface, melts, destroys and the ricochet process takes place. In the case of moist porous sand the produced crater is stretched in the direction of impact. The analysis of modeling results indicates to the presence of volatile easy-vaporized chemical compounds in the cometary surface. The hypothesis that the cometary surface consists of only ice does not agree with experimental and computational data on the forming and spreading of impact ejecta. [1] http://deepimpact.jpl.nasa.gov/home/index.html [2] M. F. A'Hearn et al, Deep Impact: Excavating Comet Tempel 1 // Science, 2005, v.310, pp. 258-264
The Fermi Pasta Ulam 'numerical experiment': history and pedagogical perspectives
NASA Astrophysics Data System (ADS)
Dauxois, Thierry; Peyrard, Michel; Ruffo, Stefano
2005-09-01
The pioneering Fermi-Pasta-Ulam (FPU) numerical experiment played a major role in the history of computer simulation because it introduced this concept for the first time. Moreover, it raised a puzzling question which was answered more than 10 years later. After an introduction to this problem, we briefly review its history and then suggest some simple numerical experiments, with the Matlab© code provided, to study various aspects of the 'FPU' problem.
Fire extinct experiments with water mist by adding additives
NASA Astrophysics Data System (ADS)
Yang, Lijun; Zhao, Jianbo
2011-12-01
The effects of fire extinguishment with water mist by adding different additives were studied. Tens of chemical substances (including alkali metal salt, dilution agent and surface active agent) were selected as additives due to their different extinct mechanisms. At first the performance of fire extinguishment with single additive was studied, then the effects of the same kinds of chemical substances under the same mass fraction were compared to study their influences on the fire extinguishment factors, including extinct time, fire temperature and oxygen concentration from which the fire extinct mechanism with additives could be concluded. Based on this the experiments were conducted to study the cooperate effect of the complexity of different additives. It indicated the relations between different firefighting mechanisms and different additives were competitive. From a large number of experiments the extinct mechanism with water mist by adding additives was concluded and an optimal compounding additive was selected.
Non-robust numerical simulations of analogue extension experiments
NASA Astrophysics Data System (ADS)
Naliboff, John; Buiter, Susanne
2016-04-01
Numerical and analogue models of lithospheric deformation provide significant insight into the tectonic processes that lead to specific structural and geophysical observations. As these two types of models contain distinct assumptions and tradeoffs, investigations drawing conclusions from both can reveal robust links between first-order processes and observations. Recent studies have focused on detailed comparisons between numerical and analogue experiments in both compressional and extensional tectonics, sometimes involving multiple lithospheric deformation codes and analogue setups. While such comparisons often show good agreement on first-order deformation styles, results frequently diverge on second-order structures, such as shear zone dip angles or spacing, and in certain cases even on first-order structures. Here, we present finite-element experiments that are designed to directly reproduce analogue "sandbox" extension experiments at the cm-scale. We use material properties and boundary conditions that are directly taken from analogue experiments and use a Drucker-Prager failure model to simulate shear zone formation in sand. We find that our numerical experiments are highly sensitive to numerous numerical parameters. For example, changes to the numerical resolution, velocity convergence parameters and elemental viscosity averaging commonly produce significant changes in first- and second-order structures accommodating deformation. The sensitivity of the numerical simulations to small parameter changes likely reflects a number of factors, including, but not limited to, high angles of internal friction assigned to sand, complex, unknown interactions between the brittle sand (used as an upper crust equivalent) and viscous silicone (lower crust), highly non-linear strain weakening processes and poor constraints on the cohesion of sand. Our numerical-analogue comparison is hampered by (a) an incomplete knowledge of the fine details of sand failure and sand
CONVECTIVE DIFFUSION FIELD MEASUREMENTS COMPARED WITH LABORATORY AND NUMERICAL EXPERIMENTS
Some of the more fundamental diffusion parameters measured in the CONDORS convective diffusion field experiment are compared with laboratory experiment and numerical modeling results by means of nondimensionalizations using convective scaling (i.e., mixing depth, z sub i, for len...
Polymer Photooxidation: An Experiment to Demonstrate the Effect of Additives.
ERIC Educational Resources Information Center
Allen, Norman S.; McKellar, John F.
1979-01-01
This undergraduate experiment shows that the inclusion of an appropriate additive can have a very marked effect on the physical properties of a polymer. The polymer used is polypropylene and the additives are 2-hydroxy-4-octyloxy-benzophenone and benzophenone. (BB)
Numerical Simulation of High Drag Reduction in a Turbulent Channel Flow with Polymer Additives
NASA Technical Reports Server (NTRS)
Dubief, Yves
2003-01-01
The addition of small amounts of long chain polymer molecules to wall-bounded flows can lead to dramatic drag reduction. Although this phenomenon has been known for about fifty years, the action of the polymers and its effect on turbulent structures are still unclear. Detailed experiments have characterized two distinct regimes (Warholic et al. 1999), which are referred to as low drag reduction (LDR) and high drag reduction (HDR). The first regime exhibits similar statistical trends as Newtonian flow: the log-law region of the mean velocity profile remains parallel to that of the Newtonian ow but its lower bound moves away from the wall and the upward shift of the log-region is a function of drag reduction, DR. Although streamwise fluctuations are increased and transverse ones are reduced, the shape of the rms velocity profiles is not qualitatively modified. At higher drag reductions, of the order of 40-50%, the ow enters the HDR regime for which the slope of the log-law is dramatically augmented and the Reynolds shear stress is small (Warholic et al. 1999; Ptasinski et al. 2001). The drag reduction is eventually bounded by a maximum drag reduction (MDR) (Virk & Mickley 1970) which is a function of the Reynolds number. While several experiments report mean velocity profiles very close to the empirical profile of Virk & Mickley (1970) for MDR conditions, the observations regarding the structure of turbulence can differ significantly. For instance, Warholic et al. (1999) measured a near-zero Reynolds shear stress, whereas a recent experiment (Ptasinski et al. 2001) shows evidence of non-negligible Reynolds stress in their MDR flow. To the knowledge of the authors, only the LDR regime has been documented in numerical simulations (Sureshkumar et al. 1997; Dimitropoulos et al. 1998; Min et al. 2001; Dubief & Lele 2001; Sibilla & Baron 2002). This paper discusses the simulation of polymer drag reduced channel ow at HDR using the FENE-P (Finite Elastic non
Numerical Simulation and Cold Modeling experiments on Centrifugal Casting
NASA Astrophysics Data System (ADS)
Keerthiprasad, Kestur Sadashivaiah; Murali, Mysore Seetharam; Mukunda, Pudukottah Gopaliengar; Majumdar, Sekhar
2011-02-01
In a centrifugal casting process, the fluid flow eventually determines the quality and characteristics of the final product. It is difficult to study the fluid behavior here because of the opaque nature of melt and mold. In the current investigation, numerical simulations of the flow field and visualization experiments on cold models have been carried out for a centrifugal casting system using horizontal molds and fluids of different viscosities to study the effect of different process variables on the flow pattern. The effects of the thickness of the cylindrical fluid annulus formed inside the mold and the effects of fluid viscosity, diameter, and rotational speed of the mold on the hollow fluid cylinder formation process have been investigated. The numerical simulation results are compared with corresponding data obtained from the cold modeling experiments. The influence of rotational speed in a real-life centrifugal casting system has also been studied using an aluminum-silicon alloy. Cylinders of different thicknesses are cast at different rotational speeds, and the flow patterns observed visually in the actual castings are found to be similar to those recorded in the corresponding cold modeling experiments. Reasonable agreement is observed between the results of numerical simulation and the results of cold modeling experiments with different fluids. The visualization study on the hollow cylinders produced in an actual centrifugal casting process also confirm the conclusions arrived at from the cold modeling experiments and numerical simulation in a qualitative sense.
Numerical Experiments In Strongly Coupled Complex (Dusty) Plasmas
NASA Astrophysics Data System (ADS)
Hou, L. J.; Ivlev A.; Hubertus M. T.; Morfill, G. E.
2010-07-01
Complex (dusty) plasma is a suspension of micron-sized charged dust particles in a weakly ionized plasma with electrons, ions, and neutral atoms or molecules. Therein, dust particles acquire a few thousand electron charges by absorbing surrounding electrons and ions, and consequently interact with each other via a dynamically screened Coulomb potential while undergoing Brownian motion due primarily to frequent collisions with the neutral molecules. When the interaction potential energy between charged dust particles significantly exceeds their kinetic energy, they become strongly coupled and can form ordered structures comprising liquid and solid states. Since the motion of charged dust particles in complex (dusty) plasmas can be directly observed in real time by using a video camera, such systems have been generally regarded as a promising model system to study many phenomena occurring in solids, liquids and other strongly-coupled systems at the kinetic level, such as phase transitions, transport processes, and collective dynamics. Complex plasma physics has now grown into a mature research field with a very broad range of interdisciplinary facets. In addition to usual experimental and theoretical study, computer simulation in complex plasma plays an important role in bridging experimental observations and theories and in understanding many interesting phenomena observed in laboratory. The present talk will focus on a class of computer simulations that are usually non-equilibrium ones with external perturbation and that mimic the real complex plasma experiments (i. e., numerical experiment). The simulation method, i. e., the so-called Brownian Dynamics methods, will be firstly reviewed and then examples, such as simulations of heat transfer and shock wave propagation, will be present.
A Comparison of Metamodeling Techniques via Numerical Experiments
NASA Technical Reports Server (NTRS)
Crespo, Luis G.; Kenny, Sean P.; Giesy, Daniel P.
2016-01-01
This paper presents a comparative analysis of a few metamodeling techniques using numerical experiments for the single input-single output case. These experiments enable comparing the models' predictions with the phenomenon they are aiming to describe as more data is made available. These techniques include (i) prediction intervals associated with a least squares parameter estimate, (ii) Bayesian credible intervals, (iii) Gaussian process models, and (iv) interval predictor models. Aspects being compared are computational complexity, accuracy (i.e., the degree to which the resulting prediction conforms to the actual Data Generating Mechanism), reliability (i.e., the probability that new observations will fall inside the predicted interval), sensitivity to outliers, extrapolation properties, ease of use, and asymptotic behavior. The numerical experiments describe typical application scenarios that challenge the underlying assumptions supporting most metamodeling techniques.
Numerical simulations of the QUELL experiment in SULTAN
Marinucci, C.
1995-03-01
The QUench Experiment on Long Length (QUELL) in the SULTAN Facility is planned to investigate the quench propagation and detection of a conductor with ITER relevant geometry and scaled performance. The objective of this study is to show the ability of QUELL to provide quench conditions relevant for ITER and to simulate the system performance, dealing in particular with the design aspects of the power supply, cryogenic system and heaters. The numerical analysis was performed with GANDALF - a 1-D code to analyze Dual Channel Cable-in-Conduit Conductors. A numerical convergence test and a comparison with another code and with analytical results have confirmed the validity of the simulations.
Estimation of stream nutrient uptake from nutrient addition experiments
Payn, Robert
2005-09-01
Nutrient uptake in streams is often quantified by determining nutrient uptake length. However, current methods for measuring nutrient uptake length are often impractical, expensive, or demonstrably incorrect. We have developed a new method to estimate ambient nutrient uptake lengths using field experiments involving several levels of nutrient addition. Data analysis involves plotting nutrient addition uptake lengths versus added concentration and extrapolating to the negative ambient concentration. This method is relatively easy, inexpensive, and based on sound theoretical development. It is more accurate than the commonly used method involving a single nutrient addition. The utility of the method is supported by field studies directly comparing our new method with isotopic tracer methods for determining uptake lengths of phosphorus, ammonium, and nitrate. Our method also provides parameters for comparing potential nutrient limitation among streams.
NASA Astrophysics Data System (ADS)
Martin-Short, R.; Edmiston, J. K.
2015-12-01
Typical hydraulic fracturing operations involve the use of a large quantity of water, which can be problematic for several reasons including possible formation (permeability) damage, disposal of waste water, and the use of precious local water resource. An alternate reservoir permeability enhancing technology not requiring water is cryogenic fracturing. This method induces controlled fracturing of rock formations by thermal shock and has potentially important applications in the geothermal and hydrocarbon industries. In this process, cryogenic fluid—such as liquid nitrogen—is injected into the subsurface, causing fracturing due to thermal gradients. These fractures may improve the formation permeability relative to that achievable by hydraulic fracturing alone. We conducted combined laboratory visualization and numerical simulations studies of thermal-shock-induced fracture initiation and propagation resulting from liquid nitrogen injection in rock and analog materials. The experiment used transparent soda-lime glass cubes to facilitate real-time visualization of fracture growth and the fracture network geometry. In this contribution, we report the effect of overall temperature difference between cryogenic fluid and solid material on the produced fracture network, by pre-heating the glass cubes to several temperatures and injecting liquid nitrogen. Temperatures are monitored at several points by thermocouple and the fracture evolution is captured visually by camera. The experiment was modeled using a customized, thermoelastic, fracture-capable numerical simulation code based on peridynamics. The performance of the numerical code was validated by the results of the laboratory experiments, and then the code was used to study the different factors affecting a cryogenic fracturing operation, including the evolution of residual stresses and constitutive relationships for material failure. In complex rock such as shale, understanding the process of cryogenic
Numerical modeling of injection experiments at The Geysers
Pruess, Karsten; Enedy, Steve
1993-01-28
Data from injection experiments in the southeast Geysers are presented that show strong interference (both negative and positive) with a neighboring production well. Conceptual and numerical models are developed that explain the negative interference (decline of production rate) in terms of heat transfer limitations and water-vapor relative permeability effects. Recovery and overrecovery following injection shut-in are attributed to boiling of injected fluid, with heat of vaporization provided by the reservoir rocks.
Numerical modeling of injection experiments at The Geysers
Pruess, K. ); Enedy, S. )
1993-01-01
Data from injection experiments in the southeast Geysers are presented that show strong interference (both negative and positive) with a neighboring production well. Conceptual and numerical models are developed that explain the negative interference (decline of production rate) in terms of heat transfer limitations and water-vapor relative permeability effects. Recovery and over-recovery following injection shut-in are attributed to boiling of injected fluid, with heat of vaporization provided by the reservoir rocks.
Effects of an additional dimension in the Young experiment
Barros, Allan Kardec
2015-09-15
The results of the Young experiment can be analyzed either by classical or Quantum Physics. The later one though leads to a more complete interpretation, based on two different patterns that appear when one works either with single or double slits. Here we show that the two patterns can be derived from a single principle, in the context of General Relativity, if one assumes an additional spatial dimension to the four known today. The found equations yield the same results as those in Quantum Mechanics.
Numerical simulation of experiments in the Giant Planet Facility
NASA Technical Reports Server (NTRS)
Green, M. J.; Davy, W. C.
1979-01-01
Utilizing a series of existing computer codes, ablation experiments in the Giant Planet Facility are numerically simulated. Of primary importance is the simulation of the low Mach number shock layer that envelops the test model. The RASLE shock-layer code, used in the Jupiter entry probe heat-shield design, is adapted to the experimental conditions. RASLE predictions for radiative and convective heat fluxes are in good agreement with calorimeter measurements. In simulating carbonaceous ablation experiments, the RASLE code is coupled directly with the CMA material response code. For the graphite models, predicted and measured recessions agree very well. Predicted recession for the carbon phenolic models is 50% higher than that measured. This is the first time codes used for the Jupiter probe design have been compared with experiments.
INEX (integrated numerical experiment) simulations of the Boeing FEL system
Tokar, R.L.; Young, L.M.; Lumpkin, A.H.; McVey, B.D.; Thode, L.E.; Bender, S.C.; Chan, K.C.D. ); Yeremian, A.D.; Dowell, D.H.; Lowrey, A.R. )
1989-01-01
The INEX (integrated numerical experiment) numerical model is applied to the 0.6 {mu}m FEL oscillator at Boeing Aerospace and Electronics Company in Seattle, WA. This system consists of a 110 MeV L-band rf linac, a beam transport line from the accelerator to the entrance of the wiggler, the 5.0 meter THUNDER variable taper wiggler, and a near concentric two mirror optical oscillator. Many aspects of the model for the electron beam accelerator and transport line agree with experimental measurements. Predictions for lasing performance are compared with data obtained in May and June 1989 using a mild tapered wiggler. We obtain good agreement with the achieved extraction efficiency, while 1D pulse simulations reproduce the observed sideband instability. 15 refs., 11 figs.
Numerical simulation of the hydrodynamic instability experiments and flow mixing
NASA Astrophysics Data System (ADS)
Bai, Jingsong; Wang, Tao; Li, Ping; Zou, Liyong; Liu, Cangli
2009-12-01
Based on the numerical methods of volume of fluid (VOF) and piecewise parabolic method (PPM) and parallel circumstance of Message Passing Interface (MPI), a parallel multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and performed to study the hydrodynamic instability and flow mixing. Firstly, the MVPPM code is verified and validated by simulating three instability cases: The first one is a Riemann problem of viscous flow on the shock tube; the second one is the hydrodynamic instability and mixing of gaseous flows under re-shocks; the third one is a half height experiment of interfacial instability, which is conducted on the AWE’s shock tube. By comparing the numerical results with experimental data, good agreement is achieved. Then the MVPPM code is applied to simulate the two cases of the interfacial instabilities of jelly models accelerated by explosion products of a gaseous explosive mixture (GEM), which are adopted in our experiments. The first is implosive dynamic interfacial instability of cylindrical symmetry and mixing. The evolving process of inner and outer interfaces, and the late distribution of mixing mass caused by Rayleigh-Taylor (RT) instability in the center of different radius are given. The second is jelly layer experiment which is initialized with one periodic perturbation with different amplitude and wave length. It reveals the complex processes of evolution of interface, and presents the displacement of front face of jelly layer, bubble head and top of spike relative to initial equilibrium position vs. time. The numerical results are in excellent agreement with that experimental images, and show that the amplitude of initial perturbations affects the evolvement of fluid mixing zone (FMZ) growth rate extremely, especially at late times.
Tsai, Pei-I; Hsu, Ching-Chi; Chen, San-Yuan; Wu, Tsung-Han; Huang, Chih-Chieh
2016-09-01
Traditional solid cages have been widely used in posterior lumbar interbody fusion (PLIF) surgery. However, solid cages significantly affect the loading mechanism of the human spine due to their extremely high structural stiffness. Previous studies proposed and investigated porous additive manufactured (AM) cages; however, their biomechanical performances were analyzed using oversimplified bone-implant numerical models. Thus, the aim of this study was to investigate the outer shape and inner porous structure of the AM cages. The outer shape of the AM cages was discovered using a simulation-based genetic algorithm; their inner porous structure was subsequently analyzed parametrically using T10-S1 multilevel spine models. Finally, six types of the AM cages, which were manufactured using selective laser melting, were tested to validate the numerical outcomes. The subsidence resistance of the optimum design was superior to the conventional cage designs. A porous AM cage with a pillar diameter of 0.4mm, a pillar angle of 40°, and a porosity of between 69% and 80% revealed better biomechanical performances. Both the numerical and experimental outcomes can help surgeons to understand the biomechanics of PLIF surgery combined with the use of AM cages. PMID:27392226
Numerical modeling of oxygen exclusion experiments of anaerobic bioventing.
Mihopoulos, Philip G; Suidan, Makram T; Sayles, Gregory D; Kaskassian, Sebastien
2002-10-01
A numerical and experimental study of transport phenomena underlying anaerobic bioventing (ABV) is presented. Understanding oxygen exclusion patterns in vadose zone environments is important in designing an ABV process for bioremediation of soil contaminated with chlorinated solvents. In particular, the establishment of an anaerobic zone of influence by nitrogen injection in the vadose zone is investigated. Oxygen exclusion experiments are performed in a pilot scale flow cell (2 x 1.1 x 0.1 m) using different venting flows and two different outflow boundary conditions (open and partially covered). Injection gas velocities are varied from 0.25 x 10(-3) to 1.0 x 10(-3) cm/s and are correlated with the ABV radius of influence. Numerical simulations are used to predict the collected experimental data. In general, reasonable agreement is found between observed and predicted oxygen concentrations. Use of impervious covers can significantly reduce the volume of forcing gas used, where an increase in oxygen exclusion efficiency is consistent with a decrease in the outflow area above the injection well. PMID:12400833
Numerical experiments on the stability of preplanetary disks
NASA Technical Reports Server (NTRS)
Cassen, P. M.; Smith, B. F.; Reynolds, R. T.; Miller, R. H.
1981-01-01
Gravitational stability of gaseous protostellar disks is relevant to theories of planetary formation. Stable gas disks favor formation of planetesimals by the accumulation of solid material; unstable disks allow the possibility of direct condensation of gaseous protoplanets. This paper presents the results of numerical experiments designed to test the stability of thin disks against large-scale, self-gravitational disruption. It is found that a disk as massive as 1 solar mass, surrounding a 1 solar mass protostar, can be stable against long-wavelength gravitational disruption if its temperature is about 300 K or greater. Stability of a cooler disk requires that it be less massive, but even at 100 K a stable disk can have an appreciable fraction (about 1/3) of a solar mass.
NASA Astrophysics Data System (ADS)
Chiggiato, Jacopo; Jarosz, Ewa; Book, Jeffrey W.; Dykes, James; Torrisi, Lucio; Poulain, Pierre-Marie; Gerin, Riccardo; Horstmann, Jochen; Beşiktepe, Şükrü
2012-01-01
During September 2008 and February 2009, the NR/V Alliance extensively sampled the waters of the Sea of Marmara within the framework of the Turkish Straits System (TSS) experiment coordinated by the NATO Undersea Research Centre. The observational effort provided an opportunity to set up realistic numerical experiments for modeling the observed variability of the Marmara Sea upper layer circulation at mesoscale resolution over the entire basin during the trial period, complementing relevant features and forcing factors revealed by numerical model results with information acquired from in situ and remote sensing datasets. Numerical model solutions from realistic runs using the Regional Ocean Modeling System (ROMS) produce a general circulation in the Sea of Marmara that is consistent with previous knowledge of the circulation drawn from past hydrographic measurements, with a westward meandering current associated with a recurrent large anticyclone. Additional idealized numerical experiments illuminate the role various dynamics play in determining the Sea of Marmara circulation and pycnocline structure. Both the wind curl and the strait flows are found to strongly influence the strength and location of the main mesoscale features. Large displacements of the pycnocline depth were observed during the sea trials. These displacements can be interpreted as storm-driven upwelling/downwelling dynamics associated with northeasterly winds; however, lateral advection associated with flow from the Straits also played a role in some displacements.
Numerical simulation of molten silicon flow; comparison with experiment
NASA Astrophysics Data System (ADS)
Kakimoto, Koichi; Nicodème, Pierre; Lecomte, Michael; Dupret, François; Crochet, Marcel J.
1991-12-01
Numerical simulation containing fluid flow, heat conduction and heat exchange by radiation has been performed using the geometry of a real Czochralski furnace for silicon single crystal growth. The flow velocity fields of molten silicon are obtained from extrapolation of the stream function, which has been newly developed using the velocity boundary layer theory. The calculated flow velocity and particle path are semi-quantitatively identical to the results obtained from X-ray radiography experiment. The calculated value of the characteristic velocity is about 10 -2 m/s. The same order of flow velocity which is obtained from the experiment has been already reported. It has also become clear from a comparison of flow velocities between experimental and calculated results that the order of the volume expansion coefficient of the molten silicon (β) is 10 -4 K -1. The flow was almost axisymmetric and steady for a specific case with low crystal and crucible rotation rates and with a shallow melt. We also found that a flow with larger azimuthal velocity component exists just beneath a crystal, while that with opposite flow direction exists near the crucible wall.
Numerical Modelling of the Deep Impact Mission Experiment
NASA Technical Reports Server (NTRS)
Wuennemann, K.; Collins, G. S.; Melosh, H. J.
2005-01-01
NASA s Deep Impact Mission (launched January 2005) will provide, for the first time ever, insights into the interior of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of a comets nucleus. Although it is usually assumed that comets consist of a very porous mixture of water ice and rock, little is known about the internal structure and in particular the constitutive material properties of a comet. It is therefore difficult to predict the dimensions of the excavated crater. Estimates of the crater size are based on laboratory experiments of impacts into various target compositions of different densities and porosities using appropriate scaling laws; they range between 10 s of meters up to 250 m in diameter [1]. The size of the crater depends mainly on the physical process(es) that govern formation: Smaller sizes are expected if (1) strength, rather than gravity, limits crater growth; and, perhaps even more crucially, if (2) internal energy losses by pore-space collapse reduce the coupling efficiency (compaction craters). To investigate the effect of pore space collapse and strength of the target we conducted a suite of numerical experiments and implemented a novel approach for modeling porosity and the compaction of pores in hydrocode calculations.
Cascade processes in stratified media: experiment and direct numerical simulation.
NASA Astrophysics Data System (ADS)
Sibgatullin, Ilias; Brouzet, Christophe; Joubaud, Sylvain; Ermanyuk, Evgeny; Dauxois, Thierry
2016-04-01
Internal gravity waves may transfer substantial part of energy in oceans and astrophysical objects, influence the background stratification, and angular momentum. Internal waves can be generated by convection in astrophysical objects, by tidal motion and interaction with orography in oceans. Internal and inertial waves obey similar system of equations. Due to very particular type of dispersive relation and the way internal waves are reflected from surfaces, in confined domains the monochromatic internal waves after sequence of reflections may form closed paths, the "wave attractors" [1]. Presently, linear theory of wave attractors is quite elaborated and a principal interest of research is focused on nonlinear regimes and unstable configurations, overturning events and mixing. We have performed direct numerical simulation of wave attractors which closely reproduces experiments [2] being carried out in Ecole Normal Superior de Lyon (ENS de Lyon). Direct numerical simulation is realized with the help of spectral element approach and code nek5000. Triadic resonance is confirmed as the first instability which appears on the most energetic ray of the attractor at sufficiently large forcing. With further increase of the forcing amplitude the daughter waves also become unstable resulting in a sophisticated cascade process which was first observed experimentally. For very high forcing amplitude interaction of focused waves with the walls results in appearance of small-scale folded structures. Their interaction with principal flow is the subject of further research. 1. Maas, L. R. M. & Lam, F.-P. A., Geometric focusing of internal waves. J. Fluid Mech, 1995,. 300, 1-41 2. Scolan, H., Ermanyuk, E., Dauxois, T., 2013, Physical Review Letters, 110, 234501
NASA Technical Reports Server (NTRS)
Smalheer, C. V.
1973-01-01
The chemistry of lubricant additives is discussed to show what the additives are chemically and what functions they perform in the lubrication of various kinds of equipment. Current theories regarding the mode of action of lubricant additives are presented. The additive groups discussed include the following: (1) detergents and dispersants, (2) corrosion inhibitors, (3) antioxidants, (4) viscosity index improvers, (5) pour point depressants, and (6) antifouling agents.
McKenna, S.A.; Rautman, C.A.
1996-08-01
A review of pertinent literature reveals techniques which may be practical for upscaling saturated hydraulic conductivity at Yucca Mountain: geometric mean, spatial averaging, inverse numerical modeling, renormalization, and a perturbation technique. Isotropic realizations of log hydraulic conductivity exhibiting various spatial correlation lengths are scaled from the point values to five discrete scales through these techniques. For the variances in log{sub 10} saturated hydraulic conductivity examined here, geometric mean, numerical inverse and renormalization adequately reproduce point scale fluxes across the modeled domains. Fastest particle velocities and dispersion measured on the point scale are not reproduced by the upscaled fields. Additional numerical experiments examine the utility of power law averaging on a geostatistical realization of a cross-section similar to the cross-sections that will be used in the 1995 groundwater travel time calculations. A literature review on scaling techniques for thermal and mechanical properties is included. 153 refs., 29 figs., 6 tabs.
Guo, Hongsheng; Liu, Fengshan; Smallwood, Gregory J.; Guelder, OEmer L.
2006-04-15
The influence of hydrogen addition to the fuel of an atmosphere pressure coflow laminar ethylene-air diffusion flame on soot formation was studied by numerical simulation. A detailed gas-phase reaction mechanism, which includes aromatic chemistry up to four rings, and complex thermal and transport properties were used. The fully coupled elliptic governing equations were solved. The interactions between soot and gas-phase chemistry were taken into account. Radiation heat transfer from CO{sub 2}, CO, H{sub 2}O, and soot was calculated using the discrete-ordinates method coupled to a statistical narrow-band-correlated K-based wide-band model. The predicted results were compared with the available experimental data and analyzed. It is indicated that the addition of hydrogen to the fuel in an ethylene-air diffusion flame suppresses soot formation through the effects of dilution and chemistry. This result is in agreement with available experiments. The simulations further suggest that the chemically inhibiting effect of hydrogen addition on soot formation is due to the decrease of hydrogen atom concentration in soot surface growth regions and higher concentration of molecular hydrogen in the lower flame region. (author)
Numerical experiments with rubble piles: equilibrium shapes and spins
NASA Astrophysics Data System (ADS)
Richardson, Derek C.; Elankumaran, Pradeep; Sanderson, Robyn E.
2005-02-01
We present numerical experiments investigating the shape and spin limits of self-gravitating "perfect" rubble piles that consist of identical, smooth, rigid, spherical particles with configurable normal coefficient of restitution and no sliding friction. Such constructs are currently employed in a variety of investigations, ranging from the formation of asteroid satellites to the dynamical properties of Saturn's densest rings. We find that, owing to cannonball stacking behavior, rubble piles can maintain non-spherical shapes without bulk spin, unlike a fluid, and can spin faster than a perfect fluid before shedding mass, consistent with the theory for the more general continuum rubble pile model (Holsapple, 2004, Icarus 172, 272-303). Rubble piles that reassemble following a catastrophic disruption reconfigure themselves to lie within stability limits predicted by the continuum theory. We also find that coarse configurations consisting of a small number of particles are more resistant to tidal disruption than fine configurations with many particles. Overall this study shows that idealized rubble piles behave qualitatively in a manner similar to certain granular materials, at least in the limit where global shape readjustments and/or mass shedding begins. The limits obtained here may provide constraints on the possible internal structure of some small Solar System bodies that have extreme shapes or are under high stress. Amalthea is presented as a case study.
Numerical experiments of fracture-induced velocity and attenuation anisotropy
NASA Astrophysics Data System (ADS)
Carcione, J. M.; Picotti, S.; Santos, J. E.
2012-12-01
Fractures are common in the Earth's crust due to different factors, for instance, tectonic stresses and natural or artificial hydraulic fracturing caused by a pressurized fluid. A dense set of fractures behaves as an effective long-wavelength anisotropic medium, leading to azimuthally varying velocity and attenuation of seismic waves. Effective in this case means that the predominant wavelength is much longer than the fracture spacing. Here, fractures are represented by surface discontinuities in the displacement u and particle velocity v as ?, where the brackets denote the discontinuity across the surface, ? is a fracture stiffness and ? is a fracture viscosity. We consider an isotropic background medium, where a set of fractures are embedded. There exists an analytical solution—with five stiffness components—for equispaced plane fractures and an homogeneous background medium. The theory predicts that the equivalent medium is transversely isotropic and viscoelastic. We then perform harmonic numerical experiments to compute the stiffness components as a function of frequency, by using a Galerkin finite-element procedure, and obtain the complex velocities of the medium as a function of frequency and propagation direction, which provide the phase velocities, energy velocities (wavefronts) and quality factors. The algorithm is tested with the analytical solution and then used to obtain the stiffness components for general heterogeneous cases, where fractal variations of the fracture compliances and background stiffnesses are considered.
Experiments and numerical investigations of wave propagation in thermal plumes
NASA Astrophysics Data System (ADS)
Laudenbach, N.
2001-12-01
In laboratory experiments thermal plumes are created by injecting hot corn syrup into a column of cold syrup. The viscosity contrast is up to a factor of 1000. Solitary waves, that propagate upwards in the plume conduit, are generated by enhancing the injection rate for a few seconds. For the measurement of the thermal structure of the plume we have implemented a method based on the deflection of a laser beam passing through the plume. Continuous scanning provides a new radial temperature profile each second, which allows detailed studies of the thermal structure of solitary waves. A PIV - (particle image volecimetry) method provides the velocity structure of the thermal plume. Measurements were taken for plume heads, conduits and propagating waves. Comparison between experimental results and numerical 2-D axisymmetric simulations shows a good agreement of the temperature profiles and velocity fields in the plume conduit and waves. Because of thermal diffusion, the conduit widens with height, while its central temperature decreases. The solitary waves start with the same temperature as the unperturbed conduit, however, we find that the temperature in the waves decreases less rapide with rising height. This can be explained by the faster upward propagation and the trapping of fluid within the soliton. If solitary waves exists in mantle plumes, this would imply that they arrive at the bottom of the lithosphere with a larger excess temperature than what the plumes normally exhibits, which could explain strong variations of melt generation with time.
Experiments and Numerical Investigations of Wave Propagation In Thermal Plumes.
NASA Astrophysics Data System (ADS)
Laudenbach, N.; Christensen, U. R.
In laboratory experiments thermal plumes are created by injecting hot corn syrup into a column of cold syrup. The viscosity contrast is up to a factor of 1000. Solitary waves, that propagate upwards in the plume conduit, are generated by enhancing the injection rate for a few seconds. For the measurement of the thermal structure of the plume we have implemented a method based on the deflection of a laser beam passing through the plume. Continuous scanning provides a new radial temperature profile each second, which allows detailed studies of the thermal structure of solitary waves. A PIV - (particle image volecimetry) method provides the velocity structure of the thermal plume. Measurements were taken for plume heads, conduits and propagating waves. Comparison between experimental results and numerical 2-D axisymmetric simulations shows a good agreement of the temperature profiles and velocity fields in the plume conduit and waves. Because of thermal diffusion, the conduit widens with height, while its central temperature decreases. The solitary waves start with the same temperature as the unperturbed conduit, however, we find that the temperature in the waves decreases less rapide with rising height. This can be explained by the faster upward propagation and the trapping of fluid within the soliton. If solitary waves exists in mantle plumes, this would imply that they arrive at the bottom of the lithosphere with a larger excess temperature than what the plumes normally exhibits. Especially for weak hotspots solitary waves could have strong influence on the variation of melt generation with time.
Rider, William; Kamm, J. R.; Tomkins, C. D.; Zoldi, C. A.; Prestridge, K. P.; Marr-Lyon, M.; Rightley, P. M.; Benjamin, R. F.
2002-01-01
We consider the detailed structures of mixing flows for Richtmyer-Meshkov experiments of Prestridge et al. [PRE 00] and Tomkins et al. [TOM 01] and examine the most recent measurements from the experimental apparatus. Numerical simulations of these experiments are performed with three different versions of high resolution finite volume Godunov methods. We compare experimental data with simulations for configurations of one and two diffuse cylinders of SF{sub 6} in air using integral measures as well as fractal analysis and continuous wavelet transforms. The details of the initial conditions have a significant effect on the computed results, especially in the case of the double cylinder. Additionally, these comparisons reveal sensitive dependence of the computed solution on the numerical method.
Convective melting in a magma chamber: theory and numerical experiment.
NASA Astrophysics Data System (ADS)
Simakin, A.
2012-04-01
conditions for the intruded magma to crystallize first and then switch to the roof melting or only crystallize were not defined. We did this in our numerical experiments in terms of the initial magma and roof rocks temperatures for particular sill size. Neglecting strong viscosity variation in the boundary layer at the melting front leads to the overestimation of the melting rate by H&S model on approximately 70% at Tm=940oC. At Tm =800oC effect of the crystals present in descending plumes compensates viscosity increase and numerical Um practically coincides with theoretical one (difference 8%). Some researchers (Huber et al., 2010) use empirical and scaling results obtained from stagnant-lid convection (Davaille and Jaupart, 1993). We find that the later model is not applicable to the melting problem since super-exponential dependence of the viscosity from temperature is valid providing full solidification below eutectic temperature Ts. "Melting temperature" at the stagnant-lid style of convection is defined by Arhenius rheological parameters and bulk melt temperature and can be less than Ts. Our numerical study was applied to the estimation of the possible time frame and efficiency of the remelting of the silicic pyroclastics by superheated rhyolites in the caldera environment (Simakin and Bindeman, 2012). Literature. 1)Davaille, A. and Jaupart, C. (1993) J. Fluid. Mech., 253: 141-166. 2) Huppert, H.E. and Sparks, R.S.J. (1988) J. Petrol., 29: 599-624. 3)Huber, C., Bachmann, O., Dufek, J. (2010) J. Volcanol. Geotherm. Res., 195: 97-105. 4)Jaupart, C. and Brandeis, G. (1986) Earth Planet. Sci. Lett. 80: 183-199. 5)Simakin, A.G. and Bindeman, I.N. (2012) Remelting in caldera and rift environments and the genesis of hot, "recycled" rhyolites. Earth Planet. Sci. Lett. (in review). 6) Snyder, D. (2000) Earth Planet. Sci. Lett. 175: 257-273.
A numerical experiment on light pollution from distant sources
NASA Astrophysics Data System (ADS)
Kocifaj, M.
2011-08-01
To predict the light pollution of the night-time sky realistically over any location or measuring point on the ground presents quite a difficult calculation task. Light pollution of the local atmosphere is caused by stray light, light loss or reflection of artificially illuminated ground objects or surfaces such as streets, advertisement boards or building interiors. Thus it depends on the size, shape, spatial distribution, radiative pattern and spectral characteristics of many neighbouring light sources. The actual state of the atmospheric environment and the orography of the surrounding terrain are also relevant. All of these factors together influence the spectral sky radiance/luminance in a complex manner. Knowledge of the directional behaviour of light pollution is especially important for the correct interpretation of astronomical observations. From a mathematical point of view, the light noise or veil luminance of a specific sky element is given by a superposition of scattered light beams. Theoretical models that simulate light pollution typically take into account all ground-based light sources, thus imposing great requirements on CPU and MEM. As shown in this paper, a contribution of distant sources to the light pollution might be essential under specific conditions of low turbidity and/or Garstang-like radiative patterns. To evaluate the convergence of the theoretical model, numerical experiments are made for different light sources, spectral bands and atmospheric conditions. It is shown that in the worst case the integration limit is approximately 100 km, but it can be significantly shortened for light sources with cosine-like radiative patterns.
A Numerical Method for Determining Diffusivity from Annealing Experiments
NASA Astrophysics Data System (ADS)
Harris-Kuhlman, K. R.; Kulcinski, G. L.
1998-12-01
Terrestrial analogs of lunar ilmenite (FeTiO3) have been implanted with solar-wind energy 4He at 4 keV and 3He at 3 keV using Plasma Source Ion Implantation (PSII). Isochronal annealing of the samples revealed thermally induced 4He evolution similar to the helium release of the Apollo 11 regoliths reported by Pepin, et. al., [1970]. These annealing experiments are analyzed with a three dimensional numerical method based on Fick's law for diffusion. An iterative method is used to calculate the diffusivity. The code uses an assumed diffusivity to calculate the amount of gas released during a temperature step. The initial depth profile of the implanted species is generated using the TRIM electronic stopping code [Ziegler, 1996]. The calculated value is compared to the measured value and a linear regression is used to calculate a new diffusivity until there is convergence within a specified tolerance level. The diffusivity as a function of temperature is then fitted to an Arrhenius equation. Analysis of results for 4 keV 4He on ilmenite shows two distinct regions of Arrehnius behavior with activation energies of 0.5 +/- 0.1 eV at emperatures below 800 deg C and 1.5 +/- 0.2 eV at temperatures from 800 deg C to 1100 deg C. Pepin, R. O., L. E. Nyquist, D. Phinney, and D. C. Black (1970) "Rare Gases in Apollo 11 Lunar Material," Proceedings of the Apollo 11 Lunar Science Conference, 2, pp. 1435-1454. Ziegler, J. P. (1996) SRIM Instruction Manual: The Stopping and Range of Ions in Matter, (Yorktown, New York: IBM - Research); based on Ziegler, J. P., J. P. Biersack and U. Littmark, The Stopping and Range of Ions in Solids, (New York: Pergamon Press, 1985).
Software reliability: Additional investigations into modeling with replicated experiments
NASA Technical Reports Server (NTRS)
Nagel, P. M.; Schotz, F. M.; Skirvan, J. A.
1984-01-01
The effects of programmer experience level, different program usage distributions, and programming languages are explored. All these factors affect performance, and some tentative relational hypotheses are presented. An analytic framework for replicated and non-replicated (traditional) software experiments is presented. A method of obtaining an upper bound on the error rate of the next error is proposed. The method was validated empirically by comparing forecasts with actual data. In all 14 cases the bound exceeded the observed parameter, albeit somewhat conservatively. Two other forecasting methods are proposed and compared to observed results. Although demonstrated relative to this framework that stages are neither independent nor exponentially distributed, empirical estimates show that the exponential assumption is nearly valid for all but the extreme tails of the distribution. Except for the dependence in the stage probabilities, Cox's model approximates to a degree what is being observed.
Floret Test, Numerical Simulations of the Dent, Comparison with Experiments
Lefrancois, A.; Cutting, J.; Gagliardi, F.; Tarver, C.; Tran, T.
2006-02-14
The Floret test has been developed as a screening test to study the performance of a small amount of HE. Numerical simulations have been performed recently using CTH. The objective of this study is to perform numerical simulations in order to better understand the shock waves interactions, involved in the dent formation. Different 3D wedge configurations have been tested using the Ignition and Growth reactive flow model for the HE receptor with Ls-Dyna.
NASA Technical Reports Server (NTRS)
Shia, Run-Lie; Ha, Yuk Lung; Wen, Jun-Shan; Yung, Yuk L.
1990-01-01
Extensive testing of the advective scheme proposed by Prather (1986) has been carried out in support of the California Institute of Technology-Jet Propulsion Laboratory two-dimensional model of the middle atmosphere. The original scheme is generalized to include higher-order moments. In addition, it is shown how well the scheme works in the presence of chemistry as well as eddy diffusion. Six types of numerical experiments including simple clock motion and pure advection in two dimensions have been investigated in detail. By comparison with analytic solutions, it is shown that the new algorithm can faithfully preserve concentration profiles, has essentially no numerical diffusion, and is superior to a typical fourth-order finite difference scheme.
Trescott, Peter C.; Pinder, George Francis; Larson, S.P.
1976-01-01
The model will simulate ground-water flow in an artesian aquifer, a water-table aquifer, or a combined artesian and water-table aquifer. The aquifer may be heterogeneous and anisotropic and have irregular boundaries. The source term in the flow equation may include well discharge, constant recharge, leakage from confining beds in which the effects of storage are considered, and evapotranspiration as a linear function of depth to water. The theoretical development includes presentation of the appropriate flow equations and derivation of the finite-difference approximations (written for a variable grid). The documentation emphasizes the numerical techniques that can be used for solving the simultaneous equations and describes the results of numerical experiments using these techniques. Of the three numerical techniques available in the model, the strongly implicit procedure, in general, requires less computer time and has fewer numerical difficulties than do the iterative alternating direction implicit procedure and line successive overrelaxation (which includes a two-dimensional correction procedure to accelerate convergence). The documentation includes a flow chart, program listing, an example simulation, and sections on designing an aquifer model and requirements for data input. It illustrates how model results can be presented on the line printer and pen plotters with a program that utilizes the graphical display software available from the Geological Survey Computer Center Division. In addition the model includes options for reading input data from a disk and writing intermediate results on a disk.
Numerical Simulation of the Perrin-Like Experiments
ERIC Educational Resources Information Center
Mazur, Zygmunt; Grech, Dariusz
2008-01-01
A simple model of the random Brownian walk of a spherical mesoscopic particle in viscous liquids is proposed. The model can be solved analytically and simulated numerically. The analytic solution gives the known Einstein-Smoluchowski diffusion law r[superscript 2] = 2Dt, where the diffusion constant D is expressed by the mass and geometry of a…
NASA Astrophysics Data System (ADS)
Jiang, C. X.; Cheng, J. P.; Li, F. C.
2015-01-01
This paper attempts to introduce a numerical simulation procedure to simulate water-entry problems influenced by turbulent drag-reducing additives in a viscous incompressible medium. Firstly we performed a numerical investigation on water-entry supercavities in water and turbulent drag-reducing solution at the impact velocity of 28.4 m/s to confirm the accuracy of the numerical method. Based on the verification, projectile entering water and turbulent drag-reducing solution at relatively high velocity of 142.7 m/s (phase transition is considered) is simulated. The cross viscosity equation was adopted to represent the shear-thinning characteristic of aqueous solution of drag-reducing additives. The configuration and dynamic characteristics of water entry supercavity, flow resistance were discussed respectively. It was obtained that the numerical simulation results are in consistence with experimental data. Numerical results show that the supercavity length in drag-reducing solution is larger than one in water and the velocity attenuates faster at high velocity than at low velocity; the influence of drag-reducing solution is more obvious at high impact velocity. Turbulent drag-reducing additives have the great potential for enhancement of supercavity.
Mechanical characterisation of Dacron graft: Experiments and numerical simulation.
Bustos, Claudio A; García-Herrera, Claudio M; Celentano, Diego J
2016-01-01
Experimental and numerical analyses focused on the mechanical characterisation of a woven Dacron vascular graft are presented. To that end, uniaxial tensile tests under different orientations have been performed to study the anisotropic behaviour of the material. These tests have been used to adjust the parameters of a hyperelastic anisotropic constitutive model which is applied to predict through numerical simulation the mechanical response of this material in the ring tensile test. The obtained results show that the model used is capable of representing adequately the nonlinear elastic region and, in particular, it captures the progressive increase of the rigidity and the anisotropy due to the stretching of the Dacron. The importance of this research lies in the possibility of predicting the graft׳s mechanical response under generalized loading such as those that occur under physiological conditions after surgical procedures. PMID:26627367
Taylor bubbles at high viscosity ratios: experiments and numerical simulations
NASA Astrophysics Data System (ADS)
Hewakandamby, Buddhika; Hasan, Abbas; Azzopardi, Barry; Xie, Zhihua; Pain, Chris; Matar, Omar
2015-11-01
The Taylor bubble is a single long bubble which nearly fills the entire cross section of a liquid-filled circular tube, often occurring in gas-liquid slug flows in many industrial applications, particularly oil and gas production. The objective of this study is to investigate the fluid dynamics of three-dimensional Taylor bubble rising in highly viscous silicone oil in a vertical pipe. An adaptive unstructured mesh modelling framework is adopted here which can modify and adapt anisotropic unstructured meshes to better represent the underlying physics of bubble rising and reduce computational effort without sacrificing accuracy. The numerical framework consists of a mixed control volume and finite element formulation, a `volume of fluid'-type method for the interface-capturing based on a compressive control volume advection method, and a force-balanced algorithm for the surface tension implementation. Experimental results for the Taylor bubble shape and rise velocity are presented, together with numerical results for the dynamics of the bubbles. A comparison of the simulation predictions with experimental data available in the literature is also presented to demonstrate the capabilities of our numerical method. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.
NASA Technical Reports Server (NTRS)
Bowman, Kenneth P.; Sacks, Jerome; Chang, Yue-Fang
1993-01-01
Methods for the design and analysis of numerical experiments that are especially useful and efficient in multidimensional parameter spaces are presented. The analysis method, which is similar to kriging in the spatial analysis literature, fits a statistical model to the output of the numerical model. The method is applied to a fully nonlinear, global, equivalent-barotropic dynamical model. The statistical model also provides estimates for the uncertainty of predicted numerical model output, which can provide guidance on where in the parameter space to conduct further experiments, if necessary. The method can provide significant improvements in the efficiency with which numerical sensitivity experiments are conducted.
Slump Flows inside Pipes: Numerical Results and Comparison with Experiments
NASA Astrophysics Data System (ADS)
Malekmohammadi, S.; Naccache, M. F.; Frigaard, I. A.; Martinez, D. M.
2008-07-01
In this work an analysis of the buoyancy-driven slumping flow inside a pipe is presented. This flow usually occurs when an oil well is sealed by a plug cementing process, where a cement plug is placed inside the pipe filled with a lower density fluid, displacing it towards the upper cylinder wall. Both the cement and the surrounding fluids have a non Newtonian behavior. The cement is viscoplastic and the surrounding fluid presents a shear thinning behavior. A numerical analysis was performed to evaluate the effects of some governing parameters on the slump length development. The conservation equations of mass and momentum were solved via a finite volume technique, using Fluent software (Ansys Inc.). The Volume of Fluid surface-tracking method was used to obtain the interface between the fluids and the slump length as a function of time. The results were obtained for different values of fluids densities differences, fluids rheology and pipe inclinations. The effects of these parameters on the interface shape and on the slump length versus time curve were analyzed. Moreover, the numerical results were compared to experimental ones, but some differences are observed, possibly due to chemical effects at the interface.
Numerical Experiments for Storm Surge Inundation in Korean Coastal Area
NASA Astrophysics Data System (ADS)
Yoon, J.; Shim, J.; Jun, K.
2012-12-01
Sea-level rising due to climate change following the global warming and the increased intensity of typhoon are magnifying inundation hazards up to the unpredictable level, resulting from the typhoon surge in Korea and other coastal states around the world. Typhoon is the most serious natural disaster in Korean coastal area. Many people died by storm surge inundation every year. And typhoon caused a lot of damage to property. Climate changes due to global warming are producing a stronger natural disaster. Coastal zones have been damaged by typhoons and accompanying storm surge. Especially, the most serious loss of life and terrible property damage caused by typhoon Maemi in 2003. The typhoon Maemi invaded Korean Peninsula leaving property loss of $ 4 Billion and killing 131 people. After then, there has been an increased interest in these coastal zone problems. If storm surges coincide with high tides, the loss of life and property damage due to high waters arc even worse. Therefore it is desirable to accurately forecast the amount water level increase. In this study, using a numerical model FVCOM(finite volume coastal circulation model, Chen et al.,2004), storm surge was simulated to examine its fluctuation characteristics for the coastal area behind Masan, Yeosu and Busan city in Korea. In the numerical model, a moving boundary condition(wet-dry treatment) was incorporated to explain wave inundation. To simulate the inundation scenario, the model grids were extended up to the area inside the lowland in application of the digital elevation data(DEM) made by precisely combining the aero-LiDAR survey data and bathymetry data for the 3 demonstration regions of Busan, Masan and Yeosu. Minimum grid of 300 m unstructured triangular mesh applied to calculate the storm surge was adopted as a grid system. And the minimum grid size of 30 m was built near Busan, Masan and Yeosu area which are the fine coastal regions and where the inundation is simulated. Numerically
Preliminary Results from Numerical Experiments on the Summer 1980 Heat Wave and Drought
NASA Technical Reports Server (NTRS)
Wolfson, N.; Atlas, R.; Sud, Y. C.
1985-01-01
During the summer of 1980, a prolonged heat wave and drought affected the United States. A preliminary set of experiments has been conducted to study the effect of varying boundary conditions on the GLA model simulation of the heat wave. Five 10-day numerical integrations with three different specifications of boundary conditions were carried out: a control experiment which utilized climatological boundary conditions, an SST experiment which utilized summer 1980 sea-surface temperatures in the North Pacific, but climatological values elsewhere, and a Soil Moisture experiment which utilized the values of Mintz-Serafini for the summer, 1980. The starting dates for the five forecasts were 11 June, 7 July, 21 July, 22 August, and 6 September of 1980. These dates were specifically chosen as days when a heat wave was already established in order to investigate the effect of soil moistures or North Pacific sea-surface temperatures on the model's ability to maintain the heat wave pattern. The experiments were evaluated in terms of the heat wave index for the South Plains, North Plains, Great Plains and the entire U.S. In addition a subjective comparison of map patterns has been performed.
Hurtado, Pablo I; Garrido, Pedro L
2010-04-01
Most systems, when pushed out of equilibrium, respond by building up currents of locally conserved observables. Understanding how microscopic dynamics determines the averages and fluctuations of these currents is one of the main open problems in nonequilibrium statistical physics. The additivity principle is a theoretical proposal that allows to compute the current distribution in many one-dimensional nonequilibrium systems. Using simulations, we validate this conjecture in a simple and general model of energy transport, both in the presence of a temperature gradient and in canonical equilibrium. In particular, we show that the current distribution displays a Gaussian regime for small current fluctuations, as prescribed by the central limit theorem, and non-Gaussian (exponential) tails for large current deviations, obeying in all cases the Gallavotti-Cohen fluctuation theorem. In order to facilitate a given current fluctuation, the system adopts a well-defined temperature profile different from that of the steady state and in accordance with the additivity hypothesis predictions. System statistics during a large current fluctuation is independent of the sign of the current, which implies that the optimal profile (as well as higher-order profiles and spatial correlations) are invariant upon current inversion. We also demonstrate that finite-time joint fluctuations of the current and the profile are well described by the additivity functional. These results suggest the additivity hypothesis as a general and powerful tool to compute current distributions in many nonequilibrium systems. PMID:20481672
Three-Dimensional Numerical Modeling of Magnetohydrodynamic Augmented Propulsion Experiment
NASA Technical Reports Server (NTRS)
Turner, M. W.; Hawk, C. W.; Litchford, R. J.
2009-01-01
Over the past several years, NASA Marshall Space Flight Center has engaged in the design and development of an experimental research facility to investigate the use of diagonalized crossed-field magnetohydrodynamic (MHD) accelerators as a possible thrust augmentation device for thermal propulsion systems. In support of this effort, a three-dimensional numerical MHD model has been developed for the purpose of analyzing and optimizing accelerator performance and to aid in understanding critical underlying physical processes and nonideal effects. This Technical Memorandum fully summarizes model development efforts and presents the results of pretest performance optimization analyses. These results indicate that the MHD accelerator should utilize a 45deg diagonalization angle with the applied current evenly distributed over the first five inlet electrode pairs. When powered at 100 A, this configuration is expected to yield a 50% global efficiency with an 80% increase in axial velocity and a 50% increase in centerline total pressure.
Zargarian, A; Esfahanian, M; Kadkhodapour, J; Ziaei-Rad, S
2016-03-01
In this paper, the effects of cell geometry and relative density on the high-cycle fatigue behavior of Titanium scaffolds produced by selective laser melting and electron beam melting techniques were numerically investigated by finite element analysis. The regular titanium lattice samples with three different unit cell geometries, namely, diamond, rhombic dodecahedron and truncated cuboctahedron, and the relative density range of 0.1-0.3 were analyzed under uniaxial cyclic compressive loading. A failure event based algorithm was employed to simulate fatigue failure in the cellular material. Stress-life approach was used to model fatigue failure of both bulk (struts) and cellular material. The predicted fatigue life and the damage pattern of all three structures were found to be in good agreement with the experimental fatigue investigations published in the literature. The results also showed that the relationship between fatigue strength and cycles to failure obeyed the power law. The coefficient of power function was shown to depend on relative density, geometry and fatigue properties of the bulk material while the exponent was only dependent on the fatigue behavior of the bulk material. The results also indicated the failure surface at an angle of 45° to the loading direction. PMID:26706539
NASA Technical Reports Server (NTRS)
Fowlis, W. W. (Editor); Davis, M. H. (Editor)
1981-01-01
The atmospheric general circulation experiment (AGCE) numerical design for Spacelab flights was studied. A spherical baroclinic flow experiment which models the large scale circulations of the Earth's atmosphere was proposed. Gravity is simulated by a radial dielectric body force. The major objective of the AGCE is to study nonlinear baroclinic wave flows in spherical geometry. Numerical models must be developed which accurately predict the basic axisymmetric states and the stability of nonlinear baroclinic wave flows. A three dimensional, fully nonlinear, numerical model and the AGCE based on the complete set of equations is required. Progress in the AGCE numerical design studies program is reported.
Analysis of Numerical Simulation Results of LIPS-200 Lifetime Experiments
NASA Astrophysics Data System (ADS)
Chen, Juanjuan; Zhang, Tianping; Geng, Hai; Jia, Yanhui; Meng, Wei; Wu, Xianming; Sun, Anbang
2016-06-01
Accelerator grid structural and electron backstreaming failures are the most important factors affecting the ion thruster's lifetime. During the thruster's operation, Charge Exchange Xenon (CEX) ions are generated from collisions between plasma and neutral atoms. Those CEX ions grid's barrel and wall frequently, which cause the failures of the grid system. In order to validate whether the 20 cm Lanzhou Ion Propulsion System (LIPS-200) satisfies China's communication satellite platform's application requirement for North-South Station Keeping (NSSK), this study analyzed the measured depth of the pit/groove on the accelerator grid's wall and aperture diameter's variation and estimated the operating lifetime of the ion thruster. Different from the previous method, in this paper, the experimental results after the 5500 h of accumulated operation of the LIPS-200 ion thruster are presented firstly. Then, based on these results, theoretical analysis and numerical calculations were firstly performed to predict the on-orbit lifetime of LIPS-200. The results obtained were more accurate to calculate the reliability and analyze the failure modes of the ion thruster. The results indicated that the predicted lifetime of LIPS-200's was about 13218.1 h which could satisfy the required lifetime requirement of 11000 h very well.
Experiments and numerical simulation of mixing under supercritical conditions
NASA Astrophysics Data System (ADS)
Schmitt, T.; Rodriguez, J.; Leyva, I. A.; Candel, S.
2012-05-01
Supercritical pressure conditions designate a situation where the working fluid pressure is above the critical point. Among these conditions, it is interesting to identify a transcritical range which corresponds to cases where the pressure is above the critical point, but the injection temperature is below the critical value. This situation is of special interest because it raises fundamental issues which have technological relevance in the analysis of flows in liquid rocket engines. This situation is here envisaged by analyzing the behavior of a nitrogen shear coaxial jet comprising an inner stream injected at temperatures close to the critical temperature and a coaxial flow at a higher temperature. Experiments are carried out both in the absence of external modulation and by imposing a large amplitude transverse acoustic field. Real gas large eddy simulations are performed for selected experiments. The combination of experiments and calculations is used to evaluate effects of injector geometry and operating parameters. Calculations retrieve what is observed experimentally when the momentum flux ratio of the outer to the inner stream J= (ρ _eu_e^2)/(ρ _iu_i^2) is varied. Results exhibit the change in flow structure and the development of a recirculation region when this parameter exceeds a critical value. The instantaneous flow patterns for different momentum flux ratios are used in a second stage to characterize the dynamical behavior of the flow in terms of power spectral density of velocity and density fluctuations. Results obtained under acoustic modulation provide insight into mixing enhancement of coaxial streams with a view of its possible consequences in high frequency combustion instabilities. It is shown in particular that the presence of strong acoustic modulations notably reduces the high density jet core length, indicating an increased mixing efficiency. This behavior is more pronounced when the jet is placed at the location of maximum transverse
Numerically Modeling Pulsed-Current, Kinked Wire Experiments
NASA Astrophysics Data System (ADS)
Filbey, Gordon; Kingman, Pat
1999-06-01
The U.S. Army Research Laboratory (ARL) has embarked on a program to provide far-term land fighting vehicles with electromagnetic armor protection. Part of this work seeks to establish robust simulations of magneto-solid-mechanics phenomena. Whether describing violent rupture of a fuse link resulting from a large current pulse or the complete disruption of a copper shaped-charge jet subjected to high current densities, the simulations must include effects of intense Lorentz body forces and rapid Ohmic heating. Material models are required that describe plasticity, flow and fracture, conductivity, and equation of state (EOS) parameters for media in solid, liquid, and vapor phases. An extended version of the Eulerian wave code CTH has been used to predict the apex motion of a V-shaped (``kinked'') copper wire 3mm in diameter during a 400 kilo-amp pulse. These predictions, utilizing available material, EOS, and conductivity data for copper and the known characteristics of an existing capacitor-bank pulsed power supply, were then used to configure an experiment. The experiments were in excellent agreement with the prior simulations. Both computational and experimental results (including electrical data and flash X-rays) will be presented.
NUMERICAL SIMULATIONS OF CONVERGING SHOCKS IN PULSED POWER DRIVEN EXPERIMENTS
R. KANZLEITER; W. ATCHISON; ET AL
2000-12-01
The final shot of the current Near Term Liner Experiment (NTLX) series occurred on September 29, 2000. Utilization of a pulsed power source with a standardized liner/target ''cartridge'' produced a uniform implosion to drive hydrodynamic experiments. Diagnostics showed that high quality data of shock propagation can be obtained from pulsed power liner drivers as in the current NTLX series. Very good agreement in calculating shock locations was obtained between the codes used to model the NTLX series, RAGE and RAVEN. RAVEN also accurately predicts liner/target impact as measured by B-Dot probes. Large differences are observed between the calculated and measured positions of converging shock waves even in simple geometrical configurations. Liner/target impact is accurately calculated and similar results are produced for shock velocities in Lucite. RAGE and RAVEN use different hydrodynamic algorithms, yet agree, this focuses current efforts on EOS issues within the outer tin target to resolve discrepancies. Further diagnostics covering shock breakout from the outer tin target and shock propagation shortly thereafter would be highly beneficial.
Numerical Simulation of Receptivity for a Transition Experiment
NASA Technical Reports Server (NTRS)
Collis, S. Scott; Joslin, R. D. (Technical Monitor)
2000-01-01
The cost of fuel to overcome turbulence induced viscous drag on a commercial airplane constitutes a significant fraction of the operating cost of an airline. Achieving laminar flow and maintaining it over a large portion of the wing can significantly reduce the viscous drag, and hence the cost. Design of such laminar-flow-control wings and their practical operation requires the ability to accurately and reliably predict the transition from laminar to turbulent flow. The transition process begins with the conversion of environmental and surface disturbances into the instability waves of the flow by a process called receptivity. The goal of the current research project has been to improve the prediction of transition through a better understanding of the physics of receptivity. The initial objective of this work was to investigate the specific stability and receptivity characteristics of a particular experimental investigation of boundary layer receptivity at NASA Langley. Some simulation results using direct solutions of the linearized Navier-Stokes equations which modeled this experiment where presented in the 1999 APS DFD meeting. However, based on these initial investigations, it became clear that to cover the vast receptivity parameter space required for a practical transition prediction tool, more efficient methods would be required. Thus, the focus of this research was shifted from modeling this particular experiment to formulating and developing new techniques that could efficiently yet accurately predict receptivity for a wide range of disturbance conditions.
Numerical controlled polishing, continued force wear and part correction experiments
Hannah, P.R.; Day, R.D.; Hatch, D.J.; McClure, E.R.
1994-09-01
This abstract reports the near completion of the first phase of this program. It is the aim of this program to provide the operator of a N/C diamond turning machine or N/C grinding machine (jig grinder) with the wear characteristics necessary to achieve uniform material removal. The second phase of this program addresses a different problem, although solving this problem is highly dependent on the results of the first phase. Diamond turned, or any lathe turned surface, exhibits regular tool marks due to the tool passing over the surface being cut. Changes in depth of cut, feed rate and work rpm will change the character of these groves, but will not eliminate them. Optical surfaces produced by this process exhibit increased scattering as the light wavelength decreases limiting their use; at least for optical purposes, to IR and some visible applications. Utilizing wear information gathered in the first part of this program we will attempt to reduce these residual tool marks by polishing. The polishing of diamond turned surfaces is not new. Diamond turned metal surfaces, especially in electroless nickel and high phosphorus nickel electroplate have been polished to improve their scatter characteristics. What we believe is unique is the use of a spherical wheel, rotating on axis and being moved over the part in a prescribed manner by numerical control. Over the past year we have made some major changes in our polishing methods and procedures. We have listed below these changes, as a refresher for the reader as to our previous procedures. These changes will be addressed in the body of the text.
Numerical Simulations of the Boundary Layer Transition Flight Experiment
NASA Technical Reports Server (NTRS)
Tang, Chun Y.; Trumble, Kerry A.; Campbell, Charles H.; Lessard, Victor R.; Wood, William A.
2010-01-01
Computational Fluid Dynamics (CFD) simulations were used to study the possible effects that the Boundary Layer Transition (BLT) Flight Experiments may have on the heating environment of the Space Shuttle during its entry to Earth. To investigate this issue, hypersonic calculations using the Data-Parallel Line Relaxation (DPLR) and Langley Aerothermodynamic Upwind Relaxation (LAURA) CFD codes were computed for a 0.75 tall protuberance at flight conditions of Mach 15 and 18. These initial results showed high surface heating on the BLT trip and the areas surrounding the protuberance. Since the predicted peak heating rates would exceed the thermal limits of the materials selected to construct the BLT trip, many changes to the geometry were attempted in order to reduce the surface heat flux. The following paper describes the various geometry revisions and the resulting heating environments predicted by the CFD codes.
Numerical study of the effect of water addition on gas explosion.
Liang, Yuntao; Zeng, Wen
2010-02-15
Through amending the SENKIN code of CHEMKIN III chemical kinetics package, a computational model of gas explosion in a constant volume bomb was built, and the detailed reaction mechanism (GRI-Mech 3.0) was adopted. The mole fraction profiles of reactants, some selected free radicals and catastrophic gases in the process of gas explosion were analyzed by this model. Furthermore, through the sensitivity analysis of the reaction mechanism of gas explosion, the dominant reactions that affect gas explosion and the formation of catastrophic gases were found out. At the same time, the inhibition mechanisms of water on gas explosion and the formation of catastrophic gases were analyzed. The results show that the induced explosion time is prolonged, and the mole fractions of reactant species such as CH(4), O(2) and catastrophic gases such as CO, CO(2) and NO are decreased as water is added to the mixed gas. With the water fraction in the mixed gas increasing, the sensitivities of the dominant reactions contributing to CH(4), CO(2) are decreased and the sensitivity coefficients of CH(4), CO and NO mole fractions are also decreased. The inhibition of gas explosion with water addition can be ascribed to the significant decrease of H, O and OH in the process of gas explosion due to the water presence. PMID:19811873
Resistivity inversion in 2-D anisotropic media: numerical experiments
NASA Astrophysics Data System (ADS)
Wiese, Timothy; Greenhalgh, Stewart; Zhou, Bing; Greenhalgh, Mark; Marescot, Laurent
2015-04-01
Many rocks and layered/fractured sequences have a clearly expressed electrical anisotropy although it is rare in practice to incorporate anisotropy into resistivity inversion. In this contribution, we present a series of 2.5-D synthetic inversion experiments for various electrode configurations and 2-D anisotropic models. We examine and compare the image reconstructions obtained using the correct anisotropic inversion code with those obtained using the false but widely used isotropic assumption. Superior reconstruction in terms of reduced data misfit, true anomaly shape and position, and anisotropic background parameters were obtained when the correct anisotropic assumption was employed for medium to high coefficients of anisotropy. However, for low coefficient values the isotropic assumption produced better-quality results. When an erroneous isotropic inversion is performed on medium to high level anisotropic data, the images are dominated by patterns of banded artefacts and high data misfits. Various pole-pole, pole-dipole and dipole-dipole data sets were investigated and evaluated for the accuracy of the inversion result. The eigenvalue spectra of the pseudo-Hessian matrix and the formal resolution matrix were also computed to determine the information content and goodness of the results. We also present a data selection strategy based on high sensitivity measurements which drastically reduces the number of data to be inverted but still produces comparable results to that of the comprehensive data set. Inversion was carried out using transversely isotropic model parameters described in two different co-ordinate frames for the conductivity tensor, namely Cartesian versus natural or eigenframe. The Cartesian frame provided a more stable inversion product. This can be simply explained from inspection of the eigenspectra of the pseudo-Hessian matrix for the two model descriptions.
Numerical simulation of the experiment of electrical explosion of aluminum foil
NASA Astrophysics Data System (ADS)
Shutov, A. V.
2015-11-01
Numerical simulation of the experiment of Korobenko et al (2007 Phys. Rev. B 75 064208) in strongly coupled plasma of aluminum have been fulfilled. The results of numerical simulation and the experiment are compared. It is established that the hydrodynamic flows in the experiment can be assumed one-dimensional. The elastic-plastic effects in the dynamics of aluminum foil are also insignificant. The focus in the modeling is devoted to the study of the dynamics of the thermodynamic states of aluminum and their spatial homogeneity. It is emphasized the strong influence of the thermal conductivity for such experiments.
NASA Technical Reports Server (NTRS)
Duque, Earl P. N.; Johnson, Wayne; vanDam, C. P.; Chao, David D.; Cortes, Regina; Yee, Karen
1999-01-01
Accurate, reliable and robust numerical predictions of wind turbine rotor power remain a challenge to the wind energy industry. The literature reports various methods that compare predictions to experiments. The methods vary from Blade Element Momentum Theory (BEM), Vortex Lattice (VL), to variants of Reynolds-averaged Navier-Stokes (RaNS). The BEM and VL methods consistently show discrepancies in predicting rotor power at higher wind speeds mainly due to inadequacies with inboard stall and stall delay models. The RaNS methodologies show promise in predicting blade stall. However, inaccurate rotor vortex wake convection, boundary layer turbulence modeling and grid resolution has limited their accuracy. In addition, the inherently unsteady stalled flow conditions become computationally expensive for even the best endowed research labs. Although numerical power predictions have been compared to experiment. The availability of good wind turbine data sufficient for code validation experimental data that has been extracted from the IEA Annex XIV download site for the NREL Combined Experiment phase II and phase IV rotor. In addition, the comparisons will show data that has been further reduced into steady wind and zero yaw conditions suitable for comparisons to "steady wind" rotor power predictions. In summary, the paper will present and discuss the capabilities and limitations of the three numerical methods and make available a database of experimental data suitable to help other numerical methods practitioners validate their own work.
Naghieh, S; Karamooz Ravari, M R; Badrossamay, M; Foroozmehr, E; Kadkhodaei, M
2016-06-01
In recent years, thanks to additive manufacturing technology, researchers have gone towards the optimization of bone scaffolds for the bone reconstruction. Bone scaffolds should have appropriate biological as well as mechanical properties in order to play a decisive role in bone healing. Since the fabrication of scaffolds is time consuming and expensive, numerical methods are often utilized to simulate their mechanical properties in order to find a nearly optimum one. Finite element analysis is one of the most common numerical methods that is used in this regard. In this paper, a parametric finite element model is developed to assess the effects of layers penetration׳s effect on inter-layer adhesion, which is reflected on the mechanical properties of bone scaffolds. To be able to validate this model, some compression test specimens as well as bone scaffolds are fabricated with biocompatible and biodegradable poly lactic acid using fused deposition modeling. All these specimens are tested in compression and their elastic modulus is obtained. Using the material parameters of the compression test specimens, the finite element analysis of the bone scaffold is performed. The obtained elastic modulus is compared with experiment indicating a good agreement. Accordingly, the proposed finite element model is able to predict the mechanical behavior of fabricated bone scaffolds accurately. In addition, the effect of post-heating of bone scaffolds on their elastic modulus is investigated. The results demonstrate that the numerically predicted elastic modulus of scaffold is closer to experimental outcomes in comparison with as-built samples. PMID:26874065
CO2 permeability of fractured cap rocks - experiments and numerical simulations (CO2Seals)
NASA Astrophysics Data System (ADS)
(Draeger), Ines Rick; Clauser, Christoph
2010-05-01
In CO2 sequestration and underground gas storage the sealing capacity of a cap rock is of paramount importance. The main question is therefore how the leakage of CO2 through fissures and faults within the cap rock may affect the CO2 sealing efficiency of low-permeable seal lithotypes. In many cases, these structures provide the main pathways for leakage of CO2. Here, we provide an overview of one part of the joint research project CO2Seals, which deals with the effect of structural features - such as tectonic faults and fissures in the overburden - on the migration of CO2 in addition to mineralogical, petrophysical, and geochemical properties of different lithotypes. The primary contribution of the entire project consists of an improvement of the present quantitative understanding of CO2 transport and retention processes and associated interactions in cap rocks between rock and CO2 or brine. To this end, we are adapting different numerical tools for simulating the relevant petrophysical and geochemical processes of CO2 in cap rocks, in close operation with: (1) large-scale CO2-percolation experiments on fractured cap rock samples; (2) permeability, gas breakthrough, and diffusion experiments; (3) measurements of the mechanical stability of cap rocks and the geochemical alterations of fault zone rock. The observed resulting changes in petrophysical properties, such as porosity, relative rock permeability (CO2 and brine), and fault permeability provide basics for the following numerical simulations. For example, first permeability tests of a marl and clay cap rock out of Cretaceous and Jurassic formations revealed gas permeability of 10-18 m2 down to 10-22 m2. In addition, first percolation experiments indicated that the influence of fault zones on the measured CO2 permeability of clays is very low. Furthermore, numerical bench-scale models are performed to provide confidence for the subsequent transfer to reservoir systems. Large-scale numerical models were created
NASA Astrophysics Data System (ADS)
Xia, Daqing; Xu, Youping
1998-06-01
In first paper of articles, the physical and calculating schemes of the water-bearing numerical model are described. The model is developed by bearing all species of hydrometeors in a conventional numerical model in which the dynamic framework of hydrostatic equilibrium is taken. The main contributions are: the mixing ratios of all species of hydrometeors are added as the prognostic variables of model, the prognostic equations of these hydrometeors are introduced, the cloud physical framework is specially designed, some technical measures are used to resolve a series of physical, mathematical and computational problems arising from water-bearing; and so on. The various problems (in such aspects as the designs of physical and calculating schemes and the composition of computational programme) which are exposed in feasibility test, in sensibility test, and especially in operational forecasting experiments are successfully resolved using a lot of technical measures having been developed from researches and tests. Finally, the operational forecasting running of the water-bearing numerical model and its forecasting system is realized stably and reliably, and the fine forecasts are obtained. All of these mentioned above will be described in second paper.
ERIC Educational Resources Information Center
Iler, H. Darrell; Brown, Amber; Landis, Amanda; Schimke, Greg; Peters, George
2014-01-01
A numerical analysis of the free radical addition polymerization system is described that provides those teaching polymer, physical, or advanced organic chemistry courses the opportunity to introduce students to numerical methods in the context of a simple but mathematically stiff chemical kinetic system. Numerical analysis can lead students to an…
NASA Astrophysics Data System (ADS)
Papafotiou, Alexandros; Ganz, Christina; Altfelder, Sven; Noell, Ursula; Neuweiler, Insa
2010-05-01
The unsaturated zone has a prominent role for groundwater resources, as it controls through flow and transport any mass exchange between atmosphere and groundwater. However, providing reliable predictions for the unsaturated zone is very demanding, as it is dominated by complex two-phase flow processes that produce high uncertainty with respect to the hydraulic properties. When modeling unsaturated flow, the typically unknown spatial distribution of hydraulic properties in the soil constitutes a primary source of uncertainty. Even if information on the exact distribution is known, additional uncertainty may stem from the non-uniqueness of the hydraulic properties, most profoundly expressed through hysteresis in the capillary pressure-saturation relationship, also known as water retention curve. In this work, we present modeling considerations for predicting an infiltration and drainage event in the unsaturated zone during a field experiment. The experiment was performed by infiltrating brilliant-blue solution while monitoring the plume movement with ERT. After the completion of infiltration (and the consequent drainage), the upper 1 meter of the soil was excavated in slices to obtain the 3D distribution of water saturation and pressure. Numerical simulations are carried out with a two-phase flow model. The results illustrate possibilities and limitations of predicting such flow processes based on the experimental information available. We demonstrate the influence and significance of hysteresis by comparing experimental findings with model runs that explicitly consider wetting and drying conditions in the experiment. Our approach allows us to identify key processes that have to be accounted for. In a feedback loop with the design of future experiments we aim at improving input specifications necessary for reliable predictive modeling of unsaturated flow.
ERIC Educational Resources Information Center
Van Atta, Robert E.; Van Atta, R. Lewis
1980-01-01
Provides a gas chromatography experiment that exercises the quantitative technique of standard addition to the analysis for a minor component, methyl salicylate, in a commercial product, "wintergreen rubbing alcohol." (CS)
NASA Astrophysics Data System (ADS)
Volk, R.; Calzavarini, E.; Verhille, G.; Lohse, D.; Mordant, N.; Pinton, J.-F.; Toschi, F.
2008-08-01
We compare experimental data and numerical simulations for the dynamics of inertial particles with finite density in turbulence. In the experiment, bubbles and solid particles are optically tracked in a turbulent flow of water using an Extended Laser Doppler Velocimetry technique. The probability density functions (PDF) of particle accelerations and their auto-correlation in time are computed. Numerical results are obtained from a direct numerical simulation in which a suspension of passive pointwise particles is tracked, with the same finite density and the same response time as in the experiment. We observe a good agreement for both the variance of acceleration and the autocorrelation time scale of the dynamics; small discrepancies on the shape of the acceleration PDF are observed. We discuss the effects induced by the finite size of the particles, not taken into account in the present numerical simulations.
Voices from the Classroom: Experiences of Teachers of Deaf Students with Additional Disabilities
ERIC Educational Resources Information Center
Musyoka, Millicent Malinda; Gentry, Mary Anne; Bartlett, James Joseph
2016-01-01
The purpose of this research is to investigate experiences of K-12 classroom teachers of deaf students with additional disabilities. Today, more deaf and hard of hearing students are identified as having additional disabilities (Bruce, DiNatale & Ford, 2008; Ewing, 2011; Gallaudet Research Institute, 2011; Jones, Jones & Ewing, 2006;…
NASA Astrophysics Data System (ADS)
Kikuchi, T.; Sakai, Y.; Komori, T.; Sato, T.; Hasegawa, J.; Horioka, K.; Takahashi, K.; Sasaki, T.; Harada, Nob
2016-05-01
Tune depression in a compact beam equipment is estimated, and numerical simulation results are compared with an experimental one for the compact beam simulator in a driver of heavy ion inertial fusion. The numerical simulation with multi-particle tracking is carried out, corresponding to the experimental condition, and the result is discussed with the experimental one. It is expected that the numerical simulation developed in this paper is useful tool to investigate the beam dynamics in the experiment with the compact beam simulator.
NASA Astrophysics Data System (ADS)
Sturz, L.; Wu, M.; Zimmermann, G.; Ludwig, A.; Ahmadein, M.
2015-06-01
Solidification benchmark experiments on columnar and equiaxed dendritic growth, as well as the columnar-equiaxed transition have been carried out under diffusion-dominated conditions for heat and mass transfer in a low-gravity environment. The system under investigation is the transparent organic alloy system Neopentylglycol-37.5wt.-%(d)Camphor, processed aboard a TEXUS sounding rocket flight. Solidifications was observed by standard optical methods in addition to measurements of the thermal fields within the sheet like experimental cells of 1 mm thickness. The dendrite tip kinetic, primary dendrite arm spacing, temporal and spatial temperature evolution, columnar tip velocity and the critical parameters at the CET have been analysed. Here we focus on a detailed comparison of the experiment “TRACE” with a 5-phase volume averaging model to validate the numerical model and to give insight into the corresponding physical mechanisms and parameters leading to CET. The results are discussed in terms of sensitivity versus numerical parameters.
Response of non-added solutes during nutrient addition experiments in streams
NASA Astrophysics Data System (ADS)
Rodriguez-Cardona, B.; Wymore, A.; Koenig, L.; Coble, A. A.; McDowell, W. H.
2015-12-01
Nutrient addition experiments, such as Tracer Additions for Spiraling Curve Characterization (TASCC), have become widely popular as a means to study nutrient uptake dynamics in stream ecosystems. However, the impact of these additions on ambient concentrations of non-added solutes is often overlooked. TASCC addition experiments are ideal for assessing interactions among solutes because it allows for the characterization of multiple solute concentrations across a broad range of added nutrient concentrations. TASCC additions also require the addition of a conservative tracer (NaCl) to track changes in conductivity during the experimental manipulation. Despite its use as a conservative tracer, chloride (Cl) and its associated sodium (Na) might change the concentrations of other ions and non-added nutrients through ion exchange or other processes. Similarly, additions of biologically active solutes might change the concentrations of other non-added solutes. These methodological issues in nutrient addition experiments have been poorly addressed in the literature. Here we examine the response of non-added solutes to pulse additions (i.e. TASCC) of NaCl plus nitrate (NO3-), ammonium, and phosphate across biomes including temperate and tropical forests, and arctic taiga. Preliminary results demonstrate that non-added solutes respond to changes in the concentration of these added nutrients. For example, concentrations of dissolved organic nitrogen (DON) in suburban headwater streams of New Hampshire both increase and decrease in response to NO3- additions, apparently due to biotic processes. Similarly, cations such as potassium, magnesium, and calcium also increase during TASCC experiments, likely due to cation exchange processes associated with Na addition. The response of non-added solutes to short-term pulses of added nutrients and tracers needs to be carefully assessed to ensure that nutrient uptake metrics are accurate, and to detect biotic interactions that may
NASA Technical Reports Server (NTRS)
Kwon, J. H.
1977-01-01
Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.
Numerical simulation studies of the LBNL heavy-ion beam combiner experiment
Fawley, W.M.; Seidl, P.; Haber, I.; Friedman, A.; Grote, D.P.
1997-01-01
Transverse beam combining is a cost-saving option employed in many designs for heavy-ion inertial fusion energy drivers. A major area of interest, both theoretically and experimentally, is the resultant transverse phase space dilution during the beam merging process. Currently, a prototype combining experiment is underway at LBNL and we have employed a variety of numerical descriptions to aid in both the initial design of the experiment data. These range from simple envelope codes to detailed 2- and 3-D PIC simulations. We compare the predictions of the different numerical models to each other and to experimental data at different longitudinal positions.
Influence of Ar addition on ozone generation in a non-thermal plasma—a numerical investigation
NASA Astrophysics Data System (ADS)
Chen, Hsin Liang; Lee, How Ming; Chen, Shiaw Huei; Wei, Ta Chin; Been Chang, Moo
2010-10-01
A numerical model based on a dielectric barrier discharge is developed in this study to investigate the influence of Ar addition on ozone generation. The simulation results show good agreement with the experimental data, confirming the validity of the numerical model. The mechanisms regarding how the Ar addition affects ozone generation are investigated with the assistance of a numerical simulation by probing into the following two questions, (1) why the ozone concentration just slightly decreases in the low specific input energy (SIE, the ratio of discharge power to gas flow rate) region even if the inlet O2 concentration is substantially decreased and (2) why the variation of the increased rate of ozone concentration with SIE (i.e. the variation in the slope of ozone concentration versus SIE) is more significant for an O2/Ar mixture plasma. As SIE is relatively low, ozone decomposition through electron-impact and radical attack reactions is less significant because of low ozone concentration and gas temperature. Therefore, the ozone concentration depends mainly on the amount of oxygen atoms generated. The simulation results indicate that the amount of oxygen atoms generated per electronvolt for Ar concentrations of 0%, 10%, 30%, 50% and 80% are 0.178, 0.174, 0.169, 0.165 and 0.166, respectively, explaining why the ozone concentration does not decrease linearly with the inlet O2 concentration in the low SIE region. On the other hand, the simulation results show that increasing Ar concentration would lead to a lower reduced field and a higher gas temperature. The former would lead to an increase in the rate constant of e + O3 → e + O + O2 while the latter would result in a decrease in the rate constant of O + O2 + M → O3 + M and an increase in that of O3 + O → 2O2. The changes in the rate constants of these reactions would have a negative effect on ozone generation, which is the rationale for the second question.
NASA Astrophysics Data System (ADS)
Mininni, P.; Dmitruk, P.; Odier, P.; Pinton, J.-F.; Plihon, N.; Verhille, G.; Volk, R.; Bourgoin, M.
2014-05-01
We analyze time series stemming from experiments and direct numerical simulations of hydrodynamic and magnetohydrodynamic turbulence. Simulations are done in periodic boxes, but with a volumetric forcing chosen to mimic the geometry of the flow in the experiments, the von Kármán swirling flow between two counterrotating impellers. Parameters in the simulations are chosen to (within computational limitations) allow comparisons between the experiments and the numerical results. Conducting fluids are considered in all cases. Two different configurations are considered: a case with a weak externally imposed magnetic field and a case with self-sustained magnetic fields. Evidence of long-term memory and 1/f noise is observed in experiments and simulations, in the case with weak magnetic field associated with the hydrodynamic behavior of the shear layer in the von Kármán flow, and in the dynamo case associated with slow magnetohydrodynamic behavior of the large-scale magnetic field.
Test particle acceleration in a numerical MHD experiment of an anemone jet
NASA Astrophysics Data System (ADS)
Rosdahl, K. J.; Galsgaard, K.
2010-02-01
Aims: To use a 3D numerical MHD experiment representing magnetic flux emerging into an open field region as a background field for tracing charged particles. The interaction between the two flux systems generates a localised current sheet where MHD reconnection takes place. We investigate how efficiently the reconnection region accelerates charged particles and what kind of energy distribution they acquire. Methods: The particle tracing is done numerically using the Guiding Center Approximation on individual data sets from the numerical MHD experiment. Results: We derive particle and implied photon distribution functions having power law forms, and look at the impact patterns of particles hitting the photosphere. We find that particles reach energies far in excess of those seen in observations of solar flares. However the structure of the impact region in the photosphere gives a good representation of the topological structure of the magnetic field. Three movies are only available in electronic form at http://www.aanda.org
NASA Astrophysics Data System (ADS)
Degtyarev, Alexander; Khramushin, Vasily
2016-02-01
The paper deals with the computer implementation of direct computational experiments in fluid mechanics, constructed on the basis of the approach developed by the authors. The proposed approach allows the use of explicit numerical scheme, which is an important condition for increasing the effciency of the algorithms developed by numerical procedures with natural parallelism. The paper examines the main objects and operations that let you manage computational experiments and monitor the status of the computation process. Special attention is given to a) realization of tensor representations of numerical schemes for direct simulation; b) realization of representation of large particles of a continuous medium motion in two coordinate systems (global and mobile); c) computing operations in the projections of coordinate systems, direct and inverse transformation in these systems. Particular attention is paid to the use of hardware and software of modern computer systems.
Tsebers, A.O.
1985-04-01
This paper is an attempt to simulate the magnetorheological behavior of a suspension of ellipsoidal ferromagnetic particles in a numerical experiment. Accuracy of the calculations used are achieved and illustrated in the paper. It is shown that the relative error in the calculation of the characteristic viscosity does not exceed 5%.
NASA Technical Reports Server (NTRS)
Shih, Hsin-Yi; Tien, James S.; Ferkul, Paul (Technical Monitor)
2001-01-01
The recently developed numerical model of concurrent-flow flame spread over thin solids has been used as a simulation tool to help the designs of a space experiment. The two-dimensional and three-dimensional, steady form of the compressible Navier-Stokes equations with chemical reactions are solved. With the coupled multi-dimensional solver of the radiative heat transfer, the model is capable of answering a number of questions regarding the experiment concept and the hardware designs. In this paper, the capabilities of the numerical model are demonstrated by providing the guidance for several experimental designing issues. The test matrix and operating conditions of the experiment are estimated through the modeling results. The three-dimensional calculations are made to simulate the flame-spreading experiment with realistic hardware configuration. The computed detailed flame structures provide the insight to the data collection. In addition, the heating load and the requirements of the product exhaust cleanup for the flow tunnel are estimated with the model. We anticipate that using this simulation tool will enable a more efficient and successful space experiment to be conducted.
NASA Astrophysics Data System (ADS)
Brune, S.; Czaplinska, D.; Piazolo, S.; Wilson, C. J. L.; Quinteros, J.
2015-12-01
Most numerical models of lithosphere deformation approximate the rheological behavior of polymineralic crust and mantle via single-phase flow laws assuming that the weakest or most abundant material controls the bulk rheology. However, previous work showed that in two phase aggregates the bulk viscosity of the dominant phase is significantly affected by second phase particles. Here we combine two unconventional approaches to quantify the relative impact of such particles on strain localisation and bulk response: (1) We run centimetre-scale numerical models of a matrix with inclusions using the elasto-visco-plastic FEM software Slim3D. Recrystallization-induced weakening processes in the matrix, i.e. grain boundary migration and nucleation, are approximated using strain-dependent viscous softening. (2) We conduct high T, constant strain rate deformation experiments with a matrix of deuterated ice (D2O) containing rigid or soft particles, i.e. calcite and graphite, respectively. Ice is a valuable rock analogue, as it replicates the microstructural and fabric changes as well as the non-Newtonian response of other anisotropic minerals, such as olivine and quartz. The laboratory experiments exhibit two types of rheological behaviour: stress partitioning between ice and particles and strain localization in rheologically softer material. To quantify the contribution of both response types, we calibrate numerical simulations with data derived from laboratory experiments. The strain rate, stress, and viscosity evolution of the numerical experiment provides insight to non-linear strain localization processes, particle motion and time-dependent stress concentrations during the deformation. We fit the parameters of the viscous softening function and thereby quantify the amount of additional weakening in the matrix of ice mixtures in comparison to pure ice, which allows to constrain softening parameters used in large-scale simulations of glacial flow and lithosphere deformation.
Finsterle, S.; Moridis, G.J.; Pruess, K.; Persoff, P.
1994-01-01
The emplacement of liquids under controlled viscosity conditions is investigated by means of numerical simulations. Design calculations are performed for a laboratory experiment on a decimeter scale, and a field experiment on a meter scale. The purpose of the laboratory experiment is to study the behavior of multiple gout plumes when injected in a porous medium. The calculations for the field trial aim at designing a grout injection test from a vertical well in order to create a grout plume of a significant extent in the subsurface.
NASA Astrophysics Data System (ADS)
Řidký, V.; Šidlof, P.; Vlček, V.
2013-04-01
The work is devoted to comparing measured data with the results of numerical simulations. As mathematical model was used mathematical model whitout turbulence for incompressible flow In the experiment was observed the behavior of designed NACA0015 airfoil in airflow. For the numerical solution was used OpenFOAM computational package, this is open-source software based on finite volume method. In the numerical solution is prescribed displacement of the airfoil, which corresponds to the experiment. The velocity at a point close to the airfoil surface is compared with the experimental data obtained from interferographic measurements of the velocity field. Numerical solution is computed on a 3D mesh composed of about 1 million ortogonal hexahedron elements. The time step is limited by the Courant number. Parallel computations are run on supercomputers of the CIV at Technical University in Prague (HAL and FOX) and on a computer cluster of the Faculty of Mechatronics of Liberec (HYDRA). Run time is fixed at five periods, the results from the fifth periods and average value for all periods are then be compared with experiment.
Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment
NASA Technical Reports Server (NTRS)
Cambier, Jean-Luc
1991-01-01
A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.
Yates, Kimberly K.; DuFore, Christopher M.; Robbins, Lisa L.
2013-01-01
Use of different approaches for manipulating seawater chemistry during ocean acidification experiments has confounded comparison of results from various experimental studies. Some of these discrepancies have been attributed to whether addition of acid (such as hydrochloric acid, HCl) or carbon dioxide (CO2) gas has been used to adjust carbonate system parameters. Experimental simulations of carbonate system parameter scenarios for the years 1766, 2007, and 2100 were performed using the carbonate speciation program CO2SYS to demonstrate the variation in seawater chemistry that can result from use of these approaches. Results showed that carbonate system parameters were 3 percent and 8 percent lower than target values in closed-system acid additions, and 1 percent and 5 percent higher in closed-system CO2 additions for the 2007 and 2100 simulations, respectively. Open-system simulations showed that carbonate system parameters can deviate by up to 52 percent to 70 percent from target values in both acid addition and CO2 addition experiments. Results from simulations for the year 2100 were applied to empirically derived equations that relate biogenic calcification to carbonate system parameters for calcifying marine organisms including coccolithophores, corals, and foraminifera. Calculated calcification rates for coccolithophores, corals, and foraminifera differed from rates at target conditions by 0.5 percent to 2.5 percent in closed-system CO2 gas additions, from 0.8 percent to 15 percent in the closed-system acid additions, from 4.8 percent to 94 percent in open-system acid additions, and from 7 percent to 142 percent in open-system CO2 additions.
Rider, William; Kamm, J. R.; Zoldi, C. A.; Tomkins, C. D.
2002-01-01
We present detailed spatial analysis comparing experimental data and numerical simulation results for Richtmyer-Meshkov instability experiments of Prestridge et al. and Tomkins et al. These experiments consist, respectively, of one and two diffuse cylinders of sulphur hexafluoride (SF{sub 6}) impulsively accelerated by a Mach 1.2 shockwave in air. The subsequent fluid evolution and mixing is driven by the deposition of baroclinic vorticity at the interface between the two fluids. Numerical simulations of these experiments are performed with three different versions of high resolution finite volume Godunov methods, including a new weighted adaptive Runge-Kutta (WARK) scheme. We quantify the nature of the mixing using using integral measures as well as fractal analysis and continuous wavelet transforms. Our investigation of the gas cylinder configurations follows the path of our earlier studies of the geometrically and dynamically more complex gas 'curtain' experiment. In those studies, we found significant discrepancies in the details of the experimentally measured mixing and the details of the numerical simulations. Here we evaluate the effects of these hydrodynamic integration techniques on the diffuse gas cylinder simulations, which we quantitatively compare with experimental data.
NASA Technical Reports Server (NTRS)
Metcalfe, R. W.; Menon, S.; Hussain, A. K. M. F.
1985-01-01
An eduction scheme has been developed in an attempt to determine the characteristics of large-scale vortical structures in a turbulent mixing layer. This analysis scheme has been applied to a set of experimental data taken in a new, larger mixing layer facility designed to minimize boundary and resonance effects. A similar scheme has been developed to apply to the results of a direct numerical simulation of a temporally growing mixing layer. A comparison of the two approaches shows important similarities in the coherent structures. The numerical simulations indicate that low levels of coherent forcing can dramatically change the evolution of the mixing layer. In the absence of such forcing, the numerical simulations and experiments show a lack of regularity in the transverse position, spacing, amplitude, shape and spanwise coherence of the large-scale vortical structures.
Kirsch, Louise P.; Cross, Emily S.
2015-01-01
The way in which we perceive others in action is biased by one's prior experience with an observed action. For example, we can have auditory, visual, or motor experience with actions we observe others perform. How action experience via 1, 2, or all 3 of these modalities shapes action perception remains unclear. Here, we combine pre- and post-training functional magnetic resonance imaging measures with a dance training manipulation to address how building experience (from auditory to audiovisual to audiovisual plus motor) with a complex action shapes subsequent action perception. Results indicate that layering experience across these 3 modalities activates a number of sensorimotor cortical regions associated with the action observation network (AON) in such a way that the more modalities through which one experiences an action, the greater the response is within these AON regions during action perception. Moreover, a correlation between left premotor activity and participants' scores for reproducing an action suggests that the better an observer can perform an observed action, the stronger the neural response is. The findings suggest that the number of modalities through which an observer experiences an action impacts AON activity additively, and that premotor cortical activity might serve as an index of embodiment during action observation. PMID:26209850
NASA Astrophysics Data System (ADS)
Spina, L.; Colucci, S.; De'Michieli Vitturi, M.; Scheu, B.; Dingwell, D. B.
2014-12-01
Numerical modeling, joined with experimental investigations, is fundamental for studying the dynamics of magmatic fluid into the conduit, where direct observations are unattainable. Furthermore, laboratory experiments can provide invaluable data to vunalidate complex multiphase codes. With the aim on unveil the essence of nucleation process, as well as the behavior of the multiphase magmatic fluid, we performed slow decompression experiments in a shock tube system. We choose silicon oil as analogue for the magmatic melt, and saturated it with Argon at 10 MPa for 72h. The slow decompression to atmospheric conditions was monitored through a high speed camera and pressure sensors, located into the experimental conduit. The experimental conditions of the decompression process have then been reproduced numerically with a compressible multiphase solver based on OpenFOAM. Numerical simulations have been performed by the OpenFOAM compressibleInterFoam solver for 2 compressible, non-isothermal immiscible fluids, using a VOF (volume of fluid) phase-fraction based interface capturing approach. The data extracted from 2D images obtained from laboratory analyses were compared to the outcome of numerical investigation, showing the capability of the model to capture the main processes studied.
NASA Astrophysics Data System (ADS)
Christensen, Amalie; Raufaste, Christophe; Misztal, Marek; Celestini, Franck; Guidi, Maria; Ellegaard, Clive; Mathiesen, Joachim
2016-03-01
Many natural fracture systems are characterized by a single length scale, which is the distance between neighboring fractures. Examples are mud cracks and columnar jointing. In columnar jointing the origin of this scale has been a long-standing issue. Here we present a comprehensive study of columnar jointing based on experiments on cooling stearic acids, numerical simulations using both discrete and finite element methods and basic analytical calculations. We show that the diameter of columnar joints is a nontrivial function of the material properties and the cooling conditions of the system. We determine the shape of this function analytically and show that it is in agreement with the experiments and the numerical simulations.
Validation with experiments on simplified numerical prediction of hybrid rocket internal ballistics
NASA Astrophysics Data System (ADS)
Funami, Yuki; Shimada, Toru
2012-11-01
In order to design hybrid rocket engines, we have developed a numerical prediction approach to the internal ballistics. The key point is its cost performance. Therefore simple but efficient models are required. Fluid phenomenon and thermal conduction phenomenon in a solid fuel should be treated time-dependently, because characteristic times of these phenomena are longer than those of other phenomena. Besides, they are solved with the energy-flux balance equation at the solid fuel surface to determine the regression rate. It is confirmed that numerical evaluation of time- and space-averaged regression rate is the same order of magnitude as that in experiments. However, the factors n in ṙ¯ = aG¯oxn differ between calculations and experiments.
Lefrancois, A; Vandersall, K; L'Eplattenier, P; Burger, M
2006-02-06
Isentropic compression experiments and numerical simulations on TATB based HE were performed respectively at Z accelerator facility from Sandia National Laboratory and at Lawrence Livermore National Laboratory in order to study the isentrope and associated Hugoniot of this HE [1]. 3D configurations have been calculated here to test the new beta version of the electromagnetism package coupled with the dynamics in Ls-Dyna and compared with the ICE Z shot 1967.
NASA Technical Reports Server (NTRS)
Pline, Alexander D.; Werner, Mark P.; Hsieh, Kwang-Chung
1991-01-01
The Surface Tension Driven Convection Experiment (STDCE) is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the United States Microgravity Laboratory-1 (USML-1) Spacelab mission planned for June, 1992. One of the components of data collected during the experiment is a video record of the flow field. This qualitative data is then quantified using an all electric, two dimensional Particle Image Velocimetry (PIV) technique called Particle Displacement Tracking (PDT), which uses a simple space domain particle tracking algorithm. Results using the ground based STDCE hardware, with a radiant flux heating mode, and the PDT system are compared to numerical solutions obtained by solving the axisymmetric Navier Stokes equations with a deformable free surface. The PDT technique is successful in producing a velocity vector field and corresponding stream function from the raw video data which satisfactorily represents the physical flow. A numerical program is used to compute the velocity field and corresponding stream function under identical conditions. Both the PDT system and numerical results were compared to a streak photograph, used as a benchmark, with good correlation.
NASA Technical Reports Server (NTRS)
Pline, Alexander D.; Wernet, Mark P.; Hsieh, Kwang-Chung
1991-01-01
The Surface Tension Driven Convection Experiment (STDCE) is a Space Transportation System flight experiment to study both transient and steady thermocapillary fluid flows aboard the United States Microgravity Laboratory-1 (USML-1) Spacelab mission planned for June, 1992. One of the components of data collected during the experiment is a video record of the flow field. This qualitative data is then quantified using an all electric, two dimensional Particle Image Velocimetry (PIV) technique called Particle Displacement Tracking (PDT), which uses a simple space domain particle tracking algorithm. Results using the ground based STDCE hardware, with a radiant flux heating mode, and the PDT system are compared to numerical solutions obtained by solving the axisymmetric Navier Stokes equations with a deformable free surface. The PDT technique is successful in producing a velocity vector field and corresponding stream function from the raw video data which satisfactorily represents the physical flow. A numerical program is used to compute the velocity field and corresponding stream function under identical conditions. Both the PDT system and numerical results were compared to a streak photograph, used as a benchmark, with good correlation.
Haque, Rubaiyet Iftekharul; Loussert, Christophe; Sergent, Michelle; Benaben, Patrick; Boddaert, Xavier
2015-01-01
Optimization of the acoustic resonant sensor requires a clear understanding of how the output responses of the sensor are affected by the variation of different factors. During this work, output responses of a capacitive acoustic transducer, such as membrane displacement, quality factor, and capacitance variation, are considered to evaluate the sensor design. The six device parameters taken into consideration are membrane radius, backplate radius, cavity height, air gap, membrane tension, and membrane thickness. The effects of factors on the output responses of the transducer are investigated using an integrated methodology that combines numerical simulation and design of experiments (DOE). A series of numerical experiments are conducted to obtain output responses for different combinations of device parameters using finite element methods (FEM). Response surface method is used to identify the significant factors and to develop the empirical models for the output responses. Finally, these results are utilized to calculate the optimum device parameters using multi-criteria optimization with desirability function. Thereafter, the validating experiments are designed and deployed using the numerical simulation to crosscheck the responses. PMID:25894937
Numerical experiments on the modulation theory for the nonlinear atomic chain
NASA Astrophysics Data System (ADS)
Dreyer, W.; Herrmann, M.
2008-02-01
Modulation theory with periodic travelling waves is a powerful, but not rigorous tool to derive a thermodynamic description for atomic chains with nearest neighbour interactions (FPU chains). This theory is sufficiently complex to deal with strong oscillations on the microscopic scale, and therefore it is capable to describe the creation of temperature and the transport of heat on a macroscopic scale. In this paper we investigate the validity of modulation theory by means of several numerical experiments. We start with a survey on the foundations of modulation theory. In particular, we discuss the hyperbolic scaling, the notion of cold data, microscopic oscillations and Young measures, periodic and modulated travelling waves, and, finally, the resulting macroscopic conservation laws. Afterwards we discuss how the validity of a macroscopic theory may be tested within numerical simulations of the microscopic dynamics. To this end we describe an approach to thermodynamic data exploration which is motivated by the theory of Young measures, and relies on mesoscopic windows in space and time. The last part is devoted to several numerical experiments including examples with periodic boundary conditions and smooth initial data, and macroscopic Riemann problems. We interpret the outcome of these experiments in the framework of thermodynamics, and end up with two conclusions. (1) There are many examples for which modulation theory provides in fact the right thermodynamic description because it can predict both the structure of the microscopic oscillations and their macroscopic evolution correctly. (2) Modulation theory will fail if the oscillations exhibit a more complicate structure.
Integrated Numerical Experiments (INEX) and the Free-Electron Laser Physical Process Code (FELPPC)
Thode, L.E.; Chan, K.C.D.; Schmitt, M.J.; McKee, J.; Ostic, J.; Elliott, C.J.; McVey, B.D.
1990-01-01
The strong coupling of subsystem elements, such as the accelerator, wiggler, and optics, greatly complicates the understanding and design of a free electron laser (FEL), even at the conceptual level. Given the requirements for high-performance FELs, the strong coupling between the laser subsystems must be included to obtain a realistic picture of the potential operational capability. To address the strong coupling character of the FEL the concept of an Integrated Numerical Experiment (INEX) was proposed. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostics. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. The INEX approach has been applied to a large number of accelerator and FEL experiments. Overall, the agreement between INEX and the experiments is very good. Despite the success of INEX, the approach is difficult to apply to trade-off and initial design studies because of the significant manpower and computational requirements. On the other hand, INEX provides a base from which realistic accelerator, wiggler, and optics models can be developed. The Free Electron Laser Physical Process Code (FELPPC) includes models developed from INEX, provides coupling between the subsystems models and incorporates application models relevant to a specific trade-off or design study.
Haque, Rubaiyet Iftekharul; Loussert, Christophe; Sergent, Michelle; Benaben, Patrick; Boddaert, Xavier
2015-01-01
Optimization of the acoustic resonant sensor requires a clear understanding of how the output responses of the sensor are affected by the variation of different factors. During this work, output responses of a capacitive acoustic transducer, such as membrane displacement, quality factor, and capacitance variation, are considered to evaluate the sensor design. The six device parameters taken into consideration are membrane radius, backplate radius, cavity height, air gap, membrane tension, and membrane thickness. The effects of factors on the output responses of the transducer are investigated using an integrated methodology that combines numerical simulation and design of experiments (DOE). A series of numerical experiments are conducted to obtain output responses for different combinations of device parameters using finite element methods (FEM). Response surface method is used to identify the significant factors and to develop the empirical models for the output responses. Finally, these results are utilized to calculate the optimum device parameters using multi-criteria optimization with desirability function. Thereafter, the validating experiments are designed and deployed using the numerical simulation to crosscheck the responses. PMID:25894937
Wu, R.S.; Xie, X.B.
1994-12-31
The theory of spatial distribution of seismic energy density in one dimensional (1D) random media derived in part 1 (Wu, 1993) is tested by numerical experiments using a full wave propagation matrix method. The geometry of numerical experiment mimics the configuration of zero-offset VSP (Vertical Seismic Profiling) along a borehole. A procedure of octave-band frequency averaging is applied to the measured data to reduce fluctuation of spatial energy distribution, so that stable estimations of medium parameters can be achieved without resorting to ensemble averaging. Results from Monte-Carlo numerical experiments for both infinite random media and finite random slabs with or without bottom reflections show good agreement for dark-to-gray (weak to intermediate scattering compared with absorption) media. When scattering is very strong (when backscattering-absorption ratio S{sub b} > 3), results from single realization fluctuate substantially. However, most the practical situations of sedimentary rocks in the crust fall into the validity region of the energy transfer theory.
Integrated Numerical Experiments (INEX) and the Free-Electron Laser Physical Process Code (FELPPC)
NASA Astrophysics Data System (ADS)
Thode, L. E.; Chan, K. C. D.; Schmitt, M. J.; McKee, J.; Ostic, J.; Elliott, C. J.; McVey, B. D.
The strong coupling of subsystem elements, such as the accelerator, wiggler, and optics, greatly complicates the understanding and design of a free electron laser (FEL), even at the conceptual level. Given the requirements for high-performance FELs, the strong coupling between the laser subsystems must be included to obtain a realistic picture of the potential operational capability. To address the strong coupling character of the FEL the concept of an Integrated Numerical Experiment (INEX) was proposed. Unique features of the INEX approach are consistency and numerical equivalence of experimental diagnostics. The equivalent numerical diagnostics mitigates the major problem of misinterpretation that often occurs when theoretical and experimental data are compared. The INEX approach has been applied to a large number of accelerator and FEL experiments. Overall, the agreement between INEX and the experiments is very good. Despite the success of INEX, the approach is difficult to apply to trade-off and initial design studies because of the significant manpower and computational requirements. On the other hand, INEX provides a base from which realistic accelerator, wiggler, and optics models can be developed. The Free Electron Laser Physical Process Code (FELPPC) includes models developed from INEX, provides coupling between the subsystems models and incorporates application models relevant to a specific trade-off or design study.
An additional condition for Bell experiments for accepting local realistic theories
NASA Astrophysics Data System (ADS)
Nagata, Koji; Nakamura, Tadao
2013-12-01
We assume that one source of two uncorrelated spin-carrying particles emits them in a state, which can be described as a spin-1/2 bipartite pure uncorrelated state. We consider a Bell-Clauser-Horne-Shimony-Holt (Bell-CHSH) experiment with two-orthogonal-settings. We propose an additional condition for the state to be reproducible by the property of local realistic theories. We use the proposed measurement theory in order to construct the additional condition (Nagata and Nakamura in Int J Theor Phys 49:162, 2010). The condition is that local measurement outcome is . Otherwise, such an experiment does not allow for the existence of local realistic theories even in the situation that all Bell-CHSH inequalities hold. Also we derive new set of Bell inequalities when local measurement outcome is.
Hartigan, P.; Carver, R.; Foster, J. M.; Rosen, P. A.; Williams, R. J. R.; Wilde, B. H.; Coker, R. F.; Hansen, J. F.; Blue, B. E.; Frank, A.
2009-11-01
Collimated supersonic flows in laboratory experiments behave in a similar manner to astrophysical jets provided that radiation, viscosity, and thermal conductivity are unimportant in the laboratory jets and that the experimental and astrophysical jets share similar dimensionless parameters such as the Mach number and the ratio of the density between the jet and the ambient medium. When these conditions apply, laboratory jets provide a means to study their astrophysical counterparts for a variety of initial conditions, arbitrary viewing angles, and different times, attributes especially helpful for interpreting astronomical images where the viewing angle and initial conditions are fixed and the time domain is limited. Experiments are also a powerful way to test numerical fluid codes in a parameter range in which the codes must perform well. In this paper, we combine images from a series of laboratory experiments of deflected supersonic jets with numerical simulations and new spectral observations of an astrophysical example, the young stellar jet HH 110. The experiments provide key insights into how deflected jets evolve in three dimensions, particularly within working surfaces where multiple subsonic shells and filaments form, and along the interface where shocked jet material penetrates into and destroys the obstacle along its path. The experiments also underscore the importance of the viewing angle in determining what an observer will see. The simulations match the experiments so well that we can use the simulated velocity maps to compare the dynamics in the experiment with those implied by the astronomical spectra. The experiments support a model where the observed shock structures in HH 110 form as a result of a pulsed driving source rather than from weak shocks that may arise in the supersonic shear layer between the Mach disk and bow shock of the jet's working surface.
NASA Technical Reports Server (NTRS)
Over, Ann P.
2001-01-01
The Combustion Module-1 (CM-1) was a large, state-of-the-art space shuttle Spacelab facility that was designed, built, and operated on STS-83 and STS-94 by a team from the NASA Glenn Research Center composed of civil servants and local support contractors (Analex and Zin Technologies). CM-1 accomplished the incredible task of providing a safe environment to support flammable and toxic gases while providing a suite of diagnostics for science measurements more extensive than any prior shuttle experiment (or anything since). Finally, CM-1 proved that multiple science investigations can be accommodated in one facility, a crucial step for Glenn's Fluids and Combustion Facility developed for the International Space Station. However, the story does not end with CM-1. In 1998, CM-2 was authorized to take the CM-1 accomplishments a big step further by completing three major steps: Converting the entire experiment to operate in a SPACEHAB module. Conducting an extensive hardware refurbishment and upgrading diagnostics (e.g., cameras, gas chromatograph, and numerous sensors). Adding a new, completely different combustion experiment.
Charrois, Theresa; Rosenthal, Meagen; Tsuyuki, Ross T.
2012-01-01
Background: Pharmacists in Alberta can apply to the Alberta College of Pharmacists in order to obtain the designation of additional prescriber. This designation uniquely allows them to initiate therapy, in addition to other medication-related activities. Our objective was to examine specific experiences of pharmacists regarding the decision to apply and the application itself, and use this information to help inform other pharmacists who are considering additional prescribing. Methods: All pharmacists involved in a randomized, controlled trial being conducted in rural Alberta who had received their additional prescribing authorization (APA) were invited to participate. Pharmacists were contacted via e-mail and asked to respond to questions regarding their experiences in applying for APA. Responses were analyzed using content analysis and the identites of all respondents were kept anonymous. Results: Fourteen pharmacists were invited to participate. Review and examination of the responses revealed 3 main themes: motivation, hurdles and outcomes. Motivation can be understood as the reasons why they applied for their APA. Hurdles include any problems encountered of a personal, environmental or professional nature. Outcomes refer to how this designation has changed their practice. Discussion: Pharmacists had to address many factors that were unexpected during the application process; however, the eventual outcome of obtaining APA was deemed beneficial, both professionally and with regard to patient care. Conclusion: The information shared from these pharmacists will help other pharmacists, regardless of jurisdiction, overcome some of the challenges associated with obtaining advanced prescribing privileges. PMID:23509485
Numerical experiments with an implicit particle filter for the shallow water equations
NASA Astrophysics Data System (ADS)
Souopgui, I.; Chorin, A. J.; Hussaini, M.
2012-12-01
The estimation of initial conditions for the shallow water equations for a given set of later data is a well known test problem for data assimilation codes. A popular approach to this problem is the variational method (4D-Var), i.e. the computation of the mode of the posterior probability density function (pdf) via the adjoint technique. Here, we improve on 4D-Var by computing the conditional mean (the minimum least square error estimator) rather than the mode (a biased estimator) and we do so with implicit sampling, a Monte Carlo (MC) importance sampling method. The idea in implicit sampling is to first search for the high-probability region of the posterior pdf and then to find samples in this region. Because the samples are concentrated in the high-probability region, fewer samples are required than with competing MC schemes. The search for the high-probability region can be implemented by a minimization that is very similar to the minimization in 4D-Var, and we make use of a 4D-Var code in our implementation. The samples are obtained by solving algebraic equations with a random right-hand-side. These equations can be solved efficiently, so that the additional cost of our approach, compared to traditional 4D-Var, is small. The long-term goal is to assimilate experimental data, obtained with the CORIOLIS turntable in Grenoble (France), to study the drift of a vortex. We present results from numerical twin experiments as a first step towards our long-term goal. We discretize the shallow water equations on a square domain (2.5m× 2.5m) using finite differences on a staggered grid of size 28× 28 and a fourth order Runge-Kutta. We assume open boundary conditions and estimate the initial state (velocities and surface height) given noisy observations of the state. We solve the optimization problem using a 4D-Var code that relies on a L-BFGS method; the random algebraic equations are solved with random maps, i.e. we look for solutions in given, but random, directions
NASA Technical Reports Server (NTRS)
Ramachandran, Narayanan
2000-01-01
Normal vibrational modes on large spacecraft are excited by crew activity, operating machinery, and other mechanical disturbances. Periodic engine burns for maintaining vehicle attitude and random impulse type disturbances also contribute to the acceleration environment of a Spacecraft. Accelerations from these vibrations (often referred to as g-jitter) are several orders of magnitude larger than the residual accelerations from atmospheric drag and gravity gradient effects. Naturally, the effects of such accelerations have been a concern to prospective experimenters wishing to take advantage of the microgravity environment offered by spacecraft operating in low Earth orbit and the topic has been studied extensively, both numerically and analytically. However, these studies have not produced a general theory that predicts the effects of multi-spectral periodic accelerations on a general class of experiments nor have they produced scaling laws that a prospective experimenter could use to assess how his/her experiment might be affected by this acceleration environment. Furthermore, there are no actual flight experimental data that correlates heat or mass transport with measurements of the periodic acceleration environment. The present investigation approaches this problem with carefully conducted terrestrial experiments and rigorous numerical modeling thereby providing comparative theoretical and experimental data. The modeling, it is hoped will provide a predictive tool that can be used for assessing experiment response to Spacecraft vibrations.
NASA Technical Reports Server (NTRS)
Marchese, Anthony J.; Dryer, Frederick L.
1997-01-01
This program supports the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies. Experimental emphasis is on the study of simple alcohols (methanol, ethanol) and alkanes (n-heptane, n-decane) as fuels with time dependent measurements of drop size, flame-stand-off, liquid-phase composition, and finally, extinction. Experiments have included bench-scale studies at Princeton, studies in the 2.2 and 5.18 drop towers at NASA-LeRC, and both the Fiber Supported Droplet Combustion (FSDC-1, FSDC-2) and the free Droplet Combustion Experiment (DCE) studies aboard the shuttle. Test matrix and data interpretation are performed through spherically-symmetric, time-dependent numerical computations which embody detailed sub-models for physical and chemical processes. The computed burning rate, flame stand-off, and extinction diameter are compared with the respective measurements for each individual experiment. In particular, the data from FSDC-1 and subsequent space-based experiments provide the opportunity to compare all three types of data simultaneously with the computed parameters. Recent numerical efforts are extending the computational tools to consider time dependent, axisymmetric 2-dimensional reactive flow situations.
NASA Astrophysics Data System (ADS)
Saxena, A. K.; Kaushik, T. C.; Gupta, Satish C.
2010-03-01
Two low energy (1.6 and 8 kJ) portable electrically exploding foil accelerators are developed for moderately high pressure shock studies at small laboratory scale. Projectile velocities up to 4.0 km/s have been measured on Kapton flyers of thickness 125 μm and diameter 8 mm, using an in-house developed Fabry-Pérot velocimeter. An asymmetric tilt of typically few milliradians has been measured in flyers using fiber optic technique. High pressure impact experiments have been carried out on tantalum, and aluminum targets up to pressures of 27 and 18 GPa, respectively. Peak particle velocities at the target-glass interface as measured by Fabry-Pérot velocimeter have been found in good agreement with the reported equation of state data. A one-dimensional hydrodynamic code based on realistic models of equation of state and electrical resistivity has been developed to numerically simulate the flyer velocity profiles. The developed numerical scheme is validated against experimental and simulation data reported in literature on such systems. Numerically computed flyer velocity profiles and final flyer velocities have been found in close agreement with the previously reported experimental results with a significant improvement over reported magnetohydrodynamic simulations. Numerical modeling of low energy systems reported here predicts flyer velocity profiles higher than experimental values, indicating possibility of further improvement to achieve higher shock pressures.
Numerical Experiments and Flow Visualization of Drag Reduction using EMHD Control
NASA Astrophysics Data System (ADS)
O'Sullivan, Peter; Biringen, Sedat
1997-11-01
Turbulent channel flow of saltwater is studied numerically with the aim of achieving drag reduction via EMHD control using the Lorentz force of flush-mounted microtiles developed by Bandyopadhyay at NUWC. Previous numerical simulations have indicated significant local spatial deviations on the time-average skin friction (± 10%) in the vicinity of the control actuators. However, there was only negligible net viscous drag reduction (<2%). In an attempt to better understand the physics of the interaction of the controlling Lorentz force with the passage of an advecting incipient burst we have performed numerous short simulations. We scanned a wide variety of archival data from our previous simulations and extracted a good sample of burst ``candidate'' initial conditions. We then advanced these flows forward in time both with and without control. We studied the resulting burst/control dynamics with respect to control strength, duration, temporal frequency, spanwise offset and streamwise spacing of the microtiles. The spatio-temporal response of the burst and sweep structures is crucial to a successful drag reduction strategy. By identifying the key parameters in the control space we aim to narrow the focus of the design problem. We will present our findings of these numerical experiments which are based on data analysis together with flow visualization.
3D flow past transonic turbine cascade SE 1050 — Experiment and numerical simulations
NASA Astrophysics Data System (ADS)
Šimurda, D.; Fürst, J.; Luxa, M.
2013-08-01
This paper is concerned with experimental and numerical research on 3D flow past prismatic turbine cascade SE1050 (known in QNET network as open test case SE1050). The primary goal was to assess the influence of the inlet velocity profile on the flow structures in the interblade channel and on the flow field parameters at the cascade exit and to compare these findings to results of numerical simulations. Investigations of 3D flow past the cascade with non-uniform inlet velocity profile were carried out both experimentally and numerically at subsonic ( M 2is = 0.8) and at transonic ( M 2is = 1.2) regime at design angle of incidence. Experimental data was obtained using a traversing device with a five-hole conical probe. Numerically, the 3D flow was simulated by open source code OpenFOAM and in-house code. Analyses of experimental data and CFD simulations have revealed the development of distinctive vortex structures resulting from non-uniform inlet velocity profile. Origin of these structures results in increased loss of kinetic energy and spanwise shift of kinetic energy loss coefficient distribution. Differences found between the subsonic and the transonic case confirm earlier findings available in the literature. Results of CFD and experiments agree reasonably well.
NASA Astrophysics Data System (ADS)
Figueroa, Aldo; Meunier, Patrice; Cuevas, Sergio; Villermaux, Emmanuel; Ramos, Eduardo
2014-01-01
We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, "The diffusive strip method for scalar mixing in two-dimensions," J. Fluid Mech. 662, 134-172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.
Figueroa, Aldo; Meunier, Patrice; Villermaux, Emmanuel; Cuevas, Sergio; Ramos, Eduardo
2014-01-15
We present a combination of experiment, theory, and modelling on laminar mixing at large Péclet number. The flow is produced by oscillating electromagnetic forces in a thin electrolytic fluid layer, leading to oscillating dipoles, quadrupoles, octopoles, and disordered flows. The numerical simulations are based on the Diffusive Strip Method (DSM) which was recently introduced (P. Meunier and E. Villermaux, “The diffusive strip method for scalar mixing in two-dimensions,” J. Fluid Mech. 662, 134–172 (2010)) to solve the advection-diffusion problem by combining Lagrangian techniques and theoretical modelling of the diffusion. Numerical simulations obtained with the DSM are in reasonable agreement with quantitative dye visualization experiments of the scalar fields. A theoretical model based on log-normal Probability Density Functions (PDFs) of stretching factors, characteristic of homogeneous turbulence in the Batchelor regime, allows to predict the PDFs of scalar in agreement with numerical and experimental results. This model also indicates that the PDFs of scalar are asymptotically close to log-normal at late stages, except for the large concentration levels which correspond to low stretching factors.
NASA Astrophysics Data System (ADS)
Champlain, A.; Matéo-Vélez, J.-C.; Roussel, J.-F.; Hess, S.; Sarrailh, P.; Murat, G.; Chardon, J.-P.; Gajan, A.
2016-01-01
Recent high-altitude observations, made by the Lunar Dust Experiment (LDEX) experiment on board LADEE orbiting the Moon, indicate that high-altitude (>10 km) dust particle densities are well correlated with interplanetary dust impacts. They show no evidence of high dust density suggested by Apollo 15 and 17 observations and possibly explained by electrostatic forces imposed by the plasma environment and photon irradiation. This paper deals with near-surface conditions below the domain of observation of LDEX where electrostatic forces could clearly be at play. The upper and lower limits of the cohesive force between dusts are obtained by comparing experiments and numerical simulations of dust charging under ultraviolet irradiation in the presence of an electric field and mechanical vibrations. It is suggested that dust ejection by electrostatic forces is made possible by microscopic-scale amplifications due to soil irregularities. At low altitude, this process may be complementary to interplanetary dust impacts.
Iceberg capsize hydrodynamics: a comparison of laboratory experiments and numerical modeling
NASA Astrophysics Data System (ADS)
Burton, J. C.; Cathles, L. M.; Correa-Legisos, S.; Ellowitz, J.; Darnell, K.; Zhang, W. W.; MacAyeal, D. R.
2013-12-01
Large icebergs are often observed to capsize in open water near fjords. During capsize, large amounts of gravitational potential energy are released which can lead to coastal tsunamis, mixing of the water column, and possibly lead to further calving at the glacier terminus. This process is rarely studied; in nature the scale and irregular timing of the events makes observations exceedingly difficult. Here we compare laboratory experiments and numerical modeling of the capsize process to better understand the coupling of the hydrodynamic forces to the solid iceberg. Although the characteristic Reynolds number is much lower for both the laboratory model and the numerical simulations, the comparison provides a starting point to quantify and identify generic features that can be estimated in the field, such as hydrodynamic pressure, water flow velocities, vertical mixing, and elastic stresses on the iceberg itself, which could lead to fracture.
Numerical prediction of the monsoon depression of 5-7 July 1979. [Monsoon Experiment (MONEX)
NASA Technical Reports Server (NTRS)
Shukla, J.; Atlas, R.; Baker, W. E.
1981-01-01
A well defined monsoon depression was used for two assimilation and forecast experiments: (1) using conventional surface and upper air data, (2) using these data plus Monex data. The data sets were assimilated and used with a general circulation model to make numerical predictions. The model, the analysis and assimilation procedure, the differences in the analyses due to different data inputs, and the differences in the numerical predictions are described. The MONEX data have a positive impact, although the differences after 24 hr are not significant. The MONEX assimilation does not agree with manual analysis location of depression center. The 2.5 x 3 deg horizontal resolution of the prediction model is too coarse. The assimilation of geopotential height data derived from satellite soundings generated gravity waves with amplitudes similar to the meteorologically significant features investigated.
Numerical experiments on the climatic sensitivity of an atmospheric hydrologic cycle
NASA Technical Reports Server (NTRS)
Roads, J. O.
1978-01-01
It is shown for an intermediate numerical model that fractional cloudiness and relative humidity decrease with increasing temperature. The fractional cloudiness decreases at a rate about 1 per deg K. This occurs in spite of an increase in the evaporation, water transport, condensation, precipitation and cloud water content with increasing temperature. These results are quite similar to those found from models with more highly parameterized clouds, notably the NCAR model. The fractional cloudiness in this model is measured by the fractional coverage of total cloud water and the fractional coverage of positive condensation, in addition to the relative humidity. It is also shown that some of the characteristics of a temperate climate can be simulated in an intermediate numerical model with periodic, antisymmetric and symmetric boundary conditions on an f plane. Intermediate models of this sort may therefore be useful to investigate general questions about the earth's hydrologic cycle on climatic space and time scales
Comparison of numerical simulations to experiments for atomization in a jet nebulizer.
Lelong, Nicolas; Vecellio, Laurent; Sommer de Gélicourt, Yann; Tanguy, Christian; Diot, Patrice; Junqua-Moullet, Alexandra
2013-01-01
The development of jet nebulizers for medical purposes is an important challenge of aerosol therapy. The performance of a nebulizer is characterized by its output rate of droplets with a diameter under 5 µm. However the optimization of this parameter through experiments has reached a plateau. The purpose of this study is to design a numerical model simulating the nebulization process and to compare it with experimental data. Such a model could provide a better understanding of the atomization process and the parameters influencing the nebulizer output. A model based on the Updraft nebulizer (Hudson) was designed with ANSYS Workbench. Boundary conditions were set with experimental data then transient 3D calculations were run on a 4 µm mesh with ANSYS Fluent. Two air flow rate (2 L/min and 8 L/min, limits of the operating range) were considered to account for different turbulence regimes. Numerical and experimental results were compared according to phenomenology and droplet size. The behavior of the liquid was compared to images acquired through shadowgraphy with a CCD Camera. Three experimental methods, laser diffractometry, phase Doppler anemometry (PDA) and shadowgraphy were used to characterize the droplet size distributions. Camera images showed similar patterns as numerical results. Droplet sizes obtained numerically are overestimated in relation to PDA and diffractometry, which only consider spherical droplets. However, at both flow rates, size distributions extracted from numerical image processing were similar to distributions obtained from shadowgraphy image processing. The simulation then provides a good understanding and prediction of the phenomena involved in the fragmentation of droplets over 10 µm. The laws of dynamics apply to droplets down to 1 µm, so we can assume the continuity of the distribution and extrapolate the results for droplets between 1 and 10 µm. So, this model could help predicting nebulizer output with defined geometrical and
Comparison of Numerical Simulations to Experiments for Atomization in a Jet Nebulizer
Lelong, Nicolas; Vecellio, Laurent; Sommer de Gélicourt, Yann; Tanguy, Christian; Diot, Patrice; Junqua-Moullet, Alexandra
2013-01-01
The development of jet nebulizers for medical purposes is an important challenge of aerosol therapy. The performance of a nebulizer is characterized by its output rate of droplets with a diameter under 5 µm. However the optimization of this parameter through experiments has reached a plateau. The purpose of this study is to design a numerical model simulating the nebulization process and to compare it with experimental data. Such a model could provide a better understanding of the atomization process and the parameters influencing the nebulizer output. A model based on the Updraft nebulizer (Hudson) was designed with ANSYS Workbench. Boundary conditions were set with experimental data then transient 3D calculations were run on a 4 µm mesh with ANSYS Fluent. Two air flow rate (2 L/min and 8 L/min, limits of the operating range) were considered to account for different turbulence regimes. Numerical and experimental results were compared according to phenomenology and droplet size. The behavior of the liquid was compared to images acquired through shadowgraphy with a CCD Camera. Three experimental methods, laser diffractometry, phase Doppler anemometry (PDA) and shadowgraphy were used to characterize the droplet size distributions. Camera images showed similar patterns as numerical results. Droplet sizes obtained numerically are overestimated in relation to PDA and diffractometry, which only consider spherical droplets. However, at both flow rates, size distributions extracted from numerical image processing were similar to distributions obtained from shadowgraphy image processing. The simulation then provides a good understanding and prediction of the phenomena involved in the fragmentation of droplets over 10 µm. The laws of dynamics apply to droplets down to 1 µm, so we can assume the continuity of the distribution and extrapolate the results for droplets between 1 and 10 µm. So, this model could help predicting nebulizer output with defined geometrical and
"Physically-based" numerical experiment to determine the dominant hillslope processes during floods?
NASA Astrophysics Data System (ADS)
Gaume, Eric; Esclaffer, Thomas; Dangla, Patrick; Payrastre, Olivier
2016-04-01
To study the dynamics of hillslope responses during flood event, a fully coupled "physically-based" model for the combined numerical simulation of surface runoff and underground flows has been developed. A particular attention has been given to the selection of appropriate numerical schemes for the modelling of both processes and of their coupling. Surprisingly, the most difficult question to solve, from a numerical point of view, was not related to the coupling of two processes with contrasted kinetics such as surface and underground flows, but to the high gradient infiltration fronts appearing in soils, source of numerical diffusion, instabilities and sometimes divergence. The model being elaborated, it has been successfully tested against results of high quality experiments conducted on a laboratory sandy slope in the early eighties, which is still considered as a reference hillslope experimental setting (Abdul & Guilham). The model appeared able to accurately simulate the pore pressure distributions observed in this 1.5 meter deep and wide laboratory hillslope, as well as its outflow hydrograph shapes and the measured respective contributions of direct runoff and groundwater to these outflow hydrographs. Based on this great success, the same model has been used to simulate the response of a theoretical 100-meter wide and 10% sloped hillslope, with a 2 meter deep pervious soil and impervious bedrock. Three rain events have been tested: a 100 millimeter rainfall event over 10 days, over 1 day or over one hour. The simulated responses are hydrologically not realistic and especially the fast component of the response, that is generally observed in the real-world and explains flood events, is almost absent of the simulated response. Thinking a little about the whole problem, the simulation results appears totally logical according to the proposed model. The simulated response, in fact a recession hydrograph, corresponds to a piston flow of a relatively uniformly
A numerical experiment on the equilibrium and stability of a rotating galactic bar
NASA Technical Reports Server (NTRS)
Miller, R. H.; Vandervoort, P. O.; Welty, D. E.; Smith, B. F.
1982-01-01
A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.
Spinel dissolution via addition of glass forming chemicals. Results of preliminary experiments
Fox, K. M.; Johnson, F. C.
2015-11-01
Increased loading of high level waste in glass can lead to crystallization within the glass. Some crystalline species, such as spinel, have no practical impact on the chemical durability of the glass, and therefore may be acceptable from both a processing and a product performance standpoint. In order to operate a melter with a controlled amount of crystallization, options must be developed for remediating an unacceptable accumulation of crystals. This report describes preliminary experiments designed to evaluate the ability to dissolve spinel crystals in simulated waste glass melts via the addition of glass forming chemicals (GFCs).
NASA Astrophysics Data System (ADS)
Dossmann, Yvan; Paci, Alexandre; Auclair, Francis; Floor, Jochem
2010-05-01
Internal tides are suggested to play a major role in the sustaining of the global oceanic circulation [1][5]. Although the exact origin of the energy conversions occurring in stratified fluids is questioned [2], it is clear that the diapycnal energy transfers provided by the energy cascade of internal gravity waves generated at tidal frequencies in regions of steep bathymetry is strongly linked to the general circulation energy balance. Therefore a precise quantification of the energy supply by internal waves is a crucial step in forecasting climate, since it improves our understanding of the underlying physical processes. We focus on an academic case of internal waves generated over an oceanic ridge in a linearly stratified fluid. In order to accurately quantify the diapycnal energy transfers caused by internal waves dynamics, we adopt a complementary approach involving both laboratory and numerical experiments. The laboratory experiments are conducted in a 4m long tank of the CNRM-GAME fluid mechanics laboratory, well known for its large stratified water flume (e.g. Knigge et al [3]). The horizontal oscillation at precisely controlled frequency of a Gaussian ridge immersed in a linearly stratified fluid generates internal gravity waves. The ridge of e-folding width 3.6 cm is 10 cm high and spans 50 cm. We use PIV and Synthetic Schlieren measurement techniques, to retrieve the high resolution velocity and stratification anomaly fields in the 2D vertical plane across the ridge. These experiments allow us to get access to real and exhaustive measurements of a wide range of internal waves regimes by varying the precisely controlled experimental parameters. To complete this work, we carry out some direct numerical simulations with the same parameters (forcing amplitude and frequency, initial stratification, boundary conditions) as the laboratory experiments. The model used is a non-hydrostatic version of the numerical model Symphonie [4]. Our purpose is not only to
Application of numerical modelling in the design of a full-scale heated Tunnel Sealing Experiment
NASA Astrophysics Data System (ADS)
Guo, R.; Chandler, N.; Martino, J.; Dixon, D.
2005-10-01
The Tunnel Sealing Experiment (TSX) was a full-scale in situ demonstration of technology for constructing nearly water tight-seals in excavations through crystalline rock deep below the surface of the earth. The experiment has been carried out at Atomic Energy of Canada Limited's (AECL's) Underground Research Laboratory near Lac du Bonnet, Canada, in support of international programs for geologic disposal of radioactive waste. The TSX, with partners from Canada, Japan, France and the United States, was carried under conditions of high pressure (up to 4 MPa) and elevated temperature (up to 85°C). Comparing numerical model predictions with eight years of data collected from approximately 900 sensors was an important component of this experiment. Model of Transport In Fractured/porous Media (MOTIF), a finite element computer program developed by AECL for simulating fully coupled or uncoupled fluid flow, solute transport and heat transport, was used to model both the ambient temperature and heated phases of the TSX. The plan to heat the water in the TSX to 85°C was developed using model predictions and a comparison of simulated results with measurements during heating of the water in the TSX to about 50°C. The three-dimensional MOTIF simulations were conducted in parallel with axisymmetric modelling using Fast Lagrangian Analysis of Continua (FLAC), which computed the heat loss from pipes that carried the heated water through the rock to and from the experiment. The numerical model was initially used to develop a plan for operation of the experiment heaters, and then subsequently used to predict temperatures and hydraulic heads in the TSX bulkhead seals and surrounding rock. Copyright
Patel, Bhavesh G.; Das, Amita; Kaw, Predhiman; Singh, Rajesh; Kumar, Ajai
2012-07-15
Two dimensional numerical modelling based on simplified hydrodynamic evolution for an expanding plasma plume (created by laser blow off) against an ambient background gas has been carried out. A comparison with experimental observations shows that these simulations capture most features of the plasma plume expansion. The plume location and other gross features are reproduced as per the experimental observation in quantitative detail. The plume shape evolution and its dependence on the ambient background gas are in good qualitative agreement with the experiment. This suggests that a simplified hydrodynamic expansion model is adequate for the description of plasma plume expansion.
Numerical experiments on short-term meteorological effects on solar variability
NASA Technical Reports Server (NTRS)
Somerville, R. C. J.; Hansen, J. E.; Stone, P. H.; Quirk, W. J.; Lacis, A. A.
1975-01-01
A set of numerical experiments was conducted to test the short-range sensitivity of a large atmospheric general circulation model to changes in solar constant and ozone amount. On the basis of the results of 12-day sets of integrations with very large variations in these parameters, it is concluded that realistic variations would produce insignificant meteorological effects. Any causal relationships between solar variability and weather, for time scales of two weeks or less, rely upon changes in parameters other than solar constant or ozone amounts, or upon mechanisms not yet incorporated in the model.
Hidalgo, M E; Ayesa, E
2001-09-01
This paper describes a mathematical tool for local identifiability analysis that can easily be applied to high-order state-space nonlinear systems and implemented in simulators with a discrete-time approach. The methodology is based on the recursive numerical evaluation of a reduced information matrix during the simulation of a calibration experiment and in the setting-up of a group of information parameters based on geometric interpretations of this matrix. As an example of application, the proposed methodology has been used in the study of an OUR batch test from the point of view of ASM No. 1 calibration. PMID:11487118
Numerical modeling and experiments of creep crack growth under cyclic loading
Brust, F.W.
1995-12-31
This paper presents a summary of some recent studies of creep crack growth under history dependent load conditions. The effect of a proper constitutive law is illustrated. Moreover, the asymptotic fields are reconsidered under cyclic creep conditions. In addition, several experiments are modeled and the behavior of integral parameters is discussed.
Numerical Experiments and Flow Visualization of Drag Reduction using EMHD Control
NASA Astrophysics Data System (ADS)
O'Sullivan, Peter; Biringen, Sedat
1998-11-01
Turbulent channel flow of saltwater is studied numerically with the aim of achieving drag reduction via EMHD control using the Lorentz force of flush-mounted microtiles similar to those developed by Bandyopadhyay at NUWC. Our prior numerical simulations indicated significant local spatial deviations on the time-average skin friction (± 10%) in the vicinity of the control actuators. However, there was negligible net viscous drag reduction (<1%). To better understand the physics of the interaction of the controlling Lorentz force with the passage of advecting burst events we have performed short-duration simulations. We implement a burst detection scheme based on that of Alfredsson and Johansson (1984) in order to locate and track the strongest Q2 events in the flow. Subsequently, we compare the evolution of these structures, both with and without EMHD control. We find that the specific designs we have studied do not succeed in reducing the primary Reynolds stress as a burst advects above an actuator. Via flow visualization of large-scale coherent structures we find that the classic hairpin vortices which populate the boundary layer are not significantly affected by the applied control - apart from a temporary spatial phase shift. We will present our findings of these numerical experiments and discuss the prospects for a more successful design and control strategy.
NASA Astrophysics Data System (ADS)
Ergin, Selma; Dobrucalı, Erinç
2014-06-01
The exhaust smoke dispersion for a generic frigate is investigated numerically through the numerical solution of the governing fluid flow, energy, species and turbulence equations. The main objective of this work is to obtain the effects of the yaw angle, velocity ratio and buoyancy on the dispersion of the exhaust smoke. The numerical method is based on the fully conserved control-volume representation of the fully elliptic Navier-Stokes equations. Turbulence is modeled using a two-equation ( k- ɛ) model. The flow visualization tests using a 1/100 scale model of the frigate in the wind tunnel were also carried out to determine the exhaust plume path and to validate the computational results. The results show that down wash phenomena occurs for the yaw angles between ψ =10° and 20°. The results with different exhaust gas temperatures show that the buoyancy effect increases with the increasing of the exhaust gas temperature. However, its effect on the plume rise is less significant in comparison with its momentum. A good agreement between the predictions and experiment results is obtained.
Mixing of a point-source indoor pollutant: Numerical predictions and comparison with experiments
Lobscheid, C.; Gadgil, A.J.
2002-01-01
In most practical estimates of indoor pollutant exposures, it is common to assume that the pollutant is uniformly and instantaneously mixed in the indoor space. It is also commonly known that this assumption is simplistic, particularly for point sources, and for short-term or localized indoor exposures. We report computational fluid dynamics (CFD) predictions of mixing time of a point-pulse release of a pollutant in an unventilated mechanically mixed isothermal room. We aimed to determine the adequacy of the standard RANS two-equation ({kappa}-{var_epsilon}) turbulence model to predict the mixing times under these conditions. The predictions were made for the twelve mixing time experiments performed by Drescher et al. (1995). We paid attention to adequate grid resolution, suppression of numerical diffusion, and careful simulation of the mechanical blowers used in the experiments. We found that the predictions are in good agreement with experimental measurements.
The sensitivity of the general circulation to Arctic Sea ice boundaries - A numerical experiment
NASA Technical Reports Server (NTRS)
Herman, G. F.; Johnson, W. T.
1978-01-01
Results are presented for a set of numerical experiments conducted with the Goddard (formerly GISS) general circulation model. The experiments were designed to test the model atmospheric response to a single fixed and specified parameter, the total ice cover in the Davis Strait, Barents Sea, East Greenland Sea, Sea of Okhotsk and Bering Sea. Margin variations are considered that are substantially smaller than those involved in ice age or ice-free Arctic simulations. Anomaly is defined as the mean of two runs corresponding to climatological maximum sea ice conditions. Model results indicate that the ice margin anomalies are capable of altering local climates in certain regions of high and middle latitudes. Possible interactions between high latitudes and subtropical regions are suggested.
AN EXPERIMENT AND NUMERICAL SIMULATION OF SLURRY-EROSION CAUSED BY BEDLOAD ON BEDROCK
NASA Astrophysics Data System (ADS)
Oikawa, Shin; Iwasaki, Toshiki; Yamaguchi, Satomi; Shimizu, Yasuyuki; Kimura, Ichiro
The bed degradation is going rapidly in the upper part of the Ishikari River recently, owing to exposed bedrock erosion caused by bedload on the bedrock. In order to estimate the erosion rate of bedrock, we performed an experiment by a circular channel flume with artificial bedrock made of plaster. We applied the erosion rate estimated by the present experiment to a numerical simulation of bed deformation on the bedrock by using horizontal 2D flow model and the bedload layer model. The result of the simulation shows that gut pattern appeared in the simulated results which is very similar to the experimental result. It was found that the bed degradation was progressed by the erosion caused by sediment transport itself.
NASA Astrophysics Data System (ADS)
Qu, W. S.; Tan, L.; Cao, S. L.; Xu, Y.; Huang, J.; Xu, Q. H.
2015-01-01
As a kind of widely used device in pipe system for pressure and flow rate regulating, the valve would experience cavitation in the case when a sharp pressure drop occurs, which will induce the energy loss, noise and vibration of pipeline system, and even operational accidents. The experiment on flow resistance coefficient of a DN600 pressure-regulating valve under operation conditions from 0% to 100% openings is conducted. Based on the RNG k-e turbulence model and the Rayleigh-Plesset cavitation equation, a set of computational model is developed to simulate the turbulent flow in the valve under operational conditions from 0% to 100% openings. The computational results of flow resistance coefficient are compared to the experimental data. And the numerical simulation is employed to predict the cavitation performance of the valve at different inlet flow conditions. The transient cavitating flow is calculated to reveal the time evolution of cavitation in the valve.
NASA Astrophysics Data System (ADS)
Klassen, I.; Hillebrand, G.; Olsen, N. R. B.; Vollmer, S.; Lehmann, B.; Nestmann, F.
2013-10-01
The prediction of cohesive sediment transport requires numerical models which include the dominant physico-chemical processes of fine sediments. Mainly in terms of simulating small scale processes, flocculation of fine particles plays an important role since aggregation processes affect the transport and settling of fine-grained particles. Flocculation algorithms used in numerical models are based on and calibrated using experimental data. A good agreement between the results of the simulation and the measurements is a prerequisite for further applications of the transport functions. In this work, the sediment transport model (SSIIM) was extended by implementing a physics-based aggregation process model based on McAnally (1999). SSIIM solves the Navier-Stokes-Equations in a three-dimensional, non-orthogonal grid using the k-ɛ turbulence model. The program calculates the suspended load with the convection-diffusion equation for the sediment concentration. Experimental data from studies in annular flumes (Hillebrand, 2008; Klassen, 2009) is used to test the flocculation algorithm. Annular flumes are commonly used as a test rig for laboratory studies on cohesive sediments since the flocculation processes are not interfered with by pumps etc. We use the experiments to model measured floc sizes, affected by aggregation processes, as well as the sediment concentration of the experiment. Within the simulation of the settling behavior, we use different formulas for calculating the settling velocity (Stokes, 1850 vs. Winterwerp, 1998) and include the fractal dimension to take into account the structure of flocs. The aim of the numerical calculations is to evaluate the flocculation algorithm by comparison with the experimental data. The results from these studies have shown, that the flocculation process and the settling behaviour are very sensitive to variations in the fractal dimension. We get the best agreement with measured data by adopting a characteristic fractal
Additional experiments relative to the shelf life of Li(Si)/FeS2 thermal batteries
NASA Astrophysics Data System (ADS)
Searcy, J. Q.; Armijo, J. R.
1985-02-01
A continuing effort to develop a new thermal battery technology based on the Li(Si)/FeS2 electrochemical couple is reported. The results relate to the long shelf life requirement for thermal batteries designed by Sandia, and include topics relevant to leakage through the hermetic seal and accelerated aging experiments with materials new to the technology. Conclusions relevant to leakage through the hermetic seal are that the maximum leak rate must not exceed 1.8 x 10(-7) w, where w is the grams of Li(Si) contained by a battery, and that a bomb type leak test can be designed that is adequate for most Li(Si)/FeS2 batteries. Conclusions relevant to long term compatibility of new materials include the following: nickel is not compatible with the iron disulfide in the cathode; the CaSi2 additive used to suppress the initial voltage transient does not react or degrade during accelerated aging experiments, but the use of that material can lead to an increase in the variability of the activated lives, especially for long life batteries; Grafoil current collectors used with the cathode do not degrade in accelerated aging experiments.
Numerical simulatin of supernova-relevant laser-driven hydro experiments on OMEGA
Leibrandt, D; Robey, H F; Edwards, M J; Braun, D G; Miles, A R; Drake, R P
2004-02-10
In ongoing experiments performed on the OMEGA laser [J. M. Soures et al., Phys. Plasmas 5, 2108 (1996)] at the University of Rochester Laboratory for Laser Energetics (LLE), nanosecond laser pulses are used to drive strong blast waves into two-layer targets. Perturbations on the interface between the two materials are unstable to the Richtmyer-Meshkov instability as a result of shock transit and the Rayleigh-Taylor instability during the deceleration-phase behind the shock front. These experiments are designed to produce a strongly shocked interface whose evolution is a scaled version of the unstable hydrogen-helium interface in core-collapse supernovae such as SN 1987A. The ultimate goal of this research is to develop an understanding of the effect of hydrodynamic instabilities and the resulting transition to turbulence on supernovae observables that remain as yet unexplained. The authors are, at present, particularly interested in the development of the Rayleigh-Taylor instability through the late nonlinear stage, the transition to turbulence, and the subsequent transport of material within the turbulent region. In this paper, the results of numerical simulations of 2D single and multimode experiments are presented. These simulations are run using the 2D Arbitrary Lagrangian Eulerian (ALE) radiation hydrodynamics code CALE [R. T. Barton, Numerical Astrophysics (Jones and Bartlett, Boston, 1985)]. The simulation results are shown to compare well with experimental radiography. A buoyancy-drag model captures the behavior of the single-mode interface, but gives only partial agreement in the multi-mode cases. The Richtmyer-Meshkov and target decompression contributions to the perturbation growth are both estimated and shown to be significant. Significant dependence of the simulation results on the material equation of state (EOS) is demonstrated, and the prospect of continuing the experiments to conclusively demonstrate the transition to turbulence is discussed.
NASA Astrophysics Data System (ADS)
Richardson, Derek C.; Blum, J.; Weinhart, T.; Schwartz, S. R.; Michel, P.; Walsh, K. J.
2012-10-01
Spacecraft images of asteroids show evidence of low-gravity granular flows. Interpretation of these flows requires numerical modeling, which in turn requires code validation at laboratory scales. We have implemented a soft-sphere discrete element method (SSDEM) for modeling granular flows in our numerical code (Schwartz et al. 2012, Granular Matter 14, 363). Here we present results from a study to calibrate our code against controlled landslide experiments in order to determine the SSDEM parameters that best match real materials, to see how changes in those parameters affect the flow, and to mimic effects such as those due to irregular particle shapes. The apparatus, designed at University of Braunschweig, is a 0.6 × 0.8 m enclosed bed with a surface comprised of 10 mm diameter glass spheres glued into precisely drilled holes in a metal plate. The exact positions and depths of each of these glued spheres are input to the simulations. The experiments consist of filling the apparatus with loose glass beads (also 10 mm diameter) up to a set depth then gradually tilting the bed to note the angle of landslide initiation and the characteristics of the resulting flow. We reproduce this procedure in simulations, which we find are quite sensitive to the adopted SSDEM parameters, e.g., rolling friction and tangential damping delay landslide onset, while higher particle elasticity gives rise to faster, shorter-duration landslides. Preliminary results show a best match to the experiments (landslide initiation around 25 degrees) when adopting low static friction and no rolling friction in the simulations, but more experiments are in process. In future work, we will perform simulations in low-gravity environments representative of asteroid surfaces. This work is supported in part by grant NNX08AM39G from the NASA Office of Space Science. This study resulted from International Team collaboration #202 sponsored by ISSI in Switzerland.
The Zombie Instability: Using Numerical Simulation to Design a Laboratory Experiment
NASA Astrophysics Data System (ADS)
Wang, Meng; Pei, Suyang; Jiang, Chung-Hsiang; Hassanzadeh, Pedram; Marcus, Philip
2014-11-01
A new type of finite amplitude-instability has been found in numerical simulations of stratified, rotating, shear flows. The instability occurs via baroclinic critical layers that create linearly unstable vortex layers, which roll-up into vortices. Under the right conditions, those vortices can form a new generation of vortices, resulting in ``vortex self-replication'' that fills the fluid with vortices. Creating this instability in a laboratory would provide further evidence for the existence of the instability, which we first found in numerical simulations of protoplanetary disks. To design a laboratory experiment we need to know how the flow parameters-- shear, rotation and stratification, etc. affect the instability. To build an experiment economically, we also need to know how the finite-amplitude trigger of the instability scales with viscosity and the size of the domain. In this talk, we summarize our findings. We present a map, in terms of the experimentally controllable parameters, that shows where the instability occurs and whether the instability creates a few isolated transient vortices, a few long-lived vortices, or long-lived, self-replicating vortices that fill the entire flow.
Kamm, J.R.; Bos, R.J.
1995-06-01
In this paper the authors discuss numerical simulations of the Non-Proliferation Experiment (NPE), which was an underground explosion conducted in September 1993 in the volcanic tuff of the Nevada Test Site. The NPE source consisted of 1.29 {times} 10{sup 6} kg of ANFO-emulsion blasting agent, with the approximate energy of 1.1 kt, emplaced 389 m beneath the surface of Rainier Mesa. The authors compare detailed numerical simulations of the NPE with data collected from that experiment, and with calculations of an equally energetic nuclear explosion in identical geology. Calculated waveforms, at ranges out to approximately 1 km, agree moderately well in the time domain with free-field data, and are in qualitative agreement with free-surface records. Comparison of computed waveforms for equally energetic chemical and nuclear sources reveals relatively minor differences beyond the immediate near-source region, with the chemical source having an {approximately}25% greater seismic moment but otherwise indistinguishable (close-in) seismic source properties. 41 refs., 67 figs., 7 tabs.
Local mechanical properties of LFT injection molded parts: Numerical simulations versus experiments
NASA Astrophysics Data System (ADS)
Desplentere, F.; Soete, K.; Bonte, H.; Debrabandere, E.
2014-05-01
In predictive engineering for polymer processes, the proper prediction of material microstructure from known processing conditions and constituent material properties is a critical step forward properly predicting bulk properties in the finished composite. Operating within the context of long-fiber thermoplastics (LFT, length < 15mm) this investigation concentrates on the prediction of the local mechanical properties of an injection molded part. To realize this, the Autodesk Simulation Moldflow Insight 2014 software has been used. In this software, a fiber breakage algorithm for the polymer flow inside the mold is available. Using well known micro mechanic formulas allow to combine the local fiber length with the local orientation into local mechanical properties. Different experiments were performed using a commercially available glass fiber filled compound to compare the measured data with the numerical simulation results. In this investigation, tensile tests and 3 point bending tests are considered. To characterize the fiber length distribution of the polymer melt entering the mold (necessary for the numerical simulations), air shots were performed. For those air shots, similar homogenization conditions were used as during the injection molding tests. The fiber length distribution is characterized using automated optical method on samples for which the matrix material is burned away. Using the appropriate settings for the different experiments, good predictions of the local mechanical properties are obtained.
NASA Astrophysics Data System (ADS)
Piazolo, Sandra; Montagnat, Maurine; Prakash, Abhishek; Borthwick, Verity; Evans, Lynn; Griera, Albert; Bons, Paul D.; Svahnberg, Henrik; Prior, David J.
2015-04-01
Understanding the influence of the pre-existing microstructure on subsequent microstructural development is pivotal for the correct interpretation of rocks and ice that stayed at high homologous temperatures over a significant period of time. The microstructural behaviour of these materials through time has an important bearing on the interpretation of characteristics such as grain size, for example, using grain size statistics to detect former high strain zones that remain at high temperatures but low stress. We present a coupled experimental and modelling approach to better understand the evolution of recrystallization characteristics as a function of deformation-annealing time paths in a material with a high viscoplastic anisotropy e.g. polycrystalline ice and magnesium alloys. Deformation microstructures such as crystal bending, subgrain boundaries, grain size variation significantly influence the deformation and annealing behaviour of crystalline material. For numerical simulations we utilize the microdynamic modelling platform, Elle (www.elle.ws), taking local microstructural evolution into account to simulate the following processes: recovery within grains, rotational recrystallization, grain boundary migration and nucleation. We first test the validity of the numerical simulations against experiments, and then use the model to interpret microstructural features in natural examples. In-situ experiments are performed on laboratory grown and deformed ice and magnesium alloy. Our natural example is a deformed then recrystallized anorthosite from SW Greenland. The presented approach can be applied to many other minerals and crystalline materials.
Optimization of CFETR CSMC cabling based on numerical modeling and experiments
NASA Astrophysics Data System (ADS)
Qin, Jinggang; Dai, Chao; Liu, Bo; Wu, Yu; Liu, Fang; Liao, Guojun; Xue, Tianjun; Wei, Zhourong; Nijhuis, Arend; Zhou, Chao; Devred, Arnaud
2015-12-01
The China Fusion Engineering Test Reactor (CFETR) is a new tokamak device, whose magnet system includes toroidal field (TF), central solenoid (CS) and poloidal field (PF) coils. The main goal is to build a fusion engineering tokamak reactor with 50-200 MW fusion power and self-sufficiency by blanket, which means that the deuterium-tritium reaction in the plasma produces neutrons and alpha particles, and the neutrons react with the lithium-containing blanket surrounding the plasma, breeding the tritium by lithium-neutron reaction. To develop the manufacturing technique for the full-size CS coil, the Central Solenoid Model Coil (CSMC) project for CFETR was launched first. A Nb3Sn conductor is to be used in the CFETR CSMC, whose design refers to the ITER CS conductor with the same short-twist-pitch cable pattern. Due to the short twist pitch and relatively low void fraction, a high compaction ratio is required during cabling and the risk of strand damage is increased significantly. Although it is impossible to avoid strand deformation for this design, it is crucial to find a way to reduce strand damage as much as possible. A numerical model was used to analyze the causes of strand damage, including variation in twist pitch length as well as different mechanical properties for copper and Nb3Sn strands. Several experiments have been performed to verify the numerical results, including cabling trials for different conditions and critical current (I c) tests on strands with/without deformation. The results show that the numerical analysis is consistent with the experiments and provides the optimal cabling conditions for the CFETR CSMC.
NASA Astrophysics Data System (ADS)
Spina, Laura; Colucci, Simone; De'Michieli Vitturi, Mattia; Scheu, Bettina; Dingwell, Donald Bruce
2015-04-01
The study of the fluid dynamics of magmatic melts into the conduit, where direct observations are unattainable, was proven to be strongly enhanced by multiparametric approaches. Among them, the coupling of numerical modeling with laboratory experiments represents a fundamental tool of investigation. Indeed, the experimental approach provide invaluable data to validate complex multiphase codes. We performed decompression experiments in a shock tube system, using pure silicon oil as a proxy for the basaltic melt. A range of viscosity comprised between 1 and 1000 Pa s was investigated. The samples were saturated with Argon for 72h at 10MPa, before being slowly decompressed to atmospheric pressure. The evolution of the analogue magmatic system was monitored through a high speed camera and pressure sensors, located into the analogue conduit. The experimental decompressions have then been reproduced numerically using a multiphase solver based on OpenFOAM framework. The original compressible multiphase Openfoam solver twoPhaseEulerFoam was extended to take into account the multicomponent nature of the fluid mixtures (liquid and gas) and the phase transition. According to the experimental conditions, the simulations were run with values of fluid viscosity ranging from 1 to 1000 Pa s. The sensitivity of the model has been tested for different values of the parameters t and D, representing respectively the relaxation time for gas exsolution and the average bubble diameter, required by the Gidaspow drag model. Valuable range of values for both parameters are provided from experimental observations, i.e. bubble nucleation time and bubble size distribution at a given pressure. The comparison of video images with the outcomes of the numerical models was performed by tracking the evolution of the gas volume fraction through time. Therefore, we were able to calibrate the parameter of the model by laboratory results, and to track the fluid dynamics of experimental decompression.
Numerical Modeling of X-ray Photoionization Experiments Driven by Z-Pinch X-rays
NASA Astrophysics Data System (ADS)
Shupe, N. C.; Cohen, D. H.; MacFarlane, J. J.
2004-12-01
We have performed an initial round of experiments at the Z-Machine at Sandia National Laboratory in an attempt to create and characterize an X-ray photoionized plasma that is analogous to those found in X-ray binaries and AGNs. The ultimate goal is to benchmark X-ray spectral modeling codes that are used to analyze Chandra and XMM data from accretion powered astrophysical objects. The initial experiments involved neon and the primary measurement made was time-integrated, back-lit X-ray absorption spectroscopy of the photoionized neon. We present numerical modeling of this experiment, including non-LTE radiation hydrodynamics and spectral synthesis results, that are in good agreement with the data. We also present scaling studies for future experiments, including sythesized time-resolved X-ray emission spectra that correspond to the high-resolution spectral data being produced by the current generation of X-ray telescopes. The authors acknowledge the support of Research Corporation grant CC5489.
Ferreira, R C; Nascimento-Junior, A B; Santos, P J P; Botter-Carvalho, M L; Pinto, T K
2015-01-15
An experiment was carried out on an estuarine mudflat to assess impacts of inorganic nutrients used to fertilize sugar-cane fields on the surrounding aquatic ecosystem, through changes in the nematode community structure. During 118 days, nine quadrats each 4m(2) were sampled six times after the beginning of fertilizer addition. The fertilizer was introduced weekly in six areas, at two different concentrations (low and high doses), and three areas were used as control. The introduction of nutrients modified key nematode community descriptors. In general, the nematodes were negatively affected over the study period. However, Comesa, Metachromadora, Metalinhomoeus, Spirinia and Terschellingia were considered tolerant, and other genera showed different degrees of sensitivity. Nutrient input also affect the availability and quality of food, changing the nematode trophic structure. The use of inorganic fertilizer should be evaluated with care because of the potential for damage to biological communities of coastal aquatic systems. PMID:25499965
NASA Astrophysics Data System (ADS)
Abd El-Halim, A. A.; Kumlung, Arunsiri
2015-01-01
Until now sandy soils can be considered as one roup having common hydrophysical problems. Therefore, a laboratory experiment was conducted to evaluate the influence of bagasse as an amendment to improve hydrophysical properties of sandy soil, through the determination of bulk density, aggregatesize distribution, total porosity, hydraulic conductivity, pore-space structure and water retention. To fulfil this objective, sandy soils were amended with bagasse at the rate of 0, 0.5, 1, 2, 3 and 4% on the dry weight basis. The study results demonstrated that the addition of bagasse to sandy soils in between 3 to 4% on the dry weight basis led to a significant decrease in bulk density, hydraulic conductivity, and rapid-drainable pores, and increase in the total porosity, water-holding pores, fine capillary pores, water retained at field capacity, wilting point, and soil available water as compared with the control treatment
Additional experiments on flowability improvements of aviation fuels at low temperatures, volume 2
NASA Technical Reports Server (NTRS)
Stockemer, F. J.; Deane, R. L.
1982-01-01
An investigation was performed to study flow improver additives and scale-model fuel heating systems for use with aviation hydrocarbon fuel at low temperatures. Test were performed in a facility that simulated the heat transfer and temperature profiles anticipated in wing fuel tanks during flight of long-range commercial aircraft. The results are presented of experiments conducted in a test tank simulating a section of an outer wing integral fuel tank approximately full-scale in height, chilled through heat exchange panels bonded to the upper and lower horizontal surfaces. A separate system heated lubricating oil externally by a controllable electric heater, to transfer heat to fuel pumped from the test tank through an oil-to-fuel heat exchanger, and to recirculate the heated fuel back to the test tank.
Direct comparisons between the results of a numerical model and laboratory experiments are presented for density-stratified flows around an isolated three-dimensional hill. The numerical model integrates the Navier-Stokes equations for incompressible stratified flow using a finit...
Loophole-free Bell test using electron spins in diamond: second experiment and additional analysis.
Hensen, B; Kalb, N; Blok, M S; Dréau, A E; Reiserer, A; Vermeulen, R F L; Schouten, R N; Markham, M; Twitchen, D J; Goodenough, K; Elkouss, D; Wehner, S; Taminiau, T H; Hanson, R
2016-01-01
The recently reported violation of a Bell inequality using entangled electronic spins in diamonds (Hensen et al., Nature 526, 682-686) provided the first loophole-free evidence against local-realist theories of nature. Here we report on data from a second Bell experiment using the same experimental setup with minor modifications. We find a violation of the CHSH-Bell inequality of 2.35 ± 0.18, in agreement with the first run, yielding an overall value of S = 2.38 ± 0.14. We calculate the resulting P-values of the second experiment and of the combined Bell tests. We provide an additional analysis of the distribution of settings choices recorded during the two tests, finding that the observed distributions are consistent with uniform settings for both tests. Finally, we analytically study the effect of particular models of random number generator (RNG) imperfection on our hypothesis test. We find that the winning probability per trial in the CHSH game can be bounded knowing only the mean of the RNG bias. This implies that our experimental result is robust for any model underlying the estimated average RNG bias, for random bits produced up to 690 ns too early by the random number generator. PMID:27509823
Loophole-free Bell test using electron spins in diamond: second experiment and additional analysis
Hensen, B.; Kalb, N.; Blok, M. S.; Dréau, A. E.; Reiserer, A.; Vermeulen, R. F. L.; Schouten, R. N.; Markham, M.; Twitchen, D. J.; Goodenough, K.; Elkouss, D.; Wehner, S.; Taminiau, T. H.; Hanson, R.
2016-01-01
The recently reported violation of a Bell inequality using entangled electronic spins in diamonds (Hensen et al., Nature 526, 682–686) provided the first loophole-free evidence against local-realist theories of nature. Here we report on data from a second Bell experiment using the same experimental setup with minor modifications. We find a violation of the CHSH-Bell inequality of 2.35 ± 0.18, in agreement with the first run, yielding an overall value of S = 2.38 ± 0.14. We calculate the resulting P-values of the second experiment and of the combined Bell tests. We provide an additional analysis of the distribution of settings choices recorded during the two tests, finding that the observed distributions are consistent with uniform settings for both tests. Finally, we analytically study the effect of particular models of random number generator (RNG) imperfection on our hypothesis test. We find that the winning probability per trial in the CHSH game can be bounded knowing only the mean of the RNG bias. This implies that our experimental result is robust for any model underlying the estimated average RNG bias, for random bits produced up to 690 ns too early by the random number generator. PMID:27509823
Piepel, Gregory F. ); Cooley, Scott K. ); Peeler, David K.; Vienna, John D. ); Edwards, Tommy B.
2002-01-01
A glass composition variation study (CVS) for high-level waste (HLW) stored in Idaho is being statistically designed and performed in phases over several years. The purpose of the CVS is to investigate and model how HLW-glass properties depend on glass composition. The resulting glass property-composition models will be used to develop desirable glass formulations and for other purposes. Phases 1 and 2 of the CVS have been completed and are briefly described. This paper focuses on the CVS Phase 3 experimental design, which was chosen to augment the Phase 1 and 2 data with additional data points, as well as to account for additional glass components not studied in Phases 1 and/or 2. In total, 16 glass components were varied in the Phase 3 experimental design. The paper describes how these Phase 3 experimental design augmentation challenges were addressed using the previous data, preliminary property-composition models, and statistical mixture experiment and optimal experimental design methods and software.
Lopez-Castroman, J; Guillaume, S; Olié, E; Jaussent, I; Baca-García, E; Courtet, P
2015-01-01
Family history of suicidal behavior and personal history of childhood abuse are reported risk factors for suicide attempts and suicide completion. We aim to quantify the additive effect of family history of suicidal behavior and different subtypes of childhood abuse on suicidal behavior. We examined a sample of 496 suicide attempters, comparing individuals with family history of suicidal behavior and personal history of childhood (physical or sexual) abuse, individuals with family history of suicidal behavior only, individuals with history of early traumatic experiences only, and individuals with none of these two risk factors with regards to suicidal features. An additive effect was found for the age at the first attempt in suicide attempters with both family history of suicidal behavior and either physical or sexual abuse. No significant interactions were found between family history of suicidal behavior and childhood trauma in relation to any characteristics of suicidal behavior. Subjects presenting family history of suicidal behavior and childhood abuse attempt suicide earlier in life than subjects with just one or none of them, particularly if they were sexually abused. Other suicidality indexes were only partially or not associated with this combination of risk factors. A careful assessment of patients with both family history of suicidal behavior and childhood abuse could help to prevent future suicide attempts, particularly in young people. PMID:25259671
NASA Astrophysics Data System (ADS)
Jeanniot, Ludovic; Kusznir, Nick; Manatschal, Gianreto; Cowie, Leanne
2013-04-01
The initiation of sea-floor spreading, during the continental breakup process, requires both the rupture of the continental crust and the initiation of decompression melting. This process results in mantle upwelling and at some point decompressional melting which creates new oceanic crust. Using numerical experiments, we investigate how the deformation mode of continental lithosphere thinning and stretching controls the rupture of continental crust and lithospheric mantle, the onset of decompression melting, their relative timing, and the circumstances under which mantle exhumation may occur. We assume that the topmost continental and ocean lithosphere, corresponding to the cooler brittle seismogenic layer, deforms by extensional faulting (pure-shear deformation) and magmatic intrusion, consistent with the observations of deformation processes occurring at slow spreading ocean ridges (Cannat, 1996). We assume that deformation beneath this topmost lithosphere layer (approximately 15-20 km thick) occurs in response to passive upwelling and thermal and melt buoyancy driven small-scale convection. We use a 2D finite element viscous flow model (FeMargin) to describe lithosphere and asthenosphere deformation. This flow field is used to advect lithosphere and asthenosphere temperature and material. The finite element model is kinematically driven by Vx for the topmost upper crust inducing passive upwelling beneath that layer. A vertical velocity Vz is defined for buoyancy enhanced upwelling as predicted by Braun et al. (2000). Melt generation is predicted by decompression melting using the parameterization and methodology of Katz et al. (2003). Numerical experiments have been used to investigate the dependency of continental crust and lithosphere rupture, decompression melt initiation, rifted margin ocean-continent transition architecture and subsidence history on the half-spreading rate Vx, buoyancy driven upwelling rate Vz, the relative contribution of these deformation
NASA Astrophysics Data System (ADS)
Poonoosamy, Jenna; Kosakowski, Georg; Van Loon, Luc R.; Mäder, Urs
2015-06-01
In the context of testing reactive transport codes and their underlying conceptual models, a simple 2D reactive transport experiment was developed. The aim was to use simple chemistry and design a reproducible and fast to conduct experiment, which is flexible enough to include several process couplings: advective-diffusive transport of solutes, effect of liquid phase density on advective transport, and kinetically controlled dissolution/precipitation reactions causing porosity changes. A small tank was filled with a reactive layer of strontium sulfate (SrSO4) of two different grain sizes, sandwiched between two layers of essentially non-reacting quartz sand (SiO2). A highly concentrated solution of barium chloride was injected to create an asymmetric flow field. Once the barium chloride reached the reactive layer, it forced the transformation of strontium sulfate into barium sulfate (BaSO4). Due to the higher molar volume of barium sulfate, its precipitation caused a decrease of porosity and lowered the permeability. Changes in the flow field were observed with help of dye tracer tests. The experiments were modelled using the reactive transport code OpenGeosys-GEM. Tests with non-reactive tracers performed prior to barium chloride injection, as well as the density-driven flow (due to the high concentration of barium chloride solution), could be well reproduced by the numerical model. To reproduce the mineral bulk transformation with time, two populations of strontium sulfate grains with different kinetic rates of dissolution were applied. However, a default porosity permeability relationship was unable to account for measured pressure changes. Post mortem analysis of the strontium sulfate reactive medium provided useful information on the chemical and structural changes occurring at the pore scale at the interface that were considered in our model to reproduce the pressure evolution with time.
Poonoosamy, Jenna; Kosakowski, Georg; Van Loon, Luc R; Mäder, Urs
2015-01-01
In the context of testing reactive transport codes and their underlying conceptual models, a simple 2D reactive transport experiment was developed. The aim was to use simple chemistry and design a reproducible and fast to conduct experiment, which is flexible enough to include several process couplings: advective-diffusive transport of solutes, effect of liquid phase density on advective transport, and kinetically controlled dissolution/precipitation reactions causing porosity changes. A small tank was filled with a reactive layer of strontium sulfate (SrSO4) of two different grain sizes, sandwiched between two layers of essentially non-reacting quartz sand (SiO2). A highly concentrated solution of barium chloride was injected to create an asymmetric flow field. Once the barium chloride reached the reactive layer, it forced the transformation of strontium sulfate into barium sulfate (BaSO4). Due to the higher molar volume of barium sulfate, its precipitation caused a decrease of porosity and lowered the permeability. Changes in the flow field were observed with help of dye tracer tests. The experiments were modelled using the reactive transport code OpenGeosys-GEM. Tests with non-reactive tracers performed prior to barium chloride injection, as well as the density-driven flow (due to the high concentration of barium chloride solution), could be well reproduced by the numerical model. To reproduce the mineral bulk transformation with time, two populations of strontium sulfate grains with different kinetic rates of dissolution were applied. However, a default porosity permeability relationship was unable to account for measured pressure changes. Post mortem analysis of the strontium sulfate reactive medium provided useful information on the chemical and structural changes occurring at the pore scale at the interface that were considered in our model to reproduce the pressure evolution with time. PMID:25805363
Tsunami-induced boulder transport - combining physical experiments and numerical modelling
NASA Astrophysics Data System (ADS)
Oetjen, Jan; Engel, Max; May, Simon Matthias; Schüttrumpf, Holger; Brueckner, Helmut; Prasad Pudasaini, Shiva
2016-04-01
Coasts are crucial areas for living, economy, recreation, transportation, and various sectors of industry. Many of them are exposed to high-energy wave events. With regard to the ongoing population growth in low-elevation coastal areas, the urgent need for developing suitable management measures, especially for hazards like tsunamis, becomes obvious. These measures require supporting tools which allow an exact estimation of impact parameters like inundation height, inundation area, and wave energy. Focussing on tsunamis, geological archives can provide essential information on frequency and magnitude on a longer time scale in order to support coastal hazard management. While fine-grained deposits may quickly be altered after deposition, multi-ton coarse clasts (boulders) may represent an information source on past tsunami events with a much higher preservation potential. Applying numerical hydrodynamic coupled boulder transport models (BTM) is a commonly used approach to analyse characteristics (e.g. wave height, flow velocity) of the corresponding tsunami. Correct computations of tsunamis and the induced boulder transport can provide essential event-specific information, including wave heights, runup and direction. Although several valuable numerical models for tsunami-induced boulder transport exist (e. g. Goto et al., 2007; Imamura et al., 2008), some important basic aspects of both tsunami hydrodynamics and corresponding boulder transport have not yet been entirely understood. Therefore, our project aims at these questions in four crucial aspects of boulder transport by a tsunami: (i) influence of sediment load, (ii) influence of complex boulder shapes other than idealized rectangular shapes, (iii) momentum transfers between multiple boulders, and (iv) influence of non-uniform bathymetries and topographies both on tsunami and boulder. The investigation of these aspects in physical experiments and the correct implementation of an advanced model is an urgent need
Determining the tube bundle streamlining critical parameters using the numerical experiment method
NASA Astrophysics Data System (ADS)
Kaplunov, S. M.; Val'es, N. G.; Samolysov, A. V.; Marchevskaya, O. A.
2015-08-01
The article is devoted to development and application of mathematical models describing the most dangerous mechanisms through which vibrations are excited in tube bundles and blunt cylindrically shaped structures, and to development of reliable calculation methods for describing these models, which would make it possible to obtain prompt data for designing and subsequent operation of the considered structural elements. For solving such problems, a comprehensive approach is required, which should be based on a combined use of numerical experiments on computers and experimental investigations on full-scale equipment. The authors have developed a procedure for numerically investigating the hydrodynamic forces arising during stalled streamlining and the tube bundle vibrations caused by these forces. The procedure is based on using the developed mathematical model describing fluid-elastic excitation of vibrations in a bundle of elastic tubes placed in external cross flow. The problem of studying fluid-elastic excitation is brought to stability analysis, which is carried out with the assumption about a linear behavior of destabilizing forces for undisturbed state of elastic tubes. A theoretical investigation of the developed mathematical model was carried out, from which the necessary and sufficient condition of system stability has been obtained in terms of system dimensionless parameters (mass, damping, and velocity). An algorithm for numerically determining the matrices of linear hydrodynamic coupling coefficients for particular tube bundles is developed. The validity of the algorithm and the computer programs developed on its basis are checked by comparing the results of test calculations with the bank of known experimental data. A procedure is proposed for determining the matrices of linear hydrodynamic coupling coefficients in bundles having a regular layout of their cross section and a large number of tubes through calculating these matrices for a relatively small
NASA Astrophysics Data System (ADS)
Nassar, M.; Ginn, T. R.; Le Borgne, T.; Schreyer, L. G.; Dentz, M.
2015-12-01
The challenge in characterizing mixing-limited reaction rates between displacing and displaced groundwater solutions (Figure 1, left panel) is the in quantifying mixing extent that is controlled by small scale heterogeneity. We describe limited numerical 2D testing of the lamella approach that focuses on the deformation of the moving front, treated as a set of linearized patches termed lamellae. We simulate flow and reactive transport in a recently characterized sample of Massillon sandstone to provide Eulerian test data for comparison with the new Lagrangian lamella-based solution. In our numerical experiments particle tracking is used to approximate the lamellar positions and deformations at any given time, and reactions are calculated on each lamella in proportion to the local scalar dissipation rate. The simulated data show the effect of small scale heterogeneity including strong shearing and local collapse or coalescence (Figure 1, right panel) of the reaction front on the global reaction rate. We propose a simple approximation to handle coalescence in the lamella-based upscaling and we test it against the simulated data.
Haeufle, D F B; Günther, M; Blickhan, R; Schmitt, S
2011-01-01
Recently, the hyperbolic Hill-type force-velocity relation was derived from basic physical components. It was shown that a contractile element CE consisting of a mechanical energy source (active element AE), a parallel damper element (PDE), and a serial element (SE) exhibits operating points with hyperbolic force-velocity dependency. In this paper, the contraction dynamics of this CE concept were analyzed in a numerical simulation of quick release experiments against different loads. A hyperbolic force-velocity relation was found. The results correspond to measurements of the contraction dynamics of a technical prototype. Deviations from the theoretical prediction could partly be explained by the low stiffness of the SE, which was modeled analog to the metal spring in the hardware prototype. The numerical model and hardware prototype together, are a proof of this CE concept and can be seen as a well-founded starting point for the development of Hill-type artificial muscles. This opens up new vistas for the technical realization of natural movements with rehabilitation devices. PMID:22275541
NASA Astrophysics Data System (ADS)
Lahaye, Noé; Paci, Alexandre; Smith, Stefan Llewellyn
2016-04-01
We examine the instability of lenticular vortices -- or lenses -- in a stratified rotating fluid. The simplest configuration is one in which the lenses overlay a deep layer and have a free surface, and this can be studied using a two-layer rotating shallow water model. We report results from laboratory experiments and high-resolution direct numerical simulations of the destabilization of vortices with constant potential vorticity, and compare these to a linear stability analysis. The stability properties of the system are governed by two parameters: the typical upper-layer potential vorticity and the size (depth) of the vortex. Good agreement is found between analytical, numerical and experimental results for the growth rate and wavenumber of the instability. The nonlinear saturation of the instability is associated with conversion from potential to kinetic energy and weak emission of gravity waves, giving rise to the formation of coherent vortex multipoles with trapped waves. The impact of flow in the lower layer is examined. In particular, it is shown that the growth rate can be strongly affected and the instability can be suppressed for certain types of weak co-rotating flow.
Numerical Design Of Experiments to Analyse the Contact Conditions in Microforming
Barbier, C.; Thibaud, S.; Picart, P.; Chambert, J.
2007-05-17
In microforming, the so-called size effects can be observed in the material flow behaviour as well as in the frictional behaviour. In order to study the frictional behaviour a preliminary numerical characterization of the surface tribology has been carried out. A numerical design of experiments (DOE) is based on cylinder upsetting tests to define the influence of surface geometric properties on the resultant force. The simulations have been performed with the finite element software LS-Dyna by using an axisymmetric model. The mechanical behaviour of the cylinder specimen was described by an elastic-plastic material law, whereas the upsetting plates were assumed to be rigid. The workpiece is considered to be a copper alloy (CuZn10). The average roughness Ra and the average mean spacing Sm have been chosen to describe surface roughness properties. The tool and workpiece surfaces have been modelled using a sinusoidal profile. The five input parameters of the DOE are the amplitude and the period of the two sinusoidal profiles and the phase displacement between them. The analysis of variance shows the statistically significant parameters or interactions.
NASA Astrophysics Data System (ADS)
Hu, Q.; Li, Y.; Pan, H. L.; Liu, J. T.; Zhuang, B. T.
2015-01-01
Vane type propellant management device (PMD) is one of the key components of the vane-type surface tension tank (STT), and its fluid orbital performance directly determines the STT's success or failure. In present paper, numerical analysis and microgravity experiment study on fluid orbital performance of a vane type PMD were carried out. By using two-phase flow model of volume of fluid (VOF), fluid flow characteristics in the tank with the vane type PMD were numerically calculated, and the rules of fluid transfer and distribution were gotten. A abbreviate model test system of the vane type PMD is established and microgravity drop tower tests were performed, then fluid management and transmission rules of the vane type PMD were obtained under microgravity environment. The analysis and tests results show that the vane type PMD has good and initiative fluid orbital management ability and meets the demands of fluid orbital extrusion in the vane type STT. The results offer valuable guidance for the design and optimization of the new generation of vane type PMD, and also provide a new approach for fluid management and control in space environment.
A numerical tool for reproducing driver behaviour: experiments and predictive simulations.
Casucci, M; Marchitto, M; Cacciabue, P C
2010-03-01
This paper presents the simulation tool called SDDRIVE (Simple Simulation of Driver performance), which is the numerical computerised implementation of the theoretical architecture describing Driver-Vehicle-Environment (DVE) interactions, contained in Cacciabue and Carsten [Cacciabue, P.C., Carsten, O. A simple model of driver behaviour to sustain design and safety assessment of automated systems in automotive environments, 2010]. Following a brief description of the basic algorithms that simulate the performance of drivers, the paper presents and discusses a set of experiments carried out in a Virtual Reality full scale simulator for validating the simulation. Then the predictive potentiality of the tool is shown by discussing two case studies of DVE interactions, performed in the presence of different driver attitudes in similar traffic conditions. PMID:19249745
Numerical experiments on the long-term morphodynamics of the Colorado River Delta
NASA Astrophysics Data System (ADS)
Montaño, Yovani; Carbajal, Noel
2008-03-01
The interaction among tidal currents, sediment transport, and long-term changes of the sea bottom in the Colorado River Delta have been investigated applying a two dimensional nonlinear hydrodynamic-numerical finite differences model. The system was forced by the dominant M 2 tidal component at the open boundary. We carried out calculations to study the morphodynamics of the actual bathymetry caused by the bedload sediment transport. To investigate the origin of actual morphological features, we performed experiments using a smoothed bathymetry, in which the islands Montague, Gore, and Pelícano were eliminated. Under the imposed tidal hydrodynamics, the results indicate that the bedload transport contributes significantly in the genesis of sandbanks and in the formation and maintaining of the Montague and Gore Islands.
Fundamental Experiments and Numerical Analyses on Heat Transfer Characteristics of a Vapor Chamber
NASA Astrophysics Data System (ADS)
Koito, Yasushi; Imura, Hideaki; Mochizuki, Masataka; Saito, Yuji; Torii, Shuichi
A vapor chamber is used as a novel heat spreader to cool high-performance MPUs (microprocessor units). The vapor chamber is placed between small heat sources and a large heat sink. This paper describes the effect of heat source size on the heat transfer characteristics of the vapor chamber. First, by the experiments, the effect of heat source size on the temperature distribution of the vapor chamber is investigated, and the validity of the mathematical model of the vapor chamber is confirmed. Secondly, by the numerical analyses, the effect of heat source size on the thermal resistances inside the vapor chamber is discussed. It is found that the heat source size greatly affects the thermal resistance of the evaporator section inside the vapor chamber. Although the thermal resistance is hardly affected by the heat generation rate and the heat flux of the heat source, it increases as the heat source becomes smaller.
Lefrancois, A.; Hare, D.; L'Eplattenier, P.; Burger, M.
2006-02-13
Isentropic compression experiments and numerical simulations on LX-04 (HMX / Viton 85/15) were performed respectively at Z accelerator facility from Sandia National Laboratory and at Lawrence Livermore National Laboratory in order to study the isentrope and associated Hugoniot of this HE. 2D and 3D configurations have been calculated here to test the new beta version of the electromagnetism package coupled with the dynamics in Ls-Dyna and compared with the ICE Z shot 1067 on LX 04. The electromagnetism module is being developed in the general-purpose explicit and implicit finite element program LS-DYNA{reg_sign} in order to perform coupled mechanical/thermal/electromagnetism simulations. The Maxwell equations are solved using a Finite Element Method (FEM) for the solid conductors coupled with a Boundary Element Method (BEM) for the surrounding air (or vacuum). More details can be read in the references.
A Numerical Experiment on the Role of Surface Shear Stress in the Generation of Sound
NASA Technical Reports Server (NTRS)
Shariff, Karim; Wang, Meng; Merriam, Marshal (Technical Monitor)
1996-01-01
The sound generated due to a localized flow over an infinite flat surface is considered. It is known that the unsteady surface pressure, while appearing in a formal solution to the Lighthill equation, does not constitute a source of sound but rather represents the effect of image quadrupoles. The question of whether a similar surface shear stress term constitutes a true source of dipole sound is less settled. Some have boldly assumed it is a true source while others have argued that, like the surface pressure, it depends on the sound field (via an acoustic boundary layer) and is therefore not a true source. A numerical experiment based on the viscous, compressible Navier-Stokes equations was undertaken to investigate the issue. A small region of a wall was oscillated tangentially. The directly computed sound field was found to to agree with an acoustic analogy based calculation which regards the surface shear as an acoustically compact dipole source of sound.
NASA Astrophysics Data System (ADS)
Timofeyev, Yuriy M.
1995-12-01
Regular, long-term, g1obalscale measurements of atmospheric minor gaseous and aerosol composition (MGAC) by means ofdifferent instruments (PHOENIX, OZONE-MIR, ISTOK-1, DOPI) are planned on board the SPECTR and PRIRODA modules of the Space Station MIR during 1995-1998. The main characteristics of these devices are given. The principal goals of the space experiments are: investigations of the spatial and temporal MGAC variations, comparisons of different space-borne atmospheric chemistry sensors and their intercalibration, validation of the space MGAC measurements using different ground-based station and aircraft data, and studies of the molecular absorption in the atmosphere aimed to enhance an accuracy of radiative transfer atmospheric models. Special attention is . devoted to radiative transfer model (line-mixing, line-shift, line-narrowing, Non-LTE effects). The numerical estimations of the errors of the MGAC vertical profile retrievals using different device data are carried out.
Numerical modeling experiments of coastal upwelling at the field of Arctic fjords.
NASA Astrophysics Data System (ADS)
Kosecki, Szymon; Dzierzbicka-Głowacka, Lidia
2016-04-01
Coastal upwelling is a well described, known phenomenon in theory. Nowadays there is more and more both environmental and modeling studies about it. Upwelling especially in the Arctic fjords is a process that strongly affects hydrodynamics and even more ecosystems. It is so important, that it brings detailed question about effects and needed wind driven forcing parameters. My modeling experiment studies were strongly different than the studies that are typically carried out using numerical models. Instead of searching for this phenomenon in modeled analysis or environmental data, I did several case scenarios simulations. For those I used statistically selected wind data measured in-stiu. The hi-resolution coastal mapping, the flexible mesh discretization method and the sigma-layered three dimensional model MIKE 3 by DHI allowed me to explore this phenomenon with very good accuracy. This studies have been done in Institute of Oceanology PAS in Sopot, as a part of Centre for Polar Studies.
Rodriguez, S; Jacob, X; Gibiat, V
2016-05-01
Using transducer arrays and appropriate emission delays allow to focus acoustic waves at a chosen location in a medium. The focusing spatial accuracy depends on the accurate knowledge of its acoustic properties. When those properties are unknown, methods based on the Time-Reversal principle allow accurate focusing. Still, these methods are either intrusive (an active source has to be introduced at the target location first), either blind (the target cannot be selected in the presence of several objects.) The purpose of the present work is to achieve non-invasive accurate focusing on a selected target using inaccurate acoustic properties for the investigated medium. Potential applications are for instance noninvasive surgery based on High Intensity Focused Ultrasound (HIFU). Numerical experiments are presented and demonstrate accurate focusing on a previously designated target located in an unknown heterogeneous medium. PMID:26890791
NASA Astrophysics Data System (ADS)
Chillara, Vamshi Krishna; Lissenden, Cliff J.
2016-01-01
Interest in using the higher harmonic generation of ultrasonic guided wave modes for nondestructive evaluation continues to grow tremendously as the understanding of nonlinear guided wave propagation has enabled further analysis. The combination of the attractive properties of guided waves with the attractive properties of higher harmonic generation provides a very unique potential for characterization of incipient damage, particularly in plate and shell structures. Guided waves can propagate relatively long distances, provide access to hidden structural components, have various displacement polarizations, and provide many opportunities for mode conversions due to their multimode character. Moreover, higher harmonic generation is sensitive to changing aspects of the microstructures such as to the dislocation density, precipitates, inclusions, and voids. We review the recent advances in the theory of nonlinear guided waves, as well as the numerical simulations and experiments that demonstrate their utility.
NASA Astrophysics Data System (ADS)
Güldemeister, N.; Moser, D.; Wünnemann, K.; Hoerth, T.; Schäfer, F.
2013-09-01
Meteorite impacts can cause environmental consequences, one of which is the generation of ground motions that may exceed the magnitude of the largest earthquakes [1]. Impacts generate shock waves that attenuate with distance until they even tually turn into seismic waves. Thus, meteorite impact may be considered as a source for seismic shaking similar to earthquakes. Seismic signals have been recorded in explosion experiments [2] and in hydrocode models of large impact events such as the Chicxulub crater [3]. To determine how much of the kinetic energy Ekin of the impactoris turned into seismic energy Eseis can be investigated experimentally (by recording the acoustic emission) or by numerical models. The ratio of Eseis/Ekin is the so called seismic efficiency k. The seismic efficiency depends on material properties (porosity) and is usually estimated to range between 10-2 and 10-6 [2,4]. In the framework of the "MEMIN" (multidisciplinary experimental and modeling impact crater research network) project a suite of hypervelocity impact experiments on a decimeter scale have been carried out [5]. We use acoustic emission (AE) technique and pressure gauges in high spatiotemporal Meteorite impacts can cause environmental consequences, one of which is the generation of ground motions that may exceed the magnitude of the largest earthquakes [1]. Impacts generate shock waves that attenuate with distance until they even tually turn into seismic waves. Thus, meteorite impact may be considered as a source for seismic shaking similar to earthquakes. Seismic signals have been recorded in explosion experiments [2] and in hydrocode models of large impact events such as the Chicxulub crater [3]. To determine how much of the kinetic energy Ekin of the impactoris turned into seismic energy Eseis can be investigated experimentally (by recording the acoustic emission) or by numerical models. The ratio of Eseis/Ekin is the so called seismic efficiency k. The seismic efficiency depends
Koch, Jeffrey A.; Bernat, Thomas P.; Collins, Gilbert W.; Hammel, Bruce A.; MacKinnon, Andrew J.; Still, Charles H.; Sater, James D.; Bittner, Donald N
2003-01-15
Targets for future laser-fusion ignition experiments will consist of a frozen deuterium-tritium ice layer adhering to the inner surface of a spherical shell, and the specifications for the inner surface quality of this ice layer are extremely demanding. We have developed a numerical raytrace model in order to validate backlit optical shadowgraphy as an ice-surface diagnostic, and we have used the code to simulate shadowgraph data obtained from mathematical ice layers having known modal imperfections. We find that backlit optical shadowgraphy is a valid diagnostic of the mode spectrum of ice-surface imperfections for mode numbers as high as 80 provided the experimental data are analyzed appropriately. We also describe alternative measurement techniques, which may be more sensitive than conventional backlit shadowgraphy.
Geng, Qian; Li, Huan; Li, Yueming
2014-05-01
Experiments were performed to investigate the vibration and acoustic response characteristics of a clamped rectangular aluminum plate in thermal environments. Modal tests were carried out to study the influence of thermal environment on natural vibration. With the increment of structural temperature, natural frequencies of the plate decrease obviously. Mode shape interchange was observed for the modes with frequencies very close to each other. The thermally induced softening effect has unequal influences on the plate along the two in-plane directions. Numerical methods were also employed to study the experimental phenomena. Calculated results indicated that the initial deflection has a great influence on the natural vibration of the heated plate. Even a slight curvature can reduce the thermally induced softening effect obviously. Dynamic response tests were carried out under acoustic and mechanical excitations, and the measured results indicate that the variation in damping determines the response amplitudes at resonant peaks in the test. PMID:24815251
NASA Astrophysics Data System (ADS)
Zhao, Xiao-Feng; Huang, Si-Xun; Du, Hua-Dong
2011-02-01
This paper puts forward possibilities of refractive index profile retrieval using field measurements at an array of radio receivers in terms of variational adjoint approach. The derivation of the adjoint model begins with the parabolic wave equation for a smooth, perfectly conducting surface and horizontal polarization conditions. To deal with the ill-posed difficulties of the inversion, the regularization ideas are introduced into the establishment of the cost function. Based on steepest descent iterations, the optimal value of refractivity could be retrieved quickly at each point over height. Numerical experiments demonstrate that the method works well for low-distance signals, while it is not accurate enough for long-distance propagations. Through curve fitting processing, however, giving a good initial refractivity profile could generally improve the inversions.
Sun, Zi-Qiao; Yang, Yang; Liu, Jing
2012-02-01
This study presented the first in vivo animal experiments of using nano-cryosurgical modality to completely freezing tumor tissues embedded with large blood vessels, which is a tough issue to tackle otherwise. Three-dimensional theoretical simulations were also performed on the complex freezing problems by considering flow and heat transfer of blood flow in large vessels. According to the experimental measurements and numerical predictions, injecting the nanoparticles with high thermal conductivity into the freezing target can significantly reduce the heating effect of blood vessel, shorten the freezing time, and enlarge the freezing range. Most importantly, the introduction of nanoparticles successfully overcomes the classical challenges in completely ablating the tumor region with large blood vessel and enhancing the freezing efficacy of cryosurgery. This investigation consolidates the practical and theoretical foundation for nano-cryosurgery which suggests a highly efficient freezing strategy for treating late stage tumor. PMID:22515090
A Numerical Simulation and Statistical Modeling of High Intensity Radiated Fields Experiment Data
NASA Technical Reports Server (NTRS)
Smith, Laura J.
2004-01-01
Tests are conducted on a quad-redundant fault tolerant flight control computer to establish upset characteristics of an avionics system in an electromagnetic field. A numerical simulation and statistical model are described in this work to analyze the open loop experiment data collected in the reverberation chamber at NASA LaRC as a part of an effort to examine the effects of electromagnetic interference on fly-by-wire aircraft control systems. By comparing thousands of simulation and model outputs, the models that best describe the data are first identified and then a systematic statistical analysis is performed on the data. All of these efforts are combined which culminate in an extrapolation of values that are in turn used to support previous efforts used in evaluating the data.
Numerical simulation of the generation of secondary electrons in the High Current Experiment
Stoltz, P.H.; Furman, M.A.; Vay, J.-L.; Molvik, A.W.; Cohen, R.H.
2003-04-01
Electron effects in the High Current Experiment (HCX) are studied via computer simulation. An approximate expression for the secondary electron yield for a potassium ion striking stainless steel is derived and compared with experimental results. This approximate expression has a peak of roughly 55 electrons at normal incidence at an ion energy of 60 MeV. Using an empirical angular dependence, the secondary electron yield is combined with a numerical simulation of the HCX ion beam dynamics to obtain an estimate for the number of secondary electrons expected per ion-wall collision in the HCX. This estimate is that approximately 150-200 electrons per ion collision may result in the HCX.
NASA Astrophysics Data System (ADS)
Yoshida, Shingo; Kato, Naoyuki; Fukuda, Jun'ichi
2016-04-01
Repeating M5-class earthquakes occurred with a regular recurrence interval on the plate boundary offshore Kamaishi before the 2011 M9 Tohoku-oki earthquake. Since this event, 11 repeating events of M5- M6 have occurred with shorter recurrence intervals than before the M9 event. We performed a forecast experiment on the Kamaishi repeaters after the Tohoku-oki earthquake based on numerical simulations assuming a large patch containing a small patch on a fault model that undergoes continuous afterslip. Simulations were conducted for various parameter values to produce time history of repeated ruptures. To forecast the Ith event, we selected several modeled sequences that are able to accurately reproduce the observed event sequence up to the ( I - 1)th event. The averages of the occurrence times and magnitudes over the selected sequences were used as ensemble forecasts. We attempted to predict the 7th event to the 11th event and discussed the obtained forecast spreads and errors.
Van den Schoor, F; Norman, F; Vandebroek, L; Verplaetsen, F; Berghmans, J
2009-05-30
In this study the auto-ignition limit of ammonia/methane/air mixtures is calculated based upon a perfectly stirred reactor model with convective heat transfer. The results of four different reaction mechanisms are compared with existing experimental data at an initial temperature of 723 K with ammonia concentrations of 0-20 mol.% and methane concentrations of 2.5-10 mol.%. It is found that the calculation of the auto-ignition limit pressure at constant temperature leads to larger relative deviations between calculated and experimental results than the calculation of the auto-ignition temperature at constant pressure. In addition to the calculations, a reaction path analysis is performed to explain the observed lowering of the auto-ignition limit of methane/air mixtures by ammonia addition. It is found that this decrease is caused by the formation of NO and NO(2), which enhance the oxidation of methane at low temperatures. PMID:18926632
Numerical Analysis of a Pulse Detonation Cross Flow Heat Load Experiment
NASA Technical Reports Server (NTRS)
Paxson, Daniel E.; Naples, Andrew .; Hoke, John L.; Schauer, Fred
2011-01-01
A comparison between experimentally measured and numerically simulated, time-averaged, point heat transfer rates in a pulse detonation (PDE) engine is presented. The comparison includes measurements and calculations for heat transfer to a cylinder in crossflow and to the tube wall itself using a novel spool design. Measurements are obtained at several locations and under several operating conditions. The measured and computed results are shown to be in substantial agreement, thereby validating the modeling approach. The model, which is based in computational fluid dynamics (CFD) is then used to interpret the results. A preheating of the incoming fuel charge is predicted, which results in increased volumetric flow and subsequent overfilling. The effect is validated with additional measurements.
Moreira, L.P.; Romao, E.C.; Vieira, L.C.A.; Ferron, G.; Sampaio, A.P.
2005-08-05
A simple bend-draw experimental device is employed to analyze the behavior of narrow strips submitted to a nearly cyclic bending deformation mode followed by a steady state drawing. In this bending-drawing experiment, the strip is firstly bent over a central bead and two lateral beads by applying a controlled holding load and then is pulled out of device throughout the bead radii by a drawing load. The apparatus is mounted in a standard tensile test machine where the holding and drawing loads are recorded with an acquisition data system. The specimen is a rectangular strip cut with 320 mm long and 7 mm wide. The longitudinal (1) and width (w) strip plastic strains are determined from two hardness marks 120 mm spaced whereas the corresponding thickness (t) strain is obtained by volume conservation. Previous experiments showed a correlation between the plastic strain ({epsilon}w/{epsilon}t)BD resulting from the bending-drawing and the Lankford R-values obtained from the uniaxial tensile test. However, previous 3D numerical simulations based upon Hill's quadratic and Ferron's yield criteria revealed a better correlation between the ({epsilon}w/{epsilon}t)BD and the stress ratio {sigma}PS/{sigma}({alpha}), where {sigma}PS stands for the plane-strain tension yield stress and {sigma}({alpha}) for the uniaxial yield stress in uniaxial tension along the drawing direction making an angle {alpha} with the rolling direction. In the present work, the behavior of an IF steel sheet is firstly evaluated by means of uniaxial tensile and drawing-bending experiments conducted at every 15 degrees with respect to the rolling direction. Afterwards, the bending-drawing experiment is investigated with the commercial finite element (FE) code ABAQUS/Standard in an attempt to assess the influence of cyclic loadings upon the bending-drawing strain-ratios.
NASA Astrophysics Data System (ADS)
Elkhoury, J. E.; Detwiler, R. L.; Serajian, V.; Bruno, M. S.
2012-12-01
Geothermal energy resources are more widespread than previously thought and have the potential for providing a significant amount of sustainable clean energy worldwide. In particular, hot permeable sedimentary formations provide many advantages over traditional geothermal recovery and enhanced geothermal systems in low permeability crystalline formations. These include: (1) eliminating the need for hydraulic fracturing, (2) significant reduction in risk for induced seismicity, (3) reducing the need for surface wastewater disposal, (4) contributing to decreases in greenhouse gases, and (5) potential use for CO2 sequestration. Advances in horizontal drilling, completion, and production technology from the oil and gas industry can now be applied to unlock these geothermal resources. Here, we present experimental results from a laboratory scale circulation system and numerical simulations aimed at quantifying the heat transfer capacity of sedimentary rocks. Our experiments consist of fluid flow through a saturated and pressurized sedimentary disc of 23-cm diameter and 3.8-cm thickness heated along its circumference at a constant temperature. Injection and production ports are 7.6-cm apart in the center of the disc. We used DI de-aired water and mineral oil as working fluids and explored temperatures from 20 to 150 oC and flow rates from 2 to 30 ml/min. We performed experiments on sandstone samples (Castlegate and Kirby) with different porosity, permeability and thermal conductivity to evaluate the effect of hydraulic and thermal properties on the heat transfer capacity of sediments. The producing fluid temperature followed an exponential form with time scale transients between 15 and 45 min. Steady state outflow temperatures varied between 60% and 95% of the set boundary temperature, higher percentages were observed for lower temperatures and flow rates. We used the flow and heat transport simulator TOUGH2 to develop a numerical model of our laboratory setting. Given
NASA Astrophysics Data System (ADS)
Jeanniot, L.; Kusznir, N. J.; Manatschal, G.
2012-12-01
During the continental breakup process, the initiation of sea-floor spreading requires both the rupture of the continental crust and the initiation of decompression melting. Using numerical experiments, we investigate how the deformation mode of continental lithosphere thinning and stretching controls the rupture of continental crust and lithospheric mantle, the onset of decompression melting and their relative timing. We use a two dimensional finite element viscous flow model to describe lithosphere and asthenosphere deformation. This flow field is used to advect lithosphere and asthenosphere material and temperature. Decompression melting is predicted using the parameterization scheme of Katz et al. (2003). Consistent with the observations of deformation processes occurring at slow spreading ocean ridges (Cannat, 1996), we assume that the topmost continental and oceanic lithosphere, corresponding to the cooler brittle seismogenic layer, deforms by extensional faulting (which we approximate to pure-shear deformation) and magmatic intrusion. Beneath this topmost lithosphere layer approximately 15-20 km thick, we assume that deformation occurs in response to passive upwelling and thermal and melt buoyancy driven small-scale convection. The relative contribution of these deformation components is parameterised by the ratio Vz/Vx, where Vx is the half spreading rate applied to the topmost lithosphere deformation and Vz is the upwelling velocity associated with the small scale convection. We use a series of numerical experiments to investigate the dependency of continental crust and lithosphere rupture, decompression melt initiation, rifted margin ocean-continent transition architecture and subsidence history on the half-spreading rate Vx, buoyancy driven upwelling rate Vz, the ratio Vz/Vx and upper lithosphere pure-shear width W. Based on the numerical experiment results we explore a polyphase evolution of deformation modes leading to continental breakup, sea
NASA Astrophysics Data System (ADS)
Hedgpeth, A.; Beilman, D.; Crow, S. E.
2014-12-01
Arctic soil organic matter (SOM) mineralization processes are fundamental to the functioning of high latitude soils in relation to nutrients, stability, and feedbacks to atmospheric CO2 and climate. The arctic permafrost zone covers 25% of the northern hemisphere and contains 1672Pg of soil carbon (C). 88% of this C currently resides in frozen soils that are vulnerable to environmental change. For instance, arctic growing seasons may be lengthened, resulting in an increase in plant productivity and rate of below ground labile C inputs as root exudates. Understanding controls on Arctic SOM dynamics requires recognition that labile C inputs have the potential to significantly affect mineralization of previously stable SOM, also known as 'priming effects'. We conducted a substrate addition incubation experiment to quantify and compare respiration in highly organic (42-48 %C) permafrost soils along a north-south transect in western Canada. Near surface soils (10-20 cm) were collected from permafrost peatland sites in the Mackenzie River Basin from 69.2-62.6°N. The surface soils are fairly young (Δ14C values > -140.0) and can be assumed to contain relatively reactive soil carbon. To assess whether addition of labile substrate alters SOM decomposition dynamics, 4.77-11.75 g of permafrost soil were spiked with 0.5 mg D-glucose g-1 soil and incubated at 5°C. A mass balance approach was used to determin substrate-induced respiration and preliminary results suggest a potential for positive priming in these C-rich soils. Baseline respiration rates from the three sites were similar (0.067-0.263 mg CO2 g-1 soil C) yet show some site-specific trends. The rate at which added substrate was utilized within these soils suggests that other factors besides temperature and soil C content are controlling substrate consumption and its effect on SOM decomposition. Microbial activity can be stimulated by substrate addition to such an extent that SOM turnover is enhanced, suggesting that
NASA Astrophysics Data System (ADS)
Zaussinger, F.; Futterer, B.; Egbers, C.
2012-12-01
Thermal convection is one important driving mechanism of flow in the earth mantle. Setting up a self-gravitating buoyancy in a spherical shell geometry is the limiting factor for laboratory experiments to analyze velocity flow structures and heat transport. The geophysical flow model 'GeoFlow II', which is located at the Columbus module on the ISS, realizes such a central gravity. Under microgravity conditions a central dielectrophoretic force field is applied to a fluid filled spherical annulus. In contrast to the first mission 'GeoFlow I' the electro-hydrodynamical volume expansion coefficient of the working fluid has a strong dependence on the temperature and leads to pattern, which are related to a strong temperature dependent viscosity of the fluid. Even though the oil's viscosity itself is temperature-dependent, too, the maximum of viscosity contrast is only up to 1.5. The optical measurement of the fluid flow is based on the Wollaston shearing interferometry, since the on orbit setup avoids the use of measurement particles. This technique leads to fringe patterns. Simulations with RESPECT and GAIAA tend to verify the experimentally observed patterns by different numerical models.
Oki, Delwyn S.
1998-01-01
-calculated freshwater-saltwater interface location for the future recharge and pumping conditions. Model results indicate that an additional 10 million gallons per day (beyond the 1995-allocated rates) of freshwater can potentially be developed from northern Oahu. Various distributions of pumping can be used to obtain the additional 10 million gallons per day of water. The quality of the water pumped will be dependent on site-specific factors and cannot be predicted on the basis of model results. If the additional 10 million gallons per day pumpage is restricted to the Kawailoa and Waialua areas, model results indicate that a regional drawdown (relative to the water-level distribution associated with the 1995-allocated pumping rates) of less than 0.6 foot can be maintained in these two areas. The additional pumping, however, would cause salinity increases in water pumped by existing deep wells. In addition, increases in salinity may occur at other wells in areas where the model indicates no significant problem with upconing.
An experiment in software reliability: Additional analyses using data from automated replications
NASA Technical Reports Server (NTRS)
Dunham, Janet R.; Lauterbach, Linda A.
1988-01-01
A study undertaken to collect software error data of laboratory quality for use in the development of credible methods for predicting the reliability of software used in life-critical applications is summarized. The software error data reported were acquired through automated repetitive run testing of three independent implementations of a launch interceptor condition module of a radar tracking problem. The results are based on 100 test applications to accumulate a sufficient sample size for error rate estimation. The data collected is used to confirm the results of two Boeing studies reported in NASA-CR-165836 Software Reliability: Repetitive Run Experimentation and Modeling, and NASA-CR-172378 Software Reliability: Additional Investigations into Modeling With Replicated Experiments, respectively. That is, the results confirm the log-linear pattern of software error rates and reject the hypothesis of equal error rates per individual fault. This rejection casts doubt on the assumption that the program's failure rate is a constant multiple of the number of residual bugs; an assumption which underlies some of the current models of software reliability. data raises new questions concerning the phenomenon of interacting faults.
Additions and Improvements to the FLASH Code for Simulating High Energy Density Physics Experiments
NASA Astrophysics Data System (ADS)
Lamb, D. Q.; Daley, C.; Dubey, A.; Fatenejad, M.; Flocke, N.; Graziani, C.; Lee, D.; Tzeferacos, P.; Weide, K.
2015-11-01
FLASH is an open source, finite-volume Eulerian, spatially adaptive radiation hydrodynamics and magnetohydrodynamics code that incorporates capabilities for a broad range of physical processes, performs well on a wide range of computer architectures, and has a broad user base. Extensive capabilities have been added to FLASH to make it an open toolset for the academic high energy density physics (HEDP) community. We summarize these capabilities, with particular emphasis on recent additions and improvements. These include advancements in the optical ray tracing laser package, with methods such as bi-cubic 2D and tri-cubic 3D interpolation of electron number density, adaptive stepping and 2nd-, 3rd-, and 4th-order Runge-Kutta integration methods. Moreover, we showcase the simulated magnetic field diagnostic capabilities of the code, including induction coils, Faraday rotation, and proton radiography. We also describe several collaborations with the National Laboratories and the academic community in which FLASH has been used to simulate HEDP experiments. This work was supported in part at the University of Chicago by the DOE NNSA ASC through the Argonne Institute for Computing in Science under field work proposal 57789; and the NSF under grant PHY-0903997.
NASA Astrophysics Data System (ADS)
Zhong, L.; Oostrom, M.; Wietsma, T. W.; Covert, M. A.
2008-10-01
Low-permeability zones are typically bypassed when remedial fluids are injected into subsurface heterogeneous aquifer systems. Therefore, contaminants in the bypassed areas may not be contacted by the amendments in the remedial fluid, which may significantly prolong remediation operations. Laboratory experiments and numerical studies have been conducted to investigate the use of a shear-thinning polymer (Xanthan gum) to improve access to low-permeability zones in heterogeneous systems. The chemicals sodium mono-phosphate and the surfactant MA-80 were used as the remedial amendments. The impact of polymer concentration, fluid injection rate, and permeability contrast in the heterogeneous systems has been studied in a series of eleven two-dimensional flow cell experiments. The Subsurface Transport over Multiple Phases (STOMP) simulator was modified to include polymer-induced shear-thinning effects. The experimental and simulation results clearly show that using the polymer leads to an enhanced delivery of remedial amendments to lower-permeability zones and an increased sweeping efficiency. An added benefit of using the polymer is the stabilization of the displacing front when density differences exist between displaced and displacing fluids. The modified STOMP simulator was able to predict the experimental observed fluid displacing behavior well and might be used to predict subsurface remediation performance when a shear-thinning fluid is used to remediate a heterogeneous system at larger scales.
Zhong, L; Oostrom, M; Wietsma, T W; Covert, M A
2008-10-23
Low-permeability zones are typically bypassed when remedial fluids are injected into subsurface heterogeneous aquifer systems. Therefore, contaminants in the bypassed areas may not be contacted by the amendments in the remedial fluid, which may significantly prolong remediation operations. Laboratory experiments and numerical studies have been conducted to investigate the use of a shear-thinning polymer (Xanthan gum) to improve access to low-permeability zones in heterogeneous systems. The chemicals sodium mono-phosphate and the surfactant MA-80 were used as the remedial amendments. The impact of polymer concentration, fluid injection rate, and permeability contrast in the heterogeneous systems has been studied in a series of eleven two-dimensional flow cell experiments. The Subsurface Transport over Multiple Phases (STOMP) simulator was modified to include polymer-induced shear-thinning effects. The experimental and simulation results clearly show that using the polymer leads to an enhanced delivery of remedial amendments to lower-permeability zones and an increased sweeping efficiency. An added benefit of using the polymer is the stabilization of the displacing front when density differences exist between displaced and displacing fluids. The modified STOMP simulator was able to predict the experimental observed fluid displacing behavior well and might be used to predict subsurface remediation performance when a shear-thinning fluid is used to remediate a heterogeneous system at larger scales. PMID:18786743
NASA Astrophysics Data System (ADS)
Tzeferacos, Petros; Fatenejad, Milad; Flocke, Norbert; Graziani, Carlo; Gregori, Gianluca; Lamb, Donald; Lee, Dongwook; Meinecke, Jena; Scopatz, Anthony; Weide, Klaus
2014-10-01
In this study we present high-resolution numerical simulations of laboratory experiments that study the turbulent amplification of magnetic fields generated by laser-driven colliding jets. The radiative magneto-hydrodynamic (MHD) simulations discussed here were performed with the FLASH code and have assisted in the analysis of the experimental results obtained from the Vulcan laser facility. In these experiments, a pair of thin Carbon foils is placed in an Argon-filled chamber and is illuminated to create counter-propagating jets. The jets carry magnetic fields generated by the Biermann battery mechanism and collide to form a highly turbulent region. The interaction is probed using a wealth of diagnostics, including induction coils that are capable of providing the field strength and directionality at a specific point in space. The latter have revealed a significant increase in the field's strength due to turbulent amplification. Our FLASH simulations have allowed us to reproduce the experimental findings and to disentangle the complex processes and dynamics involved in the colliding flows. This work was supported in part at the University of Chicago by DOE NNSA ASC.
Numerical Calculations of 3-D High-Lift Flows and Comparison with Experiment
NASA Technical Reports Server (NTRS)
Compton, William B, III
2015-01-01
Solutions were obtained with the Navier-Stokes CFD code TLNS3D to predict the flow about the NASA Trapezoidal Wing, a high-lift wing composed of three elements: the main-wing element, a deployed leading-edge slat, and a deployed trailing-edge flap. Turbulence was modeled by the Spalart-Allmaras one-equation turbulence model. One case with massive separation was repeated using Menter's two-equation SST (Menter's Shear Stress Transport) k-omega turbulence model in an attempt to improve the agreement with experiment. The investigation was conducted at a free stream Mach number of 0.2, and at angles of attack ranging from 10.004 degrees to 34.858 degrees. The Reynolds number based on the mean aerodynamic chord of the wing was 4.3 x 10 (sup 6). Compared to experiment, the numerical procedure predicted the surface pressures very well at angles of attack in the linear range of the lift. However, computed maximum lift was 5% low. Drag was mainly under predicted. The procedure correctly predicted several well-known trends and features of high-lift flows, such as off-body separation. The two turbulence models yielded significantly different solutions for the repeated case.
Numerical simulation of mud erosion rate in sand-mud alternate layer and comparison with experiment
NASA Astrophysics Data System (ADS)
Yoshida, T.; Yamaguchi, T.; Oyama, H.; Sato, T.
2015-12-01
For gas production from methane hydrates in sand-mud alternate layers, depressurization method is expected as feasible. After methane hydrate is dissociated, gas and water flow in pore space. There is a concern about the erosion of mud surface and it may result in flow blockage that disturbs the gas production. As a part of a Japanese National hydrate research program (MH21, funded by METI), we developed a numerical simulation of water-induced mud erosion in pore-scale sand-mud domains to model such mud erosion. The size of which is of the order of 100 micro meter. Water flow is simulated using a lattice Boltzmann method (LBM) and mud surface is treated as solid boundary with arbitrary shape, which changes with time. Periodic boundary condition is adopted at the domain boundaries, except for the surface of mud layers and the upper side. Shear stress acting on the mud surface is calculated using a momentum-exchange method. Mud layer is eroded when the shear stress exceeds a threshold coined a critical shear stress. In this study, we compared the simulated mud erosion rate with experimental data acquired from an experiment using artificial sand-mud core. As a result, the simulated erosion rate agrees well with that of the experiment.
Laser beam propagation through bulk nonlinear media: Numerical simulation and experiment
NASA Astrophysics Data System (ADS)
Kovsh, Dmitriy I.
This dissertation describes our efforts in modeling the propagation of high intensity laser pulses through optical systems consisting of one or multiple nonlinear elements. These nonlinear elements can be up to 103 times thicker than the depth of focus of the laser beam, so that the beam size changes drastically within the medium. The set of computer codes developed are organized in a software package (NLO_BPM). The ultrafast nonlinearities of the bound-electronic n2 and two-photon absorption as well as time dependent excited-state, free-carrier and thermal nonlinearities are included in the codes for modeling propagation of picosecond to nanosecond pulses and pulse trains. Various cylindrically symmetric spatial distributions of the input beam are modeled. We use the cylindrical symmetry typical of laser outputs to reduce the CPU and memory requirements making modeling a real- time task on PC's. The hydrodynamic equations describing the rarefaction of the medium due to heating and electrostriction are solved in the transient regime to determine refractive index changes on a nanosecond time scale. This effect can be simplified in some cases by an approximation that assumes an instantaneous expansion. We also find that the index change obtained from the photo-acoustic equation overshoots its steady-state value once the ratio between the pulse width and the acoustic transit time is greater than unity. We numerically study the sensitivity of the closed- aperture Z-scan experiment to nonlinear refraction for various input beam profiles. If the beam has a ring structure with a minimum (or zero) on axis in the far field, the sensitivity of Z-scan measurements can be increased by up to one order of magnitude. The linear propagation module integrated with the nonlinear beam propagation codes allows the simulation of typical experiments such as Z-scan and optical limiting experiments. We have used these codes to model the performance of optical limiters. We study two of the
NASA Astrophysics Data System (ADS)
Bilardello, D.
2014-12-01
Understanding depositional remanent magnetizations (DRMs) bears implications on interpreting paleomagnetic and paleointensity records extracted from sedimentary rocks. Laboratory deposition experiments have yielded DRMs with shallow remanent inclinations and revealed a field dependence of the magnetization (M), which is orders of magnitude lower than the saturation remanence. To investigate these observations further, experiments involving differently shaped particles were performed. Spherical particles confirmed the field dependence of both the inclination error and M and the fact that the DRM acquired experimentally is lower than saturation. A sediment concentration dependence of the inclination error was observed, indicating a dependance of the inclination error on the sediment load/burial depth or the sedimentation rate. Other outcome was the certainty that spherical particles alone can lead to substantial inclination shallowing. Numerical simulations of settling spherical particles indicated that DRM should be ~10 times lower than the saturation remanence and predicted that rolling of the grains on the sediment surface and particle interactions during settling can produce a substantial shallowing of the inclination and lowering of the remanence, bringing the simulations in close agreement to the experimental results. Experiments involving platy particles, instead allowed interesting comparisons and gave insight into the behavior of differently shaped particles, for instance yielding smaller amounts of shallowing than spheres, in contrast to general belief. Viewing DRM as an anisotropic process allows fitting the experimental results with tensors (kDRM). The ratios of kvertical over khorizontal are in good agreement to the ratios of M obtained in vertical over horizontal experimental fields, which should be equivalent to the widely used inclination shallowing factor f. Experimental results were highly repeatabile, however not always as repeatable for both M and
Single Droplet Combustion of Decane in Microgravity: Experiments and Numerical Modeling
NASA Technical Reports Server (NTRS)
Dietrich, D. L.; Struk, P. M.; Ikegam, M.; Xu, G.
2004-01-01
This paper presents experimental data on single droplet combustion of decane in microgravity and compares the results to a numerical model. The primary independent experiment variables are the ambient pressure and oxygen mole fraction, pressure, droplet size (over a relatively small range) and ignition energy. The droplet history (D(sup 2) history) is non-linear with the burning rate constant increasing throughout the test. The average burning rate constant, consistent with classical theory, increased with increasing ambient oxygen mole fraction and was nearly independent of pressure, initial droplet size and ignition energy. The flame typically increased in size initially, and then decreased in size, in response to the shrinking droplet. The flame standoff increased linearly for the majority of the droplet lifetime. The flame surrounding the droplet extinguished at a finite droplet size at lower ambient pressures and an oxygen mole fraction of 0.15. The extinction droplet size increased with decreasing pressure. The model is transient and assumes spherical symmetry, constant thermo-physical properties (specific heat, thermal conductivity and species Lewis number) and single step chemistry. The model includes gas-phase radiative loss and a spherically symmetric, transient liquid phase. The model accurately predicts the droplet and flame histories of the experiments. Good agreement requires that the ignition in the experiment be reasonably approximated in the model and that the model accurately predict the pre-ignition vaporization of the droplet. The model does not accurately predict the dependence of extinction droplet diameter on pressure, a result of the simplified chemistry in the model. The transient flame behavior suggests the potential importance of fuel vapor accumulation. The model results, however, show that the fractional mass consumption rate of fuel in the flame relative to fuel vaporized is close to 1.0 for all but the lowest ambient oxygen mole
Characterization of stony soils' hydraulic conductivity using laboratory and numerical experiments
NASA Astrophysics Data System (ADS)
Beckers, Eléonore; Pichault, Mathieu; Pansak, Wanwisa; Degré, Aurore; Garré, Sarah
2016-08-01
Determining soil hydraulic properties is of major concern in various fields of study. Although stony soils are widespread across the globe, most studies deal with gravel-free soils, so that the literature describing the impact of stones on the hydraulic conductivity of a soil is still rather scarce. Most frequently, models characterizing the saturated hydraulic conductivity of stony soils assume that the only effect of rock fragments is to reduce the volume available for water flow, and therefore they predict a decrease in hydraulic conductivity with an increasing stoniness. The objective of this study is to assess the effect of rock fragments on the saturated and unsaturated hydraulic conductivity. This was done by means of laboratory experiments and numerical simulations involving different amounts and types of coarse fragments. We compared our results with values predicted by the aforementioned predictive models. Our study suggests that it might be ill-founded to consider that stones only reduce the volume available for water flow. We pointed out several factors of the saturated hydraulic conductivity of stony soils that are not considered by these models. On the one hand, the shape and the size of inclusions may substantially affect the hydraulic conductivity. On the other hand, laboratory experiments show that an increasing stone content can counteract and even overcome the effect of a reduced volume in some cases: we observed an increase in saturated hydraulic conductivity with volume of inclusions. These differences are mainly important near to saturation. However, comparison of results from predictive models and our experiments in unsaturated conditions shows that models and data agree on a decrease in hydraulic conductivity with stone content, even though the experimental conditions did not allow testing for stone contents higher than 20 %.
NASA Technical Reports Server (NTRS)
1974-01-01
The design and rationale of an advanced labeled release experiment based on single addition of soil and multiple sequential additions of media into each of four test chambers are outlined. The feasibility for multiple addition tests was established and various details of the methodology were studied. The four chamber battery of tests include: (1) determination of the effect of various atmospheric gases and selection of that gas which produces an optimum response; (2) determination of the effect of incubation temperature and selection of the optimum temperature for performing Martian biochemical tests; (3) sterile soil is dosed with a battery of C-14 labeled substrates and subjected to experimental temperature range; and (4) determination of the possible inhibitory effects of water on Martian organisms is performed initially by dosing with 0.01 ml and 0.5 ml of medium, respectively. A series of specifically labeled substrates are then added to obtain patterns in metabolic 14CO2 (C-14)O2 evolution.
Additive Manufacturing, Design, Testing, and Fabrication: A Full Engineering Experience at JSC
NASA Technical Reports Server (NTRS)
Zusack, Steven
2016-01-01
I worked on several projects this term. While most projects involved additive manufacturing, I was also involved with two design projects, two testing projects, and a fabrication project. The primary mentor for these was Richard Hagen. Secondary mentors were Hai Nguyen, Khadijah Shariff, and fabrication training from James Brown. Overall, my experience at JSC has been successful and what I have learned will continue to help me in my engineering education and profession long after I leave. My 3D printing projects ranged from less than a 1 cubic centimeter to about 1 cubic foot and involved several printers using different printing technologies. It was exciting to become familiar with printing technologies such as industrial grade FDM (Fused Deposition Modeling), the relatively new SLA (Stereolithography), and PolyJet. My primary duty with the FDM printers was to model parts that came in from various sources to print effectively and efficiently. Using methods my mentor taught me and the Stratasys Insight software, I was able to minimize imperfections, hasten build time, improve strength for specific forces (tensile, shear, etc...), and reduce likelihood of a print-failure. Also using FDM, I learned how to repair a part after it was printed. This is done by using a special kind of glue that chemically melts the two faces of plastic parts together to form a fused interface. My first goal with SLA technology was to bring the printer back to operational readiness. In becoming familiar with the Pegasus SLA printer, I researched the leveling, laser settings, and different vats to hold liquid material. With this research, I was successfully able to bring the Pegasus back online and have successfully printed multiple sample parts as well as functional parts. My experience with PolyJet technology has been focused on an understanding of the abilities/limits, costs, and the maintenance for daily use. Still upcoming will be experience with using a composite printer that uses FDM
Nitrate removal in stream ecosystems measured by 15N addition experiments: Total uptake
Mulholland, Patrick J; Hall, Robert; Tank, Jennifer; Sobota, Daniel; O'Brien, Jon; Webster, Jackson; Valett, H. Maurice; Dodds, Walter; Poole, Geoff; Peterson, Chris G.; Meyer, Judy; McDowell, William; Johnson, Sherri; Hamilton, Stephen; Gregory, Stanley; Grimm, Nancy; Dahm, Cliff; Cooper, Lee W; Ashkenas, Linda; Thomas, Suzanne; Sheibley, Rich; Potter, Jody; Niederlehner, Bobbie; Johnson, Laura; Helton, Ashley; Crenshaw, Chelsea; Burgin, Amy; Bernot, Melody; Beaulieu, Jake; Arango, Clay
2009-01-01
We measured uptake length of {sup 15}NO{sub 3}{sup -} in 72 streams in eight regions across the United States and Puerto Rico to develop quantitative predictive models on controls of NO{sub 3}{sup -} uptake length. As part of the Lotic Intersite Nitrogen Experiment II project, we chose nine streams in each region corresponding to natural (reference), suburban-urban, and agricultural land uses. Study streams spanned a range of human land use to maximize variation in NO{sub 3}{sup -} concentration, geomorphology, and metabolism. We tested a causal model predicting controls on NO{sub 3}{sup -} uptake length using structural equation modeling. The model included concomitant measurements of ecosystem metabolism, hydraulic parameters, and nitrogen concentration. We compared this structural equation model to multiple regression models which included additional biotic, catchment, and riparian variables. The structural equation model explained 79% of the variation in log uptake length (S{sub Wtot}). Uptake length increased with specific discharge (Q/w) and increasing NO{sub 3}{sup -} concentrations, showing a loss in removal efficiency in streams with high NO{sub 3}{sup -} concentration. Uptake lengths shortened with increasing gross primary production, suggesting autotrophic assimilation dominated NO{sub 3}{sup -} removal. The fraction of catchment area as agriculture and suburban-urban land use weakly predicted NO{sub 3}{sup -} uptake in bivariate regression, and did improve prediction in a set of multiple regression models. Adding land use to the structural equation model showed that land use indirectly affected NO{sub 3}{sup -} uptake lengths via directly increasing both gross primary production and NO{sub 3}{sup -} concentration. Gross primary production shortened S{sub Wtot}, while increasing NO{sub 3}{sup -} lengthened S{sub Wtot} resulting in no net effect of land use on NO{sub 3}{sup -} removal.
Nitrate removal in stream ecosystems measured by 15N addition experiments: Total uptake
Hall, R.O., Jr.; Tank, J.L.; Sobota, D.J.; Mulholland, P.J.; O'Brien, J. M.; Dodds, W.K.; Webster, J.R.; Valett, H.M.; Poole, G.C.; Peterson, B.J.; Meyer, J.L.; McDowell, W.H.; Johnson, S.L.; Hamilton, S.K.; Grimm, N. B.; Gregory, S.V.; Dahm, Clifford N.; Cooper, L.W.; Ashkenas, L.R.; Thomas, S.M.; Sheibley, R.W.; Potter, J.D.; Niederlehner, B.R.; Johnson, L.T.; Helton, A.M.; Crenshaw, C.M.; Burgin, A.J.; Bernot, M.J.; Beaulieu, J.J.; Arangob, C.P.
2009-01-01
We measured uptake length of 15NO-3 in 72 streams in eight regions across the United States and Puerto Rico to develop quantitative predictive models on controls of NO-3 uptake length. As part of the Lotic Intersite Nitrogen eXperiment II project, we chose nine streams in each region corresponding to natural (reference), suburban-urban, and agricultural land uses. Study streams spanned a range of human land use to maximize variation in NO-3 concentration, geomorphology, and metabolism. We tested a causal model predicting controls on NO-3 uptake length using structural equation modeling. The model included concomitant measurements of ecosystem metabolism, hydraulic parameters, and nitrogen concentration. We compared this structural equation model to multiple regression models which included additional biotic, catchment, and riparian variables. The structural equation model explained 79% of the variation in log uptake length (S Wtot). Uptake length increased with specific discharge (Q/w) and increasing NO-3 concentrations, showing a loss in removal efficiency in streams with high NO-3 concentration. Uptake lengths shortened with increasing gross primary production, suggesting autotrophic assimilation dominated NO-3 removal. The fraction of catchment area as agriculture and suburban-urban land use weakly predicted NO-3 uptake in bivariate regression, and did improve prediction in a set of multiple regression models. Adding land use to the structural equation model showed that land use indirectly affected NO-3 uptake lengths via directly increasing both gross primary production and NO-3 concentration. Gross primary production shortened SWtot, while increasing NO-3 lengthened SWtot resulting in no net effect of land use on NO- 3 removal. ?? 2009.
Numerical simulations of impacts involving porous bodies. II. Comparison with laboratory experiments
NASA Astrophysics Data System (ADS)
Jutzi, Martin; Michel, Patrick; Hiraoka, Kensuke; Nakamura, Akiko M.; Benz, Willy
2009-06-01
In this paper, we compare the outcome of high-velocity impact experiments on porous targets, composed of pumice, with the results of simulations by a 3D SPH hydrocode in which a porosity model has been implemented. The different populations of small bodies of our Solar System are believed to be composed, at least partially, of objects with a high degree of porosity. To describe the fragmentation of such porous objects, a different model is needed than that used for non-porous bodies. In the case of porous bodies, the impact process is not only driven by the presence of cracks which propagate when a stress threshold is reached, it is also influenced by the crushing of pores and compaction. Such processes can greatly affect the whole body's response to an impact. Therefore, another physical model is necessary to improve our understanding of the collisional process involving porous bodies. Such a model has been developed recently and introduced successfully in a 3D SPH hydrocode [Jutzi, M., Benz, W., Michel, P., 2008. Icarus 198, 242-255]. Basic tests have been performed which already showed that it is implemented in a consistent way and that theoretical solutions are well reproduced. However, its full validation requires that it is also capable of reproducing the results of real laboratory impact experiments. Here we present simulations of laboratory experiments on pumice targets for which several of the main material properties have been measured. We show that using the measured material properties and keeping the remaining free parameters fixed, our numerical model is able to reproduce the outcome of these experiments carried out under different impact conditions. This first complete validation of our model, which will be tested for other porous materials in the future, allows us to start addressing problems at larger scale related to small bodies of our Solar System, such as collisions in the Kuiper Belt or the formation of a family by the disruption of a porous
Zhong, Lirong; Oostrom, Martinus; Wietsma, Thomas W.; Covert, Matthew A.
2008-07-29
Abstract Heterogeneity is often encountered in subsurface contamination characterization and remediation. Low-permeability zones are typically bypassed when remedial fluids are injected into subsurface heterogeneous aquifer systems. Therefore, contaminants in the bypassed areas may not be contacted by the amendments in the remedial fluid, which may significantly prolong the remediation operations. Laboratory experiments and numerical studies have been conducted to develop the Mobility-Controlled Flood (MCF) technology for subsurface remediation and to demonstrate the capability of this technology in enhancing the remedial amendments delivery to the lower permeability zones in heterogeneous systems. Xanthan gum, a bio-polymer, was used to modify the viscosity of the amendment-containing remedial solutions. Sodium mono-phosphate and surfactant were the remedial amendment used in this work. The enhanced delivery of the amendments was demonstrated in two-dimensional (2-D) flow cell experiments, packed with heterogeneous systems. The impact of polymer concentration, fluid injection rate, and permeability contract in the heterogeneous systems has been studied. The Subsurface Transport over Multiple Phases (STOMP) simulator was modified to include polymer-induced shear thinning effects. Shear rates of polymer solutions were computed from pore-water velocities using a relationship proposed in the literature. Viscosity data were subsequently obtained from empirical viscosity-shear rate relationships derived from laboratory data. The experimental and simulation results clearly show that the MCF technology is capable of enhancing the delivery of remedial amendments to subsurface lower permeability zones. The enhanced delivery significantly improved the NAPL removal from these zones and the sweeping efficiency on a heterogeneous system was remarkably increased when a polymer fluid was applied. MCF technology is also able to stabilize the fluid displacing front when there is a
Analysis of the Source Physics Experiment SPE4 Prime Using State-Of Parallel Numerical Tools.
NASA Astrophysics Data System (ADS)
Vorobiev, O.; Ezzedine, S. M.; Antoun, T.; Glenn, L.
2015-12-01
This work describes a methodology used for large scale modeling of wave propagation from underground chemical explosions conducted at the Nevada National Security Site (NNSS) fractured granitic rock. We show that the discrete natures of rock masses as well as the spatial variability of the fabric of rock properties are very important to understand ground motions induced by underground explosions. In order to build a credible conceptual model of the subsurface we integrated the geological, geomechanical and geophysical characterizations conducted during recent test at the NNSS as well as historical data from the characterization during the underground nuclear test conducted at the NNSS. Because detailed site characterization is limited, expensive and, in some instances, impossible we have numerically investigated the effects of the characterization gaps on the overall response of the system. We performed several computational studies to identify the key important geologic features specific to fractured media mainly the joints characterized at the NNSS. We have also explored common key features to both geological environments such as saturation and topography and assess which characteristics affect the most the ground motion in the near-field and in the far-field. Stochastic representation of these features based on the field characterizations has been implemented into LLNL's Geodyn-L hydrocode. Simulations were used to guide site characterization efforts in order to provide the essential data to the modeling community. We validate our computational results by comparing the measured and computed ground motion at various ranges for the recently executed SPE4 prime experiment. We have also conducted a comparative study between SPE4 prime and previous experiments SPE1 and SPE3 to assess similarities and differences and draw conclusions on designing SPE5.
NUMERICAL EXPERIMENTS ON THE TWO-STEP EMERGENCE OF TWISTED MAGNETIC FLUX TUBES IN THE SUN
Toriumi, S.; Yokoyama, T.
2011-07-10
We present the new results of the two-dimensional numerical experiments on the cross-sectional evolution of a twisted magnetic flux tube rising from the deeper solar convection zone (-20,000 km) to the corona through the surface. The initial depth is 10 times deeper than most of the previous calculations focusing on the flux emergence from the uppermost convection zone. We find that the evolution is illustrated by the following two-step process. The initial tube rises due to its buoyancy, subject to aerodynamic drag due to the external flow. Because of the azimuthal component of the magnetic field, the tube maintains its coherency and does not deform to become a vortex roll pair. When the flux tube approaches the photosphere and expands sufficiently, the plasma on the rising tube accumulates to suppress the tube's emergence. Therefore, the flux decelerates and extends horizontally beneath the surface. This new finding owes to our large-scale simulation, which simultaneously calculates the dynamics within the interior as well as above the surface. As the magnetic pressure gradient increases around the surface, magnetic buoyancy instability is triggered locally and, as a result, the flux rises further into the solar corona. We also find that the deceleration occurs at a higher altitude than assumed in our previous experiment using magnetic flux sheets. By conducting parametric studies, we investigate the conditions for the two-step emergence of the rising flux tube: field strength {approx}> 1.5 x 10{sup 4} G and the twist {approx}> 5.0 x 10{sup -4} km{sup -1} at -20,000 km depth.
NASA Astrophysics Data System (ADS)
Jankovic, I.
2002-05-01
Flow and transport in porous formations are analyzed using numerical simulations. Hydraulic conductivity is treated as a spatial random function characterized by a probability density function and a two-point covariance function. Simulations are performed for a multi-indicator conductivity structure developed by Gedeon Dagan (personal communication). This conductivity structure contains inhomogeneities (inclusions) of elliptical and ellipsoidal geometry that are embedded in a homogeneous background. By varying the distribution of sizes and conductivities of inclusions, any probability density function and two-point covariance may be reproduced. The multi-indicator structure is selected since it yields simple approximate transport solutions (Aldo Fiori, personal communication) and accurate numerical solutions (based on the Analytic Element Method). The dispersion is examined for two conceptual models. Both models are based on the multi-indicator conductivity structure. The first model is designed to examine dispersion in aquifers with continuously varying conductivity. The inclusions in this model cover as much area/volume of the porous formation as possible. The second model is designed for aquifers that contain clay/sand/gravel lenses embedded in otherwise homogeneous background. The dispersion in both aquifer types is simulated numerically. Simulation results are compared to those obtained using simple approximate solutions. In order to infer transport statistics that are representative of an infinite domain using the numerical experiments, the inclusions are placed in a domain that was shaped as a large ellipse (2D) and a large spheroid (3D) that were submerged in an unbounded homogeneous medium. On a large scale, the large body of inclusions behaves like a single large inhomogeneity. The analytic solution for a uniform flow past the single inhomogeneity of such geometry yields uniform velocity inside the domain. The velocity differs from that at infinity and
Early electromagnetic waves from earthquake rupturing: II. validation and numerical experiments
NASA Astrophysics Data System (ADS)
Gao, Yongxin; Chen, Xiaofei; Hu, Hengshan; Zhang, Jie
2013-03-01
We validate the branch-cut integration (BCI) technique presented in the companion paper. The numerical result shows that the early electromagnetic (EM) signal calculated by the BCI method is in good agreement with that in the full waveform calculated by the real-axis integration method. We further find that to calculate the early EM signal only the integrals along the vertical branch cuts that are around the k0 and kem branch points are needed, whereas neither the integrals along the vertical branch cuts around the Pf(P), S and Ps branch points nor the residues of the poles are necessary. We conduct numerical experiments to analyse the early EM signal generated by an earthquake in a porous half-space, including its component analysis, sensitivities to the conductivity, viscosity and recording depth, and radiation pattern. The component analysis shows that the total early EM (emTotal) signal is not a single wave but a combination of three kinds of EM waves, namely, the direct emd wave radiated from the source, the reflected emr waves converted from the emd wave, the direct P and S waves at the free surface and the critically refracted EM0 waves which are also converted from the emd wave, the direct P and S waves at the free surface. Three pulses, namely, the emd-, P- and S-converted pulses are identified in the emTotal signal according to their different arrival times. The emd-converted pulse arrives immediately after the occurrence of the earthquake and it is a combination of the emd wave, the emr and EM0 waves converted from the emd wave at the free surface. The P-converted (S-converted, respectively) pulse owns an arrival time approximately equal to that spent by the P wave (S wave, respectively) travelling from the hypocentre to the epicentre, and it is a combination of the emr and EM0 waves converted from the P wave (S wave, respectively) at the free surface. The P-converted pulse is usually weaker than the S- and emd-converted pulses in the electrogram and is
NASA Astrophysics Data System (ADS)
Hilbe, Michael; Anselmetti, Flavio S.
2015-02-01
Previous studies of the sediments of Lake Lucerne have shown that massive subaqueous mass movements affecting unconsolidated sediments on lateral slopes are a common process in this lake, and, in view of historical reports describing damaging waves on the lake, it was suggested that tsunamis generated by mass movements represent a considerable natural hazard on the lakeshores. Newly performed numerical simulations combining two-dimensional, depth-averaged models for mass-movement propagation and for tsunami generation, propagation and inundation reproduce a number of reported tsunami effects. Four analysed mass-movement scenarios—three based on documented slope failures involving volumes of 5.5 to 20.8 × 106 m3—show peak wave heights of several metres and maximum runup of 6 to >10 m in the directly affected basins, while effects in neighbouring basins are less drastic. The tsunamis cause large-scale inundation over distances of several hundred metres on flat alluvial plains close to the mass-movement source areas. Basins at the ends of the lake experience regular water-level oscillations with characteristic periods of several minutes. The vulnerability of potentially affected areas has increased dramatically since the times of the damaging historical events, recommending a thorough evaluation of the hazard.
Dynamics of Soil Water Evaporation during Soil Drying: Laboratory Experiment and Numerical Analysis
Han, Jiangbo; Zhou, Zhifang
2013-01-01
Laboratory and numerical experiments were conducted to investigate the evolution of soil water evaporation during a continuous drying event. Simulated soil water contents and temperatures by the calibrated model well reproduced measured values at different depths. Results show that the evaporative drying process could be divided into three stages, beginning with a relatively high evaporation rate during stage 1, followed by a lower rate during transient stage and stage 2, and finally maintaining a very low and constant rate during stage 3. The condensation zone was located immediately below the evaporation zone in the profile. Both peaks of evaporation and condensation rate increased rapidly during stage 1 and transition stage, decreased during stage 2, and maintained constant during stage 3. The width of evaporation zone kept a continuous increase during stages 1 and 2 and maintained a nearly constant value of 0.68 cm during stage 3. When the evaporation zone totally moved into the subsurface, a dry surface layer (DSL) formed above the evaporation zone at the end of stage 2. The width of DSL also presented a continuous increase during stage 2 and kept a constant value of 0.71 cm during stage 3. PMID:24489492
NASA Astrophysics Data System (ADS)
Chern, J.; Tao, W.; Shen, B.
2011-12-01
The Madden-Julian oscillation (MJO) is the dominant component of intraseasonal variability in the tropic. It interacts and influences a wide range of weather and climate phenomena across different temporal and spatial scales. Despite the important role the MJO plays in the weather and climate system, past multi-model MJO intercomparison studies have shown that current global general circulation models (GCMs) still have considerable shortcomings in representing and forecasting this phenomenon. To improve representation of MJO and tropical convective cloud systems in global model, an Multiscale Modeling Framework (MMF) in which a cloud-resolving model takes the place of the sing-column cumulus parameterization used in convectional GCMs has been successfully developed at NAAS Goddard (Tao et al. 2009). To evaluate and improve the ability of this modeling system in representation and prediction of the MJO, several numerical hindcast experiments of a few selected MJO events during YOTC have been carried out. The ability of the model to simulate the MJO events is examined using diagnostic and skill metrics developed by the CLIVAR MJO Working Group Project as well as comparisons with a high-resolution global mesoscale model simulations, satellite observations, and analysis dataset. Several key variables associated with the MJO are investigated, including precipitation, outgoing longwave radiation, large-scale circulation, surface latent heat flux, low-level moisture convergence, vertical structure of moisture and hydrometers, and vertical diabatic heating profiles to gain insight of cloud processes associated with the MJO events.
NASA Astrophysics Data System (ADS)
Tirel, C.; Gautier, P.; van Hinsbergen, D.; Wortel, R.
2007-12-01
The Cycladic extensional province (Greece) contains classical examples of metamorphic core complexes (MCCs), where exhumation was accommodated along multiple interfering and/or sequentially developed syn- and antithetic extensional detachment zones. Previous studies on the development of MCCs did not take into account the possible interference between multiple and closely spaced MCCs. In the present study, we have performed new lithosphere-scale experiments in which the deformation is not a priori localized so as to explore the conditions of the development of several MCCs in a direction parallel to extension. In a narrow range of conditions, MCCs are closely spaced, interfere with each other, and develop in sequence. From a comparison between numerical results and geological observations, we find that the Cyclades metamorphic core complexes are in good agreement with the model in terms of Moho geometry and depth, kinematic and structural history, timing and duration of core complex formation and metamorphic history. We infer that, for Cycladic MCC-type to develop, an initial crustal thickness prior to the onset of post-orogenic extension between 40 and 44 km, a boundary velocity close to 2 cm/yr and an initial thermal lithospheric thickness of about 60 km are required. The latter may be explained by a significant heating due to delamination of subducting continental crust or vigorous small-scale thermal convection.
Numerical experiments of magnetic reconnection in the solar flare and CME current sheet
NASA Astrophysics Data System (ADS)
Mei, Zhixing; Lin, Jun; Shen, Chengcai
2012-07-01
Magnetic reconnection plays a critical role in the energy conversion in the solar eruption. This paper performs a set of MHD experiments for the magnetic reconnection process in a current sheet formed in a disrupting magnetic configuration. The eruption results from the loss of equilibrium in the magnetic configuration that includes a current-carrying flux rope, which is used to model the filament floating in the corona. In order to study the fine structure and micro process inside the current sheet (CS), the mesh refinement technology is used to depress the numerical diffusion. A uniform physical diffusion is applied and results in a Lundquist number S=10^4 in the vicinity of CS. Because of the advantage of the foregoing setting, some features appear with high resolution, including plasmoids due to the tearing mode and the plasmoid instabilities, turbulence regions, and the slow mode shocks. Inside CS, magnetic reconnection goes through the Sweet-Parker and the fractal fashions, and eventually, it displays a time-dependent Petschek pattern. Our results seem to support the concept of fractal reconnection suggested by Shibata et al. (1995) and Shibata & Tanuma (2001). And our results suggest that the CS evolves through a Sweet-Parker reconnection prior to the fast reconnection stage. For the first time, the detailed features and/or fine structures inside the CME/flare CS in the eruption were investigated in this work.
NASA Astrophysics Data System (ADS)
Mei, Z.; Shen, C.; Wu, N.; Lin, J.; Murphy, N. A.; Roussev, I. I.
2012-10-01
Magnetic reconnection plays a critical role in energy conversion during solar eruptions. This paper presents a set of magnetohydrodynamic experiments for the magnetic reconnection process in a current sheet (CS) formed in the wake of the rising flux rope. The eruption results from the loss of equilibrium in a magnetic configuration that includes a current-carrying flux rope, representing a pre-existing filament. In order to study the fine structure and micro processes inside the CS, mesh refinement is used to reduce the numerical diffusion. We start with a uniform, explicitly defined resistivity which results in a Lundquist number S = 104 in the vicinity of CS. The use of mesh refinement allows the simulation to capture high-resolution features such as plasmoids from the tearing mode and plasmoid instability regions of turbulence and slow-mode shocks. Inside the CS, magnetic reconnection goes through the Sweet-Parker and the fractal stages, and eventually displays a time-dependent Petschek pattern. Our results support the concept of fractal reconnection suggested by Shibata et al. and Shibata & Tanuma, and also suggest that the CS evolves through Sweet-Parker reconnection prior to the fast reconnection stage. For the first time, the detailed features and/or fine structures inside the coronal mass ejection/flare CS in the eruption were investigated in this work.
Numerical and laboratory experiments on the dynamics of plume-ridge interaction. Progress report
Kincaid, C.; Gable, C.W.
1995-09-01
Mantle plumes and passive upwelling beneath ridges are the two dominant modes of mantle transport and thermal/chemical fluxing between the Earth`s deep interior and surface. While plumes and ridges independently contribute to crustal accretion, they also interact and the dispersion of plumes within the upper mantle is strongly modulated by mid-ocean ridges. The simplest mode of interaction, with the plume centered on the ridge, has been well documented and modeled. The remaining question is how plumes and ridges interact when the plume is located off-axis; it has been suggested that a pipeline-like flow from the off-axis plume to the ridge axis at the base of the rigid lithosphere may develop. Mid-ocean ridges migrating away from hot mantle plumes can be affected by plume discharges over long times and ridge migration distances. Salient feature of this model is that off-axis plumes communicate with the ridge through a channel resulting from the refraction and dispersion of an axi-symmetric plume conduit along the base of the sloping lithosphere. To test the dynamics of this model, a series of numerical and laboratory dynamic experiments on the problem of a fixed ridge and an off-axis buoyant upwelling were conducted. Results are discussed.
Numerical simulations of recent proton acceleration experiments with sub-100 TW laser systems
NASA Astrophysics Data System (ADS)
Sinigardi, Stefano
2016-09-01
Recent experiments carried out at the Italian National Research Center, National Optics Institute Department in Pisa, are showing interesting results regarding maximum proton energies achievable with sub-100 TW laser systems. While laser systems are being continuously upgraded in laboratories around the world, at the same time a new trend on stabilizing and making ion acceleration results reproducible is growing in importance. Almost all applications require a beam with fixed performance, so that the energy spectrum and the total charge exhibit moderate shot to shot variations. This result is surely far from being achieved, but many paths are being explored in order to reach it. Some of the reasons for this variability come from fluctuations in laser intensity and focusing, due to optics instability. Other variation sources come from small differences in the target structure. The target structure can vary substantially, when it is impacted by the main pulse, due to the prepulse duration and intensity, the shape of the main pulse and the total energy deposited. In order to qualitatively describe the prepulse effect, we will present a two dimensional parametric scan of its relevant parameters. A single case is also analyzed with a full three dimensional simulation, obtaining reasonable agreement between the numerical and the experimental energy spectrum.
What can large-scale magnetohydrodynamic numerical experiments tell us about coronal heating?
Peter, H
2015-05-28
The upper atmosphere of the Sun is governed by the complex structure of the magnetic field. This controls the heating of the coronal plasma to over a million kelvin. Numerical experiments in the form of three-dimensional magnetohydrodynamic simulations are used to investigate the intimate interaction between magnetic field and plasma. These models allow one to synthesize the coronal emission just as it would be observed by real solar instrumentation. Large-scale models encompassing a whole active region form evolving coronal loops with properties similar to those seen in extreme ultraviolet light from the Sun, and reproduce a number of average observed quantities. This suggests that the spatial and temporal distributions of the heating as well as the energy distribution of individual heat deposition events in the model are a good representation of the real Sun. This provides evidence that the braiding of fieldlines through magneto-convective motions in the photosphere is a good concept to heat the upper atmosphere of the Sun. PMID:25897097
Hansen, C.; Victor, B.; Morgan, K.; Hossack, A.; Sutherland, D.; Jarboe, T.; Nelson, B. A.; Marklin, G.
2015-05-15
We present application of three scalar metrics derived from the Biorthogonal Decomposition (BD) technique to evaluate the level of agreement between macroscopic plasma dynamics in different data sets. BD decomposes large data sets, as produced by distributed diagnostic arrays, into principal mode structures without assumptions on spatial or temporal structure. These metrics have been applied to validation of the Hall-MHD model using experimental data from the Helicity Injected Torus with Steady Inductive helicity injection experiment. Each metric provides a measure of correlation between mode structures extracted from experimental data and simulations for an array of 192 surface-mounted magnetic probes. Numerical validation studies have been performed using the NIMROD code, where the injectors are modeled as boundary conditions on the flux conserver, and the PSI-TET code, where the entire plasma volume is treated. Initial results from a comprehensive validation study of high performance operation with different injector frequencies are presented, illustrating application of the BD method. Using a simplified (constant, uniform density and temperature) Hall-MHD model, simulation results agree with experimental observation for two of the three defined metrics when the injectors are driven with a frequency of 14.5 kHz.
Dynamics of soil water evaporation during soil drying: laboratory experiment and numerical analysis.
Han, Jiangbo; Zhou, Zhifang
2013-01-01
Laboratory and numerical experiments were conducted to investigate the evolution of soil water evaporation during a continuous drying event. Simulated soil water contents and temperatures by the calibrated model well reproduced measured values at different depths. Results show that the evaporative drying process could be divided into three stages, beginning with a relatively high evaporation rate during stage 1, followed by a lower rate during transient stage and stage 2, and finally maintaining a very low and constant rate during stage 3. The condensation zone was located immediately below the evaporation zone in the profile. Both peaks of evaporation and condensation rate increased rapidly during stage 1 and transition stage, decreased during stage 2, and maintained constant during stage 3. The width of evaporation zone kept a continuous increase during stages 1 and 2 and maintained a nearly constant value of 0.68 cm during stage 3. When the evaporation zone totally moved into the subsurface, a dry surface layer (DSL) formed above the evaporation zone at the end of stage 2. The width of DSL also presented a continuous increase during stage 2 and kept a constant value of 0.71 cm during stage 3. PMID:24489492
On the generation of sound by turbulent convection. I - A numerical experiment. [in solar interior
NASA Technical Reports Server (NTRS)
Bogdan, Thomas J.; Cattaneo, Fausto; Malagoli, Andrea
1993-01-01
Motivated by the problem of the origin of the solar p-modes, we study the generation of acoustic waves by turbulent convection. Our approach uses the results of high-resolution 3D simulations as the experimental basis for our investigation. The numerical experiment describes the evolution of a horizontally periodic layer of vigorously convecting fluid. The sound is measured by a procedure, based on a suitable linearization of the equations of compressible convection that allows the amplitude of the acoustic field to be determined. Through this procedure we identify unambiguously some 400 acoustic modes. The total energy of the acoustic field is found to be a fraction of a percent of the kinetic energy of the convection. The amplitudes of the observed modes depend weakly on (horizontal) wavenumber but strongly on frequency. The line widths of the observed modes typically exceed the natural linewidths of the modes as inferred from linear theory. This broadening appears to be related to the (stochastic) interaction between the modes and the underlying turbulence which causes abrupt, episodic events during which the phase coherence of the modes is lost.
NASA Astrophysics Data System (ADS)
Beucher, Romain; Huismans, Ritske S.
2016-04-01
Extension of the continental lithosphere can lead to the formation of a wide range of rifted margins styles with contrasting tectonic and geomorphological characteristics. It is now understood that many of these characteristics depend on the manner extension is distributed depending on (among others factors) rheology, structural inheritance, thermal structure and surface processes. The relative importance and the possible interactions of these controlling factors is still largely unknown. Here we investigate the feedbacks between tectonics and the transfers of material at the surface resulting from erosion, transport, and sedimentation. We use large-scale (1200 x 600 km) and high-resolution (~1km) numerical experiments coupling a 2D upper-mantle-scale thermo-mechanical model with a plan-form 2D surface processes model (SPM). We test the sensitivity of the coupled models to varying crust-lithosphere rheology and erosional efficiency ranging from no-erosion to very efficient erosion. We discuss how fast, when and how the topography of the continents evolves and how it can be compared to actual passive margins escarpment morphologies. We show that although tectonics is the main factor controlling the rift geometry, transfers of masses at the surface affect the timing of faulting and the initiation of sea-floor spreading. We discuss how such models may help to understand the evolution of high-elevated passive margins around the world.
Numerical Experiments on Acoustic Reciprocity in Compressible Potential Flows in Ducts
NASA Astrophysics Data System (ADS)
EVERSMAN, W.
2001-09-01
A reciprocity theorem for the scattering matrix for the propagation of acoustic modes in a duct with acoustically hard walls or with acoustically absorbing walls has been given in a companion publication. It was found that for a source at a specified end of the duct, suitably scaled reflection matrices in direct and reverse flow have a reciprocal relationship. Scaled transmission matrices obtained for direct flow and reversed flow with simultaneous switching of source location from one end to the other also have a reciprocal relationship. A reverse flow theorem for the equivalent one-dimensional propagation model, which is a good approximation to the three-dimensional model at low frequencies, was also obtained. In this case, using reciprocity and acoustic power conservation arguments it is additionally found that the acoustic power transmission coefficient is the same for a source at either end of the duct for a given flow direction. This result leads to an invariance theorem which relates acoustic power propagated due to sources of equal pressure amplitude at the two ends of the duct. A numerical verification of these reciprocal relationships is given here for propagation in axially symmetric (circular and annular) ducts with multi-modal propagation and at low frequencies when a one-dimensional model is appropriate.
The Additional-Mass Effect of Plates as Determined by Experiments
NASA Technical Reports Server (NTRS)
Gracey, William
1941-01-01
The apparent increase in the inertia properties of a body moving in a fluid medium has been called the additional-mass effect. This report presents a resume of test procedures and results of experimental determinations of the additional-mass effect of flat plates. In addition to data obtained from various foreign sources and from a NACA investigation in 1933, the results of tests recently conducted by the National Advisory Committee for Aeronautics are included.
NASA Astrophysics Data System (ADS)
Beyer, C.; Ballarini, E.; Bauer, R.; Griebler, C.; Bauer, S.
2011-12-01
The biodegradation of oxidizable hydrocarbon contaminants in the subsurface requires the presence of compatible microbial communities as well as sufficient amounts of electron acceptors and nutrients. In this context, transverse mixing, driven by dispersion and diffusion, is one of the main mechanisms governing the availability of dissolved electron acceptors at a hydrocarbon plume fringe. Aerobic and anaerobic biodegradation of hydrocarbons limited by transverse mixing has been studied experimentally in 2D bench-scale flow-through tanks, filled with a saturated porous medium. Flow of groundwater through the tanks was induced by pumping water at one side through injection ports, and simultaneously extracting water at the other side of the tank. An ethylbenzene plume was established by injection through the central inlet port. A mixture of unlabeled and fully deuterium-labeled isotopomers was used in order to investigate the spatial distribution of degradation processes via monitoring of compound-specific stable isotope fractionation. In the first phase of the experiment, aerobic biodegradation was studied. For this purpose, the tank was recharged with water containing oxygen as a dissolved electron acceptor and the aerobic strain Pseudomonas putida F1 was inoculated. Later, nitrate was added to the recharge water as an additional electron acceptor and the denitrifying strain Aromatoleum aromaticum EbN1 was amended to study competitive aerobic/anaerobic biodegradation. A numerical reactive transport model of the experiment was set up for a model based interpretation of the observed degradation patterns. In a sensitivity analysis, the influence of the relevant hydrodynamic parameters on the observable distributions of ethylbenzene isotopomers, oxygen and nitrate was studied. Subsequent model calibration allowed for a good agreement with ethylbenzene concentrations measured at the tank outlet ports as well as oxygen concentrations, which were measured at several
A numerical tracing experiment of water erosion at the plot scale
NASA Astrophysics Data System (ADS)
Nord, G.; Esteves, M.
2009-04-01
A numerical experiment of water erosion was set up on a plot of 100m2 (20m long by 5m wide) using PSEM_2D (Plot Soil Erosion Model). Various environmental factors were tested: six topographies combining different slopes and microtopography, four extreme rainfall events representative of the Temperate and the Mediterranean climates, four different textured soils, and two conditions of upstream injected discharge. 192 simulations were run. The results enabled to study the processes that control erosion and sediment yield at the plot scale. A notable contribution of this study was the possibility to trace numerically eroded sediment. Sediment originated from five successive zones of 20m2 (4m long by 5m wide) was labelled in a different way and traced to assess re-deposition within the plot and exportation at the outlet of the plot. Over the 100m2 plot, the contribution of rainfall erosion processes was dominant for most of the simulations without upstream injected discharge. However, in certain conditions of runoff discharge, slope, and critical shear stress, flow erosion was activated in the downer part of the plot, becoming the major source of sediment and producing substantial sediment yield. For the simulations with upstream injected discharge, flow erosion processes were generally the primary source of sediment over the whole plot and sediment yield was very high. Sediment that left the 100m2 plot originated mostly from the downer part of the slope. In the case of low slope and break in slope, sediment came exclusively from the toeslope pointing out the efficiency of the break in slope to reduce significantly sediment yield. In the case of topographies with microrelief, an increase in slope caused an increase in both sediment yield and maximum travel distance of sediment. In the case of fine sediment, plane surfaces, and upstream injected discharge, the contribution of sediment from upper parts of the plot and the maximum travel distance of sediment were enhanced
NASA Astrophysics Data System (ADS)
Sotiropoulos, F.; Kang, S.; Chamorro, L. P.; Hill, C.
2011-12-01
experimentally in the St. Anthony Falls Laboratory Main Channel. The experiments and simulations are compared with each other and shown to be in very good agreement both in terms of the mean flow and the turbulence statistics. The results are analyzed to study the structure of turbulence in the wake of the turbine and also identify the effects of turbulent fluctuations in the approach flow on the power produced by the turbine. Overall our results make a strong case that high-resolution numerical modeling, validated with detailed laboratory measurements, is a viable tool for assessing and optimizing the performance of MHK devices.
Raghavan, Narendran; Dehoff, Ryan; Pannala, Sreekanth; Simunovic, Srdjan; Kirka, Michael; Turner, John; Carlson, Neil; Babu, Sudarsanam S.
2016-04-26
The fabrication of 3-D parts from CAD models by additive manufacturing (AM) is a disruptive technology that is transforming the metal manufacturing industry. The correlation between solidification microstructure and mechanical properties has been well understood in the casting and welding processes over the years. This paper focuses on extending these principles to additive manufacturing to understand the transient phenomena of repeated melting and solidification during electron beam powder melting process to achieve site-specific microstructure control within a fabricated component. In this paper, we have developed a novel melt scan strategy for electron beam melting of nickel-base superalloy (Inconel 718) andmore » also analyzed 3-D heat transfer conditions using a parallel numerical solidification code (Truchas) developed at Los Alamos National Laboratory. The spatial and temporal variations of temperature gradient (G) and growth velocity (R) at the liquid-solid interface of the melt pool were calculated as a function of electron beam parameters. By manipulating the relative number of voxels that lie in the columnar or equiaxed region, the crystallographic texture of the components can be controlled to an extent. The analysis of the parameters provided optimum processing conditions that will result in columnar to equiaxed transition (CET) during the solidification. Furthermore, the results from the numerical simulations were validated by experimental processing and characterization thereby proving the potential of additive manufacturing process to achieve site-specific crystallographic texture control within a fabricated component.« less
NASA Astrophysics Data System (ADS)
Li, Y.; Ma, X.; Su, N.
2013-12-01
The movement of water and solute into and through the vadose zone is, in essence, an issue of immiscible displacement in pore-space network of a soil. Therefore, multiphase flow and transport in porous media, referring to three medium: air, water, and the solute, pose one of the largest unresolved challenges for porous medium fluid seepage. However, this phenomenon has always been largely neglected. It is expected that a reliable analysis model of the multi-phase flow in soil can truly reflect the process of natural movement about the infiltration, which is impossible to be observed directly. In such cases, geophysical applications of the nuclear magnetic resonance (NMR) provides the opportunity to measure the water movements into soils directly over a large scale from tiny pore to regional scale, accordingly enable it available both on the laboratory and on the field. In addition, the NMR provides useful information about the pore space properties. In this study, we proposed both laboratory and field experiments to measure the multi-phase flow parameters, together with optimize the model in computer programming based on the fractional partial differential equations (fPDE). In addition, we establish, for the first time, an infiltration model including solute flowing with water, which has huge influence on agriculture and soil environment pollution. Afterwards, with data collected from experiments, we simulate the model and analyze the spatial variability of parameters. Simulations are also conducted according to the model to evaluate the effects of airflow on water infiltration and other effects such as solute and absorption. It has significant meaning to oxygen irrigation aiming to higher crop yield, and shed more light into the dam slope stability. In summary, our framework is a first-time model added in solute to have a mathematic analysis with the fPDE and more instructive to agriculture activities.
NASA Astrophysics Data System (ADS)
Sawyer, A. H.; Cardenas, M. B.; Buttles, J. L.
2010-12-01
The flow of river water around wood debris creates pressure gradients along the riverbed that drive a large zone of river-groundwater mixing, or hyporheic exchange. River water downwells into the riverbed upstream of a channel-spanning log and upwells downstream. Using laboratory flume experiments and coupled CFD-groundwater flow simulations, we develop a predictive relationship for the pressure wave created by a channel-spanning log and resulting hyporheic flow. Amplitude of the pressure wave (and thus hyporheic exchange) increases with channel Froude number and blockage ratio (log diameter:channel flow depth). This relationship can be used to predict losses or gains in river-groundwater connectivity due to removal or addition of instream wood debris. We tested this relationship in a second-order stream in Valles Caldera National Preserve (NM). Log additions created alternating zones of upwelling and downwelling in a reach that was previously losing throughout. Though losing conditions limited the magnitude of exchange, the fundamental shift in exchange patterns from losing to nested flows could influence habitat complexity and transport of nutrients and energy.
NASA Astrophysics Data System (ADS)
Chen, Alexander Y.; Beloborodov, Andrei M.
2014-11-01
We present the first self-consistent global simulations of pulsar magnetospheres with operating e ± discharge. We focus on the simple configuration of an aligned or anti-aligned rotator. The star is spun up from a zero (vacuum) state to a high angular velocity, and we follow the coupled evolution of its external electromagnetic field and plasma particles using the "particle-in-cell" method. A plasma magnetosphere begins to form through the extraction of particles from the star; these particles are accelerated by the rotation-induced electric field, producing curvature radiation and igniting e ± discharge. We follow the system evolution for several revolution periods, longer than required to reach a quasi-steady state. Our numerical experiment puts to test previous ideas for the plasma flow and gaps in the pulsar magnetosphere. We first consider rotators capable of producing pairs out to the light cylinder through photon-photon collisions. We find that their magnetospheres are similar to the previously obtained force-free solutions with a Y-shaped current sheet. The magnetosphere continually ejects e ± pairs and ions. Pair creation is sustained by a strong electric field along the current sheet. We observe powerful curvature and synchrotron emission from the current sheet, consistent with Fermi observations of gamma-ray pulsars. We then study pulsars that can only create pairs in the strong-field region near the neutron star, well inside the light cylinder. We find that both aligned and anti-aligned rotators relax to the "dead" state with suppressed pair creation and electric currents, regardless of the discharge voltage.
Sediment pulses in mountain rivers: 2. Comparison between experiments and numerical predictions
NASA Astrophysics Data System (ADS)
Cui, Yantao; Parker, Gary; Pizzuto, James; Lisle, Thomas E.
2003-09-01
Mountain rivers in particular are prone to sediment input in the form of pulses rather than a more continuous supply. These pulses often enter in the form of landslides from adjacent hillslopes or debris flows from steeper tributaries. The activities of humans such as timber harvesting, road building, and urban development can increase the frequency of sediment pulses. The question as to how mountain rivers accommodate pulses of sediment thus becomes of practical as well as academic significance. In part 1 [, 2003], the results of three laboratory experiments on sediment pulses are reported. It was found there that the pulses were eliminated from the flume predominantly by dispersion of the topographic high. Significant translation was observed only when the pulse material was substantially finer than the ambient load in the river. Here the laboratory data are used to test a numerical model originally devised for predicting the evolution of sediment pulses in field-scale gravel bed streams. The model successfully reproduces the predominantly dispersive deformation of the experimental pulses. Rates of dispersion are generally underestimated, largely because bed load transport rates are underestimated by the transport equation used in the model. The model reproduces the experimental data best when the pulse is significantly coarser than the ambient sediment. In this case, the model successfully predicts the formation and downstream progradation of a delta that formed in the backwater zone of the pulse in run 3. The performance of the model is less successful when the pulse is composed primarily of sand. This is likely because the bed load equation used in the study is specifically designed for gravel. When the model is adapted to conditions characteristic of large, sand bed rivers with low Froude numbers, it predicts substantial translation of pulses as well as dispersion.