NASA Astrophysics Data System (ADS)
Yang, Wenming; An, Hui; Amin, Maghbouli; Li, Jing
2014-11-01
A 3-dimensional computational fluid dynamics modeling is conducted on a direct injection diesel engine fueled by biodiesel using multi-dimensional software KIVA4 coupled with CHEMKIN. To accurately predict the oxidation of saturated and unsaturated agents of the biodiesel fuel, a multicomponent advanced combustion model consisting of 69 species and 204 reactions combined with detailed oxidation pathways of methyl decenoate (C11H22O2), methyl-9-decenoate (C11H20O2) and n-heptane (C7H16) is employed in this work. In order to better represent the real fuel properties, the detailed chemical and thermo-physical properties of biodiesel such as vapor pressure, latent heat of vaporization, liquid viscosity and surface tension were calculated and compiled into the KIVA4 fuel library. The nitrogen monoxide (NO) and carbon monoxide (CO) formation mechanisms were also embedded. After validating the numerical simulation model by comparing the in-cylinder pressure and heat release rate curves with experimental results, further studies have been carried out to investigate the effect of combustion chamber design on flow field, subsequently on the combustion process and performance of diesel engine fueled by biodiesel. Research has also been done to investigate the impact of fuel injector location on the performance and emissions formation of diesel engine.
3-Dimensional Topographic Models for the Classroom
NASA Technical Reports Server (NTRS)
Keller, J. W.; Roark, J. H.; Sakimoto, S. E. H.; Stockman, S.; Frey, H. V.
2003-01-01
We have recently undertaken a program to develop educational tools using 3-dimensional solid models of digital elevation data acquired by the Mars Orbital Laser Altimeter (MOLA) for Mars as well as a variety of sources for elevation data of the Earth. This work is made possible by the use of rapid prototyping technology to construct solid 3-Dimensional models of science data. We recently acquired rapid prototyping machine that builds 3-dimensional models in extruded plastic. While the machine was acquired to assist in the design and development of scientific instruments and hardware, it is also fully capable of producing models of spacecraft remote sensing data. We have demonstrated this by using Mars Orbiter Laser Altimeter (MOLA) topographic data and Earth based topographic data to produce extruded plastic topographic models which are visually appealing and instantly engage those who handle them.
Incorporating 3-dimensional models in online articles
Cevidanes, Lucia H. S.; Ruellasa, Antonio C. O.; Jomier, Julien; Nguyen, Tung; Pieper, Steve; Budin, Francois; Styner, Martin; Paniagua, Beatriz
2015-01-01
Introduction The aims of this article were to introduce the capability to view and interact with 3-dimensional (3D) surface models in online publications, and to describe how to prepare surface models for such online 3D visualizations. Methods Three-dimensional image analysis methods include image acquisition, construction of surface models, registration in a common coordinate system, visualization of overlays, and quantification of changes. Cone-beam computed tomography scans were acquired as volumetric images that can be visualized as 3D projected images or used to construct polygonal meshes or surfaces of specific anatomic structures of interest. The anatomic structures of interest in the scans can be labeled with color (3D volumetric label maps), and then the scans are registered in a common coordinate system using a target region as the reference. The registered 3D volumetric label maps can be saved in .obj, .ply, .stl, or .vtk file formats and used for overlays, quantification of differences in each of the 3 planes of space, or color-coded graphic displays of 3D surface distances. Results All registered 3D surface models in this study were saved in .vtk file format and loaded in the Elsevier 3D viewer. In this study, we describe possible ways to visualize the surface models constructed from cone-beam computed tomography images using 2D and 3D figures. The 3D surface models are available in the article’s online version for viewing and downloading using the reader’s software of choice. These 3D graphic displays are represented in the print version as 2D snapshots. Overlays and color-coded distance maps can be displayed using the reader’s software of choice, allowing graphic assessment of the location and direction of changes or morphologic differences relative to the structure of reference. The interpretation of 3D overlays and quantitative color-coded maps requires basic knowledge of 3D image analysis. Conclusions When submitting manuscripts, authors can
NASA Astrophysics Data System (ADS)
Yoshida, Hiroyuki; Misawa, Takeharu; Takase, Kazuyuki
Two-fluid model can simulate two-phase flow by computational cost less than detailed two-phase flow simulation method such as interface tracking method or particle interaction method. Therefore, two-fluid model is useful for thermal hydraulic analysis in large-scale domain such as a rod bundle. Japan Atomic Energy Agency (JAEA) develops three dimensional two-fluid model analysis code ACE-3D that adopts boundary fitted coordinate system in order to simulate complex shape flow channel. In this paper, boiling two-phase flow analysis in a tight-lattice rod bundle was performed by ACE-3D code. The parallel computation using 126 CPUs was applied to this analysis. In the results, the void fraction, which distributes in outermost region of rod bundle, is lower than that in center region of rod bundle. The tendency of void fraction distribution agreed with the measurement results by neutron radiography qualitatively. To evaluate effects of two-phase flow model used in ACE-3D code, numerical simulation of boiling two-phase in tight-lattice rod bundle with no lift force model was also performed. From the comparison of calculated results, it was concluded that the effects of lift force model were not so large for overall void fraction distribution of tight-lattice rod bundle. However, the lift force model is important for local void fraction distribution of fuel bundles.
Water uptake by a maize root system - An explicit numerical 3-dimensional simulation.
NASA Astrophysics Data System (ADS)
Leitner, Daniel; Schnepf, Andrea; Klepsch, Sabine; Roose, Tiina
2010-05-01
Water is one of the most important resources for plant growth and function. An accurate modelling of the unsaturated flow is not only substantial to predict water uptake but also important to describe nutrient movement regarding water saturation and transport. In this work we present a model for water uptake. The model includes the simultaneous flow of water inside the soil and inside the root network. Water saturation in the soil volume is described by the Richards equation. Water flow inside the roots' xylem is calculated using the Poiseuille law for water flow in a cylindrical tube. The water saturation in the soil as well as water uptake of the root system is calculated numerically in three dimensions. We study water uptake of a maize plant in a confined pot under different supply scenarios. The main improvement of our approach is that the root surfaces act as spatial boundaries of the soil volume. Therefore water influx into the root is described by a surface flux instead of a volume flux, which is commonly given by an effective sink term. For the numerical computation we use the following software: The 3-dimensional maize root architecture is created by a root growth model based on L-Systems (Leitner et al 2009). A mesh of the surrounding soil volume is created using the meshing software DistMesh (Persson & Strang 2004). Using this mesh the partial differential equations are solved with the finite element method using Comsol Multiphysics 3.5a. Modelling results are related to accepted water uptake models from literature (Clausnitzer & Hopmans 1994, Roose & Fowler 2004, Javaux et al 2007). This new approach has several advantages. By considering the individual roots it is possible to analyse the influence of overlapping depletion zones due to inter root competition. Furthermore, such simulations can be used to estimate the influence of simplifying assumptions that are made in the development of effective models. The model can be easily combined with a nutrient
Development and Validation of a 3-Dimensional CFB Furnace Model
NASA Astrophysics Data System (ADS)
Vepsäläinen, Arl; Myöhänen, Karl; Hyppäneni, Timo; Leino, Timo; Tourunen, Antti
At Foster Wheeler, a three-dimensional CFB furnace model is essential part of knowledge development of CFB furnace process regarding solid mixing, combustion, emission formation and heat transfer. Results of laboratory and pilot scale phenomenon research are utilized in development of sub-models. Analyses of field-test results in industrial-scale CFB boilers including furnace profile measurements are simultaneously carried out with development of 3-dimensional process modeling, which provides a chain of knowledge that is utilized as feedback for phenomenon research. Knowledge gathered by model validation studies and up-to-date parameter databases are utilized in performance prediction and design development of CFB boiler furnaces. This paper reports recent development steps related to modeling of combustion and formation of char and volatiles of various fuel types in CFB conditions. Also a new model for predicting the formation of nitrogen oxides is presented. Validation of mixing and combustion parameters for solids and gases are based on test balances at several large-scale CFB boilers combusting coal, peat and bio-fuels. Field-tests including lateral and vertical furnace profile measurements and characterization of solid materials provides a window for characterization of fuel specific mixing and combustion behavior in CFB furnace at different loads and operation conditions. Measured horizontal gas profiles are projection of balance between fuel mixing and reactions at lower part of furnace and are used together with both lateral temperature profiles at bed and upper parts of furnace for determination of solid mixing and combustion model parameters. Modeling of char and volatile based formation of NO profiles is followed by analysis of oxidizing and reducing regions formed due lower furnace design and mixing characteristics of fuel and combustion airs effecting to formation ofNO furnace profile by reduction and volatile-nitrogen reactions. This paper presents
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
NASA Technical Reports Server (NTRS)
Bravo, Ramiro H.; Chen, Ching-Jen
1992-01-01
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
3-Dimensional modeling of protein structures distinguishes closely related phytoplasmas
Technology Transfer Automated Retrieval System (TEKTRAN)
Phytoplasmas (formerly mycoplasmalike organisms, MLOs) are cell wall-less bacteria that inhabit phloem tissue of plants and are transmitted from plant-to-plant by phloem-feeding insects. Numerous diseases affecting hundreds of plant species in many botanical families are attributed to infections by...
An Explicit 3-Dimensional Model for Reactive Transport of Nitrogen in Tile Drained Fields
NASA Astrophysics Data System (ADS)
Hill, D. J.; Valocchi, A. J.; Hudson, R. J.
2001-12-01
Recently, there has been increased interest in nitrate contamination of groundwater in the Midwest because of its link to surface water eutrophication, especially in the Gulf of Mexico. The vast majority of this nitrate is the product of biologically mediated transformation of fertilizers containing ammonia in the vadose zone of agricultural fields. For this reason, it is imperative that mathematical models, which can serve as useful tools to evaluate both the impact of agricultural fertilizer applications and nutrient-reducing management practices, are able to specifically address transport in the vadose zone. The development of a 3-dimensional explicit numerical model to simulate the movement and transformation of nitrogen species through the subsurface on the scale of an individual farm plot will be presented. At this scale, nitrogen fate and transport is controlled by a complex coupling among hydrologic, agricultural and biogeochemical processes. The nitrogen model is a component of a larger modeling effort that focuses upon conditions typical of those found in agricultural fields in Illinois. These conditions include non-uniform, multi-dimensional, transient flow in both saturated and unsaturated zones, geometrically complex networks of tile drains, coupled surface-subsurface-tile flow, and dynamic levels of dissolved oxygen in the soil profile. The advection-dispersion-reaction equation is solved using an operator-splitting approach, which is a flexible and straightforward strategy. Advection is modeled using a total variation diminishing scheme, dispersion is modeled using an alternating direction explicit method, and reactions are modeled using rate law equations. The model's stability and accuracy will be discussed, and test problems will be presented.
3-dimensional orthodontics visualization system with dental study models and orthopantomograms
NASA Astrophysics Data System (ADS)
Zhang, Hua; Ong, S. H.; Foong, K. W. C.; Dhar, T.
2005-04-01
The aim of this study is to develop a system that provides 3-dimensional visualization of orthodontic treatments. Dental plaster models and corresponding orthopantomogram (dental panoramic tomogram) are first digitized and fed into the system. A semi-auto segmentation technique is applied to the plaster models to detect the dental arches, tooth interstices and gum margins, which are used to extract individual crown models. 3-dimensional representation of roots, generated by deforming generic tooth models with orthopantomogram using radial basis functions, is attached to corresponding crowns to enable visualization of complete teeth. An optional algorithm to close the gaps between deformed roots and actual crowns by using multi-quadratic radial basis functions is also presented, which is capable of generating smooth mesh representation of complete 3-dimensional teeth. User interface is carefully designed to achieve a flexible system with as much user friendliness as possible. Manual calibration and correction is possible throughout the data processing steps to compensate occasional misbehaviors of automatic procedures. By allowing the users to move and re-arrange individual teeth (with their roots) on a full dentition, this orthodontic visualization system provides an easy and accurate way of simulation and planning of orthodontic treatment. Its capability of presenting 3-dimensional root information with only study models and orthopantomogram is especially useful for patients who do not undergo CT scanning, which is not a routine procedure in most orthodontic cases.
Creating 3-dimensional Models of the Photosphere using the SIR Code
NASA Astrophysics Data System (ADS)
Thonhofer, S.; Utz, D.; Jurčák, J.; Pauritsch, J.; Hanslmeier, A.; Lemmerer, B.
A high-resolution 3-dimensional model of the photospheric magnetic field is essential for the investigation of magnetic features such as sunspots, pores or smaller elements like single flux tubes seen as magnetic bright points. The SIR code is an advanced inversion code that retrieves physical quantities, e.g. magnetic field, from Stokes profiles. Based on this code, we developed a program for automated inversion of Hinode SOT/SP data and for storing these results in 3-dimensional data cubes in the form of fits files. We obtained models of the temperature, magnetic field strength, magnetic field angles and LOS-velocity in a region of the quiet sun. We will give a first discussion of those parameters in regards of small scale magnetic fields and what we can obtain and learn in the future.
MT3D: a 3 dimensional magnetotelluric modeling program (user's guide and documentation for Rev. 1)
Nutter, C.; Wannamaker, P.E.
1980-11-01
MT3D.REV1 is a non-interactive computer program written in FORTRAN to do 3-dimensional magnetotelluric modeling. A 3-D volume integral equation has been adapted to simulate the MT response of a 3D body in the earth. An integro-difference scheme has been incorporated to increase the accuracy. This is a user's guide for MT3D.REV1 on the University of Utah Research Institute's (UURI) PRIME 400 computer operating under PRIMOS IV, Rev. 17.
NASA Astrophysics Data System (ADS)
Kushner, Mark J.; Grapperhaus, Michael J.
1996-10-01
Inductively Coupled Plasma (ICP) reactors have the potential for scaling to large area substrates while maintaining azimuthal symmetry or side-to-side uniformity across the wafer. Asymmetric etch properties in these devices have been attributed to transmission line properties of the coil, internal structures (such as wafer clamps) and non-uniform gas injection or pumping. To investigate the origins of asymmetric etch properties, a 3-dimensional hybrid model has been developed. The hybrid model contains electromagnetic, electric circuit, electron energy equation, and fluid modules. Continuity and momentum equations are solved in the fluid module along with Poisson's equation. We will discuss results for ion and radical flux uniformity to the substrate while varying the transmission line characteristics of the coil, symmetry of gas inlets/pumping, and internal structures. Comparisons will be made to expermental measurements of etch rates. ^*Work supported by SRC, NSF, ARPA/AFOSR and LAM Research.
Using 3-dimensional printing to create presurgical models for endodontic surgery.
Bahcall, James K
2014-09-01
Advances in endodontic surgery--from both a technological and procedural perspective-have been significant over the last 18 years. Although these technologies and procedural enhancements have significantly improved endodontic surgical treatment outcomes, there is still an ongoing challenge of overcoming the limitations of interpreting preoperative 2-dimensional (2-D) radiographic representation of a 3-dimensional (3-D) in vivo surgical field. Cone-beam Computed Tomography (CBCT) has helped to address this issue by providing a 3-D enhancement of the 2-D radiograph. The next logical step to further improve a presurgical case 3-D assessment is to create a surgical model from the CBCT scan. The purpose of this article is to introduce 3-D printing of CBCT scans for creating presurgical models for endodontic surgery. PMID:25197746
NASA Astrophysics Data System (ADS)
Gumus, Kutalmis; Erkaya, Halil
2013-04-01
In Terrestrial laser scanning (TLS) applications, it is necessary to take into consideration the conditions that affect the scanning process, especially the general characteristics of the laser scanner, geometric properties of the scanned object (shape, size, etc.), and its spatial location in the environment. Three dimensional models obtained with TLS, allow determining the geometric features and relevant magnitudes of the scanned object in an indirect way. In order to compare the spatial location and geometric accuracy of the 3-dimensional model created by Terrestrial laser scanning, it is necessary to use measurement tools that give more precise results than TLS. Geometric comparisons are performed by analyzing the differences between the distances, the angles between surfaces and the measured values taken from cross-sections between the data from the 3-dimensional model created with TLS and the values measured by other measurement devices The performance of the scanners, the size and shape of the scanned objects are tested using reference objects the sizes of which are determined with high precision. In this study, the important points to consider when choosing reference objects were highlighted. The steps up to processing the point clouds collected by scanning, regularizing these points and modeling in 3 dimensions was presented visually. In order to test the geometric correctness of the models obtained by Terrestrial laser scanners, sample objects with simple geometric shapes such as cubes, rectangular prisms and cylinders that are made of concrete were used as reference models. Three dimensional models were generated by scanning these reference models with Trimble Mensi GS 100. The dimension of the 3D model that is created from point clouds was compared with the precisely measured dimensions of the reference objects. For this purpose, horizontal and vertical cross-sections were taken from the reference objects and generated 3D models and the proximity of
Siler, Drew L; Faulds, James E; Mayhew, Brett
2013-04-16
Geothermal systems in the Great Basin, USA, are controlled by a variety of fault intersection and fault interaction areas. Understanding the specific geometry of the structures most conducive to broad-scale geothermal circulation is crucial to both the mitigation of the costs of geothermal exploration (especially drilling) and to the identification of geothermal systems that have no surface expression (blind systems). 3-dimensional geologic modeling is a tool that can elucidate the specific stratigraphic intervals and structural geometries that host geothermal reservoirs. Astor Pass, NV USA lies just beyond the northern extent of the dextral Pyramid Lake fault zone near the boundary between two distinct structural domains, the Walker Lane and the Basin and Range, and exhibits characteristics of each setting. Both northwest-striking, left-stepping dextral faults of the Walker Lane and kinematically linked northerly striking normal faults associated with the Basin and Range are present. Previous studies at Astor Pass identified a blind geothermal system controlled by the intersection of west-northwest and north-northwest striking dextral-normal faults. Wells drilled into the southwestern quadrant of the fault intersection yielded 94°C fluids, with geothermometers suggesting a maximum reservoir temperature of 130°C. A 3-dimensional model was constructed based on detailed geologic maps and cross-sections, 2-dimensional seismic data, and petrologic analysis of the cuttings from three wells in order to further constrain the structural setting. The model reveals the specific geometry of the fault interaction area at a level of detail beyond what geologic maps and cross-sections can provide.
NASA Astrophysics Data System (ADS)
Metz, P. D.
A FORTRAN computer program called GROCS (GRound Coupled Systems) has been developed to study 3-dimensional underground heat flow. Features include the use of up to 30 finite elements or blocks of Earth which interact via finite difference heat flow equations and a subprogram which sets realistic time and depth dependent boundary conditions. No explicit consideration of mositure movement or freezing is given. GROCS has been used to model the thermal behavior of buried solar heat storage tanks (with and without insulation) and serpentine pipe fields for solar heat pump space conditioning systems. The program is available independently or in a form compatible with specially written TRNSYS component TYPE subroutines. The approach taken in the design of GROCS, the mathematics contained and the program architecture, are described. Then, the operation of the stand-alone version is explained. Finally, the validity of GROCS is discussed.
A 3-dimensional DTI MRI-based model of GBM growth and response to radiation therapy.
Hathout, Leith; Patel, Vishal; Wen, Patrick
2016-09-01
Glioblastoma (GBM) is both the most common and the most aggressive intra-axial brain tumor, with a notoriously poor prognosis. To improve this prognosis, it is necessary to understand the dynamics of GBM growth, response to treatment and recurrence. The present study presents a mathematical diffusion-proliferation model of GBM growth and response to radiation therapy based on diffusion tensor (DTI) MRI imaging. This represents an important advance because it allows 3-dimensional tumor modeling in the anatomical context of the brain. Specifically, tumor infiltration is guided by the direction of the white matter tracts along which glioma cells infiltrate. This provides the potential to model different tumor growth patterns based on location within the brain, and to simulate the tumor's response to different radiation therapy regimens. Tumor infiltration across the corpus callosum is simulated in biologically accurate time frames. The response to radiation therapy, including changes in cell density gradients and how these compare across different radiation fractionation protocols, can be rendered. Also, the model can estimate the amount of subthreshold tumor which has extended beyond the visible MR imaging margins. When combined with the ability of being able to estimate the biological parameters of invasiveness and proliferation of a particular GBM from serial MRI scans, it is shown that the model has potential to simulate realistic tumor growth, response and recurrence patterns in individual patients. To the best of our knowledge, this is the first presentation of a DTI-based GBM growth and radiation therapy treatment model. PMID:27572745
Gálvez, Jorge A; Gralewski, Kevin; McAndrew, Christine; Rehman, Mohamed A; Chang, Benjamin; Levin, L Scott
2016-03-01
Children are not typically considered for hand transplantation for various reasons, including the difficulty of finding an appropriate donor. Matching donor-recipient hands and forearms based on size is critically important. If the donor's hands are too large, the recipient may not be able to move the fingers effectively. Conversely, if the donor's hands are too small, the appearance may not be appropriate. We present an 8-year-old child evaluated for a bilateral hand transplant following bilateral amputation. The recipient forearms and model hands were modeled from computed tomography imaging studies and replicated as anatomic models with a 3-dimensional printer. We modified the scale of the printed hand to produce 3 proportions, 80%, 100% and 120%. The transplant team used the anatomical models during evaluation of a donor for appropriate match based on size. The donor's hand size matched the 100%-scale anatomical model hand and the transplant team was activated. In addition to assisting in appropriate donor selection by the transplant team, the 100%-scale anatomical model hand was used to create molds for prosthetic hands for the donor. PMID:26810827
NASA Astrophysics Data System (ADS)
Zamora, A.; Gutierrez, A. E.; Velasco, A. A.
2014-12-01
2- and 3-Dimensional models obtained from the inversion of geophysical data are widely used to represent the structural composition of the Earth and to constrain independent models obtained from other geological data (e.g. core samples, seismic surveys, etc.). However, inverse modeling of gravity data presents a very unstable and ill-posed mathematical problem, given that solutions are non-unique and small changes in parameters (position and density contrast of an anomalous body) can highly impact the resulting model. Through the implementation of an interior-point method constrained optimization technique, we improve the 2-D and 3-D models of Earth structures representing known density contrasts mapping anomalous bodies in uniform regions and boundaries between layers in layered environments. The proposed techniques are applied to synthetic data and gravitational data obtained from the Rio Grande Rift and the Cooper Flat Mine region located in Sierra County, New Mexico. Specifically, we improve the 2- and 3-D Earth models by getting rid of unacceptable solutions (those that do not satisfy the required constraints or are geologically unfeasible) given the reduction of the solution space.
The Effectiveness of an Interactive 3-Dimensional Computer Graphics Model for Medical Education
Konishi, Takeshi; Tamura, Yoko; Moriguchi, Hiroki
2012-01-01
Background Medical students often have difficulty achieving a conceptual understanding of 3-dimensional (3D) anatomy, such as bone alignment, muscles, and complex movements, from 2-dimensional (2D) images. To this end, animated and interactive 3-dimensional computer graphics (3DCG) can provide better visual information to users. In medical fields, research on the advantages of 3DCG in medical education is relatively new. Objective To determine the educational effectiveness of interactive 3DCG. Methods We divided 100 participants (27 men, mean (SD) age 17.9 (0.6) years, and 73 women, mean (SD) age 18.1 (1.1) years) from the Health Sciences University of Mongolia (HSUM) into 3DCG (n = 50) and textbook-only (control) (n = 50) groups. The control group used a textbook and 2D images, while the 3DCG group was trained to use the interactive 3DCG shoulder model in addition to a textbook. We conducted a questionnaire survey via an encrypted satellite network between HSUM and Tokushima University. The questionnaire was scored on a 5-point Likert scale from strongly disagree (score 1) to strongly agree (score 5). Results Interactive 3DCG was effective in undergraduate medical education. Specifically, there was a significant difference in mean (SD) scores between the 3DCG and control groups in their response to questionnaire items regarding content (4.26 (0.69) vs 3.85 (0.68), P = .001) and teaching methods (4.33 (0.65) vs 3.74 (0.79), P < .001), but no significant difference in the Web category. Participants also provided meaningful comments on the advantages of interactive 3DCG. Conclusions Interactive 3DCG materials have positive effects on medical education when properly integrated into conventional education. In particular, our results suggest that interactive 3DCG is more efficient than textbooks alone in medical education and can motivate students to understand complex anatomical structures. PMID:23611759
3-Dimensional Modeling of Capacitively and Inductively Coupled Plasma Etching Systems
NASA Astrophysics Data System (ADS)
Rauf, Shahid
2008-10-01
Low temperature plasmas are widely used for thin film etching during micro and nano-electronic device fabrication. Fluid and hybrid plasma models were developed 15-20 years ago to understand the fundamentals of these plasmas and plasma etching. These models have significantly evolved since then, and are now a major tool used for new plasma hardware design and problem resolution. Plasma etching is a complex physical phenomenon, where inter-coupled plasma, electromagnetic, fluid dynamics, and thermal effects all have a major influence. The next frontier in the evolution of fluid-based plasma models is where these models are able to self-consistently treat the inter-coupling of plasma physics with fluid dynamics, electromagnetics, heat transfer and magnetostatics. We describe one such model in this paper and illustrate its use in solving engineering problems of interest for next generation plasma etcher design. Our 3-dimensional plasma model includes the full set of Maxwell equations, transport equations for all charged and neutral species in the plasma, the Navier-Stokes equation for fluid flow, and Kirchhoff's equations for the lumped external circuit. This model also includes Monte Carlo based kinetic models for secondary electrons and stochastic heating, and can take account of plasma chemistry. This modeling formalism allows us to self-consistently treat the dynamics in commercial inductively and capacitively coupled plasma etching reactors with realistic plasma chemistries, magnetic fields, and reactor geometries. We are also able to investigate the influence of the distributed electromagnetic circuit at very high frequencies (VHF) on the plasma dynamics. The model is used to assess the impact of azimuthal asymmetries in plasma reactor design (e.g., off-center pump, 3D magnetic field, slit valve, flow restrictor) on plasma characteristics at frequencies from 2 -- 180 MHz. With Jason Kenney, Ankur Agarwal, Ajit Balakrishna, Kallol Bera, and Ken Collins.
Ku, L.P.; Kolibal, J.G.; Liew, S.L.
1985-09-01
The computational models of the TFTR constructed for the radiation transport analysis for the Q approx. 1 demonstration are summarized and reviewed. These models can be characterized by the dimensionality required to describe the geometry, and by the numerical methods of solving the transport equation. Results obtained with these models in the test cell are compared and discussed.
Andriani, Frank; Garfield, Jackie; Fusenig, Norbert E; Garlick, Jonathan A
2004-01-20
We have developed novel 3-dimensional in vitro and in vivo tissue models that mimic premalignant disease of human stratified epithelium in order to analyze the stromal contribution of extracellular matrix and basement membrane proteins to the progression of intraepithelial neoplasia. Three-dimensional, organotypic cultures were grown either on a de-epidermalized human dermis with pre-existing basement membrane components on its surface (AlloDerm), on a Type I collagen gel that lacked basement membrane proteins or on polycarbonate membranes coated with purified extracellular matrix proteins. When tumor cells (HaCaT-II4) were mixed with normal keratinocytes (4:1/normals:HaCaT-II4), tumor cells selectively attached, persisted and proliferated at the dermal-epidermal interface in vitro and generated dysplastic tissues when transplanted to nude mice only when grown in the presence of the AlloDerm substrate. This stromal interface was permissive for tumor cell attachment due to the rapid assembly of structured basement membrane. When tumor cells were mixed with normal keratinocytes and grown on polycarbonate membranes coated with individual extracellular matrix or basement membrane components, selective attachment and significant intraepithelial expansion occurred only on laminin 1 and Type IV collagen-coated membranes. This preferential adhesion of tumor cells restricted the synthesis of laminin 5 to basal cells where it was deposited in a polarized distribution. Western blot analysis revealed that tumor cell attachment was not due to differences in the synthesis or processing of laminin 5. Thus, intraepithelial progression towards premalignant disease is dependent on the selective adhesion of cells with malignant potential to basement membrane proteins that provide a permissive template for their persistence and expansion. PMID:14648700
Inouye, Joshua M; Lin, Kant Y; Perry, Jamie L; Blemker, Silvia S
2016-02-01
The convexity of the dorsal surface of the velum is critical for normal velopharyngeal (VP) function and is largely attributed to the levator veli palatini (LVP) and musculus uvulae (MU). Studies have correlated a concave or flat nasal velar surface to symptoms of VP dysfunction including hypernasality and nasal air emission. In the context of surgical repair of cleft palates, the MU has been given relatively little attention in the literature compared with the larger LVP. A greater understanding of the mechanics of the MU will provide insight into understanding the influence of a dysmorphic MU, as seen in cleft palate, as it relates to VP function. The purpose of this study was to quantify the contributions of the MU to VP closure in a computational model. We created a novel 3-dimensional (3D) finite element model of the VP mechanism from magnetic resonance imaging data collected from an individual with healthy noncleft VP anatomy. The model components included the velum, posterior pharyngeal wall (PPW), LVP, and MU. Simulations were based on the muscle and soft tissue mechanical properties from the literature. We found that, similar to previous hypotheses, the MU acts as (i) a space-occupying structure and (ii) a velar extensor. As a space-occupying structure, the MU helps to nearly triple the midline VP contact length. As a velar extensor, the MU acting alone without the LVP decreases the VP distance 62%. Furthermore, activation of the MU decreases the LVP activation required for closure almost 3-fold, from 20% (without MU) to 8% (with MU). Our study suggests that any possible salvaging and anatomical reconstruction of viable MU tissue in a cleft patient may improve VP closure due to its mechanical function. In the absence or dysfunction of MU tissue, implantation of autologous or engineered tissues at the velar midline, as a possible substitute for the MU, may produce a geometric convexity more favorable to VP closure. In the future, more complex models will
NASA Astrophysics Data System (ADS)
Fitzenz, D. D.; Miller, S. A.
2001-12-01
We present preliminary results from a 3-dimensional fault interaction model, with the fault system specified by the geometry and tectonics of the San Andreas Fault (SAF) system. We use the forward model for earthquake generation on interacting faults of Fitzenz and Miller [2001] that incorporates the analytical solutions of Okada [85,92], GPS-constrained tectonic loading, creep compaction and frictional dilatancy [Sleep and Blanpied, 1994, Sleep, 1995], and undrained poro-elasticity. The model fault system is centered at the Big Bend, and includes three large strike-slip faults (each discretized into multiple subfaults); 1) a 300km, right-lateral segment of the SAF to the North, 2) a 200km-long left-lateral segment of the Garlock fault to the East, and 3) a 100km-long right-lateral segment of the SAF to the South. In the initial configuration, three shallow-dipping faults are also included that correspond to the thrust belt sub-parallel to the SAF. Tectonic loading is decomposed into basal shear drag parallel to the plate boundary with a 35mm yr-1 plate velocity, and East-West compression approximated by a vertical dislocation surface applied at the far-field boundary resulting in fault-normal compression rates in the model space about 4mm yr-1. Our aim is to study the long-term seismicity characteristics, tectonic evolution, and fault interaction of this system. We find that overpressured faults through creep compaction are a necessary consequence of the tectonic loading, specifically where high normal stress acts on long straight fault segments. The optimal orientation of thrust faults is a function of the strike-slip behavior, and therefore results in a complex stress state in the elastic body. This stress state is then used to generate new fault surfaces, and preliminary results of dynamically generated faults will also be presented. Our long-term aim is to target measurable properties in or around fault zones, (e.g. pore pressures, hydrofractures, seismicity
3-Dimensional Marine CSEM Modeling by Employing TDFEM with Parallel Solvers
NASA Astrophysics Data System (ADS)
Wu, X.; Yang, T.
2013-12-01
In this paper, parallel fulfillment is developed for forward modeling of the 3-Dimensional controlled source electromagnetic (CSEM) by using time-domain finite element method (TDFEM). Recently, a greater attention rises on research of hydrocarbon (HC) reservoir detection mechanism in the seabed. Since China has vast ocean resources, seeking hydrocarbon reservoirs become significant in the national economy. However, traditional methods of seismic exploration shown a crucial obstacle to detect hydrocarbon reservoirs in the seabed with a complex structure, due to relatively high acquisition costs and high-risking exploration. In addition, the development of EM simulations typically requires both a deep knowledge of the computational electromagnetics (CEM) and a proper use of sophisticated techniques and tools from computer science. However, the complexity of large-scale EM simulations often requires large memory because of a large amount of data, or solution time to address problems concerning matrix solvers, function transforms, optimization, etc. The objective of this paper is to present parallelized implementation of the time-domain finite element method for analysis of three-dimensional (3D) marine controlled source electromagnetic problems. Firstly, we established a three-dimensional basic background model according to the seismic data, then electromagnetic simulation of marine CSEM was carried out by using time-domain finite element method, which works on a MPI (Message Passing Interface) platform with exact orientation to allow fast detecting of hydrocarbons targets in ocean environment. To speed up the calculation process, SuperLU of an MPI (Message Passing Interface) version called SuperLU_DIST is employed in this approach. Regarding the representation of three-dimension seabed terrain with sense of reality, the region is discretized into an unstructured mesh rather than a uniform one in order to reduce the number of unknowns. Moreover, high-order Whitney
Werner, Heron; Rolo, Liliam Cristine; Araujo Júnior, Edward; Dos Santos, Jorge Roberto Lopes
2014-03-01
Technological innovations accompanying advances in medicine have given rise to the possibility of obtaining better-defined fetal images that assist in medical diagnosis and contribute toward genetic counseling offered to parents during the prenatal period. In this article, we show our innovative experience of diagnosing fetal malformations through correlating 3-dimensional ultrasonography, magnetic resonance imaging, and computed tomography, which are accurate techniques for fetal assessment, with a fetal image reconstruction technique to create physical fetal models. PMID:24901782
NASA Astrophysics Data System (ADS)
Majumder, Sambit; Majumder, Abhik; Bhaumik, Swapan
2016-07-01
The present microelectronics market demands devices with high power dissipation capabilities having enhanced cooling per unit area. The drive for miniaturizing the devices to even micro level dimensions is shooting up the applied heat flux on such devices, resulting in complexity in heat transfer and cooling management. In this paper, a method of CPU processor cooling is introduced where active and passive cooling techniques are incorporated simultaneously. A heat sink consisting of fins is designed, where water flows internally through the mini-channel fins and air flows externally. Three dimensional numerical simulations are performed for large set of Reynolds number in laminar region using finite volume method for both developing flows. The dimensions of mini-channel fins are varied for several aspect ratios such as 1, 1.33, 2 and 4. Constant temperature (T) boundary condition is applied at heat sink base. Channel fluid temperature, pressure drop are analyzed to obtain best cooling option in the present study. It has been observed that as the aspect ratio of the channel decreases Nusselt number decreases while pressure drop increases. However, Nusselt number increases with increase in Reynolds number.
Maschio, Federico; Pandya, Mirali; Olszewski, Raphael
2016-01-01
BACKGROUND The objective of this study was to investigate the accuracy of 3-dimensional (3D) plastic (ABS) models generated using a low-cost 3D fused deposition modelling printer. MATERIAL AND METHODS Two human dry mandibles were scanned with a cone beam computed tomography (CBCT) Accuitomo device. Preprocessing consisted of 3D reconstruction with Maxilim software and STL file repair with Netfabb software. Then, the data were used to print 2 plastic replicas with a low-cost 3D fused deposition modeling printer (Up plus 2®). Two independent observers performed the identification of 26 anatomic landmarks on the 4 mandibles (2 dry and 2 replicas) with a 3D measuring arm. Each observer repeated the identifications 20 times. The comparison between the dry and plastic mandibles was based on 13 distances: 8 distances less than 12 mm and 5 distances greater than 12 mm. RESULTS The mean absolute difference (MAD) was 0.37 mm, and the mean dimensional error (MDE) was 3.76%. The MDE decreased to 0.93% for distances greater than 12 mm CONCLUSIONS Plastic models generated using the low-cost 3D printer UPplus2® provide dimensional accuracies comparable to other well-established rapid prototyping technologies. Validated low-cost 3D printers could represent a step toward the better accessibility of rapid prototyping technologies in the medical field. PMID:27003456
Maschio, Federico; Pandya, Mirali; Olszewski, Raphael
2016-01-01
Background The objective of this study was to investigate the accuracy of 3-dimensional (3D) plastic (ABS) models generated using a low-cost 3D fused deposition modelling printer. Material/Methods Two human dry mandibles were scanned with a cone beam computed tomography (CBCT) Accuitomo device. Preprocessing consisted of 3D reconstruction with Maxilim software and STL file repair with Netfabb software. Then, the data were used to print 2 plastic replicas with a low-cost 3D fused deposition modeling printer (Up plus 2®). Two independent observers performed the identification of 26 anatomic landmarks on the 4 mandibles (2 dry and 2 replicas) with a 3D measuring arm. Each observer repeated the identifications 20 times. The comparison between the dry and plastic mandibles was based on 13 distances: 8 distances less than 12 mm and 5 distances greater than 12 mm. Results The mean absolute difference (MAD) was 0.37 mm, and the mean dimensional error (MDE) was 3.76%. The MDE decreased to 0.93% for distances greater than 12 mm. Conclusions Plastic models generated using the low-cost 3D printer UPplus2® provide dimensional accuracies comparable to other well-established rapid prototyping technologies. Validated low-cost 3D printers could represent a step toward the better accessibility of rapid prototyping technologies in the medical field. PMID:27003456
In vitro 3-dimensional tumor model for radiosensitivity of HPV positive OSCC cell lines
Zhang, Mei; Rose, Barbara; Lee, C Soon; Hong, Angela M
2015-01-01
The incidence of oropharyngeal squamous cell carcinoma (OSCC) is increasing due to the rising prevalence of human papillomavirus (HPV) positive OSCC. HPV positive OSCC is associated with better outcomes than HPV negative OSCC. Our aim was to explore the possibility that this favorable prognosis is due to the enhanced radiosensitivity of HPV positive OSCC. HPV positive OSCC cell lines were generated from the primary OSCCs of 2 patients, and corresponding HPV positive cell lines generated from nodal metastases following xenografting in nude mice. Monolayer and 3 dimensional (3D) culture techniques were used to compare the radiosensitivity of HPV positive lines with that of 2 HPV negative OSCC lines. Clonogenic and protein assays were used to measure survival post radiation. Radiation induced cell cycle changes were studied using flow cytometry. In both monolayer and 3D culture, HPV positive cells exhibited a heterogeneous appearance whereas HPV negative cells tended to be homogeneous. After irradiation, HPV positive cells had a lower survival in clonogenic assays and lower total protein levels in 3D cultures than HPV negative cells. Irradiated HPV positive cells showed a high proportion of cells in G1/S phase, increased apoptosis, an increased proliferation rate, and an inability to form 3D tumor clumps. In conclusion, HPV positive OSCC cells are more radiosensitive than HPV negative OSCC cells in vitro, supporting a more radiosensitive nature of HPV positive OSCC. PMID:26046692
The 3-dimensional, 4-channel model of human visual sensitivity to grayscale scrambles.
Silva, Andrew E; Chubb, Charles
2014-08-01
Previous research supports the claim that human vision has three dimensions of sensitivity to grayscale scrambles (textures composed of randomly scrambled mixtures of different grayscales). However, the preattentive mechanisms (called here "field-capture channels") that confer this sensitivity remain obscure. The current experiments sought to characterize the specific field-capture channels that confer this sensitivity using a task in which the participant is required to detect the location of a small patch of one type of grayscale scramble in an extended background of another type. Analysis of the results supports the existence of four field-capture channels: (1) the (previously characterized) "blackshot" channel, sharply tuned to the blackest grayscales; (2) a (previously unknown) "gray-tuned" field-capture channel whose sensitivity is zero for black rising sharply to maximum sensitivity for grayscales slightly darker than mid-gray then decreasing to half-height for brighter grayscales; (3) an "up-ramped" channel whose sensitivity is zero for black, increases linearly with increasing grayscale reaching a maximum near white; (4) a (complementary) "down-ramped" channel whose sensitivity is maximal for black, decreases linearly reaching a minimum near white. The sensitivity functions of field-capture channels (3) and (4) are linearly dependent; thus, these four field-capture channels collectively confer sensitivity to a 3-dimensional space of histogram variations. PMID:24932891
Toward Scientific Numerical Modeling
NASA Technical Reports Server (NTRS)
Kleb, Bil
2007-01-01
Ultimately, scientific numerical models need quantified output uncertainties so that modeling can evolve to better match reality. Documenting model input uncertainties and verifying that numerical models are translated into code correctly, however, are necessary first steps toward that goal. Without known input parameter uncertainties, model sensitivities are all one can determine, and without code verification, output uncertainties are simply not reliable. To address these two shortcomings, two proposals are offered: (1) an unobtrusive mechanism to document input parameter uncertainties in situ and (2) an adaptation of the Scientific Method to numerical model development and deployment. Because these two steps require changes in the computational simulation community to bear fruit, they are presented in terms of the Beckhard-Harris-Gleicher change model.
Schmidt, Marianne; Scholz, Claus-Juergen; Polednik, Christine; Roller, Jeanette
2016-04-01
In the present study a panel of 12 head and neck cancer (HNSCC) cell lines were tested for spheroid formation. Since the size and morphology of spheroids is dependent on both cell adhesion and proliferation in the 3-dimensional (3D) context, morphology of HNSCC spheroids was related to expression of E-cadherin and the proliferation marker Ki67. In HNSCC cell lines the formation of tight regular spheroids was dependent on distinct E-cadherin expression levels in monolayer cultures, usually resulting in upregulation following aggregation into 3D structures. Cell lines expressing only low levels of E-cadherin in monolayers produced only loose cell clusters, frequently decreasing E-cadherin expression further upon aggregation. In these cell lines no epidermal growth factor receptor (EGFR) upregulation occurred and proliferation generally decreased in spheroids/aggregates independent of E-cadherin expression. In a second approach a global gene expression analysis of the larynx carcinoma cell line HLaC78 monolayer and the corresponding spheroids was performed. A global upregulation of gene expression in HLaC78 spheroids was related to genes involved in cell adhesion, cell junctions and cytochrome P450-mediated metabolism of xenobiotics. Downregulation was associated with genes controlling cell cycle, DNA-replication and DNA mismatch repair. Analyzing the expression of selected genes of each functional group in monolayer and spheroid cultures of all 12 cell lines revealed evidence for common gene expression shifts in genes controlling cell junctions, cell adhesion, cell cycle and DNA replication as well as genes involved in the cytochrome P450-mediated metabolism of xenobiotics. PMID:26797047
Goodall, Nicola; Kisiswa, Lilian; Prashar, Ankush; Faulkner, Stuart; Tokarczuk, Paweł; Singh, Krish; Erichsen, Jonathan T; Guggenheim, Jez; Halfter, Willi; Wride, Michael A
2009-10-01
Magnetic resonance imaging (MRI) is a powerful tool for generating 3-dimensional structural and functional image data. MRI has already proven valuable in creating atlases of mouse and quail development. Here, we have exploited high resolution MRI to determine the parameters necessary to acquire images of the chick embryo eye. Using a 9.4 Tesla (400 MHz) high field ultra-shielded and refrigerated magnet (Bruker), MRI was carried out on paraformaldehyde-fixed chick embryos or heads at E4, E6, E8, and E10. Image data were processed using established and custom packages (MRICro, ImageJ, ParaVision, Bruker and mri3dX). Voxel dimensions ranged from 62.5 microm to 117.2 microm. We subsequently used the images obtained from the MRI data in order to make precise measurements of chick embryo eye surface area, volume and axial length from E4 to E10. MRI was validated for accurate sizing of ocular tissue features by direct comparison with previously published literature. Furthermore, we demonstrate the utility of high resolution MRI for making accurate measurements of morphological changes due to experimental manipulation of chick eye development, thereby facilitating a better understanding of the effects on chick embryo eye development and growth of such manipulations. Chondroitin sulphate or heparin were microinjected into the vitreous cavity of the right eyes of each of 3 embryos at E5. At E10, embryos were fixed and various eye parameters (volume, surface area, axial length and equatorial diameter) were determined using MRI and normalised with respect to the un-injected left eyes. Statistically significant alterations in eye volume (p < 0.05; increases with chondroitin sulphate and decreases with heparin) and changes in vitreous homogeneity were observed in embryos following microinjection of glycosaminoglycans. Furthermore, in the heparin-injected eyes, significant disturbances at the vitreo-retinal boundary were observed as well as retinal folding and detachment
Application of 3-Dimensional Printing Technology to Construct an Eye Model for Fundus Viewing Study
Li, Xinhua; Gao, Zhishan; Yuan, Dongqing; Liu, Qinghuai
2014-01-01
Objective To construct a life-sized eye model using the three-dimensional (3D) printing technology for fundus viewing study of the viewing system. Methods We devised our schematic model eye based on Navarro's eye and redesigned some parameters because of the change of the corneal material and the implantation of intraocular lenses (IOLs). Optical performance of our schematic model eye was compared with Navarro's schematic eye and other two reported physical model eyes using the ZEMAX optical design software. With computer aided design (CAD) software, we designed the 3D digital model of the main structure of the physical model eye, which was used for three-dimensional (3D) printing. Together with the main printed structure, polymethyl methacrylate(PMMA) aspherical cornea, variable iris, and IOLs were assembled to a physical eye model. Angle scale bars were glued from posterior to periphery of the retina. Then we fabricated other three physical models with different states of ammetropia. Optical parameters of these physical eye models were measured to verify the 3D printing accuracy. Results In on-axis calculations, our schematic model eye possessed similar size of spot diagram compared with Navarro's and Bakaraju's model eye, much smaller than Arianpour's model eye. Moreover, the spherical aberration of our schematic eye was much less than other three model eyes. While in off- axis simulation, it possessed a bit higher coma and similar astigmatism, field curvature and distortion. The MTF curves showed that all the model eyes diminished in resolution with increasing field of view, and the diminished tendency of resolution of our physical eye model was similar to the Navarro's eye. The measured parameters of our eye models with different status of ametropia were in line with the theoretical value. Conclusions The schematic eye model we designed can well simulate the optical performance of the human eye, and the fabricated physical one can be used as a tool in fundus
Development of a liquid jet model for implementation in a 3-dimensional Eularian analysis tool
NASA Astrophysics Data System (ADS)
Buschman, Francis X., III
The ability to model the thermal behavior of a nuclear reactor is of utmost importance to the reactor designer. Condensation is an important phenomenon when modeling a reactor system's response to a Loss Of Coolant Accident (LOCA). Condensation is even more important with the use of passive safety systems which rely on condensation heat transfer for long term cooling. The increasing use of condensation heat transfer, including condensation on jets of water, in safety systems puts added pressure to correctly model this phenomenon with thermal-hydraulic system and sub-channel analysis codes. In this work, a stand alone module with which to simulate condensation on a liquid jet was developed and then implemented within a reactor vessel analysis code to improve that code's handling of jet condensation. It is shown that the developed liquid jet model vastly improves the ability of COBRA-TF to model condensation on turbulent liquid jets. The stand alone jet model and the coupled liquid jet COBRA-TF have been compared to experimental data. Jet condensation heat transfer experiments by Celata et al. with a variety of jet diameters, velocities, and subcooling were utilized to evaluate the models. A sensitivity study on the effects of noncondensables on jet condensation was also carried out using the stand alone jet model.
[Rapid 3-Dimensional Models of Cerebral Aneurysm for Emergency Surgical Clipping].
Konno, Takehiko; Mashiko, Toshihiro; Oguma, Hirofumi; Kaneko, Naoki; Otani, Keisuke; Watanabe, Eiju
2016-08-01
We developed a method for manufacturing solid models of cerebral aneurysms, with a shorter printing time than that involved in conventional methods, using a compact 3D printer with acrylonitrile-butadiene-styrene(ABS)resin. We further investigated the application and utility of this printing system in emergency clipping surgery. A total of 16 patients diagnosed with acute subarachnoid hemorrhage resulting from cerebral aneurysm rupture were enrolled in the present study. Emergency clipping was performed on the day of hospitalization. Digital Imaging and Communication in Medicine(DICOM)data obtained from computed tomography angiography(CTA)scans were edited and converted to stereolithography(STL)file formats, followed by the production of 3D models of the cerebral aneurysm by using the 3D printer. The mean time from hospitalization to the commencement of surgery was 242 min, whereas the mean time required for manufacturing the 3D model was 67 min. The average cost of each 3D model was 194 Japanese Yen. The time required for manufacturing the 3D models shortened to approximately 1 hour with increasing experience of producing 3D models. Favorable impressions for the use of the 3D models in clipping were reported by almost all neurosurgeons included in this study. Although 3D printing is often considered to involve huge costs and long manufacturing time, the method used in the present study requires shorter time and lower costs than conventional methods for manufacturing 3D cerebral aneurysm models, thus making it suitable for use in emergency clipping. PMID:27506842
Visualization of the 3-dimensional flow around a model with the aid of a laser knife
NASA Technical Reports Server (NTRS)
Borovoy, V. Y.; Ivanov, V. V.; Orlov, A. A.; Kharchenko, V. N.
1984-01-01
A method for visualizing the three-dimensional flow around models of various shapes in a wind tunnel at a Mach number of 5 is described. A laser provides a planar light flux such that any plane through the model can be selectively illuminated. The shape of shock waves and separation regions is then determined by the intensity of light scattered by soot particles in the flow.
High fidelity 3-dimensional models of beam-electron cloud interactions in circular accelerators
NASA Astrophysics Data System (ADS)
Feiz Zarrin Ghalam, Ali
Electron cloud is a low-density electron profile created inside the vacuum chamber of circular machines with positively charged beams. Electron cloud limits the peak current of the beam and degrades the beams' quality through luminosity degradation, emittance growth and head to tail or bunch to bunch instability. The adverse effects of electron cloud on long-term beam dynamics becomes more and more important as the beams go to higher and higher energies. This problem has become a major concern in many future circular machines design like the Large Hadron Collider (LHC) under construction at European Center for Nuclear Research (CERN). Due to the importance of the problem several simulation models have been developed to model long-term beam-electron cloud interaction. These models are based on "single kick approximation" where the electron cloud is assumed to be concentrated at one thin slab around the ring. While this model is efficient in terms of computational costs, it does not reflect the real physical situation as the forces from electron cloud to the beam are non-linear contrary to this model's assumption. To address the existing codes limitation, in this thesis a new model is developed to continuously model the beam-electron cloud interaction. The code is derived from a 3-D parallel Particle-In-Cell (PIC) model (QuickPIC) originally used for plasma wakefield acceleration research. To make the original model fit into circular machines environment, betatron and synchrotron equations of motions have been added to the code, also the effect of chromaticity, lattice structure have been included. QuickPIC is then benchmarked against one of the codes developed based on single kick approximation (HEAD-TAIL) for the transverse spot size of the beam in CERN-LHC. The growth predicted by QuickPIC is less than the one predicted by HEAD-TAIL. The code is then used to investigate the effect of electron cloud image charges on the long-term beam dynamics, particularly on the
Remanent magnetization and 3-dimensional density model of the Kentucky anomaly region
NASA Technical Reports Server (NTRS)
Mayhew, M. A.; Estes, R. H.; Myers, D. M.
1984-01-01
A three-dimensional model of the Kentucky body was developed to fit surface gravity and long wavelength aeromagnetic data. Magnetization and density parameters for the model are much like those of Mayhew et al (1982). The magnetic anomaly due to the model at satellite altitude is shown to be much too small by itself to account for the anomaly measured by Magsat. It is demonstrated that the source region for the satellite anomaly is considerably more extensive than the Kentucky body sensu stricto. The extended source region is modeled first using prismatic model sources and then using dipole array sources. Magnetization directions for the source region found by inversion of various combinations of scalar and vector data are found to be close to the main field direction, implying the lack of a strong remanent component. It is shown by simulation that in a case (such as this) where the geometry of the source is known, if a strong remanent component is present its direction is readily detectable, but by scalar data as readily as vector data.
NASA Astrophysics Data System (ADS)
Griffioen, J.; van der Grift, B.; Maas, D.; van den Brink, C.; Zaadnoordijk, J. W.
2003-04-01
Groundwater is contaminated at the regional scale by agricultural activities and atmospheric deposition. A 3-D transport model was set-up for a phreatic drinking water winning, where the groundwater composition was monitored accurately. The winning is situated at an area with unconsolidated Pleistocene deposits. The land use is nature and agriculture. Annual mass-balances were determined using a wide range of historic data. The modelling approach for the unsaturated zone was either simple box models (Cl, NO_3 and SO_4) or 1-D transport modelling using HYDRUS (Cd). The modelling approach for the saturated zone used a multiple solute version of MT3D, where denitrification associated with pyrite oxidation and sorption of Cd were included. The solute transport calculations were performed for the period 1950--2030. The results obtained for the year 2000 were used as input concentration for the period 2000--2030. A comparison between the calculated and the measured concentrations of groundwater abstracted for Cl, NO_3 and SO_4 yields the following. First, the input at the surface is rather well estimated. Second, the redox reactivity of the first two aquifers is negligible around the winning, which is confirmed by respiration experiments using anaerobically sampled aquifer sediments. The reactivity of the third aquifer, which is a marine deposit and lies at least 30 meters below surface, is considerable. The discrepancies between modelled and measured output are explained by lack of knowledge about the subsurface reactivity and/or wrong estimates of surface loading and leaching from the unsaturated zone. The patterns for other hydrogeochemical variables such as Ca, HCO_3 may further constrain this lack of knowledge. The results for Cd indicate that Cd becomes strongly retarded, despite the low reactivity of the sandy sediments. The winning is rather insensitive to Cd contamination (but the surface water drainage network is not). Two major uncertainties for input of Cd
Sugimoto, Yoshihisa; Tanaka, Masato; Nakahara, Ryuichi; Misawa, Haruo; Kunisada, Toshiyuki; Ozaki, Toshifumi
2012-01-01
An 11 year-old girl had 66 degrees of kyphosis in the thoracolumbar junction. For the purpose of planning for kyphosis correction, we created a 3-D, full-scale model of the spine and consulted spinal navigation. Three-dimensional models are generally used as tactile guides to verify the surgical approach and portray the anatomic relations specific to a given patient. We performed posterior fusion from Th10 to L3, and vertebral column resection of Th12 and L1. Screw entry points, directions, lengths and diameters were determined by reference to navigation. Both tools were useful in the bone resection. We could easily detect the posterior element to be resected using the 3D model. During the anterior bony resection, navigation helped us to check the disc level and anterior wall of the vertebrae, which were otherwise difficult to detect due to their depth in the surgical field. Thus, the combination of navigation and 3D models helped us to safely perform surgery for a patient with complex spinal deformity. PMID:23254585
2015-01-01
Summary This paper introduces a quasi-3-dimensional (Q3D) viscoelastic model and software tool for use in atomic force microscopy (AFM) simulations. The model is based on a 2-dimensional array of standard linear solid (SLS) model elements. The well-known 1-dimensional SLS model is a textbook example in viscoelastic theory but is relatively new in AFM simulation. It is the simplest model that offers a qualitatively correct description of the most fundamental viscoelastic behaviors, namely stress relaxation and creep. However, this simple model does not reflect the correct curvature in the repulsive portion of the force curve, so its application in the quantitative interpretation of AFM experiments is relatively limited. In the proposed Q3D model the use of an array of SLS elements leads to force curves that have the typical upward curvature in the repulsive region, while still offering a very low computational cost. Furthermore, the use of a multidimensional model allows for the study of AFM tips having non-ideal geometries, which can be extremely useful in practice. Examples of typical force curves are provided for single- and multifrequency tapping-mode imaging, for both of which the force curves exhibit the expected features. Finally, a software tool to simulate amplitude and phase spectroscopy curves is provided, which can be easily modified to implement other controls schemes in order to aid in the interpretation of AFM experiments. PMID:26734515
Solares, Santiago D.
2015-11-26
This study introduces a quasi-3-dimensional (Q3D) viscoelastic model and software tool for use in atomic force microscopy (AFM) simulations. The model is based on a 2-dimensional array of standard linear solid (SLS) model elements. The well-known 1-dimensional SLS model is a textbook example in viscoelastic theory but is relatively new in AFM simulation. It is the simplest model that offers a qualitatively correct description of the most fundamental viscoelastic behaviors, namely stress relaxation and creep. However, this simple model does not reflect the correct curvature in the repulsive portion of the force curve, so its application in the quantitative interpretationmore » of AFM experiments is relatively limited. In the proposed Q3D model the use of an array of SLS elements leads to force curves that have the typical upward curvature in the repulsive region, while still offering a very low computational cost. Furthermore, the use of a multidimensional model allows for the study of AFM tips having non-ideal geometries, which can be extremely useful in practice. Examples of typical force curves are provided for single- and multifrequency tappingmode imaging, for both of which the force curves exhibit the expected features. Lastly, a software tool to simulate amplitude and phase spectroscopy curves is provided, which can be easily modified to implement other controls schemes in order to aid in the interpretation of AFM experiments.« less
Solares, Santiago D.
2015-11-26
This study introduces a quasi-3-dimensional (Q3D) viscoelastic model and software tool for use in atomic force microscopy (AFM) simulations. The model is based on a 2-dimensional array of standard linear solid (SLS) model elements. The well-known 1-dimensional SLS model is a textbook example in viscoelastic theory but is relatively new in AFM simulation. It is the simplest model that offers a qualitatively correct description of the most fundamental viscoelastic behaviors, namely stress relaxation and creep. However, this simple model does not reflect the correct curvature in the repulsive portion of the force curve, so its application in the quantitative interpretation of AFM experiments is relatively limited. In the proposed Q3D model the use of an array of SLS elements leads to force curves that have the typical upward curvature in the repulsive region, while still offering a very low computational cost. Furthermore, the use of a multidimensional model allows for the study of AFM tips having non-ideal geometries, which can be extremely useful in practice. Examples of typical force curves are provided for single- and multifrequency tappingmode imaging, for both of which the force curves exhibit the expected features. Lastly, a software tool to simulate amplitude and phase spectroscopy curves is provided, which can be easily modified to implement other controls schemes in order to aid in the interpretation of AFM experiments.
NASA Astrophysics Data System (ADS)
Garay, M. J.; Diner, D. J.; Martonchik, J. V.; Davis, A. B.
2011-12-01
Knowledge of the detailed 3-dimensional structure of clouds and atmospheric aerosols is vital for correctly modeling their radiative effects and interpreting optical remote sensing measurements of scattered sunlight. We will describe a set of new observations made by the Multiangle SpectroPolarimetric Imager (MSPI) from the ground and from the NASA ER-2 aircraft. MSPI is being developed and tested at JPL as a payload for the preliminary Aerosol-Cloud-Ecosystems (PACE) satellite mission, which is expected to fly near the end of the decade. MSPI builds upon experience gained from the Multi-angle Imaging SpectroRadiometer (MISR) currently orbiting on NASA's Terra satellite. Ground-MSPI and Air-MSPI are two prototype cameras operating in the ultraviolet (UV) to the visible/near-infrared (VNIR) range mounted on gimbals that acquire imagery in a pushbroom fashion, including polarization in selected spectral bands with demonstrated high polarimetric accuracy (0.5% uncertainty in degree of linear polarization). The spatial resolution of Ground-MSPI is 1 m for objects at a distance of 3 km. From the operational altitude of the ER-2, Air-MSPI has a ground resolution of approximately 10 m at nadir. This resolution, coupled with good calibration and high polarimetric performance means that MSPI can be used to derive radiatively important parameters of aerosols and clouds using intensity and polarization information together. As part of the effort for developing retrieval algorithms for the instrument, we have employed an extremely flexible 3-dimensional vector radiative transfer code. We will show example imagery from both MSPI cameras and describe how these scenes are modeled using this code. We will also discuss some of the important unknowns and limitations of this observational approach.
An approximate single fluid 3-dimensional magnetohydrodynamic equilibrium model with toroidal flow
NASA Astrophysics Data System (ADS)
Cooper, W. A.; Hirshman, S. P.; Chapman, I. T.; Brunetti, D.; Faustin, J. M.; Graves, J. P.; Pfefferlé, D.; Raghunathan, M.; Sauter, O.; Tran, T. M.; Aiba, N.
2014-09-01
An approximate model for a single fluid three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium with pure isothermal toroidal flow with imposed nested magnetic flux surfaces is proposed. It recovers the rigorous toroidal rotation equilibrium description in the axisymmetric limit. The approximation is valid under conditions of nearly rigid or vanishing toroidal rotation in regions with significant 3D deformation of the equilibrium flux surfaces. Bifurcated helical core equilibrium simulations of long-lived modes in the MAST device demonstrate that the magnetic structure is only weakly affected by the flow but that the 3D pressure distortion is important. The pressure is displaced away from the major axis and therefore is not as noticeably helically deformed as the toroidal magnetic flux under the subsonic flow conditions measured in the experiment. The model invoked fails to predict any significant screening by toroidal plasma rotation of resonant magnetic perturbations in MAST free boundary computations.
3-DIMENSIONAL Geological Mapping and Modeling Activities at the Geological Survey of Norway
NASA Astrophysics Data System (ADS)
Jarna, A.; Bang-Kittilsen, A.; Haase, C.; Henderson, I. H. C.; Høgaas, F.; Iversen, S.; Seither, A.
2015-10-01
Geology and all geological structures are three-dimensional in space. Geology can be easily shown as four-dimensional when time is considered. Therefore GIS, databases, and 3D visualization software are common tools used by geoscientists to view, analyse, create models, interpret and communicate geological data. The NGU (Geological Survey of Norway) is the national institution for the study of bedrock, mineral resources, surficial deposits and groundwater and marine geology. The interest in 3D mapping and modelling has been reflected by the increase of number of groups and researches dealing with 3D in geology within NGU. This paper highlights 3D geological modelling techniques and the usage of these tools in bedrock, geophysics, urban and groundwater studies at NGU, same as visualisation of 3D online. The examples show use of a wide range of data, methods, software and an increased focus on interpretation and communication of geology in 3D. The goal is to gradually expand the geospatial data infrastructure to include 3D data at the same level as 2D.
Pazera, Pawel; Zorkun, Berna; Katsaros, Christos; Ludwig, Björn
2015-01-01
Objectives To test the applicability, accuracy, precision, and reproducibility of various 3D superimposition techniques for radiographic data, transformed to triangulated surface data. Methods Five superimposition techniques (3P: three-point registration; AC: anterior cranial base; AC + F: anterior cranial base + foramen magnum; BZ: both zygomatic arches; 1Z: one zygomatic arch) were tested using eight pairs of pre-existing CT data (pre- and post-treatment). These were obtained from non-growing orthodontic patients treated with rapid maxillary expansion. All datasets were superimposed by three operators independently, who repeated the whole procedure one month later. Accuracy was assessed by the distance (D) between superimposed datasets on three form-stable anatomical areas, located on the anterior cranial base and the foramen magnum. Precision and reproducibility were assessed using the distances between models at four specific landmarks. Non parametric multivariate models and Bland-Altman difference plots were used for analyses. Results There was no difference among operators or between time points on the accuracy of each superimposition technique (p>0.05). The AC + F technique was the most accurate (D<0.17 mm), as expected, followed by AC and BZ superimpositions that presented similar level of accuracy (D<0.5 mm). 3P and 1Z were the least accurate superimpositions (0.79
A 3-dimensional in vitro model of epithelioid granulomas induced by high aspect ratio nanomaterials
2011-01-01
Background The most common causes of granulomatous inflammation are persistent pathogens and poorly-degradable irritating materials. A characteristic pathological reaction to intratracheal instillation, pharyngeal aspiration, or inhalation of carbon nanotubes is formation of epithelioid granulomas accompanied by interstitial fibrosis in the lungs. In the mesothelium, a similar response is induced by high aspect ratio nanomaterials, including asbestos fibers, following intraperitoneal injection. This asbestos-like behaviour of some engineered nanomaterials is a concern for their potential adverse health effects in the lungs and mesothelium. We hypothesize that high aspect ratio nanomaterials will induce epithelioid granulomas in nonadherent macrophages in 3D cultures. Results Carbon black particles (Printex 90) and crocidolite asbestos fibers were used as well-characterized reference materials and compared with three commercial samples of multiwalled carbon nanotubes (MWCNTs). Doses were identified in 2D and 3D cultures in order to minimize acute toxicity and to reflect realistic occupational exposures in humans and in previous inhalation studies in rodents. Under serum-free conditions, exposure of nonadherent primary murine bone marrow-derived macrophages to 0.5 μg/ml (0.38 μg/cm2) of crocidolite asbestos fibers or MWCNTs, but not carbon black, induced macrophage differentiation into epithelioid cells and formation of stable aggregates with the characteristic morphology of granulomas. Formation of multinucleated giant cells was also induced by asbestos fibers or MWCNTs in this 3D in vitro model. After 7-14 days, macrophages exposed to high aspect ratio nanomaterials co-expressed proinflammatory (M1) as well as profibrotic (M2) phenotypic markers. Conclusions Induction of epithelioid granulomas appears to correlate with high aspect ratio and complex 3D structure of carbon nanotubes, not with their iron content or surface area. This model offers a time- and cost
NASA Astrophysics Data System (ADS)
Kobayashi, H.; Yang, W.; Ichii, K.
2015-12-01
Global simulation of canopy scale sun-induced chlorophyll fluorescence with a 3 dimensional radiative transfer modelHideki Kobayashi, Wei Yang, and Kazuhito IchiiDepartment of Environmental Geochemical Cycle Research, Japan Agency for Marine-Earth Science and Technology3173-25, Showa-machi, Kanazawa-ku, Yokohama, Japan.Plant canopy scale sun-induced chlorophyll fluorescence (SIF) can be observed from satellites, such as Greenhouse gases Observation Satellite (GOSAT), Orbiting Carbon Observatory-2 (OCO-2), and Global Ozone Monitoring Experiment-2 (GOME-2), using Fraunhofer lines in the near infrared spectral domain [1]. SIF is used to infer photosynthetic capacity of plant canopy [2]. However, it is not well understoond how the leaf-level SIF emission contributes to the top of canopy directional SIF because SIFs observed by the satellites use the near infrared spectral domain where the multiple scatterings among leaves are not negligible. It is necessary to quantify the fraction of emission for each satellite observation angle. Absorbed photosynthetically active radiation of sunlit leaves are 100 times higher than that of shaded leaves. Thus, contribution of sunlit and shaded leaves to canopy scale directional SIF emission should also be quantified. Here, we show the results of global simulation of SIF using a 3 dimensional radiative transfer simulation with MODIS atmospheric (aerosol optical thickness) and land (land cover and leaf area index) products and a forest landscape data sets prepared for each land cover category. The results are compared with satellite-based SIF (e.g. GOME-2) and the gross primary production empirically estimated by FLUXNET and remote sensing data.
Numerical model of electromagnetic scattering off a subterranean 3-dimensional dielectric
Dease, C.G.; Didwall, E.M.
1983-08-01
As part of the effort to develop On-Site Inspection (OSI) techniques for verification of compliance to a Comprehensive Test Ban Treaty (CTBT), a computer code was developed to predict the interaction of an electromagnetic (EM) wave with an underground cavity. Results from the code were used to evaluate the use of surface electromagnetic exploration techniques for detection of underground cavities or rubble-filled regions characteristic of underground nuclear explosions.
Bhave, Madhura Satish; Hassanbhai, Ammar Mansoor; Anand, Padmaja; Luo, Kathy Qian; Teoh, Swee Hin
2015-01-01
Traditional cancer treatments, such as chemotherapy and radiation therapy continue to have limited efficacy due to tumor hypoxia. While bacterial cancer therapy has the potential to overcome this problem, it comes with the risk of toxicity and infection. To circumvent these issues, this paper investigates the anti-tumor effects of non-viable bacterial derivatives of Clostridium sporogenes. These non-viable derivatives are heat-inactivated C. sporogenes bacteria (IB) and the secreted bacterial proteins in culture media, known as conditioned media (CM). In this project, the effects of IB and CM on CT26 and HCT116 colorectal cancer cells were examined on a 2-Dimensional (2D) and 3-Dimensional (3D) platform. IB significantly inhibited cell proliferation of CT26 to 6.3% of the control in 72 hours for the 2D monolayer culture. In the 3D spheroid culture, cell proliferation of HCT116 spheroids notably dropped to 26.2%. Similarly the CM also remarkably reduced the cell-proliferation of the CT26 cells to 2.4% and 20% in the 2D and 3D models, respectively. Interestingly the effect of boiled conditioned media (BCM) on the cells in the 3D model was less inhibitory than that of CM. Thus, the inhibitive effect of inactivated C. sporogenes and its conditioned media on colorectal cancer cells is established. PMID:26507312
Bhave, Madhura Satish; Hassanbhai, Ammar Mansoor; Anand, Padmaja; Luo, Kathy Qian; Teoh, Swee Hin
2015-01-01
Traditional cancer treatments, such as chemotherapy and radiation therapy continue to have limited efficacy due to tumor hypoxia. While bacterial cancer therapy has the potential to overcome this problem, it comes with the risk of toxicity and infection. To circumvent these issues, this paper investigates the anti-tumor effects of non-viable bacterial derivatives of Clostridium sporogenes. These non-viable derivatives are heat-inactivated C. sporogenes bacteria (IB) and the secreted bacterial proteins in culture media, known as conditioned media (CM). In this project, the effects of IB and CM on CT26 and HCT116 colorectal cancer cells were examined on a 2-Dimensional (2D) and 3-Dimensional (3D) platform. IB significantly inhibited cell proliferation of CT26 to 6.3% of the control in 72 hours for the 2D monolayer culture. In the 3D spheroid culture, cell proliferation of HCT116 spheroids notably dropped to 26.2%. Similarly the CM also remarkably reduced the cell-proliferation of the CT26 cells to 2.4% and 20% in the 2D and 3D models, respectively. Interestingly the effect of boiled conditioned media (BCM) on the cells in the 3D model was less inhibitory than that of CM. Thus, the inhibitive effect of inactivated C. sporogenes and its conditioned media on colorectal cancer cells is established. PMID:26507312
This report presents a three-dimensional finite-element numerical model designed to simulate chemical transport in subsurface systems with temperature effect taken into account. The three-dimensional model is developed to provide (1) a tool of application, with which one is able...
NASA Astrophysics Data System (ADS)
Lu, Junqing; Keiter, Eric R.; Kushner, Mark J.
1998-10-01
Inductively Coupled Plasmas (ICPs) are being used for a variety of deposition processes for microelectronics fabrication. Of particular concern in scaling these devices to large areas is maintaining azimuthal symmetry of the reactant fluxes. Sources of nonuniformity may be physical (e.g., gas injection and side pumping) or electromagnetic (e.g., transmission line effects in the antennas). In this paper, a 3-dimensional plasma equipment model, HPEM-3D,(M. J. Kushner, J. Appl. Phys. v.82, 5312 (1997).) is used to investigate physical and electromagentic sources of azimuthal nonuniformities in deposition tools. An ionized metal physical vapor deposition (IMPVD) system will be investigated where transmission line effects in the coils produce an asymmetric plasma density. Long mean free path transport for sputtered neutrals and tensor conducitivities have been added to HPEM-3D to address this system. Since the coil generated ion flux drifts back to the target to sputter low ionization potential metal atoms, the asymmetry is reinforced by rapid ionization of the metal atoms.
Numerical modeling of nonintrusive inspection systems
Hall, J.; Morgan, J.; Sale, K.
1992-12-01
A wide variety of nonintrusive inspection systems have been proposed in the past several years for the detection of hidden contraband in airline luggage and shipping containers. The majority of these proposed techniques depend on the interaction of radiation with matter to produce a signature specific to the contraband of interest, whether drugs or explosives. In the authors` role as diagnostic specialists in the Underground Test Program over the past forty years, L-Division of the Lawrence Livermore National Laboratory has developed a technique expertise in the combined numerical and experimental modeling of these types of system. Based on their experience, they are convinced that detailed numerical modeling provides a much more accurate estimate of the actual performance of complex experiments than simple analytical modeling. Furthermore, the construction of detailed numerical prototypes allows experimenters to explore the entire region of parameter space available to them before committing their ideas to hardware. This sort of systematic analysis has often led to improved experimental designs and reductions in fielding costs. L-Division has developed an extensive suite of computer codes to model proposed experiments and possible background interactions. These codes allow one to simulate complex radiation sources, model 3-dimensional system geometries with {open_quotes}real world{close_quotes} complexity, specify detailed elemental distributions, and predict the response of almost any type of detector. In this work several examples are presented illustrating the use of these codes in modeling experimental systems at LLNL and their potential usefulness in evaluating nonintrusive inspection systems is discussed.
Pashazadeh, Saeid; Sharifi, Mohsen
2009-01-01
Existing 3-dimensional acoustic target tracking methods that use wired/wireless networked sensor nodes to track targets based on four sensing coverage do not always compute the feasible spatio-temporal information of target objects. To investigate this discrepancy in a formal setting, we propose a geometric model of the target tracking problem alongside its equivalent geometric dual model that is easier to solve. We then study and prove some properties of dual model by exploiting its relationship with algebra. Based on these properties, we propose a four coverage axis line method based on four sensing coverage and prove that four sensing coverage always yields two dual correct answers; usually one of them is infeasible. By showing that the feasible answer can be only sometimes identified by using a simple time test method such as the one proposed by ourselves, we prove that four sensing coverage fails to always yield the feasible spatio-temporal information of a target object. We further prove that five sensing coverage always gives the feasible position of a target object under certain conditions that are discussed in this paper. We propose three extensions to four coverage axis line method, namely, five coverage extent point method, five coverage extended axis lines method, and five coverage redundant axis lines method. Computation and time complexities of all four proposed methods are equal in the worst cases as well as on average being equal to Θ(1) each. Proposed methods and proved facts about capabilities of sensing coverage degree in this paper can be used in all other methods of acoustic target tracking like Bayesian filtering methods. PMID:22423198
A 3-Dimensional Model of Water-Bearing Sequences in the Dominguez Gap Region, Long Beach, California
Ponti, Daniel J.; Ehman, Kenneth D.; Edwards, Brian D.; Tinsley, John C., III; Hildenbrand, Thomas; Hillhouse, John W.; Hanson, Randall T.; McDougall, Kristen; Powell, Charles L.; Wan, Elmira; Land, Michael; Mahan, Shannon; Sarna-Wojcicki, Andrei M.
2007-01-01
A 3-dimensional computer model of the Quaternary sequence stratigraphy in the Dominguez gap region of Long Beach, California has been developed to provide a robust chronostratigraphic framework for hydrologic and tectonic studies. The model consists of 13 layers within a 16.5 by 16.1 km (10.25 by 10 mile) square area and extends downward to an altitude of -900 meters (-2952.76 feet). Ten sequences of late Pliocene to Holocene age are identified and correlated within the model. Primary data to build the model comes from five reference core holes, extensive high-resolution seismic data obtained in San Pedro Bay, and logs from several hundred water and oil wells drilled in the region. The model is best constrained in the vicinity of the Dominguez gap seawater intrusion barrier where a dense network of subsurface data exist. The resultant stratigraphic framework and geologic structure differs significantly from what has been proposed in earlier studies. An important new discovery from this approach is the recognition of ongoing tectonic deformation throughout nearly all of Quaternary time that has impacted the geometry and character of the sequences. Anticlinal folding along a NW-SE trend, probably associated with Quaternary reactivation of the Wilmington anticline, has uplifted and thinned deposits along the fold crest, which intersects the Dominguez gap seawater barrier near Pacific Coast Highway. A W-NW trending fault system that approximately parallels the fold crest has also been identified. This fault progressively displaces all but the youngest sequences down to the north and serves as the southern termination of the classic Silverado aquifer. Uplift and erosion of fining-upward paralic sequences along the crest of the young fold has removed or thinned many of the fine-grained beds that serve to protect the underlying Silverado aquifer from seawater contaminated shallow groundwater. As a result of this process, the potential exists for vertical migration of
Prideaux, Andrew R.; Song, Hong; Hobbs, Robert F.; He, Bin; Frey, Eric C.; Ladenson, Paul W.; Wahl, Richard L.; Sgouros, George
2010-01-01
Phantom-based and patient-specific imaging-based dosimetry methodologies have traditionally yielded mean organ-absorbed doses or spatial dose distributions over tumors and normal organs. In this work, radiobiologic modeling is introduced to convert the spatial distribution of absorbed dose into biologically effective dose and equivalent uniform dose parameters. The methodology is illustrated using data from a thyroid cancer patient treated with radioiodine. Methods Three registered SPECT/CT scans were used to generate 3-dimensional images of radionuclide kinetics (clearance rate) and cumulated activity. The cumulated activity image and corresponding CT scan were provided as input into an EGSnrc-based Monte Carlo calculation: The cumulated activity image was used to define the distribution of decays, and an attenuation image derived from CT was used to define the corresponding spatial tissue density and composition distribution. The rate images were used to convert the spatial absorbed dose distribution to a biologically effective dose distribution, which was then used to estimate a single equivalent uniform dose for segmented volumes of interest. Equivalent uniform dose was also calculated from the absorbed dose distribution directly. Results We validate the method using simple models; compare the dose-volume histogram with a previously analyzed clinical case; and give the mean absorbed dose, mean biologically effective dose, and equivalent uniform dose for an illustrative case of a pediatric thyroid cancer patient with diffuse lung metastases. The mean absorbed dose, mean biologically effective dose, and equivalent uniform dose for the tumor were 57.7, 58.5, and 25.0 Gy, respectively. Corresponding values for normal lung tissue were 9.5, 9.8, and 8.3 Gy, respectively. Conclusion The analysis demonstrates the impact of radiobiologic modeling on response prediction. The 57% reduction in the equivalent dose value for the tumor reflects a high level of dose
Chrysostomou, P P; Lodish, M B; Turkbey, E B; Papadakis, G Z; Stratakis, C A
2016-04-01
Primary pigmented nodular adrenocortical disease (PPNAD) is a rare type of bilateral adrenal hyperplasia leading to hypercortisolemia. Adrenal nodularity is often appreciable with computed tomography (CT); however, accurate radiologic characterization of adrenal size in PPNAD has not been studied well. We used 3-dimensional (3D) volumetric analysis to characterize and compare adrenal size in PPNAD patients, with and without Cushing's syndrome (CS). Patients diagnosed with PPNAD and their family members with known mutations in PRKAR1A were screened. CT scans were used to create 3D models of each adrenal. Criteria for biochemical diagnosis of CS included loss of diurnal variation and/or elevated midnight cortisol levels, and paradoxical increase in urinary free cortisol and/or urinary 17-hydroxysteroids after dexamethasone administration. Forty-five patients with PPNAD (24 females, 27.8±17.6 years) and 8 controls (19±3 years) were evaluated. 3D volumetric modeling of adrenal glands was performed in all. Thirty-eight patients out of 45 (84.4%) had CS. Their mean adrenal volume was 8.1 cc±4.1, 7.2 cc±4.5 (p=0.643) for non-CS, and 8.0cc±1.6 for controls. Mean values were corrected for body surface area; 4.7 cc/kg/m(2)±2.2 for CS, and 3.9 cc/kg/m(2)±1.3 for non-CS (p=0.189). Adrenal volume and midnight cortisol in both groups was positively correlated, r=0.35, p=0.03. We conclude that adrenal volume measured by 3D CT in patients with PPNAD and CS was similar to those without CS, confirming empirical CT imaging-based observations. However, the association between adrenal volume and midnight cortisol levels may be used as a marker of who among patients with PPNAD may develop CS, something that routine CT cannot do. PMID:27065461
The use of TOUGH2 for the LBL/USGS 3-dimensional site-scale model of Yucca Mountain, Nevada
Bodvarsson, G.; Chen, G.; Haukwa, C.; Kwicklis, E.
1995-12-31
The three-dimensional site-scale numerical model o the unsaturated zone at Yucca Mountain is under continuous development and calibration through a collaborative effort between Lawrence Berkeley Laboratory (LBL) and the United States Geological Survey (USGS). The site-scale model covers an area of about 30 km{sup 2} and is bounded by major fault zones to the west (Solitario Canyon Fault), east (Bow Ridge Fault) and perhaps to the north by an unconfirmed fault (Yucca Wash Fault). The model consists of about 5,000 grid blocks (elements) with nearly 20,000 connections between them; the grid was designed to represent the most prevalent geological and hydro-geological features of the site including major faults, and layering and bedding of the hydro-geological units. Submodels are used to investigate specific hypotheses and their importance before incorporation into the three-dimensional site-scale model. The primary objectives of the three-dimensional site-scale model are to: (1) quantify moisture, gas and heat flows in the ambient conditions at Yucca Mountain, (2) help in guiding the site-characterization effort (primarily by USGS) in terms of additional data needs and to identify regions of the mountain where sufficient data have been collected, and (3) provide a reliable model of Yucca Mountain that is validated by repeated predictions of conditions in new boreboles and the ESF and has therefore the confidence of the public and scientific community. The computer code TOUGH2 developed by K. Pruess at LBL was used along with the three-dimensional site-scale model to generate these results. In this paper, we also describe the three-dimensional site-scale model emphasizing the numerical grid development, and then show some results in terms of moisture, gas and heat flow.
Numerical modeling of Hall thruster
Chable, S.; Rogier, F.
2005-05-16
A stationary plasma thruster is numerically studied using different levels. An one dimensional modeling is first analyzed and compared with experimental results. A simplified model of oscillations thruster is proposed and used to control the amplitude of oscillations. A two dimensional numerical method is discussed and applied to the computation of the flow in the exhaust.
FERRARIO, VIRGILIO F.; SFORZA, CHIARELLA; SCHMITZ, JOHANNES H.; CIUSA, VERONICA; COLOMBO, ANNA
2000-01-01
A 3-dimensional computerised system with landmark representation of the soft-tissue facial surface allows noninvasive and fast quantitative study of facial growth. The aims of the present investigation were (1) to provide reference data for selected dimensions of lips (linear distances and ratios, vermilion area, volume); (2) to quantify the relevant growth changes; and (3) to evaluate sex differences in growth patterns. The 3-dimensional coordinates of 6 soft-tissue landmarks on the lips were obtained by an optoelectronic instrument in a mixed longitudinal and cross-sectional study (2023 examinations in 1348 healthy subjects between 6 y of age and young adulthood). From the landmarks, several linear distances (mouth width, total vermilion height, total lip height, upper lip height), the vermilion height-to-mouth width ratio, some areas (vermilion of the upper lip, vermilion of the lower lip, total vermilion) and volumes (upper lip volume, lower lip volume, total lip volume) were calculated and averaged for age and sex. Male values were compared with female values by means of Student's t test. Within each age group all lip dimensions (distances, areas, volumes) were significantly larger in boys than in girls (P < 0.05), with some exceptions in the first age groups and coinciding with the earlier female growth spurt, whereas the vermilion height-to-mouth width ratio did not show a corresponding sexual dimorphism. Linear distances in girls had almost reached adult dimensions in the 13–14 y age group, while in boys a large increase was still to occur. The attainment of adult dimensions was faster in the upper than in the lower lip, especially in girls. The method used in the present investigation allowed the noninvasive evaluation of a large sample of nonpatient subjects, leading to the definition of 3-dimensional normative data. Data collected in the present study could represent a data base for the quantitative description of human lip morphology from childhood to
NASA Technical Reports Server (NTRS)
Lee, S. S.; Sengupta, S.; Nwadike, E. V.; Sinha, S. K.
1982-01-01
The six-volume report: describes the theory of a three dimensional (3-D) mathematical thermal discharge model and a related one dimensional (1-D) model, includes model verification at two sites, and provides a separate user's manual for each model. The 3-D model has two forms: free surface and rigid lid. The former, verified at Anclote Anchorage (FL), allows a free air/water interface and is suited for significant surface wave heights compared to mean water depth; e.g., estuaries and coastal regions. The latter, verified at Lake Keowee (SC), is suited for small surface wave heights compared to depth (e.g., natural or man-made inland lakes) because surface elevation has been removed as a parameter. These models allow computation of time-dependent velocity and temperature fields for given initial conditions and time-varying boundary conditions. The free-surface model also provides surface height variations with time.
NASA Technical Reports Server (NTRS)
Lee, S. S.; Sengupta, S.; Tuann, S. Y.; Lee, C. R.
1982-01-01
The six-volume report: describes the theory of a three-dimensional (3-D) mathematical thermal discharge model and a related one-dimensional (1-D) model, includes model verification at two sites, and provides a separate user's manual for each model. The 3-D model has two forms: free surface and rigid lid. The former, verified at Anclote Anchorage (FL), allows a free air/water interface and is suited for significant surface wave heights compared to mean water depth; e.g., estuaries and coastal regions. The latter, verified at Lake Keowee (SC), is suited for small surface wave heights compared to depth. These models allow computation of time-dependent velocity and temperature fields for given initial conditions and time-varying boundary conditions.
Numerical models of galactic dynamos
NASA Astrophysics Data System (ADS)
Elstner, Detlef
The state of the art for dynamo models in spiral galaxies is reviewed. The comparison of numerical models with special properties of observed magnetic fields yields constraints for the turbulent diffusivity and the α-effect. The derivation of the turbulence parameters from the vertical structure of the interstellar medium gives quite reasonable values for modelling the regular magnetic fields in galaxies with an α2Ω-dynamo. Considering the differences of the turbulence between spiral arms and interarm regions, the observed interarm magnetic fields are recovered in the numerical models due to the special properties of the α2Ω-dynamo.
NASA Astrophysics Data System (ADS)
Király, E.; Székely, B.; Bata, T.; Lócsi, L.; Karátson, D.
2009-04-01
The almost global availability of medium- and high-resolution Digital Terrain Models (DTMs) paved the way of new approaches in volcanic geomorphology. The increasing importance of understanding of surface processes that act during the degradation of volcanic edifices also mean a demand for geometric modeling of their surface, in order to derive parameters from the topography that are suitable for further analysis. Our study area, the San Francisco Volcanic Field (SFVF), is a ca. 4500 km2-large volcanic region situated around the San Francisco stratovolcano at Flagstaff, Arizona (USA) that hosts some 600 scoria and lava domes, numerous lava flows with extensive volcanic ash deposits. Because of the wide range in size and age, as well as contrasting degradation of these volcanic features, several authors have analysed them in the last decades to derive general rules of their lowering. Morphometric parameters were determined that were expected to be suitable to fulfill this requirement. In his pioneering work, Wood (1980a,b) considered 40 scoria cones, while almost two decades later Hooper and Sheridan (1998) included 237 features in their study. Their manual morphometric analyses were based on topographic maps that are time consuming, therefore their limited scope can now be extended with the availability of digital data. In the initial phase of our project more than 300 cones were analysed using the classic approach (height of the cone, width of the cone and crater, etc.). Additionally the slope histogram were analysed in order to classify the cones into different evolutionary categories. These analyses led to the selection of a few volcanoes, that entered in the next processing phase. Firstly the derivation of parameters in two-dimensional approach were carried out. Horizontal and vertical cross sections were extracted from the DTM, and the resulting planar curves were analysed via parameter estimation. The horizontal planar outlines were approached with circles
Numerical Modelling of Gelating Aerosols
Babovsky, Hans
2008-09-01
The numerical simulation of the gel phase transition of an aerosol system is an interesting and demanding task. Here, we follow an approach first discussed in [6, 8] which turns out as a useful numerical tool. We investigate several improvements and generalizations. In the center of interest are coagulation diffusion systems, where the aerosol dynamics is supplemented with diffusive spreading in physical space. This leads to a variety of scenarios (depending on the coagulation kernel and the diffusion model) for the spatial evolution of the gelation area.
Adaptive Numerical Algorithms in Space Weather Modeling
NASA Technical Reports Server (NTRS)
Toth, Gabor; vanderHolst, Bart; Sokolov, Igor V.; DeZeeuw, Darren; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Nakib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav
2010-01-01
Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical
Adaptive numerical algorithms in space weather modeling
NASA Astrophysics Data System (ADS)
Tóth, Gábor; van der Holst, Bart; Sokolov, Igor V.; De Zeeuw, Darren L.; Gombosi, Tamas I.; Fang, Fang; Manchester, Ward B.; Meng, Xing; Najib, Dalal; Powell, Kenneth G.; Stout, Quentin F.; Glocer, Alex; Ma, Ying-Juan; Opher, Merav
2012-02-01
Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit
A 3-dimensional ray-trace model for predicting the performance of flashlamp-pumped laser amplifiers
Jancaitis, K.S.; Haney, S.W.; Munro, D.H.; Le Touze, G.; Cabourdin, O.
1997-02-13
We have developed a fully three-dimensional model for the performance of flashlamp pumped laser amplifiers. The model uses a reverse ray-trace technique to calculate the pumping of the laser glass by the flashlamp radiation. We have discovered several different methods by which we can speed up the calculation of the gain profile in a amplifier. The model predicts the energy-storage performance of the Beamlet amplifiers to better than 5%. This model will be used in the optimization of the National Ignition Facility (NIF) amplifier design.
Lambros, Maria P.; DeSalvo, Michael K.; Moreno, Jonathan; Mulamalla, Hari Chandana; Kondapalli, Lavanya
2015-01-01
Cancer patients who receive radiation are often afflicted by oral mucositis, a debilitating disease, characterized by mouth sores and difficulty in swallowing. Oftentimes, cancer patients afflicted with mucositis must stop life-saving therapies. Thus it is very important to prevent mucositis before it develops. Using a validated organotypic model of human oral mucosa, a 3-dimensional cell culture model of human oral keratinocytes, it has been shown that a mixture (NAC–QYD) of N-acetyl cysteine (NAC) and a traditional Chinese medicine, Qingre Liyan decoction (QYD), prevented radiation damage (Lambros et al., 2014). Here we provide detailed methods and analysis of microarray data for non-irradiated and irradiated human oral mucosal tissue with and without pretreatment with NAC, QYD and NAC-QYD. The microarray data been deposited in Gene Expression Omnibus (GEO): GSE62397. These data can be used to further elucidate the mechanisms of irradiation damage in oral mucosa and its prevention. PMID:26697327
Lambros, Maria P; DeSalvo, Michael K; Moreno, Jonathan; Mulamalla, Hari Chandana; Kondapalli, Lavanya
2015-12-01
Cancer patients who receive radiation are often afflicted by oral mucositis, a debilitating disease, characterized by mouth sores and difficulty in swallowing. Oftentimes, cancer patients afflicted with mucositis must stop life-saving therapies. Thus it is very important to prevent mucositis before it develops. Using a validated organotypic model of human oral mucosa, a 3-dimensional cell culture model of human oral keratinocytes, it has been shown that a mixture (NAC-QYD) of N-acetyl cysteine (NAC) and a traditional Chinese medicine, Qingre Liyan decoction (QYD), prevented radiation damage (Lambros et al., 2014). Here we provide detailed methods and analysis of microarray data for non-irradiated and irradiated human oral mucosal tissue with and without pretreatment with NAC, QYD and NAC-QYD. The microarray data been deposited in Gene Expression Omnibus (GEO): GSE62397. These data can be used to further elucidate the mechanisms of irradiation damage in oral mucosa and its prevention. PMID:26697327
NASA Technical Reports Server (NTRS)
Gibson, S. G.
1983-01-01
A system of computer programs was developed to model general three dimensional surfaces. Surfaces are modeled as sets of parametric bicubic patches. There are also capabilities to transform coordinates, to compute mesh/surface intersection normals, and to format input data for a transonic potential flow analysis. A graphical display of surface models and intersection normals is available. There are additional capabilities to regulate point spacing on input curves and to compute surface/surface intersection curves. Input and output data formats are described; detailed suggestions are given for user input. Instructions for execution are given, and examples are shown.
Mirsch, Johanna; Tommasino, Francesco; Frohns, Antonia; Conrad, Sandro; Durante, Marco; Scholz, Michael; Friedrich, Thomas; Löbrich, Markus
2015-10-01
Charged particles are increasingly used in cancer radiotherapy and contribute significantly to the natural radiation risk. The difference in the biological effects of high-energy charged particles compared with X-rays or γ-rays is determined largely by the spatial distribution of their energy deposition events. Part of the energy is deposited in a densely ionizing manner in the inner part of the track, with the remainder spread out more sparsely over the outer track region. Our knowledge about the dose distribution is derived solely from modeling approaches and physical measurements in inorganic material. Here we exploited the exceptional sensitivity of γH2AX foci technology and quantified the spatial distribution of DNA lesions induced by charged particles in a mouse model tissue. We observed that charged particles damage tissue nonhomogenously, with single cells receiving high doses and many other cells exposed to isolated damage resulting from high-energy secondary electrons. Using calibration experiments, we transformed the 3D lesion distribution into a dose distribution and compared it with predictions from modeling approaches. We obtained a radial dose distribution with sub-micrometer resolution that decreased with increasing distance to the particle path following a 1/r2 dependency. The analysis further revealed the existence of a background dose at larger distances from the particle path arising from overlapping dose deposition events from independent particles. Our study provides, to our knowledge, the first quantification of the spatial dose distribution of charged particles in biologically relevant material, and will serve as a benchmark for biophysical models that predict the biological effects of these particles. PMID:26392532
Mirsch, Johanna; Tommasino, Francesco; Frohns, Antonia; Conrad, Sandro; Durante, Marco; Scholz, Michael; Friedrich, Thomas; Löbrich, Markus
2015-01-01
Charged particles are increasingly used in cancer radiotherapy and contribute significantly to the natural radiation risk. The difference in the biological effects of high-energy charged particles compared with X-rays or γ-rays is determined largely by the spatial distribution of their energy deposition events. Part of the energy is deposited in a densely ionizing manner in the inner part of the track, with the remainder spread out more sparsely over the outer track region. Our knowledge about the dose distribution is derived solely from modeling approaches and physical measurements in inorganic material. Here we exploited the exceptional sensitivity of γH2AX foci technology and quantified the spatial distribution of DNA lesions induced by charged particles in a mouse model tissue. We observed that charged particles damage tissue nonhomogenously, with single cells receiving high doses and many other cells exposed to isolated damage resulting from high-energy secondary electrons. Using calibration experiments, we transformed the 3D lesion distribution into a dose distribution and compared it with predictions from modeling approaches. We obtained a radial dose distribution with sub-micrometer resolution that decreased with increasing distance to the particle path following a 1/r2 dependency. The analysis further revealed the existence of a background dose at larger distances from the particle path arising from overlapping dose deposition events from independent particles. Our study provides, to our knowledge, the first quantification of the spatial dose distribution of charged particles in biologically relevant material, and will serve as a benchmark for biophysical models that predict the biological effects of these particles. PMID:26392532
NASA Astrophysics Data System (ADS)
Bahlake, Ahmad; Farivar, Foad; Dabir, Bahram
2016-07-01
In this paper a 3-dimensional modeling of simultaneous stripping of carbon dioxide (CO2) and hydrogen sulfide (H2S) from water using hollow fiber membrane made of polyvinylidene fluoride is developed. The water, containing CO2 and H2S enters to the membrane as feed. At the same time, pure nitrogen flow in the shell side of a shell and tube hollow fiber as the solvent. In the previous methods of modeling hollow fiber membranes just one of the membranes was modeled and the results expand to whole shell and tube system. In this research the whole hollow fiber shell and tube module is modeled to reduce the errors. Simulation results showed that increasing the velocity of solvent flow and decreasing the velocity of the feed are leads to increase in the system yield. However the effect of the feed velocity on the process is likely more than the influence of changing the velocity of the gaseous solvent. In addition H2S stripping has higher yield in comparison with CO2 stripping. This model is compared to the previous modeling methods and shows that the new model is more accurate. Finally, the effect of feed temperature is studied using response surface method and the operating conditions of feed temperature, feed velocity, and solvent velocity is optimized according to synergistic effects. Simulation results show that, in the optimum operating conditions the removal percentage of H2S and CO2 are 27 and 21 % respectively.
NASA Astrophysics Data System (ADS)
Bahlake, Ahmad; Farivar, Foad; Dabir, Bahram
2015-08-01
In this paper a 3-dimensional modeling of simultaneous stripping of carbon dioxide (CO2) and hydrogen sulfide (H2S) from water using hollow fiber membrane made of polyvinylidene fluoride is developed. The water, containing CO2 and H2S enters to the membrane as feed. At the same time, pure nitrogen flow in the shell side of a shell and tube hollow fiber as the solvent. In the previous methods of modeling hollow fiber membranes just one of the membranes was modeled and the results expand to whole shell and tube system. In this research the whole hollow fiber shell and tube module is modeled to reduce the errors. Simulation results showed that increasing the velocity of solvent flow and decreasing the velocity of the feed are leads to increase in the system yield. However the effect of the feed velocity on the process is likely more than the influence of changing the velocity of the gaseous solvent. In addition H2S stripping has higher yield in comparison with CO2 stripping. This model is compared to the previous modeling methods and shows that the new model is more accurate. Finally, the effect of feed temperature is studied using response surface method and the operating conditions of feed temperature, feed velocity, and solvent velocity is optimized according to synergistic effects. Simulation results show that, in the optimum operating conditions the removal percentage of H2S and CO2 are 27 and 21 % respectively.
Comprehensive numerical modelling of tokamaks
Cohen, R.H.; Cohen, B.I.; Dubois, P.F.
1991-01-03
We outline a plan for the development of a comprehensive numerical model of tokamaks. The model would consist of a suite of independent, communicating packages describing the various aspects of tokamak performance (core and edge transport coefficients and profiles, heating, fueling, magnetic configuration, etc.) as well as extensive diagnostics. These codes, which may run on different computers, would be flexibly linked by a user-friendly shell which would allow run-time specification of packages and generation of pre- and post-processing functions, including workstation-based visualization of output. One package in particular, the calculation of core transport coefficients via gyrokinetic particle simulation, will become practical on the scale required for comprehensive modelling only with the advent of teraFLOP computers. Incremental effort at LLNL would be focused on gyrokinetic simulation and development of the shell.
Numerical models of complex diapirs
NASA Astrophysics Data System (ADS)
Podladchikov, Yu.; Talbot, C.; Poliakov, A. N. B.
1993-12-01
Numerically modelled diapirs that rise into overburdens with viscous rheology produce a large variety of shapes. This work uses the finite-element method to study the development of diapirs that rise towards a surface on which a diapir-induced topography creeps flat or disperses ("erodes") at different rates. Slow erosion leads to diapirs with "mushroom" shapes, moderate erosion rate to "wine glass" diapirs and fast erosion to "beer glass"- and "column"-shaped diapirs. The introduction of a low-viscosity layer at the top of the overburden causes diapirs to develop into structures resembling a "Napoleon hat". These spread lateral sheets.
Numerical methods used in fusion science numerical modeling
NASA Astrophysics Data System (ADS)
Yagi, M.
2015-04-01
The dynamics of burning plasma is very complicated physics, which is dominated by multi-scale and multi-physics phenomena. To understand such phenomena, numerical simulations are indispensable. Fundamentals of numerical methods used in fusion science numerical modeling are briefly discussed in this paper. In addition, the parallelization technique such as open multi processing (OpenMP) and message passing interface (MPI) parallel programing are introduced and the loop-level parallelization is shown as an example.
Tam, Tze-wai; Ang, Put O
2009-07-21
A 3-dimensional individual-based model, the ReefModel, was developed to simulate the dynamical structure of coral reef community using object-oriented techniques. Interactions among functional groups of reef organisms were simulated in the model. The behaviours of these organisms were described with simple mechanistic rules that were derived from their general behaviours (e.g. growing habits, competitive mechanisms, response to physical disturbance) observed in natural coral reef communities. The model was implemented to explore the effects of physical disturbance on the dynamical structure of a 3-coral community that was characterized with three functional coral groups: tabular coral, foliaceous coral and massive coral. Simulation results suggest that (i) the integration of physical disturbance and differential responses (disturbance sensitivity and growing habit) of corals plays an important role in structuring coral communities; (ii) diversity of coral communities can be maximal under intermediate level of acute physical disturbance; (iii) multimodality exists in the final states and dynamic regimes of individual coral group as well as coral community structure, which results from the influence of small random spatial events occurring during the interactions among the corals in the community, under acute and repeated physical disturbances. These results suggest that alternative stable states and catastrophic regime shifts may exist in a coral community under unstable physical environment. PMID:19306887
Numerical Modeling of Ocean Circulation
NASA Astrophysics Data System (ADS)
Miller, Robert N.
2007-01-01
The modelling of ocean circulation is important not only for its own sake, but also in terms of the prediction of weather patterns and the effects of climate change. This book introduces the basic computational techniques necessary for all models of the ocean and atmosphere, and the conditions they must satisfy. It describes the workings of ocean models, the problems that must be solved in their construction, and how to evaluate computational results. Major emphasis is placed on examining ocean models critically, and determining what they do well and what they do poorly. Numerical analysis is introduced as needed, and exercises are included to illustrate major points. Developed from notes for a course taught in physical oceanography at the College of Oceanic and Atmospheric Sciences at Oregon State University, this book is ideal for graduate students of oceanography, geophysics, climatology and atmospheric science, and researchers in oceanography and atmospheric science. Features examples and critical examination of ocean modelling and results Demonstrates the strengths and weaknesses of different approaches Includes exercises to illustrate major points and supplement mathematical and physical details
Numerical Modeling of Nanoelectronic Devices
NASA Technical Reports Server (NTRS)
Klimeck, Gerhard; Oyafuso, Fabiano; Bowen, R. Chris; Boykin, Timothy
2003-01-01
Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.
This report presents a three-dimensional finite-element numerical model designed to simulate chemical transport in subsurface systems with temperature effect taken into account. The three-dimensional model is developed to provide (1) a tool of application, with which one is able ...
NUMERICAL MODELS FOR PREDICTING WATERSHED ACIDIFICATION
Three numerical models of watershed acidification, including the MAGIC II, ETD, and ILWAS models, are reviewed, and a comparative study is made of the specific process formulations that are incorporated in the models to represent hydrological, geochemical, and biogeochemical proc...
Numerical modeling of odorant uptake in the rat nasal cavity.
Yang, Geoffrey C; Scherer, Peter W; Zhao, Kai; Mozell, Maxwell M
2007-03-01
An anatomically accurate 3-dimensional numerical model of the right rat nasal cavity was developed and used to compute low, medium, and high flow rate inspiratory and expiratory mucosal odorant uptake (imposed patterning) for 3 odorants with different mucus solubilities. The computed surface mass flux distributions were compared with anatomic receptor gene expression zones identified in the literature. In general, simulations predicted that odorants that were highly soluble in mucus were absorbed dorsally and medially, corresponding roughly to receptors from one of the gene expression zones. Insoluble odorants tended to be absorbed more peripherally in the rat olfactory region corresponding to the other 2 zones. These findings also agreed in general with the electroolfactogram measurements and the voltage-sensitive dye measurements reported in the literature. This numerical approach is the first to predict detailed odorant flux information across the olfactory mucosa in the rat nasal cavity during inspiratory and expiratory flow and to relate it to anatomic olfactory receptor location, physiological function, and biochemical experiment. This numerical technique can allow us to separate the contributions of imposed and inherent patterning mechanisms on the rat olfactory mucosa. PMID:17220517
NEXT-GENERATION NUMERICAL MODELING: INCORPORATING ELASTICITY, ANISOTROPY AND ATTENUATION
S. LARSEN; ET AL
2001-03-01
A new effort has been initiated between the Department of Energy (DOE) and the Society of Exploration Geophysicists (SEG) to investigate what features the next generation of numerical seismic models should contain that will best address current technical problems encountered during exploration in increasingly complex geologies. This collaborative work is focused on designing and building these new models, generating synthetic seismic data through simulated surveys of various geometries, and using these data to test and validate new and improved seismic imaging algorithms. The new models will be both 2- and 3-dimensional and will include complex velocity structures as well as anisotropy and attenuation. Considerable attention is being focused on multi-component acoustic and elastic effects because it is now widely recognized that converted phases could play a vital role in improving the quality of seismic images. An existing, validated 3-D elastic modeling code is being used to generate the synthetic data. Preliminary elastic modeling results using this code are presented here along with a description of the proposed new models that will be built and tested.
Hydroelectric structures studies using 3-dimensional methods
Harrell, T.R.; Jones, G.V.; Toner, C.K. )
1989-01-01
Deterioration and degradation of aged, hydroelectric project structures can significantly affect the operation and safety of a project. In many cases, hydroelectric headworks (in particular) have complicated geometrical configurations, loading patterns and hence, stress conditions. An accurate study of such structures can be performed using 3-dimensional computer models. 3-D computer models can be used for both stability evaluation and for finite element stress analysis. Computer aided engineering processes facilitate the use of 3-D methods in both pre-processing and post-processing of data. Two actual project examples are used to emphasize the authors' points.
3-dimensional Oil Drift Simulations
NASA Astrophysics Data System (ADS)
Wettre, C.; Reistad, M.; Hjøllo, B.Å.
Simulation of oil drift has been an ongoing activity at the Norwegian Meteorological Institute since the 1970's. The Marine Forecasting Centre provides a 24-hour service for the Norwegian Pollution Control Authority and the oil companies operating in the Norwegian sector. The response time is 30 minutes. From 2002 the service is extended to simulation of oil drift from oil spills in deep water, using the DeepBlow model developed by SINTEF Applied Chemistry. The oil drift model can be applied both for instantaneous and continuous releases. The changes in the mass of oil and emulsion as a result of evaporation and emulsion are computed. For oil spill at deep water, hydrate formation and gas dissolution are taken into account. The properties of the oil depend on the oil type, and in the present version 64 different types of oil can be simulated. For accurate oil drift simulations it is important to have the best possible data on the atmospheric and oceanic conditions. The oil drift simulations at the Norwegian Meteorological Institute are always based on the most updated data from numerical models of the atmosphere and the ocean. The drift of the surface oil is computed from the vectorial sum of the surface current from the ocean model and the wave induced Stokes drift computed from wave energy spectra from the wave prediction model. In the new model the current distribution with depth is taken into account when calculating the drift of the dispersed oil droplets. Salinity and temperature profiles from the ocean model are needed in the DeepBlow model. The result of the oil drift simulations can be plotted on sea charts used for navigation, either as trajectory plots or particle plots showing the situation at a given time. The results can also be sent as data files to be included in the user's own GIS system.
Numerical wind speed simulation model
Ramsdell, J.V.; Athey, G.F.; Ballinger, M.Y.
1981-09-01
A relatively simple stochastic model for simulating wind speed time series that can be used as an alternative to time series from representative locations is described in this report. The model incorporates systematic seasonal variation of the mean wind, its standard deviation, and the correlation speeds. It also incorporates systematic diurnal variation of the mean speed and standard deviation. To demonstrate the model capabilities, simulations were made using model parameters derived from data collected at the Hanford Meteorology Station, and results of analysis of simulated and actual data were compared.
Survey of numerical electrostimulation models.
Reilly, J Patrick
2016-06-21
This paper evaluates results of a survey of electrostimulation models of myelinated nerve. Participants were asked to determine thresholds of excitation for 18 cases involving different characteristics of the neuron, the stimulation waveform, and the electrode arrangement. Responses were received from 7 investigators using 10 models. Excitation thresholds differed significantly among these models. For example, with a 2 ms monophasic stimulus pulse and an electrode/fiber distance of 1 cm, thresholds from the least to greatest value differed by a factor of 8.3; with a 5 μs pulse, thresholds differed by the factor 3.8. Significant differences in reported simulations point to the need for experimental validation. Additional efforts are needed to develop computational models for unmyelinated C-fibers, A-delta fibers, CNS neurons, and CNS Synapses. PMID:27223870
Survey of numerical electrostimulation models
NASA Astrophysics Data System (ADS)
Reilly, J. Patrick
2016-06-01
This paper evaluates results of a survey of electrostimulation models of myelinated nerve. Participants were asked to determine thresholds of excitation for 18 cases involving different characteristics of the neuron, the stimulation waveform, and the electrode arrangement. Responses were received from 7 investigators using 10 models. Excitation thresholds differed significantly among these models. For example, with a 2 ms monophasic stimulus pulse and an electrode/fiber distance of 1 cm, thresholds from the least to greatest value differed by a factor of 8.3; with a 5 μs pulse, thresholds differed by the factor 3.8. Significant differences in reported simulations point to the need for experimental validation. Additional efforts are needed to develop computational models for unmyelinated C-fibers, A-delta fibers, CNS neurons, and CNS Synapses.
Waste glass melter numerical and physical modeling
Eyler, L.L.; Peters, R.D.; Lessor, D.L.; Lowery, P.S.; Elliott, M.L.
1991-10-01
Results of physical and numerical simulation modeling of high-level liquid waste vitrification melters are presented. Physical modeling uses simulant fluids in laboratory testing. Visualization results provide insight into convective melt flow patterns from which information is derived to support performance estimation of operating melters and data to support numerical simulation. Numerical simulation results of several melter configurations are presented. These are in support of programs to evaluate melter operation characteristics and performance. Included are investigations into power skewing and alternating current electric field phase angle in a dual electrode pair reference design and bi-modal convective stability in an advanced design. 9 refs., 9 figs., 1 tab.
Numerical modeling for underground nuclear test monitoring
NASA Astrophysics Data System (ADS)
Taylor, Steven R.; Kamm, James R.
The symposium for Numerical Modeling for Underground Nuclear Test Monitoring was held March 23-25 in Durango, Colo. Funded by the DOE Office of Arms Control and Nonproliferation (OACN) and hosted by the Source Region Program at Los Alamos National Laboratory (LANL), the meetings's purpose was to discuss the state-of-the-art in numerical simulations of nuclear explosion phenomenology with applications to test-ban monitoring. In particular, we wished to focus on the uniqueness of model fits to data, the measurement and characterization of material response models, advanced modeling techniques, and applications of modeling to monitoring problems.The concept for the meeting arose through discussions with Marv Denny, who was on assignment at Department of Energy Headquarters from Lawrence Livermore National Laboratory (LLNL). In these conversations, the following question was discussed: how are numerical modeling techniques being used to understand the effects of explosion- source phenomenology on test-ban treaty monitoring? Numerical studies are becoming increasingly important in the evaluation of capabilities for proliferation monitoring; this trend has accelerated with the curtailment of the nuclear testing program. During these discussions, the issue of the uniqueness and limitations of numerical models arose. It was decided to address these questions by convening a group of experts to present and discuss the problems associated with modeling of close-in data from explosions.
Ferrofluids: Modeling, numerical analysis, and scientific computation
NASA Astrophysics Data System (ADS)
Tomas, Ignacio
This dissertation presents some developments in the Numerical Analysis of Partial Differential Equations (PDEs) describing the behavior of ferrofluids. The most widely accepted PDE model for ferrofluids is the Micropolar model proposed by R.E. Rosensweig. The Micropolar Navier-Stokes Equations (MNSE) is a subsystem of PDEs within the Rosensweig model. Being a simplified version of the much bigger system of PDEs proposed by Rosensweig, the MNSE are a natural starting point of this thesis. The MNSE couple linear velocity u, angular velocity w, and pressure p. We propose and analyze a first-order semi-implicit fully-discrete scheme for the MNSE, which decouples the computation of the linear and angular velocities, is unconditionally stable and delivers optimal convergence rates under assumptions analogous to those used for the Navier-Stokes equations. Moving onto the much more complex Rosensweig's model, we provide a definition (approximation) for the effective magnetizing field h, and explain the assumptions behind this definition. Unlike previous definitions available in the literature, this new definition is able to accommodate the effect of external magnetic fields. Using this definition we setup the system of PDEs coupling linear velocity u, pressure p, angular velocity w, magnetization m, and magnetic potential ϕ We show that this system is energy-stable and devise a numerical scheme that mimics the same stability property. We prove that solutions of the numerical scheme always exist and, under certain simplifying assumptions, that the discrete solutions converge. A notable outcome of the analysis of the numerical scheme for the Rosensweig's model is the choice of finite element spaces that allow the construction of an energy-stable scheme. Finally, with the lessons learned from Rosensweig's model, we develop a diffuse-interface model describing the behavior of two-phase ferrofluid flows and present an energy-stable numerical scheme for this model. For a
Mathematical modeling of electrocardiograms: a numerical study.
Boulakia, Muriel; Cazeau, Serge; Fernández, Miguel A; Gerbeau, Jean-Frédéric; Zemzemi, Nejib
2010-03-01
This paper deals with the numerical simulation of electrocardiograms (ECG). Our aim is to devise a mathematical model, based on partial differential equations, which is able to provide realistic 12-lead ECGs. The main ingredients of this model are classical: the bidomain equations coupled to a phenomenological ionic model in the heart, and a generalized Laplace equation in the torso. The obtention of realistic ECGs relies on other important features--including heart-torso transmission conditions, anisotropy, cell heterogeneity and His bundle modeling--that are discussed in detail. The numerical implementation is based on state-of-the-art numerical methods: domain decomposition techniques and second order semi-implicit time marching schemes, offering a good compromise between accuracy, stability and efficiency. The numerical ECGs obtained with this approach show correct amplitudes, shapes and polarities, in all the 12 standard leads. The relevance of every modeling choice is carefully discussed and the numerical ECG sensitivity to the model parameters investigated. PMID:20033779
Numerical flow modeling of power plant windboxes
LaRose, J.A.; Hopkins, M.W.
1995-12-31
Numerical flow modeling has become an increasingly important design and analysis tool for improving the air distribution to power plant burners. Uniform air distribution allows the burners to perform as designed to achieve the lowest possible emissions and best fuel burn-out. Modifications can be made internal to the existing windbox to improve the burner-to-burner and burner peripheral air distributions. These modifications can include turning vanes, flow splitters, perforated plate, and burner shrouding. Numerical modeling allows the analysis of design trade-offs between adding flow resistance, fan power, and windbox modification construction cost. Numerical modeling has advantages over physical modeling in that actual geometric scales and air temperatures are used. Advantages over a field data based study include the ability to quickly and cheaply analyze a variety of design options without actually modifying the windbox, and the availability of significantly more data with which to interpret the results. Costs to perform a numerical study are generally one-half to one-third of the cost to perform a physical flow model and can be one-forth of the cost to perform a field study. The continued development of affordable, high speed, large memory workstations and reliable, commercially available computation fluid dynamics (CFD) software allows practical analyses of power plant windboxes. This paper discusses (1) the impact of air distribution on burner performance, (2) the methodology used to perform numerical flow modeling of power plant windboxes, and (3) the results from several windbox analyses including available post-modification observations.
Numerical noise in ocean and estuarine models
Walters, R.; Carey, G.F.
1984-01-01
Approximate methods for solving the shallow water equations may lead to solutions exhibiting large fictitious, numerically-induced oscillations. The analysis of the discrete dispersion relation and modal solutions of small wavelengths provides a powerful technique for assessing the sensitivity of alternative numerical schemes to irregular data which may lead to such oscillatory numerical noise. For those schemes where phase speed vanishes at a finite wavenumber or there are multiple roots for wavenumber, oscillation modes can exist which are uncoupled from the dynamics of the problem. The discrete modal analysis approach is used here to identify two classes of spurious oscillation modes associated respectively with the two different asymptotic limits corresponding to estuarine and large scale ocean models. The analysis provides further insight into recent numerical results for models which include large spatial scales and Coriolis acceleration. ?? 1984.
Numerical Modeling of Ablation Heat Transfer
NASA Technical Reports Server (NTRS)
Ewing, Mark E.; Laker, Travis S.; Walker, David T.
2013-01-01
A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.
Numerical Based Linear Model for Dipole Magnets
Li,Y.; Krinsky, S.; Rehak, M.
2009-05-04
In this paper, we discuss an algorithm for constructing a numerical linear optics model for dipole magnets from a 3D field map. The difference between the numerical model and K. Brown's analytic approach is investigated and clarified. It was found that the optics distortion due to the dipoles' fringe focusing must be properly taken into account to accurately determine the chromaticities. In NSLS-II, there are normal dipoles with 35-mm gap and dipoles for infrared sources with 90-mm gap. This linear model of the dipole magnets is applied to the NSLS-II lattice design to match optics parameters between the DBA cells having dipoles with different gaps.
Numerical FEM modeling in dental implantology
NASA Astrophysics Data System (ADS)
Roateşi, Iulia; Roateşi, Simona
2016-06-01
This paper is devoted to a numerical approach of the stress and displacement calculation of a system made up of dental implant, ceramic crown and surrounding bone. This is the simulation of a clinical situation involving both biological - the bone tissue, and non-biological - the implant and the crown, materials. On the other hand this problem deals with quite fine technical structure details - the threads, tapers, etc with a great impact in masticatory force transmission. Modeling the contact between the implant and the bone tissue is important to a proper bone-implant interface model and implant design. The authors proposed a three-dimensional numerical model to assess the biomechanical behaviour of this complex structure in order to evaluate its stability by determining the risk zones. A comparison between this numerical analysis and clinical cases is performed and a good agreement is obtained.
Development of a numerical model for vehicle-bridge interaction analysis of railway bridges
NASA Astrophysics Data System (ADS)
Kim, Hee Ju; Cho, Eun Sang; Ham, Jun Su; Park, Ki Tae; Kim, Tae Heon
2016-04-01
In the field of civil engineering, analyzing dynamic response was main concern for a long time. These analysis methods can be divided into moving load analysis method and moving mass analysis method, and formulating each an equation of motion has recently been studied after dividing vehicles and bridges. In this study, the numerical method is presented, which can consider the various train types and can solve the equations of motion for a vehicle-bridge interaction analysis by non-iteration procedure through formulating the coupled equations for motion. Also, 3 dimensional accurate numerical models was developed by KTX-vehicle in order to analyze dynamic response characteristics. The equations of motion for the conventional trains are derived, and the numerical models of the conventional trains are idealized by a set of linear springs and dashpots with 18 degrees of freedom. The bridge models are simplified by the 3 dimensional space frame element which is based on the Euler-Bernoulli theory. The rail irregularities of vertical and lateral directions are generated by PSD functions of the Federal Railroad Administration (FRA).
Numerical Modelling of Ground Penetrating Radar Antennas
NASA Astrophysics Data System (ADS)
Giannakis, Iraklis; Giannopoulos, Antonios; Pajewski, Lara
2014-05-01
Numerical methods are needed in order to solve Maxwell's equations in complicated and realistic problems. Over the years a number of numerical methods have been developed to do so. Amongst them the most popular are the finite element, finite difference implicit techniques, frequency domain solution of Helmontz equation, the method of moments, transmission line matrix method. However, the finite-difference time-domain method (FDTD) is considered to be one of the most attractive choice basically because of its simplicity, speed and accuracy. FDTD first introduced in 1966 by Kane Yee. Since then, FDTD has been established and developed to be a very rigorous and well defined numerical method for solving Maxwell's equations. The order characteristics, accuracy and limitations are rigorously and mathematically defined. This makes FDTD reliable and easy to use. Numerical modelling of Ground Penetrating Radar (GPR) is a very useful tool which can be used in order to give us insight into the scattering mechanisms and can also be used as an alternative approach to aid data interpretation. Numerical modelling has been used in a wide range of GPR applications including archeology, geophysics, forensic, landmine detection etc. In engineering, some applications of numerical modelling include the estimation of the effectiveness of GPR to detect voids in bridges, to detect metal bars in concrete, to estimate shielding effectiveness etc. The main challenges in numerical modelling of GPR for engineering applications are A) the implementation of the dielectric properties of the media (soils, concrete etc.) in a realistic way, B) the implementation of the geometry of the media (soils inhomogeneities, rough surface, vegetation, concrete features like fractures and rock fragments etc.) and C) the detailed modelling of the antenna units. The main focus of this work (which is part of the COST Action TU1208) is the accurate and realistic implementation of GPR antenna units into the FDTD
Numerical Modeling of Ocean Acoustic Wavefields
NASA Astrophysics Data System (ADS)
Tappert, Frederick
1997-08-01
The U.S. Navy requires real-time ``acoustic performance prediction'' models in order to optimize sonar tactics in naval combat situations. The need for numerical models that solve the acoustic wave equation in realistic ocean environments is being met by a collaborative effort between university researchers, industrial contractors, and navy laboratory workers. This paper discusses one particularly successful numerical model, called the PE/SSF model, that was originally developed by the author. Here PE stands for Parabolic Equation, a good approximation to the elliptic Helmholtz equation; and SSF stands for the Split-Step Fourier algorithm, a highly efficient marching algorithm for solving parabolic type equations. These techniques are analyzed, and examples are displayed of ocean acoustic wavefields generated by the PE/SSF model.
3-dimensional bioprinting for tissue engineering applications.
Gu, Bon Kang; Choi, Dong Jin; Park, Sang Jun; Kim, Min Sup; Kang, Chang Mo; Kim, Chun-Ho
2016-01-01
The 3-dimensional (3D) printing technologies, referred to as additive manufacturing (AM) or rapid prototyping (RP), have acquired reputation over the past few years for art, architectural modeling, lightweight machines, and tissue engineering applications. Among these applications, tissue engineering field using 3D printing has attracted the attention from many researchers. 3D bioprinting has an advantage in the manufacture of a scaffold for tissue engineering applications, because of rapid-fabrication, high-precision, and customized-production, etc. In this review, we will introduce the principles and the current state of the 3D bioprinting methods. Focusing on some of studies that are being current application for biomedical and tissue engineering fields using printed 3D scaffolds. PMID:27114828
Numerical modeling tools for chemical vapor deposition
NASA Technical Reports Server (NTRS)
Jasinski, Thomas J.; Childs, Edward P.
1992-01-01
Development of general numerical simulation tools for chemical vapor deposition (CVD) was the objective of this study. Physical models of important CVD phenomena were developed and implemented into the commercial computational fluid dynamics software FLUENT. The resulting software can address general geometries as well as the most important phenomena occurring with CVD reactors: fluid flow patterns, temperature and chemical species distribution, gas phase and surface deposition. The physical models are documented which are available and examples are provided of CVD simulation capabilities.
Teleportation of a 3-dimensional GHZ State
NASA Astrophysics Data System (ADS)
Cao, Hai-Jing; Wang, Huai-Sheng; Li, Peng-Fei; Song, He-Shan
2012-05-01
The process of teleportation of a completely unknown 3-dimensional GHZ state is considered. Three maximally entangled 3-dimensional Bell states function as quantum channel in the scheme. This teleportation scheme can be directly generalized to teleport an unknown d-dimensional GHZ state.
Numerical comparison of strong Langmuir turbulence models
NASA Technical Reports Server (NTRS)
Shen, Mei-Mei; Nicholson, D. R.
1987-01-01
Two models of Langmuir turbulence, the nonlinear Schroedinger equation and the Zakharov equations, are solved numerically for an initial value problem in which the electric field evolves from an almost flat initial condition via the modulational instability and finally saturates into a set of solitons. The two models agree well with each other only when the initial dimensionless electric field has an amplitude less than unity. An analytic soliton gas model consisting of equal-amplitude, randomly spaced, zero-speed solitons is remarkably good at reproducing the time-averaged Fourier spectra in both cases.
Numerical Modelling Of Pumpkin Balloon Instability
NASA Astrophysics Data System (ADS)
Wakefield, D.
Tensys have been involved in the numerical formfinding and load analysis of architectural stressed membrane structures for 15 years. They have recently broadened this range of activities into the `lighter than air' field with significant involvement in aerostat and heavy-lift hybrid airship design. Since early 2004 they have been investigating pumpkin balloon instability on behalf of the NASA ULDB programme. These studies are undertaken using inTENS, an in-house finite element program suite based upon the Dynamic Relaxation solution method and developed especially for the non-linear analysis and patterning of membrane structures. The paper describes the current state of an investigation that started with a numerical simulation of the lobed cylinder problem first studied by Calladine. The influence of material properties and local geometric deformation on stability is demonstrated. A number of models of complete pumpkin balloons have then been established, including a 64-gore balloon with geometry based upon Julian Nott's Endeavour. This latter clefted dramatically upon initial inflation, a phenomenon that has been reproduced in the numerical model. Ongoing investigations include the introduction of membrane contact modelling into inTENS and correlation studies with the series of large-scale ULDB models currently in preparation.
Numerical modeling of flow through orifice meters
NASA Astrophysics Data System (ADS)
Sheikholesiami, M. Z.; Patel, B. R.
1988-03-01
Numerical modeling is performed for turbulent flow through orifice meters using Creare's computer program FLUENT. FLUENT solves the time averaged Navier-Stokes equations in 2-D and 3-D Cartesian or cylindrical coordinates. Turbulence is simulated using a two equation k-epsilon or algebraic stress turbulence model. It is shown that an 80 x 60 grid distribution is sufficient to resolve the flow field around the orifice. The variations in discharge coefficient are studied as a result of variation in beta ratio, Reynolds number, upstream and downstream boundary conditions, pipe surface roughness, and upstream swirl. The effects of beta ratio and Reynolds number on the discharge coefficient are shown to be similar to the experimental data. It is also shown that the surface roughness can increase the discharge coefficient by about 0.7 percent for the range of roughness heights encountered in practice. The numerical modeling approach would be most effective if it is combined with a systematic experimental program that can supply the necessary boundary conditions. It is recommended that numerical modeling be used for the study of other flow meters.
A Numerical Model for Atomtronic Circuit Analysis
Chow, Weng W.; Straatsma, Cameron J. E.; Anderson, Dana Z.
2015-07-16
A model for studying atomtronic devices and circuits based on finite-temperature Bose-condensed gases is presented. The approach involves numerically solving equations of motion for atomic populations and coherences, derived using the Bose-Hubbard Hamiltonian and the Heisenberg picture. The resulting cluster expansion is truncated at a level giving balance between physics rigor and numerical demand mitigation. This approach allows parametric studies involving time scales that cover both the rapid population dynamics relevant to nonequilibrium state evolution, as well as the much longer time durations typical for reaching steady-state device operation. This model is demonstrated by studying the evolution of a Bose-condensed gas in the presence of atom injection and extraction in a double-well potential. In this configuration phase locking between condensates in each well of the potential is readily observed, and its influence on the evolution of the system is studied.
Numerical model for atomtronic circuit analysis
NASA Astrophysics Data System (ADS)
Chow, Weng W.; Straatsma, Cameron J. E.; Anderson, Dana Z.
2015-07-01
A model for studying atomtronic devices and circuits based on finite-temperature Bose-condensed gases is presented. The approach involves numerically solving equations of motion for atomic populations and coherences, derived using the Bose-Hubbard Hamiltonian and the Heisenberg picture. The resulting cluster expansion is truncated at a level giving balance between physics rigor and numerical demand mitigation. This approach allows parametric studies involving time scales that cover both the rapid population dynamics relevant to nonequilibrium state evolution, as well as the much longer time durations typical for reaching steady-state device operation. The model is demonstrated by studying the evolution of a Bose-condensed gas in the presence of atom injection and extraction in a double-well potential. In this configuration phase locking between condensates in each well of the potential is readily observed, and its influence on the evolution of the system is studied.
Lattice Boltzmann model for numerical relativity.
Ilseven, E; Mendoza, M
2016-02-01
In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems. PMID:26986435
Lattice Boltzmann model for numerical relativity
NASA Astrophysics Data System (ADS)
Ilseven, E.; Mendoza, M.
2016-02-01
In the Z4 formulation, Einstein equations are written as a set of flux conservative first-order hyperbolic equations that resemble fluid dynamics equations. Based on this formulation, we construct a lattice Boltzmann model for numerical relativity and validate it with well-established tests, also known as "apples with apples." Furthermore, we find that by increasing the relaxation time, we gain stability at the cost of losing accuracy, and by decreasing the lattice spacings while keeping a constant numerical diffusivity, the accuracy and stability of our simulations improve. Finally, in order to show the potential of our approach, a linear scaling law for parallelization with respect to number of CPU cores is demonstrated. Our model represents the first step in using lattice kinetic theory to solve gravitational problems.
Infrared radiation parameterizations in numerical climate models
NASA Technical Reports Server (NTRS)
Chou, Ming-Dah; Kratz, David P.; Ridgway, William
1991-01-01
This study presents various approaches to parameterizing the broadband transmission functions for utilization in numerical climate models. One-parameter scaling is applied to approximate a nonhomogeneous path with an equivalent homogeneous path, and the diffuse transmittances are either interpolated from precomputed tables or fit by analytical functions. Two-parameter scaling is applied to parameterizing the carbon dioxide and ozone transmission functions in both the lower and middle atmosphere. Parameterizations are given for the nitrous oxide and methane diffuse transmission functions.
Numerical modeling of the acoustic guitar
NASA Astrophysics Data System (ADS)
Chaigne, Antoine; Derveaux, Grégoire; Joly, Patrick; Bécache, Eliane
2003-10-01
An interactive DVD has been created, based on a numerical model of the acoustic guitar. In a first chapter, the retained physical model is described and illustrated, from the pluck to the 3D radiation field. The second chapter is devoted to the presentation of the numerical tools used for solving the equations of the model. Numerical simulations of plate vibrations and radiated sound pressure are shown in the third chapter. A number of simulated sounds are presented and analyzed in the fourth chapter. In addition, the DVD includes a discussion between a guitar maker, an acoustician, a guitar player and a mathematician. This discussion is entitled ``towards a common language.'' Its aim is to show the interest of simulations with respect to complementary professional approaches of the instrument. This DVD received the Henri Poincaré Prize from the 8th Research Film Festival of Nancy (June 2003), sponsored by the CNRS, in the category ``Documents for the scientific community and illustrations of the research for teaching purpose.''
Avoiding numerical pitfalls in social force models
NASA Astrophysics Data System (ADS)
Köster, Gerta; Treml, Franz; Gödel, Marion
2013-06-01
The social force model of Helbing and Molnár is one of the best known approaches to simulate pedestrian motion, a collective phenomenon with nonlinear dynamics. It is based on the idea that the Newtonian laws of motion mostly carry over to pedestrian motion so that human trajectories can be computed by solving a set of ordinary differential equations for velocity and acceleration. The beauty and simplicity of this ansatz are strong reasons for its wide spread. However, the numerical implementation is not without pitfalls. Oscillations, collisions, and instabilities occur even for very small step sizes. Classic solution ideas from molecular dynamics do not apply to the problem because the system is not Hamiltonian despite its source of inspiration. Looking at the model through the eyes of a mathematician, however, we realize that the right hand side of the differential equation is nondifferentiable and even discontinuous at critical locations. This produces undesirable behavior in the exact solution and, at best, severe loss of accuracy in efficient numerical schemes even in short range simulations. We suggest a very simple mollified version of the social force model that conserves the desired dynamic properties of the original many-body system but elegantly and cost efficiently resolves several of the issues concerning stability and numerical resolution.
Avoiding numerical pitfalls in social force models.
Köster, Gerta; Treml, Franz; Gödel, Marion
2013-06-01
The social force model of Helbing and Molnár is one of the best known approaches to simulate pedestrian motion, a collective phenomenon with nonlinear dynamics. It is based on the idea that the Newtonian laws of motion mostly carry over to pedestrian motion so that human trajectories can be computed by solving a set of ordinary differential equations for velocity and acceleration. The beauty and simplicity of this ansatz are strong reasons for its wide spread. However, the numerical implementation is not without pitfalls. Oscillations, collisions, and instabilities occur even for very small step sizes. Classic solution ideas from molecular dynamics do not apply to the problem because the system is not Hamiltonian despite its source of inspiration. Looking at the model through the eyes of a mathematician, however, we realize that the right hand side of the differential equation is nondifferentiable and even discontinuous at critical locations. This produces undesirable behavior in the exact solution and, at best, severe loss of accuracy in efficient numerical schemes even in short range simulations. We suggest a very simple mollified version of the social force model that conserves the desired dynamic properties of the original many-body system but elegantly and cost efficiently resolves several of the issues concerning stability and numerical resolution. PMID:23848804
Numerical modeling of explosions for nuclear monitoring
NASA Astrophysics Data System (ADS)
Stevens, J. L.
2014-12-01
Monitoring the Earth for underground nuclear explosions requires a detailed understanding of the explosion source. In this context, "source" refers to the source of seismic waves, and it is generated by the complex nonlinear near-source motion that accompanies the nuclear explosion. In particular, nuclear monitoring requires understanding the transition from the hydrodynamic to elastic regimes, and propagation of waveforms from the source to stations at distances of hundreds to thousands of kilometers. In the transition region, shear strength is critically important, as are changes in shear strength as the shock wave propagates. Numerical modeling using 1D spherically symmetric, 2D axisymmetric and full 3D calculations provides important insights into the seismic source and the waveforms it generates. Important considerations for numerical modeling include emplacement conditions (tamped or in a cavity), source type (chemical or nuclear), material models for strength and strength reduction, and geologic conditions including topography and tectonic stresses in the source region. In addition to calculating the near source ground motion, we propagate the near source solution to regional and teleseismic distances where the observations of seismic signals from nuclear explosions are made. The objectives of nuclear monitoring are detection of seismic events (earthquakes, quarry blasts and other sources in addition to nuclear explosions), accurate location of these events, discrimination of nuclear explosions from other types of sources, and estimation of nuclear explosion yield. Numerical modeling is particularly important for discrimination and yield estimation. Numerical modeling is used to understand unexpected anomalies that occur, such as the large surface waves generated by the three North Korean nuclear tests, which may have been caused by a difference in tectonic stress state between North Korea and other test sites. Another important issue that can be addressed
Advanced Numerical Model for Irradiated Concrete
Giorla, Alain B.
2015-03-01
In this report, we establish a numerical model for concrete exposed to irradiation to address these three critical points. The model accounts for creep in the cement paste and its coupling with damage, temperature and relative humidity. The shift in failure mode with the loading rate is also properly represented. The numerical model for creep has been validated and calibrated against different experiments in the literature [Wittmann, 1970, Le Roy, 1995]. Results from a simplified model are shown to showcase the ability of numerical homogenization to simulate irradiation effects in concrete. In future works, the complete model will be applied to the analysis of the irradiation experiments of Elleuch et al. [1972] and Kelly et al. [1969]. This requires a careful examination of the experimental environmental conditions as in both cases certain critical information are missing, including the relative humidity history. A sensitivity analysis will be conducted to provide lower and upper bounds of the concrete expansion under irradiation, and check if the scatter in the simulated results matches the one found in experiments. The numerical and experimental results will be compared in terms of expansion and loss of mechanical stiffness and strength. Both effects should be captured accordingly by the model to validate it. Once the model has been validated on these two experiments, it can be applied to simulate concrete from nuclear power plants. To do so, the materials used in these concrete must be as well characterized as possible. The main parameters required are the mechanical properties of each constituent in the concrete (aggregates, cement paste), namely the elastic modulus, the creep properties, the tensile and compressive strength, the thermal expansion coefficient, and the drying shrinkage. These can be either measured experimentally, estimated from the initial composition in the case of cement paste, or back-calculated from mechanical tests on concrete. If some
Multiscale numerical modeling of levee breach processes
NASA Astrophysics Data System (ADS)
Kees, C. E.; Farthing, M. W.; Akkerman, I.; Bazilevs, Y.
2010-12-01
One of the dominant failure modes of levees during flood and storm surge events is erosion-based breach formation due to high velocity flow over the back (land-side) slope. Modeling the breaching process numerically is challenging due to both physical and geometric complexity that develops and evolves during the overtopping event. The surface water flows are aerated and sediment-laden mixtures in the supercritical and turbulent regimes. The air/water free surface may undergo perturbations on the same order as the depth or even topological change (breaking). Likewise the soil/fluid interface is characterized by evolving headcuts, which are essentially moving discontinuities in the soil surface elevation. The most widely used models of levee breaching are nevertheless based on depth-integrated models of flow, sediment transport, and bed morphology. In this work our objective is to explore models with less restrictive modeling assumptions, which have become computationally tractable due to advances in both numerical methods and high-performance computing hardware. In particular, we present formulations of fully three-dimensional flow, transport, and morphological evolution for overtopping and breaching processes and apply recently developed finite element and level set methods to solve the governing equations for relevant test problems.
Numerical simulations and modeling of turbulent combustion
NASA Astrophysics Data System (ADS)
Cuenot, B.
Turbulent combustion is the basic physical phenomenon responsible for efficient energy release by any internal combustion engine. However it is accompanied by other undesirable phenomena such as noise, pollutant species emission or damaging instabilities that may even lead to the system desctruction. It is then crucial to control this phenomenon, to understand all its mecanisms and to master it in industrial systems. For long time turbulent combustion has been explored only through theory and experiment. But the rapid increase of computers power during the last years has allowed an important development of numerical simulation, that has become today an essential tool for research and technical design. Direct numerical simulation has then allowed to rapidly progress in the knowledge of turbulent flame structures, leading to new modelisations for steady averaged simulations. Recently large eddy simulation has made a new step forward by refining the description of complex and unsteady flames. The main problem that arises when performing numerical simulation of turbulent combustion is linked to the description of the flame front. Being very thin, it can not however be reduced to a simple interface as it is the location of intense chemical transformation and of strong variations of thermodynamical quantities. Capturing the internal structure of a zone with a thickness of the order of 0.1 mm in a computation with a mesh step 10 times larger being impossible, it is necessary to model the turbulent flame. Models depend on the chemical structure of the flame, on the ambiant turbulence, on the combustion regime (flamelets, distributed combustion, etc.) and on the reactants injection mode (premixed or not). One finds then a large class of models, from the most simple algebraic model with a one-step chemical kinetics, to the most complex model involving probablity density functions, cross-correlations and multiple-step or fully complex chemical kinetics.
Numerical modelling of swirling diffusive flames
NASA Astrophysics Data System (ADS)
Parra-Santos, Teresa; Perez, Ruben; Szasz, Robert Z.; Gutkowski, Artur N.; Castro, Francisco
2016-03-01
Computational Fluid Dynamics has been used to study the mixing and combustion of two confined jets whose setup and operating conditions are those of the benchmark of Roback and Johnson. Numerical model solves 3D transient Navier Stokes for turbulent and reactive flows. Averaged velocity profiles using RNG swirl dominated k-epsilon model have been validated with experimental measurements from other sources for the non reactive case. The combustion model is Probability Density Function. Bearing in mind the annular jet has swirl number over 0.5, a vortex breakdown appears in the axis of the burner. Besides, the sudden expansion with a ratio of 2 in diameter between nozzle exits and the test chamber produces the boundary layer separation with the corresponding torus shape recirculation. Contrasting the mixing and combustion models, the last one produces the reduction of the vortex breakdown.
Numerical modeling of flowing soft materials
NASA Astrophysics Data System (ADS)
Toschi, Federico; Benzi, Roberto; Bernaschi, Massimo; Perlekar, Prasad; Sbragaglia, Mauro; Succi, Sauro
2012-11-01
The structural properties of soft-flowing and non-ergodic materials, such as emulsions, foams and gels shares similarities with the three basic states of matter (solid, liquid and gas). The macroscopic properties are characterized by non-standard features such as non-Newtonian rheology, long-time relaxation, caging effects, enhanced viscosity, structural arrest, hysteresis, dynamic disorder, aging and related phenomena. Large scale non-homogeneities can develop, even under simple shear conditions, by means of the formation of macroscopic bands of widely different viscosities (``shear banding'' phenomena). We employ a numerical model based on the Lattice Boltzmann method to perform numerical simulations of soft-matter under flowing conditions. Results of 3d simulations are presented and compared to previous 2d investigations.
Kohyama, Hiroaki
2008-07-01
We construct the phase diagram of the quark-antiquark and diquark condensates at finite temperature and density in the 2+1 dimensional (3D) two flavor massless Gross-Neveu (GN) model with the 4-component quarks. In contrast to the case of the 2-component quarks, there appears the coexisting phase of the quark-antiquark and diquark condensates. This is the crucial difference between the 2-component and 4-component quark cases in the 3D GN model. The coexisting phase is also seen in the 4D Nambu Jona-Lasinio model. Then we see that the 3D GN model with the 4-component quarks bears closer resemblance to the 4D Nambu Jona-Lasinio model.
Kobayashi, Kazuyoshi; Imagama, Shiro; Muramoto, Akio; Ito, Zenya; Ando, Kei; Yagi, Hideki; Hida, Tetsuro; Ito, Kenyu; Ishikawa, Yoshimoto; Tsushima, Mikito; Ishiguro, Naoki
2015-01-01
ABSTRACT In severe spinal deformity, pain and neurological disorder may be caused by spinal cord compression. Surgery for spinal reconstruction is desirable, but may be difficult in a case with severe deformity. Here, we show the utility of a 3D NaCl (salt) model in preoperative planning of anterior reconstruction using a rib strut in a 49-year-old male patient with cervicothoracic degenerative spondylosis. We performed surgery in two stages: a posterior approach with decompression and posterior instrumentation with a pedicle screw; followed by a second operation using an anterior approach, for which we created a 3D NaCl model including the cervicothoracic lesion, spinal deformity, and ribs for anterior reconstruction. The 3D NaCl model was easily scraped compared with a conventional plaster model and was useful for planning of resection and identification of a suitable rib for grafting in a preoperative simulation. Surgery was performed successfully with reference to the 3D NaCl model. We conclude that preoperative simulation with a 3D NaCl model contributes to performance of anterior reconstruction using a rib strut in a case of cervicothoracic deformity. PMID:26412901
Posttraumatic Orbital Emphysema: A Numerical Model
Skorek, Andrzej; Kłosowski, Paweł; Plichta, Łukasz; Zmuda Trzebiatowski, Marcin; Lemski, Paweł
2014-01-01
Orbital emphysema is a common symptom accompanying orbital fracture. The pathomechanism is still not recognized and the usually assumed cause, elevated pressure in the upper airways connected with sneezing or coughing, does not always contribute to the occurrence of this type of fracture. Observations based on the finite model (simulating blowout type fracture) of the deformations of the inferior orbital wall after a strike in its lower rim. Authors created a computer numeric model of the orbit with specified features—thickness and resilience modulus. During simulation an evenly spread 14400 N force was applied to the nodular points in the inferior rim (the maximal value not causing cracking of the outer rim, but only ruptures in the inferior wall). The observation was made from 1 · 10−3 to 1 · 10−2 second after a strike. Right after a strike dislocations of the inferior orbital wall toward the maxillary sinus were observed. Afterwards a retrograde wave of the dislocation of the inferior wall toward the orbit was noticed. Overall dislocation amplitude reached about 6 mm. Based on a numeric model of the orbit submitted to a strike in the inferior wall an existence of a retrograde shock wave causing orbital emphysema has been found. PMID:25309749
Posttraumatic orbital emphysema: a numerical model.
Skorek, Andrzej; Kłosowski, Paweł; Plichta, Lukasz; Raczyńska, Dorota; Zmuda Trzebiatowski, Marcin; Lemski, Paweł
2014-01-01
Orbital emphysema is a common symptom accompanying orbital fracture. The pathomechanism is still not recognized and the usually assumed cause, elevated pressure in the upper airways connected with sneezing or coughing, does not always contribute to the occurrence of this type of fracture. Observations based on the finite model (simulating blowout type fracture) of the deformations of the inferior orbital wall after a strike in its lower rim. Authors created a computer numeric model of the orbit with specified features-thickness and resilience modulus. During simulation an evenly spread 14400 N force was applied to the nodular points in the inferior rim (the maximal value not causing cracking of the outer rim, but only ruptures in the inferior wall). The observation was made from 1 · 10(-3) to 1 · 10(-2) second after a strike. Right after a strike dislocations of the inferior orbital wall toward the maxillary sinus were observed. Afterwards a retrograde wave of the dislocation of the inferior wall toward the orbit was noticed. Overall dislocation amplitude reached about 6 mm. Based on a numeric model of the orbit submitted to a strike in the inferior wall an existence of a retrograde shock wave causing orbital emphysema has been found. PMID:25309749
Numerical modeling of vertical cavity semiconductor lasers
Chow, W.W.; Hadley, G.R.
1996-08-01
A vertical cavity surface emitting laser (VCSEL) is a diode laser whose optical cavity is formed by growing or depositing DBR mirror stacks that sandwich an active gain region. The resulting short cavity supports lasing into a single longitudinal mode normal to the wafer, making these devices ideal for a multitude of applications, ranging from high-speed communication to high-power sources (from 2D arrays). This report describes the development of a numerical VCSEL model, whose goal is to both further their understanding of these complex devices and provide a tool for accurate design and data analysis.
NASA Technical Reports Server (NTRS)
Kubrynski, Krzysztof
1991-01-01
A subcritical panel method applied to flow analysis and aerodynamic design of complex aircraft configurations is presented. The analysis method is based on linearized, compressible, subsonic flow equations and indirect Dirichlet boundary conditions. Quadratic dipol and linear source distribution on flat panels are applied. In the case of aerodynamic design, the geometry which minimizes differences between design and actual pressure distribution is found iteratively, using numerical optimization technique. Geometry modifications are modeled by surface transpiration concept. Constraints in respect to resulting geometry can be specified. A number of complex 3-dimensional design examples are presented. The software is adopted to personal computers, and as result an unexpected low cost of computations is obtained.
NASA Astrophysics Data System (ADS)
Bachche, Shivaji; Oka, Koichi
2013-06-01
This paper presents the comparative study of various color space models to determine the suitable color space model for detection of green sweet peppers. The images were captured by using CCD cameras and infrared cameras and processed by using Halcon image processing software. The LED ring around the camera neck was used as an artificial lighting to enhance the feature parameters. For color images, CieLab, YIQ, YUV, HSI and HSV whereas for infrared images, grayscale color space models were selected for image processing. In case of color images, HSV color space model was found more significant with high percentage of green sweet pepper detection followed by HSI color space model as both provides information in terms of hue/lightness/chroma or hue/lightness/saturation which are often more relevant to discriminate the fruit from image at specific threshold value. The overlapped fruits or fruits covered by leaves can be detected in better way by using HSV color space model as the reflection feature from fruits had higher histogram than reflection feature from leaves. The IR 80 optical filter failed to distinguish fruits from images as filter blocks useful information on features. Computation of 3D coordinates of recognized green sweet peppers was also conducted in which Halcon image processing software provides location and orientation of the fruits accurately. The depth accuracy of Z axis was examined in which 500 to 600 mm distance between cameras and fruits was found significant to compute the depth distance precisely when distance between two cameras maintained to 100 mm.
Numerical modeling for dilute and dense sprays
NASA Technical Reports Server (NTRS)
Chen, C. P.; Kim, Y. M.; Shang, H. M.; Ziebarth, J. P.; Wang, T. S.
1992-01-01
We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.
Cardiothoracic Applications of 3-dimensional Printing.
Giannopoulos, Andreas A; Steigner, Michael L; George, Elizabeth; Barile, Maria; Hunsaker, Andetta R; Rybicki, Frank J; Mitsouras, Dimitris
2016-09-01
Medical 3-dimensional (3D) printing is emerging as a clinically relevant imaging tool in directing preoperative and intraoperative planning in many surgical specialties and will therefore likely lead to interdisciplinary collaboration between engineers, radiologists, and surgeons. Data from standard imaging modalities such as computed tomography, magnetic resonance imaging, echocardiography, and rotational angiography can be used to fabricate life-sized models of human anatomy and pathology, as well as patient-specific implants and surgical guides. Cardiovascular 3D-printed models can improve diagnosis and allow for advanced preoperative planning. The majority of applications reported involve congenital heart diseases and valvular and great vessels pathologies. Printed models are suitable for planning both surgical and minimally invasive procedures. Added value has been reported toward improving outcomes, minimizing perioperative risk, and developing new procedures such as transcatheter mitral valve replacements. Similarly, thoracic surgeons are using 3D printing to assess invasion of vital structures by tumors and to assist in diagnosis and treatment of upper and lower airway diseases. Anatomic models enable surgeons to assimilate information more quickly than image review, choose the optimal surgical approach, and achieve surgery in a shorter time. Patient-specific 3D-printed implants are beginning to appear and may have significant impact on cosmetic and life-saving procedures in the future. In summary, cardiothoracic 3D printing is rapidly evolving and may be a potential game-changer for surgeons. The imager who is equipped with the tools to apply this new imaging science to cardiothoracic care is thus ideally positioned to innovate in this new emerging imaging modality. PMID:27149367
Convecting reference frames and invariant numerical models
NASA Astrophysics Data System (ADS)
Bihlo, Alexander; Nave, Jean-Christophe
2014-09-01
In the recent paper by Bernardini et al. [1] the discrepancy in the performance of finite difference and spectral models for simulations of flows with a preferential direction of propagation was studied. In a simplified investigation carried out using the viscous Burgers equation the authors attributed the poorer numerical results of finite difference models to a violation of Galilean invariance in the discretization and propose to carry out the computations in a reference frame moving with the bulk velocity of the flow. Here we further discuss this problem and relate it to known results on invariant discretization schemes. Non-invariant and invariant finite difference discretizations of Burgers equation are proposed and compared with the discretization using the remedy proposed by Bernardini et al.
Hwang, Minki; Song, Jun-Seop; Lee, Young-Seon; Joung, Boyoung; Pak, Hui-Nam
2016-01-01
Background We previously reported that stable rotors were observed in in-silico human atrial fibrillation (AF) models, and were well represented by dominant frequency (DF). We explored the spatiotemporal stability of DF sites in 3D-AF models imported from patient CT images of the left atrium (LA). Methods We integrated 3-D CT images of the LA obtained from ten patients with persistent AF (male 80%, 61.8 ± 13.5 years old) into an in-silico AF model. After induction, we obtained 6 seconds of AF simulation data for DF analyses in 30 second intervals (T1–T9). The LA was divided into ten sections. Spatiotemporal changes and variations in the temporal consistency of DF were evaluated at each section of the LA. The high DF area was defined as the area with the highest 10% DF. Results 1. There was no spatial consistency in the high DF distribution at each LA section during T1–T9 except in one patient (p = 0.027). 2. Coefficients of variation for the high DF area were highly different among the ten LA sections (p < 0.001), and they were significantly higher in the four pulmonary vein (PV) areas, the LA appendage, and the peri-mitral area than in the other LA sections (p < 0.001). 3. When we conducted virtual ablation of 10%, 15%, and 20% of the highest DF areas (n = 270 cases), AF was changed to atrial tachycardia (AT) or terminated at a rate of 40%, 57%, and 76%, respectively. Conclusions Spatiotemporal consistency of the DF area was observed in 10% of AF patients, and high DF areas were temporally variable. Virtual ablation of DF is moderately effective in AF termination and AF changing into AT. PMID:27459377
Numerical modelling of morphodynamics—Vilaine Estuary
NASA Astrophysics Data System (ADS)
Vested, Hans Jacob; Tessier, Caroline; Christensen, Bo Brahtz; Goubert, Evelyne
2013-04-01
The main objective of this paper is to develop a method to simulate long-term morphodynamics of estuaries dominated by fine sediments, which are subject to both tidal flow and meteorologically induced variations in freshwater run-off and wave conditions. The method is tested on the Vilaine Estuary located in South Brittany, France. The estuary is subject to a meso-macrotidal regime. The semi-diurnal tidal range varies from around 2.5 to 5 m at neap and spring, respectively. The freshwater input is controlled by a dam located approximately 8 km from the mouth of the estuary. Sediments are characterised as mostly fines, but more sandy areas are also found. The morphology of the estuary is highly influenced by the dam. It is very dynamic and changes in a complicated manner with the run-off from the dam, the tide and the wave forcing at the mouth of the estuary. Extensive hydrodynamic and sediment field data have been collected in the past and provide a solid scientific basis for studying the estuary. Based on a conceptual understanding of the morphodynamics, a numerical morphological model with coupled hydrodynamic, surface wave and sediment transport models is formulated. The numerical models are calibrated to reproduce sediment concentrations, tidal flat altimetry and overall sediment fluxes. Scaling factors are applied to a reference year to form quasi-realistic hydrodynamic forcing and river run-off, which allow for the simulations to be extended to other years. The simulation results are compared with observed bathymetric changes in the estuary during the period 1998-2005. The models and scaling factors are applied to predict the morphological development over a time scale of up to 10 years. The influence of the initial conditions and the sequence of external hydrodynamic forcing, with respect to the morphodynamic response of the estuary, are discussed.
DANA: distributed numerical and adaptive modelling framework.
Rougier, Nicolas P; Fix, Jérémy
2012-01-01
DANA is a python framework ( http://dana.loria.fr ) whose computational paradigm is grounded on the notion of a unit that is essentially a set of time dependent values varying under the influence of other units via adaptive weighted connections. The evolution of a unit's value are defined by a set of differential equations expressed in standard mathematical notation which greatly ease their definition. The units are organized into groups that form a model. Each unit can be connected to any other unit (including itself) using a weighted connection. The DANA framework offers a set of core objects needed to design and run such models. The modeler only has to define the equations of a unit as well as the equations governing the training of the connections. The simulation is completely transparent to the modeler and is handled by DANA. This allows DANA to be used for a wide range of numerical and distributed models as long as they fit the proposed framework (e.g. cellular automata, reaction-diffusion system, decentralized neural networks, recurrent neural networks, kernel-based image processing, etc.). PMID:22994650
Dutta, Debargh K; Potnis, Pushya A; Rhodes, Kelly; Wood, Steven C
2015-01-01
Multinucleate giant cells (MGCs) are formed by the fusion of 5 to 15 monocytes or macrophages. MGCs can be generated by hip implants at the site where the metal surface of the device is in close contact with tissue. MGCs play a critical role in the inflammatory processes associated with adverse events such as aseptic loosening of the prosthetic joints and bone degeneration process called osteolysis. Upon interaction with metal wear particles, endothelial cells upregulate pro-inflammatory cytokines and other factors that enhance a localized immune response. However, the role of endothelial cells in the generation of MGCs has not been completely investigated. We developed a three-dimensional peripheral tissue-equivalent model (PTE) consisting of collagen gel, supporting a monolayer of endothelial cells and human peripheral blood mononuclear cells (PBMCs) on top, which mimics peripheral tissue under normal physiological conditions. The cultures were incubated for 14 days with Cobalt chromium alloy (CoCr ASTM F75, 1-5 micron) wear particles. PBMC were allowed to transit the endothelium and harvested cells were analyzed for MGC generation via flow cytometry. An increase in forward scatter (cell size) and in the propidium iodide (PI) uptake (DNA intercalating dye) was used to identify MGCs. Our results show that endothelial cells induce the generation of MGCs to a level 4 fold higher in 3-dimentional PTE system as compared to traditional 2-dimensional culture plates. Further characterization of MGCs showed upregulated expression of tartrate resistant alkaline phosphatase (TRAP) and dendritic cell specific transmembrane protein, (DC-STAMP), which are markers of bone degrading cells called osteoclasts. In sum, we have established a robust and relevant model to examine MGC and osteoclast formation in a tissue like environment using flow cytometry and RT-PCR. With endothelial cells help, we observed a consistent generation of metal wear particle- induced MGCs, which heralds
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
Numerical Modeling of Ocular Dysfunction in Space
NASA Technical Reports Server (NTRS)
Nelson, Emily S.; Mulugeta, Lealem; Vera, J.; Myers, J. G.; Raykin, J.; Feola, A. J.; Gleason, R.; Samuels, B.; Ethier, C. R.
2014-01-01
Upon introduction to microgravity, the near-loss of hydrostatic pressure causes a marked cephalic (headward) shift of fluid in an astronaut's body. The fluid shift, along with other factors of spaceflight, induces a cascade of interdependent physiological responses which occur at varying time scales. Long-duration missions carry an increased risk for the development of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, a spectrum of ophthalmic changes including posterior globe flattening, choroidal folds, distension of the optic nerve sheath, kinking of the optic nerve and potentially permanent degradation of visual function. In the cases of VIIP found to date, the initial onset of symptoms occurred after several weeks to several months of spaceflight, by which time the gross bodily fluid distribution is well established. We are developing a suite of numerical models to simulate the effects of fluid shift on the cardiovascular, central nervous and ocular systems. These models calculate the modified mean volumes, flow rates and pressures that are characteristic of the altered quasi-homeostatic state in microgravity, including intracranial and intraocular pressures. The results of the lumped models provide initial and boundary data to a 3D finite element biomechanics simulation of the globe, optic nerve head and retrobulbar subarachnoid space. The integrated set of models will be used to investigate the evolution of the biomechanical stress state in the ocular tissues due to long-term exposure to microgravity.
Numerical modeling of transport barrier formation
Tokar, Mikhail Z.
2010-04-01
In diverse media the characteristics of mass and heat transfer may undergo spontaneous and abrupt changes in time and space. This can lead to the formation of regions with strongly reduced transport, so called transport barriers (TB). The presence of interfaces between regions with qualitatively and quantitatively different transport characteristics impose severe requirements to methods and numerical schemes used by solving of transport equations. In particular the assumptions made in standard methods about the solution behavior by representing its derivatives fail in points where the transport changes abruptly. The situation is complicated further by the fact that neither the formation time nor the positions of interfaces are known a priori. A numerical approach, operating reliably under such conditions, is proposed. It is based on the introduction of a new dependent variable related to the variation after one time step of the original one integrated over the volume. In the vicinity of any grid knot the resulting differential equation is approximated by a second order ordinary differential equation with constant coefficients. Exact analytical solutions of these equations are conjugated between knots by demanding the continuity of the total solution and its first derivative. As an example the heat transfer in media with heat conductivity decreasing abruptly when the temperature e-folding length exceeds a critical value is considered. The formation of TB both at a heating power above the critical level and caused with radiation energy losses non-linearly dependent on the temperature is modeled.
Transient Numerical Modeling of Catalytic Channels
NASA Technical Reports Server (NTRS)
Struk, Peter M.; Dietrich, Daniel L.; Miller, Fletcher J.; T'ien, James S.
2007-01-01
This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the
Modeling Biodegradation and Reactive Transport: Analytical and Numerical Models
Sun, Y; Glascoe, L
2005-06-09
The computational modeling of the biodegradation of contaminated groundwater systems accounting for biochemical reactions coupled to contaminant transport is a valuable tool for both the field engineer/planner with limited computational resources and the expert computational researcher less constrained by time and computer power. There exists several analytical and numerical computer models that have been and are being developed to cover the practical needs put forth by users to fulfill this spectrum of computational demands. Generally, analytical models provide rapid and convenient screening tools running on very limited computational power, while numerical models can provide more detailed information with consequent requirements of greater computational time and effort. While these analytical and numerical computer models can provide accurate and adequate information to produce defensible remediation strategies, decisions based on inadequate modeling output or on over-analysis can have costly and risky consequences. In this chapter we consider both analytical and numerical modeling approaches to biodegradation and reactive transport. Both approaches are discussed and analyzed in terms of achieving bioremediation goals, recognizing that there is always a tradeoff between computational cost and the resolution of simulated systems.
Numerical linearized MHD model of flapping oscillations
NASA Astrophysics Data System (ADS)
Korovinskiy, D. B.; Ivanov, I. B.; Semenov, V. S.; Erkaev, N. V.; Kiehas, S. A.
2016-06-01
Kink-like magnetotail flapping oscillations in a Harris-like current sheet with earthward growing normal magnetic field component Bz are studied by means of time-dependent 2D linearized MHD numerical simulations. The dispersion relation and two-dimensional eigenfunctions are obtained. The results are compared with analytical estimates of the double-gradient model, which are found to be reliable for configurations with small Bz up to values ˜ 0.05 of the lobe magnetic field. Coupled with previous results, present simulations confirm that the earthward/tailward growth direction of the Bz component acts as a switch between stable/unstable regimes of the flapping mode, while the mode dispersion curve is the same in both cases. It is confirmed that flapping oscillations may be triggered by a simple Gaussian initial perturbation of the Vz velocity.
Two numerical models for landslide dynamic analysis
NASA Astrophysics Data System (ADS)
Hungr, Oldrich; McDougall, Scott
2009-05-01
Two microcomputer-based numerical models (Dynamic ANalysis (DAN) and three-dimensional model DAN (DAN3D)) have been developed and extensively used for analysis of landslide runout, specifically for the purposes of practical landslide hazard and risk assessment. The theoretical basis of both models is a system of depth-averaged governing equations derived from the principles of continuum mechanics. Original features developed specifically during this work include: an open rheological kernel; explicit use of tangential strain to determine the tangential stress state within the flowing sheet, which is both more realistic and beneficial to the stability of the model; orientation of principal tangential stresses parallel with the direction of motion; inclusion of the centripetal forces corresponding to the true curvature of the path in the motion direction and; the use of very simple and highly efficient free surface interpolation methods. Both models yield similar results when applied to the same sets of input data. Both algorithms are designed to work within the semi-empirical framework of the "equivalent fluid" approach. This approach requires selection of material rheology and calibration of input parameters through back-analysis of real events. Although approximate, it facilitates simple and efficient operation while accounting for the most important characteristics of extremely rapid landslides. The two models have been verified against several controlled laboratory experiments with known physical basis. A large number of back-analyses of real landslides of various types have also been carried out. One example is presented. Calibration patterns are emerging, which give a promise of predictive capability.
Numerical modelling of ion transport in flames
NASA Astrophysics Data System (ADS)
Han, Jie; Belhi, Memdouh; Bisetti, Fabrizio; Mani Sarathy, S.
2015-11-01
This paper presents a modelling framework to compute the diffusivity and mobility of ions in flames. The (n, 6, 4) interaction potential is adopted to model collisions between neutral and charged species. All required parameters in the potential are related to the polarizability of the species pair via semi-empirical formulas, which are derived using the most recently published data or best estimates. The resulting framework permits computation of the transport coefficients of any ion found in a hydrocarbon flame. The accuracy of the proposed method is evaluated by comparing its predictions with experimental data on the mobility of selected ions in single-component neutral gases. Based on this analysis, the value of a model constant available in the literature is modified in order to improve the model's predictions. The newly determined ion transport coefficients are used as part of a previously developed numerical approach to compute the distribution of charged species in a freely propagating premixed lean CH4/O2 flame. Since a significant scatter of polarizability data exists in the literature, the effects of changes in polarizability on ion transport properties and the spatial distribution of ions in flames are explored. Our analysis shows that changes in polarizability propagate with decreasing effect from binary transport coefficients to species number densities. We conclude that the chosen polarizability value has a limited effect on the ion distribution in freely propagating flames. We expect that the modelling framework proposed here will benefit future efforts in modelling the effect of external voltages on flames. Supplemental data for this article can be accessed at http://dx.doi.org/10.1080/13647830.2015.1090018.
Numerical Modeling of Suspension HVOF Spray
NASA Astrophysics Data System (ADS)
Jadidi, M.; Moghtadernejad, S.; Dolatabadi, A.
2016-02-01
A three-dimensional two-way coupled Eulerian-Lagrangian scheme is used to simulate suspension high-velocity oxy-fuel spraying process. The mass, momentum, energy, and species equations are solved together with the realizable k-ɛ turbulence model to simulate the gas phase. Suspension is assumed to be a mixture of solid particles [mullite powder (3Al2O3·2SiO2)], ethanol, and ethylene glycol. The process involves premixed combustion of oxygen-propylene, and non-premixed combustion of oxygen-ethanol and oxygen-ethylene glycol. One-step global reaction is used for each mentioned reaction together with eddy dissipation model to compute the reaction rate. To simulate the droplet breakup, Taylor Analogy Breakup model is applied. After the completion of droplet breakup, and solvent evaporation/combustion, the solid suspended particles are tracked through the domain to determine the characteristics of the coating particles. Numerical simulations are validated against the experimental results in the literature for the same operating conditions. Seven or possibly eight shock diamonds are captured outside the nozzle. In addition, a good agreement between the predicted particle temperature, velocity, and diameter, and the experiment is obtained. It is shown that as the standoff distance increases, the particle temperature and velocity reduce. Furthermore, a correlation is proposed to determine the spray cross-sectional diameter and estimate the particle trajectories as a function of standoff distance.
Foehn wind detection using numerical modelling
NASA Astrophysics Data System (ADS)
Irimescu, A.; Caian, M.
2010-09-01
In Romania, foehn is a short-lived atmospheric phenomenon, of a low to average intensity, not always highlighted by weather station observations. When such situations occur additional data are resorted to, rendering a continuous, aggregate image, in comparison to the punctual information yielded by weather stations. This paper aims to describe how foehn is detected in northern Oltenia (the Inner Carpathian-Balkan Curvature), using numerical modelling. Results generated by the RegCM3 Regional Climatic Model thus represent an undisputed tool, their most important advantage being the 10-km spatial resolution. The presence of foehn in northern Oltenia and its climatic peculiarities have been disclosed through the analysis in time and space of the meteorological elements specific to the phenomenon (air temperature, wind speed and direction etc) over a 40-year interval (1961-2000). The paper presents a new methodology that can be used to estimate the probability of production and the foehn characteristics (intensity, duration etc.). Interpretation of the RegCM3 model results has led to the statistical analysis of foehn occurrences within the studied area during the cold season (December, January and February). The resulted climatology, with fine resolution, can be used in foehn forecast of predictability.
Numerical modeling of bubble dynamics in magmas
NASA Astrophysics Data System (ADS)
Huber, Christian; Su, Yanqing; Parmigiani, Andrea
2014-05-01
Understanding the complex non-linear physics that governs volcanic eruptions is contingent on our ability to characterize the dynamics of bubbles and its effect on the ascending magma. The exsolution and migration of bubbles has also a great impact on the heat and mass transport in and out of magma bodies stored at shallow depths in the crust. Multiphase systems like magmas are by definition heterogeneous at small scales. Although mixture theory or homogenization methods are convenient to represent multiphase systems as a homogeneous equivalent media, these approaches do not inform us on possible feedbacks at the pore-scale and can be significantly misleading. In this presentation, we discuss the development and application of bubble-scale multiphase flow modeling to address the following questions : How do bubbles impact heat and mass transport in magma chambers ? How efficient are chemical exchanges between the melt and bubbles during magma decompression? What is the role of hydrodynamic interactions on the deformation of bubbles while the magma is sheared? Addressing these questions requires powerful numerical methods that accurately model the balance between viscous, capillary and pressure stresses. We discuss how these bubble-scale models can provide important constraints on the dynamics of magmas stored at shallow depth or ascending to the surface during an eruption.
Benchmarking numerical freeze/thaw models
NASA Astrophysics Data System (ADS)
Rühaak, Wolfram; Anbergen, Hauke; Molson, John; Grenier, Christophe; Sass, Ingo
2015-04-01
The modeling of freezing and thawing of water in porous media is of increasing interest, and for which very different application areas exist. For instance, the modeling of permafrost regression with respect to climate change issues is one area, while others include geotechnical applications in tunneling and for borehole heat exchangers which operate at temperatures below the freezing point. The modeling of these processes requires the solution of a coupled non-linear system of partial differential equations for flow and heat transport in space and time. Different code implementations have been developed in the past. Analytical solutions exist only for simple cases. Consequently, an interest has arisen in benchmarking different codes with analytical solutions, experiments and purely numerical results, similar to the long-standing DECOVALEX and the more recent "Geothermal Code Comparison" activities. The name for this freezing/ thawing benchmark consortium is INTERFROST. In addition to the well-known so-called Lunardini solution for a 1D case (case T1), two different 2D problems will be presented, one which represents melting of a frozen inclusion (case TH2) and another which represents the growth or thaw of permafrost around a talik (case TH3). These talik regions are important for controlling groundwater movement within a mainly frozen ground. First results of the different benchmark results will be shown and discussed.
Numerical modeling of polar mesocyclones generation mechanisms
NASA Astrophysics Data System (ADS)
Sergeev, Dennis; Stepanenko, Victor
2013-04-01
parameters, lateral boundary conditions are varied in the typically observed range. The approach is fully nonlinear: we use a three-dimensional non-hydrostatic mesoscale model NH3D_MPI [1] coupled with one-dimensional water body model LAKE. A key method used in the present study is the analysis of eddy kinetic and available potential energy budgets. References 1. Mikushin, D.N., and Stepanenko, V.M., The implementation of regional atmospheric model numerical algorithms for CBEA-based clusters. Lecture Notes in Computer Science, Parallel Processing and Applied Mathematics, 2010, vol. 6067, p. 525-534. 2. Rasmussen, E., and Turner, J. (eds), Polar Lows: Mesoscale Weather Systems in the Polar Regions. Cambridge: Cambridge University Press, 2003, 612 pp. 3. Yanase, W., and Niino, H., Dependence of Polar Low Development on Baroclinicity and Physical Processes: An Idealized High-Resolution Experiment, J. Atmos. Sci., 2006, vol. 64, p. 3044-3067.
Numerical geometry of map and model assessment.
Wriggers, Willy; He, Jing
2015-11-01
We are describing best practices and assessment strategies for the atomic interpretation of cryo-electron microscopy (cryo-EM) maps. Multiscale numerical geometry strategies in the Situs package and in secondary structure detection software are currently evolving due to the recent increases in cryo-EM resolution. Criteria that aim to predict the accuracy of fitted atomic models at low (worse than 8Å) and medium (4-8 Å) resolutions remain challenging. However, a high level of confidence in atomic models can be achieved by combining such criteria. The observed errors are due to map-model discrepancies and due to the effect of imperfect global docking strategies. Extending the earlier motion capture approach developed for flexible fitting, we use simulated fiducials (pseudoatoms) at varying levels of coarse-graining to track the local drift of structural features. We compare three tracking approaches: naïve vector quantization, a smoothly deformable model, and a tessellation of the structure into rigid Voronoi cells, which are fitted using a multi-fragment refinement approach. The lowest error is an upper bound for the (small) discrepancy between the crystal structure and the EM map due to different conditions in their structure determination. When internal features such as secondary structures are visible in medium-resolution EM maps, it is possible to extend the idea of point-based fiducials to more complex geometric representations such as helical axes, strands, and skeletons. We propose quantitative strategies to assess map-model pairs when such secondary structure patterns are prominent. PMID:26416532
Numerical Models of Ophiolite Genesis and Obduction
NASA Astrophysics Data System (ADS)
Guilmette, C.; Beaumont, C.; Jamieson, R.
2013-12-01
Ophiolites are relics of oceanic lithosphere tectonically emplaced in continental settings. They are diagnostic features of continental suture zones, where they mark past plate boundaries. Even after having been studied for more than 40 years, the mechanisms involved in the genesis and subsequent obduction of ophiolites over continental margins are still debated. We present the results of 2D thermal-mechanical numerical models that successfully reproduce characteristics of natural examples like the Semail, Bay of Islands, Yarlung-Zangbo, and Coast Range ophiolites. The numerical models are upper mantle scale and use pressure-, temperature- and strain-dependent viscous-plastic rheologies. Both divergent and convergent velocity boundary conditions are used and tectonic boundary forces are monitored. The models start with the rifting of a stable continent, followed by development of an ocean ridge and accretion of oceanic lithosphere at a total rate of 3 cm/y. Once a specified ocean size/age is achieved, the velocity boundary conditions are reversed leading to convergence and the spontaneous inception of a suduction zone at the mid-ocean ridge. We present results for models including different ages of oceans (40 to 90 Ma) and different convergence velocities (5 to 15 cm/y). The interaction between the lower plate passive margin and the oceanic upper plate results in 5 different tectonic styles. These differ mainly by the presence or absence of oceanic spreading in the upper plate (back-arc basin), leading to supra-subduction zone ophiolites vs. MORB-type, and by the behaviour of the oceanic slab, e.g., slab rollback vs. breakoff. The evolution of effective slab pull is interpreted to be the major control on the resulting tectonic style. Low effective slab pull models (young oceans and fast convergence rates) fail to obduct an ophiolite. Strong effective slab pull models (old oceans and lower convergence rates) result in subduction zone retreat and spontaneous oceanic
3-dimensional imaging at nanometer resolutions
Werner, James H.; Goodwin, Peter M.; Shreve, Andrew P.
2010-03-09
An apparatus and method for enabling precise, 3-dimensional, photoactivation localization microscopy (PALM) using selective, two-photon activation of fluorophores in a single z-slice of a sample in cooperation with time-gated imaging for reducing the background radiation from other image planes to levels suitable for single-molecule detection and spatial location, are described.
Numerical models of wind-driven circulation in lakes
Cheng, R.T.; Powell, T.M.; Dillon, T.M.
1976-01-01
The state-of-the-art of numerical modelling of large-scale wind-driven circulation in lakes is presented. The governing equations which describe this motion are discussed along with the appropriate numerical techniques necessary to solve them in lakes. The numerical models are categorized into three large primary groups: the layered models, the Ekman-type models, and the other three-dimensional models. Discussions and comparison of models are given and future research directions are suggested. ?? 1976.
Numerical modeling of subaqueous sand dune morphodynamics
NASA Astrophysics Data System (ADS)
Doré, Arnaud; Bonneton, Philippe; Marieu, Vincent; Garlan, Thierry
2016-03-01
The morphodynamic evolution of subaqueous sand dunes is investigated, using a 2-D Reynolds-averaged Navier-Stokes numerical model. A laboratory experiment where dunes are generated under stationary unidirectional flow conditions is used as a reference case. The model reproduces the evolution of the erodible bed until a state of equilibrium is reached. In particular, the simulation exhibits the different stages of the bed evolution, e.g., the incipient ripple generation, the nonlinear bed form growing phase, and the dune field equilibrium phase. The results show good agreement in terms of dune geometrical dimensions and time to equilibrium. After the emergence of the first ripple field, the bed growth is driven by cascading merging sequences between bed forms of different heights. A sequence extracted from the simulation shows how the downstream bed form is first eroded before merging with the upstream bed form. Superimposed bed forms emerge on the dune stoss sides during the simulation. An analysis of the results shows that they emerge downstream of a slight deflection on the dune profile. The deflection arises due to a modification of the sediment flux gradient consecutive to a reduction in the turbulence relaxation length while the upstream bed form height decreases. As they migrate, superimposed bed forms grow on the dune stoss side and eventually provoke the degeneration of the dune crest. Cascading merging sequences and superimposed bed forms dynamics both influence the dune field evolution and size and therefore play a fundamental role in the dune field self-organization process.
Numerical modeling of fluidic flow meters
NASA Astrophysics Data System (ADS)
Choudhury, D.; Patel, B. R.
1992-05-01
The transient fluid flow in fluidic flow meters has been modeled using Creare.x's flow modeling computer program FLUENT/BFC that solves the Navier-Stokes equations in general curvilinear coordinates. The numerical predictions of fluid flow in a fluidic flow meter have been compared with the available experimental results for a particular design, termed the PC-4 design. Overall flow structures such as main jet bending, and primary and secondary vortices predicted by FLUENT/BFC are in excellent agreement with flow visualization results. The oscillation frequencies of the PC-4 design have been predicted for a range of flow rates encompassing laminar and turbulent flow and the results are in good agreement with experiments. The details of the flow field predictions reveal that an important factor that determines the onset of oscillations in the fluidic flow meter is the feedback jet momentum relative to the main jet momentum. The insights provided by the analysis of the PC-4 fluidic flow meter design have led to an improved design. The improved design has sustained oscillations at lower flow rates compared with the PC-4 design and has a larger rangeability.
Understanding Etna flank instability through numerical models
NASA Astrophysics Data System (ADS)
Apuani, Tiziana; Corazzato, Claudia; Merri, Andrea; Tibaldi, Alessandro
2013-02-01
As many active volcanoes, Mount Etna shows clear evidence of flank instability, and different mechanisms were suggested to explain this flank dynamics, based on the recorded deformation pattern and character. Shallow and deep deformations, mainly associated with both eruptive and seismic events, are concentrated along recognised fracture and fault systems, mobilising the eastern and south-eastern flank of the volcano. Several interacting causes were postulated to control the phenomenon, including gravity force, magma ascent along the feeding system, and a very complex local and/or regional tectonic activity. Nevertheless, the complexity of such dynamics is still an open subject of research and being the volcano flanks heavily urbanised, the comprehension of the gravitative dynamics is a major issue for public safety and civil protection. The present research explores the effects of the main geological features (in particular the role of the subetnean clays, interposed between the Apennine-Maghrebian flysch and the volcanic products) and the role of weakness zones, identified by fracture and fault systems, on the slope instability process. The effects of magma intrusions are also investigated. The problem is addressed by integrating field data, laboratory tests and numerical modelling. A bi- and tri-dimensional stress-strain analysis was performed by a finite difference numerical code (FLAC and FLAC3D), mainly aimed at evaluating the relationship among geological features, volcano-tectonic structures and magmatic activity in controlling the deformation processes. The analyses are well supported by dedicated structural-mechanical field surveys, which allowed to estimate the rock mass strength and deformability parameters. To take into account the uncertainties which inevitably occur in a so complicated model, many efforts were done in performing a sensitivity analysis along a WNW-ESE section crossing the volcano summit and the Valle del Bove depression. This was
Numerical model of circumpolar Antarctic ice shelves
Johnson, R.C.
1985-01-01
Extensive floating ice shelves in the Antarctic have been proposed to explain the discrepancies between Pleistocene high sea levels shown by dated coral reefs and coeval low sea levels inferred from glacial ice volumes calculated from oxygen isotope ratios in deep sea cores. A numerical model using the floating shelf creep analysis of Weertman (1957) has provided a plausible basis for the acceptance of such shelves. Shelf outer limits were set at 55/sup 0/S in East Antarctica and 58/sup 0/S in West Antarctica, based in part on diatom-deficient deep sea sediments deposited prior to the Holocene. Precipitation varied from 10 gm cm/sup -2/yr/sup -1/ at 75/sup 0/S to 80 gm cm/sup -2/yr/sup -1/ at 55/sup 0/S. Mean air temperatures varied from -35/sup 0/C at the 75/sup 0/S coast to -17/sup 0/C at the outer limits. Isotope ratios were those of present Antarctic precipitation at corresponding model shelf temperatures. In the calculation, a steady state is assumed. Integration begins at the coast with summation over successive years as creep and continental ice discharge move the integration element to the outer limits. The oceanic oxygen isotope ratio change required by the discrepancies in the record is 0.40 to 0.50 ppmil. Using the flow law constant of 4.2 and a creep activation energy of 134 kjoules mol/sup -1/, the resulting change is 0.44 ppmil. Difference results reflect the uncertainties associated with the critical creep constants used in the modeling. Nevertheless, the results suggest that a quantity of Antarctic shelf ice comparable to ice volumes in major Northern glacial areas existed at times during the Pleistocene.
Numerical Modelling of Mesoscale Atmospheric Dispersion.
NASA Astrophysics Data System (ADS)
Moran, Michael D.
Mesoscale atmospheric dispersion is more complicated than smaller-scale dispersion because the mean wind field can no longer be considered steady or horizontally homogeneous over mesoscale time and space scales. Wind shear also plays a more important role on the mesoscale, and horizontal dispersion can be enhanced and even dominated by vertical wind shear through either the simultaneous or delayed interaction of horizontal differential advection and vertical mixing over one or two diurnal periods. The CSU mesoscale atmospheric dispersion modelling system has been used in this study to simulate the transport and diffusion of a perfluorocarbon gas for episodic releases made during two North American mesoscale dispersion field experiments, the 1980 Great Plains tracer experiment and the 1983 Cross-Appalachian Tracer Experiment (CAPTEX). Ground -level and elevated tracer concentrations were measured out to distances of 600 km from the source in the first experiment and 1100 km in the second. The physiography of the two experimental domains was quite different, permitting isolation and examination of the roles of terrain forcing and differential advection in mesoscale atmospheric dispersion. Suites of numerical experiments of increasing complexity were carried out for both case studies. The experiments differed in the realism of their representation of both the synoptic-scale flow and the underlying terrain. The Great Plains nocturnal low-level jet played an important role in the first case while temporal changes in the synoptic -scale flow were very significant in the second case. The contributions of differential advection and mesoscale deformation to mesoscale dispersion dominated those of small-scale turbulent diffusion for both cases, and Pasquill's (1962) delayed-shear-enhancement mechanism for lateral dispersion was found to be particularly important. This study was also the first quantitative evaluation of the CSU mesoscale dispersion modelling system with
The 3-dimensional cellular automata for HIV infection
NASA Astrophysics Data System (ADS)
Mo, Youbin; Ren, Bin; Yang, Wencao; Shuai, Jianwei
2014-04-01
The HIV infection dynamics is discussed in detail with a 3-dimensional cellular automata model in this paper. The model can reproduce the three-phase development, i.e., the acute period, the asymptotic period and the AIDS period, observed in the HIV-infected patients in a clinic. We show that the 3D HIV model performs a better robustness on the model parameters than the 2D cellular automata. Furthermore, we reveal that the occurrence of a perpetual source to successively generate infectious waves to spread to the whole system drives the model from the asymptotic state to the AIDS state.
Precise numerical modeling of next generation multimode fiber based links
NASA Astrophysics Data System (ADS)
Maksymiuk, L.; Stepniak, G.
2015-12-01
In order to numerically model modern multimode fiber based links we are required to take into account modal and chromatic dispersion, profile dispersion and spectral dependent coupling. In this paper we propose a complete numerical model which not only is precise but also versatile. Additionally to the detailed mathematical description of the model we provide also a bunch of numerical calculations performed with the use of the model.
Biochemical Applications Of 3-Dimensional Fluorescence Spectrometry
NASA Astrophysics Data System (ADS)
Leiner, Marc J.; Wolfbeis, Otto S.
1988-06-01
We investigated the 3-dimensional fluorescence of complex mixtures of bioloquids such as human serum, serum ultrafiltrate, human urine, and human plasma low density lipoproteins. The total fluorescence of human serum can be divided into a few peaks. When comparing fluorescence topograms of sera, from normal and cancerous subjects, we found significant differences in tryptophan fluorescence. Although the total fluorescence of human urine can be resolved into 3-5 distinct peaks, some of them. do not result from single fluorescent urinary metabolites, but rather from. several species having similar spectral properties. Human plasma, low density lipoproteins possess a native fluorescence that changes when submitted to in-vitro autoxidation. The 3-dimensional fluorescence demonstrated the presence of 7 fluorophores in the lipid domain, and 6 fluorophores in the protein. dovain- The above results demonstrated that 3-dimensional fluorescence can resolve the spectral properties of complex ,lxtures much better than other methods. Moreover, other parameters than excitation and emission wavelength and intensity (for instance fluorescence lifetime, polarization, or quenchability) may be exploited to give a multidl,ensio,a1 matrix, that is unique for each sample. Consequently, 3-dimensio:Hhal fluorescence as such, or in combination with separation techniques is therefore considered to have the potential of becoming a useful new H.ethod in clinical chemistry and analytical biochemistry.
Numerical models for high beta magnetohydrodynamic flow
Brackbill, J.U.
1987-01-01
The fundamentals of numerical magnetohydrodynamics for highly conducting, high-beta plasmas are outlined. The discussions emphasize the physical properties of the flow, and how elementary concepts in numerical analysis can be applied to the construction of finite difference approximations that capture these features. The linear and nonlinear stability of explicit and implicit differencing in time is examined, the origin and effect of numerical diffusion in the calculation of convective transport is described, and a technique for maintaining solenoidality in the magnetic field is developed. Many of the points are illustrated by numerical examples. The techniques described are applicable to the time-dependent, high-beta flows normally encountered in magnetically confined plasmas, plasma switches, and space and astrophysical plasmas. 40 refs.
Large-signal numerical and analytical HBT models
Teeter, D.A.; East, J.R.; Mains, R.K.; Haddad, G.I. )
1993-05-01
Several large-signal HBT models are investigated in this paper to determine their usefulness at millimeter-wave frequencies. The most detailed model involves numerically solving moments of the Boltzmann Transport Equation. A description of the numerical model is given along with several simulated results. The numerical model is then used to evaluate two analytical HBT models, the conventional Gummel-Poon model and a modified Ebers-Moll model. It is found that the commonly used Gummel-Poon model exhibits poor agreement with numerical and experimental data at millimeter-wave frequencies due to neglect of transit-time delays. Improved agreement between measured and modeled data result by including transit-time effects in an Ebers-Moll model. This simple model has direct application to millimeter-wave power amplifier and oscillator design. Several measured results are presented to help verify the simple model.
Validation of Numerical Shallow Water Models for Tidal Lagoons
Eliason, D.; Bourgeois, A.
1999-11-01
An analytical solution is presented for the case of a stratified, tidally forced lagoon. This solution, especially its energetics, is useful for the validation of numerical shallow water models under stratified, tidally forced conditions. The utility of the analytical solution for validation is demonstrated for a simple finite difference numerical model. A comparison is presented of the energetics of the numerical and analytical solutions in terms of the convergence of model results to the analytical solution with increasing spatial and temporal resolution.
Advanced in turbulence physics and modeling by direct numerical simulations
NASA Technical Reports Server (NTRS)
Reynolds, W. C.
1987-01-01
The advent of direct numerical simulations of turbulence has opened avenues for research on turbulence physics and turbulence modeling. Direct numerical simulation provides values for anything that the scientist or modeler would like to know about the flow. An overview of some recent advances in the physical understanding of turbulence and in turbulence modeling obtained through such simulations is presented.
Technology Transfer Automated Retrieval System (TEKTRAN)
When Lagrangian stochastic models for turbulent dispersion are applied to complex flows, some type of ad hoc intervention is almost always necessary to eliminate unphysical behavior in the numerical solution. This paper discusses numerical considerations when solving the Langevin-based particle velo...
Conceptual and Numerical Models for UZ Flow and Transport
H. Liu
2000-03-03
The purpose of this Analysis/Model Report (AMR) is to document the conceptual and numerical models used for modeling of unsaturated zone (UZ) fluid (water and air) flow and solute transport processes. This is in accordance with ''AMR Development Plan for U0030 Conceptual and Numerical Models for Unsaturated Zone (UZ) Flow and Transport Processes, Rev 00''. The conceptual and numerical modeling approaches described in this AMR are used for models of UZ flow and transport in fractured, unsaturated rock under ambient and thermal conditions, which are documented in separate AMRs. This AMR supports the UZ Flow and Transport Process Model Report (PMR), the Near Field Environment PMR, and the following models: Calibrated Properties Model; UZ Flow Models and Submodels; Mountain-Scale Coupled Processes Model; Thermal-Hydrologic-Chemical (THC) Seepage Model; Drift Scale Test (DST) THC Model; Seepage Model for Performance Assessment (PA); and UZ Radionuclide Transport Models.
NUMERICAL MODELING OF FINE SEDIMENT PHYSICAL PROCESSES.
Schoellhamer, David H.
1985-01-01
Fine sediment in channels, rivers, estuaries, and coastal waters undergo several physical processes including flocculation, floc disruption, deposition, bed consolidation, and resuspension. This paper presents a conceptual model and reviews mathematical models of these physical processes. Several general fine sediment models that simulate some of these processes are reviewed. These general models do not directly simulate flocculation and floc disruption, but the conceptual model and existing functions are shown to adequately model these two processes for one set of laboratory data.
Numerical Modeling in Geodynamics: Success, Failure and Perspective
NASA Astrophysics Data System (ADS)
Ismail-Zadeh, A.
2005-12-01
A real success in numerical modeling of dynamics of the Earth can be achieved only by multidisciplinary research teams of experts in geodynamics, applied and pure mathematics, and computer science. The success in numerical modeling is based on the following basic, but simple, rules. (i) People need simplicity most, but they understand intricacies best (B. Pasternak, writer). Start from a simple numerical model, which describes basic physical laws by a set of mathematical equations, and move then to a complex model. Never start from a complex model, because you cannot understand the contribution of each term of the equations to the modeled geophysical phenomenon. (ii) Study the numerical methods behind your computer code. Otherwise it becomes difficult to distinguish true and erroneous solutions to the geodynamic problem, especially when your problem is complex enough. (iii) Test your model versus analytical and asymptotic solutions, simple 2D and 3D model examples. Develop benchmark analysis of different numerical codes and compare numerical results with laboratory experiments. Remember that the numerical tool you employ is not perfect, and there are small bugs in every computer code. Therefore the testing is the most important part of your numerical modeling. (iv) Prove (if possible) or learn relevant statements concerning the existence, uniqueness and stability of the solution to the mathematical and discrete problems. Otherwise you can solve an improperly-posed problem, and the results of the modeling will be far from the true solution of your model problem. (v) Try to analyze numerical models of a geological phenomenon using as less as possible tuning model variables. Already two tuning variables give enough possibilities to constrain your model well enough with respect to observations. The data fitting sometimes is quite attractive and can take you far from a principal aim of your numerical modeling: to understand geophysical phenomena. (vi) If the number of
Numerical bifurcation analysis of immunological models with time delays
NASA Astrophysics Data System (ADS)
Luzyanina, Tatyana; Roose, Dirk; Bocharov, Gennady
2005-12-01
In recent years, a large number of mathematical models that are described by delay differential equations (DDEs) have appeared in the life sciences. To analyze the models' dynamics, numerical methods are necessary, since analytical studies can only give limited results. In turn, the availability of efficient numerical methods and software packages encourages the use of time delays in mathematical modelling, which may lead to more realistic models. We outline recently developed numerical methods for bifurcation analysis of DDEs and illustrate the use of these methods in the analysis of a mathematical model of human hepatitis B virus infection.
Numerically Controlled Machining Of Wind-Tunnel Models
NASA Technical Reports Server (NTRS)
Kovtun, John B.
1990-01-01
New procedure for dynamic models and parts for wind-tunnel tests or radio-controlled flight tests constructed. Involves use of single-phase numerical control (NC) technique to produce highly-accurate, symmetrical models in less time.
Numerical modeling of waveguide heated microwave plasmas
Venkateswaran, S.; Schwer, D.A.; Merkle, C.L.
1993-12-01
Waveguide-heated microwave plasmas for space propulsion applications are analyzed by a two-dimensional numerical solution of the combined Navier-Stokes and Maxwell equations. Two waveguide configurations -- one purely transmitting and the other with a reflecting end wall -- are considered. Plasma stability and absorption characteristics are studied and contrasted with the characteristic of resonant cavity heated plasmas. In addition, preliminary estimates of the overall efficiency and the thrust and specific impulse of the propulsion system are also made. The computational results are used to explain experimental trends and to better understand the working of these devices.
Software Simplifies the Sharing of Numerical Models
NASA Technical Reports Server (NTRS)
2014-01-01
To ease the sharing of climate models with university students, Goddard Space Flight Center awarded SBIR funding to Reston, Virginia-based Parabon Computation Inc., a company that specializes in cloud computing. The firm developed a software program capable of running climate models over the Internet, and also created an online environment for people to collaborate on developing such models.
Material model library for explicit numerical codes
Hofmann, R.; Dial, B.W.
1982-08-01
A material model logic structure has been developed which is useful for most explicit finite-difference and explicit finite-element Lagrange computer codes. This structure has been implemented and tested in the STEALTH codes to provide an example for researchers who wish to implement it in generically similar codes. In parallel with these models, material parameter libraries have been created for the implemented models for materials which are often needed in DoD applications.
Numerical MHD codes for modeling astrophysical flows
NASA Astrophysics Data System (ADS)
Koldoba, A. V.; Ustyugova, G. V.; Lii, P. S.; Comins, M. L.; Dyda, S.; Romanova, M. M.; Lovelace, R. V. E.
2016-05-01
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy conservation. The code has been implemented on several different coordinate systems: 2.5D axisymmetric cylindrical coordinates, 2D Cartesian coordinates, 2D plane polar coordinates, and fully 3D cylindrical coordinates. Viscosity and diffusivity are implemented in the code to control the accretion rate in the disk and the rate of penetration of the disk matter through the magnetic field lines. The code has been utilized for the numerical investigations of a number of different astrophysical problems, several examples of which are shown.
Cumulus clouds - Numerical models, observations and entrainment
NASA Technical Reports Server (NTRS)
Simpson, J.
1983-01-01
The first computer simulation of the organization phase of a buoyant atmospheric thermal is described. Although crude, it showed the spontaneous development of a rounded tight-gradient 'cap' and internal vortical circulation. The complexities involved in these 'field of motion' models in part motivated the development of entity models, based upon laboratory thermals. These one-dimensional models and their uses with observations are briefly described as well as their limitations. Finally, an application of Schlesinger's three-dimensional model to a GATE cumulus situation clarifies many apparently conflicting observations and postulates, thereby raising further challenging questions to be addressed jointly by the more sophisticated measuring and modeling tools available in the 1980's.
Numerical modeling of laser thermal propulsion flows
NASA Technical Reports Server (NTRS)
Mccay, T. D.; Thoenes, J.
1984-01-01
An review of the problems associated with modeling laser thermal propulsion flows, a synopsis of the status of such models, and the attributes of a successful model are presented. The continuous gaseous hydrogen laser-supported combustion wave (LSCW) thruster, for which a high-energy laser system (preferably space-based) should exist by the time the propulsion technology is developed, is considered in particular. The model proposed by Raizer (1970) is based on the assumptions of one-dimensional flow at constant pressure with heat conduction as the principal heat transfer mechanism. Consideration is given to subsequent models which account for radiative transfer into the ambient gas; provide a two-dimensional generalization of Raizer's analysis for the subsonic propagation of laser sparks in air; include the effect of forward plasma radiation in a one-dimensional model; and attempt a time-dependent (elliptic) solution of the full Navier-Stokes equations for the flow in a simple axisymmetric thruster. Attention is also given to thruster and nozzle flow models and thermodynamic and transport properties.
Paleoclimate validation of a numerical climate model
Schelling, F.J.; Church, H.W.; Zak, B.D.; Thompson, S.L.
1994-04-01
An analysis planned to validate regional climate model results for a past climate state at Yucca Mountain, Nevada, against paleoclimate evidence for the period is described. This analysis, which will use the GENESIS model of global climate nested with the RegCM2 regional climate model, is part of a larger study for DOE`s Yucca Mountain Site Characterization Project that is evaluating the impacts of long term future climate change on performance of the potential high level nuclear waste repository at Yucca Mountain. The planned analysis and anticipated results are presented.
Numerical modeling of fresh concrete flow through porous medium
NASA Astrophysics Data System (ADS)
Kolařík, F.; Patzák, B.; Zeman, J.
2016-06-01
The paper focuses on a numerical modeling of a non-Newtonian fluid flow in a porous domain. It presents combination of a homogenization approach to obtain permeability from the underlying micro-structure with coupling of a Stokes and Darcy flow through the interface on the macro level. As a numerical method we employed the Finite Element method. The results obtained from the homogenization approach are validated against fully resolved solution computed by direct numerical simulation.
Numerical modelling of fracture in human arteries.
Ferrara, A; Pandolfi, A
2008-10-01
We present 3D finite element models of atherosclerotic arteries, used to investigate the influence of the geometry and tissue properties on the plaque rupture caused by overexpansion. We adopted a geometry reconstructed from a contiguous set of in vitro magnetic resonance images of a damaged artery. The artery wall is divided in three layers (adventitia, media and intima) and is discretized into tetrahedral finite elements. The artery material is described with a hyperelastic two-fiber anisotropic model proposed by Holzapfel et al. 2000. A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elasticity 61(1):1-48, while the plaque is assumed to be transversely isotropic. Cracks induced by mechanical actions are represented through cohesive surfaces, and are allowed to develop along solid elements boundaries only. Fractures are explicitly introduced in the discretized model at the locations where the tensile strength of the material is reached. PMID:19230149
Numerical Calculation of Model Rocket Trajectories.
ERIC Educational Resources Information Center
Keeports, David
1990-01-01
Discussed is the use of model rocketry to teach the principles of Newtonian mechanics. Included are forces involved; calculations for vertical launches; two-dimensional trajectories; and variations in mass, drag, and launch angle. (CW)
On numerical modeling of one-dimensional geothermal histories
Haugerud, R.A.
1989-01-01
Numerical models of one-dimensional geothermal histories are one way of understanding the relations between tectonics and transient thermal structure in the crust. Such models can be powerful tools for interpreting geochronologic and thermobarometric data. A flexible program to calculate these models on a microcomputer is available and examples of its use are presented. Potential problems with this approach include the simplifying assumptions that are made, limitations of the numerical techniques, and the neglect of convective heat transfer. ?? 1989.
NASA Astrophysics Data System (ADS)
Vasquez, D. A.; Swift, J. N.; Tan, S.; Darrah, T. H.
2013-12-01
The integration of precise geochemical analyses with quantitative engineering modeling into an interactive GIS system allows for a sophisticated and efficient method of reservoir engineering and characterization. Geographic Information Systems (GIS) is utilized as an advanced technique for oil field reservoir analysis by combining field engineering and geological/geochemical spatial datasets with the available systematic modeling and mapping methods to integrate the information into a spatially correlated first-hand approach in defining surface and subsurface characteristics. Three key methods of analysis include: 1) Geostatistical modeling to create a static and volumetric 3-dimensional representation of the geological body, 2) Numerical modeling to develop a dynamic and interactive 2-dimensional model of fluid flow across the reservoir and 3) Noble gas geochemistry to further define the physical conditions, components and history of the geologic system. Results thus far include using engineering algorithms for interpolating electrical well log properties across the field (spontaneous potential, resistivity) yielding a highly accurate and high-resolution 3D model of rock properties. Results so far also include using numerical finite difference methods (crank-nicholson) to solve for equations describing the distribution of pressure across field yielding a 2D simulation model of fluid flow across reservoir. Ongoing noble gas geochemistry results will also include determination of the source, thermal maturity and the extent/style of fluid migration (connectivity, continuity and directionality). Future work will include developing an inverse engineering algorithm to model for permeability, porosity and water saturation.This combination of new and efficient technological and analytical capabilities is geared to provide a better understanding of the field geology and hydrocarbon dynamics system with applications to determine the presence of hydrocarbon pay zones (or
Modelling asteroid brightness variations. I - Numerical methods
NASA Technical Reports Server (NTRS)
Karttunen, H.
1989-01-01
A method for generating lightcurves of asteroid models is presented. The effects of the shape of the asteroid and the scattering law of a surface element are distinctly separable, being described by chosen functions that can easily be changed. The shape is specified by means of two functions that yield the length of the radius vector and the normal vector of the surface at a given point. The general shape must be convex, but spherical concavities producing macroscopic shadowing can also be modeled.
Mathematical and Numerical Modeling of Turbulent Flows.
Vedovoto, João M; Serfaty, Ricardo; Da Silveira Neto, Aristeu
2015-01-01
The present work is devoted to the development and implementation of a computational framework to perform numerical simulations of low Mach number turbulent flows over complex geometries. The algorithm under consideration is based on a classical predictor-corrector time integration scheme that employs a projection method for the momentum equations. The domain decomposition strategy is adopted for distributed computing, displaying very satisfactory levels of speed-up and efficiency. The Immersed Boundary Methodology is used to characterize the presence of a complex geometry. Such method demands two separate grids: An Eulerian, where the transport equations are solved with a Finite Volume, second order discretization and a Lagrangian domain, represented by a non-structured shell grid representing the immersed geometry. The in-house code developed was fully verified by the Method of Manufactured Solutions, in both Eulerian and Lagrangian domains. The capabilities of the resulting computational framework are illustrated on four distinct cases: a turbulent jet, the Poiseuille flow, as a matter of validation of the implemented Immersed Boundary methodology, the flow over a sphere covering a wide range of Reynolds numbers, and finally, with the intention of demonstrating the applicability of Large Eddy Simulations - LES - in an industrial problem, the turbulent flow inside an industrial fan. PMID:26131642
Saturn's North Polar Hexagon Numerical Modeling Results
NASA Astrophysics Data System (ADS)
Morales-Juberias, R.; Sayanagi, K. M.; Dowling, T. E.
2008-12-01
In 1980, Voyager images revealed the presence of a circumpolar wave at 78 degrees planetographic latitude in the northern hemisphere of Saturn. It was notable for having a dominant planetary wavenumber-six zonal mode, and for being stationary with respect to Saturn's Kilometric Radiation rotation rate measured by Voyager. The center of this hexagonal feature was coincident with the center of a sharp eastward jet with a peak speed of 100 ms-1 and it had a meridional width of about 4 degrees. This hexagonal feature was confirmed in 1991 through ground-based observations, and it was observed again in 2006 with the Cassini VIMS instrument. The latest observations highlight the longevity of the hexagon and suggest that it extends at least several bars deep into the atmosphere. We use the Explicit Planetary Isentropic Code (EPIC) to perform high-resolution numerical simulations of this unique feature. We show that a wavenumber six instability mode arises naturally from initially barotropic jets when seeded with weak random turbulence. We also discuss the properties of the wave activity on the background vertical stability, zonal wind, planetary rotation rate and adjacent vortices. Computational resources were provided by the New Mexico Computing Applications Center and New Mexico Institute of Mining and Technology and the Comparative Planetology Laboratory at the University of Louisville.
Numerical model of crater lake eruptions
NASA Astrophysics Data System (ADS)
Morrissey, M.; Gisler, G.; Weaver, R.; Gittings, M.
2010-12-01
We present results from a numerical investigation of subaqueous eruptions involving superheated steam released through a lake mimicking the volcanic setting at Mt. Ruapehu. The simulations were conducted using an adaptive mesh, multi-material, hydrodynamics code with thermal conduction SAGE, (Simple Adaptive Grid Eulerian). Parameters investigated include eruption pressure, lake level and mass of superheated vapor. The simulations produced a spectrum of eruption styles from vapor cavities to radial jets that resulted in hazards that ranged from small-scale waves to high amplitude surges that reached and cascaded over the edge of the crater rim. There was an overall tendency for lake surface activity to increase (including wave amplitude) with increasing mass of superheated vapor and eruption pressure. Surface waves were induced by the formation and collapse of a gas cavity. The collapse of the cavity is considered to play a major role in the characteristic features observed during a subaqueous eruption. The additional mass of superheated vapor produced a larger cavity that displaced a larger area of the lake surface resulting in fast moving surges upon the collapse of the cavity. High lake levels (>90 m) appear to suppress the development of explosive jetting activity when eruption pressures are <10 MPa. At very large eruption pressures (>10 MPa), vertical jets and radial ejections of steam and water can occur in water depths >90 m. Less explosive eruption styles can produce hazardous events such as lahars by the outward movement of surface waves over the crater rim.
Evaluation of wave runup predictions from numerical and parametric models
Stockdon, Hilary F.; Thompson, David M.; Plant, Nathaniel G.; Long, Joseph W.
2014-01-01
Wave runup during storms is a primary driver of coastal evolution, including shoreline and dune erosion and barrier island overwash. Runup and its components, setup and swash, can be predicted from a parameterized model that was developed by comparing runup observations to offshore wave height, wave period, and local beach slope. Because observations during extreme storms are often unavailable, a numerical model is used to simulate the storm-driven runup to compare to the parameterized model and then develop an approach to improve the accuracy of the parameterization. Numerically simulated and parameterized runup were compared to observations to evaluate model accuracies. The analysis demonstrated that setup was accurately predicted by both the parameterized model and numerical simulations. Infragravity swash heights were most accurately predicted by the parameterized model. The numerical model suffered from bias and gain errors that depended on whether a one-dimensional or two-dimensional spatial domain was used. Nonetheless, all of the predictions were significantly correlated to the observations, implying that the systematic errors can be corrected. The numerical simulations did not resolve the incident-band swash motions, as expected, and the parameterized model performed best at predicting incident-band swash heights. An assimilated prediction using a weighted average of the parameterized model and the numerical simulations resulted in a reduction in prediction error variance. Finally, the numerical simulations were extended to include storm conditions that have not been previously observed. These results indicated that the parameterized predictions of setup may need modification for extreme conditions; numerical simulations can be used to extend the validity of the parameterized predictions of infragravity swash; and numerical simulations systematically underpredict incident swash, which is relatively unimportant under extreme conditions.
Numerical Poisson-Boltzmann Model for Continuum Membrane Systems.
Botello-Smith, Wesley M; Liu, Xingping; Cai, Qin; Li, Zhilin; Zhao, Hongkai; Luo, Ray
2013-01-01
Membrane protein systems are important computational research topics due to their roles in rational drug design. In this study, we developed a continuum membrane model utilizing a level set formulation under the numerical Poisson-Boltzmann framework within the AMBER molecular mechanics suite for applications such as protein-ligand binding affinity and docking pose predictions. Two numerical solvers were adapted for periodic systems to alleviate possible edge effects. Validation on systems ranging from organic molecules to membrane proteins up to 200 residues, demonstrated good numerical properties. This lays foundations for sophisticated models with variable dielectric treatments and second-order accurate modeling of solvation interactions. PMID:23439886
Numerical model for the Programmable Multirole Furnace (PMZF)
NASA Technical Reports Server (NTRS)
Kassemi, M.; Panzarella, C. H.; Destro-Sidik, K. E.; Krolikowski, C. R.; Licht, B. W.
1993-01-01
The present account of the Programmable Multizone Furnace numerical model uses various examples to illustrate the ways in which the model serves as an optimization, test, prediction, and visualization tool; a numerical PID-control algorithm obtains the desired sample temperature distributions and allows the model to solve an inverse heat transfer problem where the desired sample temperature profile is the input and the required heater power distribution is the output of numerical simulations. Parametric studies show how the total power consumption of the furnace is affected by such design variables as the conductivity.
Numerical Modeling of Ophthalmic Response to Space
NASA Technical Reports Server (NTRS)
Nelson, E. S.; Myers, J. G.; Mulugeta, L.; Vera, J.; Raykin, J.; Feola, A.; Gleason, R.; Samuels, B.; Ethier, C. R.
2015-01-01
To investigate ophthalmic changes in spaceflight, we would like to predict the impact of blood dysregulation and elevated intracranial pressure (ICP) on Intraocular Pressure (IOP). Unlike other physiological systems, there are very few lumped parameter models of the eye. The eye model described here is novel in its inclusion of the human choroid and retrobulbar subarachnoid space (rSAS), which are key elements in investigating the impact of increased ICP and ocular blood volume. Some ingenuity was required in modeling the blood and rSAS compartments due to the lack of quantitative data on essential hydrodynamic quantities, such as net choroidal volume and blood flowrate, inlet and exit pressures, and material properties, such as compliances between compartments.
Mathematical and Numerical Analyses of Peridynamics for Multiscale Materials Modeling
Gunzburger, Max
2015-02-17
We have treated the modeling, analysis, numerical analysis, and algorithmic development for nonlocal models of diffusion and mechanics. Variational formulations were developed and finite element methods were developed based on those formulations for both steady state and time dependent problems. Obstacle problems and optimization problems for the nonlocal models were also treated and connections made with fractional derivative models.
NASA Astrophysics Data System (ADS)
Bailey, Brian N.
2016-07-01
When Lagrangian stochastic models for turbulent dispersion are applied to complex atmospheric flows, some type of ad hoc intervention is almost always necessary to eliminate unphysical behaviour in the numerical solution. Here we discuss numerical strategies for solving the non-linear Langevin-based particle velocity evolution equation that eliminate such unphysical behaviour in both Reynolds-averaged and large-eddy simulation applications. Extremely large or `rogue' particle velocities are caused when the numerical integration scheme becomes unstable. Such instabilities can be eliminated by using a sufficiently small integration timestep, or in cases where the required timestep is unrealistically small, an unconditionally stable implicit integration scheme can be used. When the generalized anisotropic turbulence model is used, it is critical that the input velocity covariance tensor be realizable, otherwise unphysical behaviour can become problematic regardless of the integration scheme or size of the timestep. A method is presented to ensure realizability, and thus eliminate such behaviour. It was also found that the numerical accuracy of the integration scheme determined the degree to which the second law of thermodynamics or `well-mixed condition' was satisfied. Perhaps more importantly, it also determined the degree to which modelled Eulerian particle velocity statistics matched the specified Eulerian distributions (which is the ultimate goal of the numerical solution). It is recommended that future models be verified by not only checking the well-mixed condition, but perhaps more importantly by checking that computed Eulerian statistics match the Eulerian statistics specified as inputs.
Numerical Modeling of Left-Handed Metamaterials
Burke, G J; Champagne, N J; Sharpe, R M
2001-11-06
The EIGER method of moments program with periodic Green's function was used to model a periodic array of strips and split-ring resonators. Left-handed propagation due to negative index of refraction is demonstrated in a frequency band. The effective material parameters versus frequency are extracted from the EIGER solution.
Three-dimensional numerical model of holographic grating formation in photopolymer materials
NASA Astrophysics Data System (ADS)
Li, Haoyu; Qi, Yue; Malallah, Ra'ed; Sheridan, John T.
2015-05-01
When the large thickness is used as the holographic storage materials, a non-ignorable problem is the light intensity attenuation in depth due to high absorptive of the dye. For this reason more completely modeling the evolutions inside the material is necessary to consider into the developed standard kinetic model. In this paper the photo-polymerization processes during the large thickness holographic grating formation are analyzed. A 3-dimensional algorithm is present by deriving the system partial differential rate equations governing each associated chemical species, and using the finite difference approximation, these equations can be solved numerically. This extended model describes the time varying behaviors of the non-uniform photo-physical and the photochemical evolutions in photopolymer materials. In this model both dye molecules consumption and light energy absorption are calculated time varyingly, and then the polymer and monomer concentrations distributions are obtained. Applying the Lorenz-Lorenz relationship, the non-uniform grating formatted in material depth, and its refractive index, which is distorted from ideal sinusoidal spatial distribution, can be more accurately predicted.
Numerical modelling of instantaneous plate tectonics
NASA Technical Reports Server (NTRS)
Minster, J. B.; Haines, E.; Jordan, T. H.; Molnar, P.
1974-01-01
Assuming lithospheric plates to be rigid, 68 spreading rates, 62 fracture zones trends, and 106 earthquake slip vectors are systematically inverted to obtain a self-consistent model of instantaneous relative motions for eleven major plates. The inverse problem is linearized and solved iteratively by a maximum-likelihood procedure. Because the uncertainties in the data are small, Gaussian statistics are shown to be adequate. The use of a linear theory permits (1) the calculation of the uncertainties in the various angular velocity vectors caused by uncertainties in the data, and (2) quantitative examination of the distribution of information within the data set. The existence of a self-consistent model satisfying all the data is strong justification of the rigid plate assumption. Slow movement between North and South America is shown to be resolvable.
Analytical and numerical modeling for flexible pipes
NASA Astrophysics Data System (ADS)
Wang, Wei; Chen, Geng
2011-12-01
The unbonded flexible pipe of eight layers, in which all the layers except the carcass layer are assumed to have isotropic properties, has been analyzed. Specifically, the carcass layer shows the orthotropic characteristics. The effective elastic moduli of the carcass layer have been developed in terms of the influence of deformation to stiffness. With consideration of the effective elastic moduli, the structure can be properly analyzed. Also the relative movements of tendons and relative displacements of wires in helical armour layer have been investigated. A three-dimensional nonlinear finite element model has been presented to predict the response of flexible pipes under axial force and torque. Further, the friction and contact of interlayer have been considered. Comparison between the finite element model and experimental results obtained in literature has been given and discussed, which might provide practical and technical support for the application of unbonded flexible pipes.
Advanced Numerical Modeling of Turbulent Atmospheric Flows
NASA Astrophysics Data System (ADS)
Kühnlein, Christian; Dörnbrack, Andreas; Gerz, Thomas
The present chapter introduces the method of computational simulation to predict and study turbulent atmospheric flows. This includes a description of the fundamental approach to computational simulation and the practical implementation using the technique of large-eddy simulation. In addition, selected contributions from IPA scientists to computational model development and various examples for applications are given. These examples include homogeneous turbulence, convective boundary layers, heated forest canopy, buoyant thermals, and large-scale flows with baroclinic wave instability.
Numerical Modelling of the Mining Induced Horizontal Displacement
NASA Astrophysics Data System (ADS)
Tajduś, Krzysztof
2015-12-01
The paper presents results of numerical calculations and modeling of mining-induced surface deformation based on Finite Element Method (FEM). Applying the numerical method discussed to calculations allows us to assume a larger number of factors, such as rock mass structure, fracture network, rock properties, etc., which essentially affect the results obtained. On the basis of an elastic transversely isotropic model, an analysis of horizontal displacement distribution and surface subsidence was carried out for two sample regions of mines. The results of numerical calculations were later compared with the measured values. Such an analysis proved that the applied numerical model properly described distribution and values of subsidence and slope of subsidence trough, though there were serious differences in the values of calculated horizontal displacement, especially in areas of far influence range. In order to improve the matching, the influence of boundary conditions of the model on the value of calculated horizontal displacement was analyzed. The results are presented in graphs.
Numerical Modelling of Embankment on Soft Clay
NASA Astrophysics Data System (ADS)
Nujid, M. M.; Taha, M. R.
2016-07-01
This paper aims to predict deformation of embankment on soft clay of Muar. The prediction performance focusing on displacement at critical fill height of 5.5 m. The study was based on reported result in 1992. With the aid of computer intelligence, the advanced constitutive soil models could be adopted to analyze the soft clay behavior. The COMSOL Multiphysics (v4.4) has been used to simulate the problem with coupled physics available in the software. The vertical displacements are in good agreement close to published result.
Numerical modelling of the TFTR ICRH antennas
Kress, M.; Ho, Y.L.; Grossmann, W.; Drobot, A. ); Batchelor, D.B.; Ryan, P.M.; Carter, M. )
1991-01-01
A general purpose 3D electromagnetic field solver code, ARGUS, is being used to analyze the TFTR Antennas. To date, the vacuum radiation patterns produced by the bay M and L antennas have been obtained and reported. Recent work has concentrated on antenna performance comparison and understanding the role of geometry on performance (e.g., the impact of end-effects on current 2D models). Additional diagnostics such as evaluation of phase velocity and strap inductance are being implemented to enhance our understanding and to better compare with measurements. 6 refs., 2 figs.
Terrane accretion: Insights from numerical modelling
NASA Astrophysics Data System (ADS)
Vogt, Katharina; Gerya, Taras
2016-04-01
The oceanic crust is not homogenous, but contains significantly thicker crust than norm, i.e. extinct arcs, spreading ridges, detached continental fragments, volcanic piles or oceanic swells. These (crustal) fragments may collide with continental crust and form accretionary complexes, contributing to its growth. We analyse this process using a thermo-mechanical computer model (i2vis) of an ocean-continent subduction zone. In this model the oceanic plate can bend spontaneously under the control of visco-plastic rheologies. It moreover incorporates effects such as mineralogical phase changes, fluid release and consumption, partial melting and melt extraction. Based on our 2-D experiments we suggest that the lithospheric buoyancy of the downgoing slab and the rheological strength of crustal material may result in a variety of accretionary processes. In addition to terrane subduction, we are able to identify three distinct modes of terrane accretion: frontal accretion, basal accretion and underplating plateaus. We show that crustal fragments may dock onto continental crust and cease subduction, be scrapped off the downgoing plate, or subduct to greater depth prior to slab break off and subsequent exhumation. Direct consequences of these processes include slab break off, subduction zone transference, structural reworking, formation of high-pressure terranes, partial melting and crustal growth.
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
Numerical analysis and modeling of atmospheric phenomena
NASA Technical Reports Server (NTRS)
Stone, Peter H.
1994-01-01
For the past 22 years Grant NGR 22-009-727 has been supporting research in the Center for Meteorology and Physical Oceanography (and its predecessors) in a wide variety of diagnostic and modeling studies of atmospheric and ocean phenomena. Professor Jule Charney was the initial Principal Investigator. Professor Peter Stone joined him as co-Principal Investigator in 1975 and became the sole Principal Investigator in 1981. During its lifetime the Grant has supported in whole or in part 11 Master's theses, 14 Ph.D. theses, and 45 papers published in refereed scientific journals. All of these theses and papers (with bibliographic references) are listed below. All but one of the theses were used to fulfill the requirements for MIT (Massachusetts Institute of Technology) degrees and are available from the MIT libraries. The one exception is F. Chen's Ph.D. thesis which was for a Harvard degree and is available from the Harvard libraries. In addition to the work described in the citations listed below, the Grant has supported Research Assistant Amy Solomon during the past two years to carry out a study of how baroclinic adjustment is affected by vertical resolution, vertical temperature structure, and dissipation. Ms. Solomon plans to use this project for her Ph.D. thesis. Support for this project will continue under NASA Grant NAG 5-2490, 'The Factors Controlling Poleward Heat Transport in Climate Models.'
Experimental, Numerical and Observational Models in Geodynamics
NASA Astrophysics Data System (ADS)
Lithgow-Bertelloni, Carolina
2015-04-01
Geodynamics, the study of the forces that drives all Earth's processes is a rich field that deeply connects all aspects of geological and geophysical studies, from surface observations of the sedimentary record to knowledge of deep Earth structure from mineral physics and seismology. In the context of the solid Earth geodynamics primarily focuses on lithosphere and mantle dynamics, while core dynamics is the purview of geomagnetism. I will focus this talk on the former, its historical context and future developments. We have known the equations of motion and mechanics for ~200 years, but only relatively recently can they be solved with enough accuracy and resolution to do geology. We have made great strides since Arthur Holmes conceptual models of mantle flow, thanks to computational and experimental advances. We can know model plate boundaries globally with resolutions in the order of a few kms and image temperature and velocity simultaneously in the laboratory in 3D and non-intrusively. We have also learned a great deal about the physics of the Earth, from composition to rheology. New theories on plate boundary rheology are paving the way for self-consistent generation of plates from mantle flow. New computational methods allow for adaptive meshing, fabric development and history, so we can study deformation and compare directly to geological observations in mountain ranges and continental rifts. We can use ever more sophisticated images of mantle structure from seismic and other geophysical data to probe the relationship between melting, flow and dynamical processes. We can reconstruct landscapes and relief, plate motions and sedimentation and ask how much the mantle has contributed to drainage reversal, sedimentation and climate change. The future of the field is ever brighter.
Seismoelectric numerical modeling on a grid
Haines, S.S.; Pride, S.R.
2006-01-01
Our finite-difference algorithm provides a new method for simulating how seismic waves in arbitrarily heterogeneous porous media generate electric fields through an electrokinetic mechanism called seismoelectric coupling. As the first step in our simulations, we calculate relative pore-fluid/grain-matrix displacement by using existing poroelastic theory. We then calculate the electric current resulting from the grain/fluid displacement by using seismoelectric coupling theory. This electrofiltration current acts as a source term in Poisson's equation, which then allows us to calculate the electric potential distribution. We can safely neglect induction effects in our simulations because the model area is within the electrostatic near field for the depth of investigation (tens to hundreds of meters) and the frequency ranges (10 Hz to 1 kHz) of interest for shallow seismoelectric surveys.We can independently calculate the electric-potential distribution for each time step in the poroelastic simulation without loss of accuracy because electro-osmotic feedback (fluid flow that is perturbed by generated electric fields) is at least 105 times smaller than flow that is driven by fluid-pressure gradients and matrix acceleration, and is therefore negligible. Our simulations demonstrate that, distinct from seismic reflections, the seismoelectric interface response from a thin layer (at least as thin as one-twentieth of the seismic wavelength) is considerably stronger than the response from a single interface. We find that the interface response amplitude decreases as the lateral extent of a layer decreases below the width of the first Fresnel zone. We conclude, on the basis of our modeling results and of field results published elsewhere, that downhole and/or crosswell survey geometries and time-lapse applications are particularly well suited to the seismoelectric method. ?? 2006 Society of Exploration Geophysicists.
Theoretical and numerical modelling of shocks in dusty plasmas
Eliasson, B.; Shukla, P.K.
2005-10-31
The formation of dust acoustic (DA) and dust ion-acoustic (DIA) shocks are are studied theoretically and numerically by means of simple-wave solutions and a comparison between fluid and kinetic model for DIA waves. A fluid model admits sharp discontinuities at the shock front while the kinetic model involves Landau-damping of the the shock front.
USER GUIDE FOR THE ENHANCED HYDRODYNAMICAL-NUMERICAL MODEL
This guide provides the documentation required for used of the Enhanced Hydrodynamical-Numerical Model on operational problems. The enhanced model is a multilayer Hansen type model extended to handle near-shore processes by including: Non-linear term extension to facilitate small...
Numerical model for learning concepts of streamflow simulation
DeLong, L.L.
1993-01-01
Numerical models are useful for demonstrating principles of open-channel flow. Such models can allow experimentation with cause-and-effect relations, testing concepts of physics and numerical techniques. Four PT is a numerical model written primarily as a teaching supplement for a course in one-dimensional stream-flow modeling. Four PT options particularly useful in training include selection of governing equations, boundary-value perturbation, and user-programmable constraint equations. The model can simulate non-trivial concepts such as flow in complex interconnected channel networks, meandering channels with variable effective flow lengths, hydraulic structures defined by unique three-parameter relations, and density-driven flow.The model is coded in FORTRAN 77, and data encapsulation is used extensively to simplify maintenance and modification and to enhance the use of Four PT modules by other programs and programmers.
Numerical modeling of elastodynamic radiation and scattering
Savic, M.; Ziolkowski, A.M.
1994-12-31
This paper presents a study on two problems: the two-dimensional distributed surface load problem, and the scattering of elastodynamic waves from fractures. The analysis is done with the aid of the finite-difference technique. If the dimensions of a surface mechanical source (vibrator or piezoelectric transducer) are not small compared to the wavelength, one should not use the point source or plane wave representation when modeling radiation from such sources. Here the authors demonstrate the solution of the uniformly distributed surface load problem using the finite-difference (FD) technique. The scattering of transient elasto-dynamic waves from a fracture whose extent is large compared with the wavelength and whose width is small compared with the wavelength and whose width is small compared with the wavelength is one of the classical problems in seismology and non-destructive testing (NDT). Many researchers have provided analytical solutions based on different approximations for the unknown field (displacement or particle velocity) scattered from an idealized half-plane or the a strip of finite extent. Again, the authors demonstrate the full wavefield solution using the finite-difference technique. The technique presented here is aimed for the interpretation of seismic data from hydraulic fracturing experiments.
Squeal noise in simple numerical brake models
NASA Astrophysics Data System (ADS)
Oberst, S.; Lai, J. C. S.
2015-09-01
Since the early 1920s, automotive disc brake squeal has caused warranty issues and customer dissatisfaction. Despite a good deal of progress achieved, predicting brake squeal propensity is as difficult as ever as not all mechanisms and interactions are known owing to their highly fugitive complex nature. In recent years, research has been focused on the prediction of unstable vibration modes by the complex eigenvalue analysis (CEA) for the mode-coupling type of instability. There has been very limited consideration given to the calculation of the acoustic radiation properties due to friction contact between a pad and a rotor. Recent analyses using a forced response analysis with harmonic contact pressure excitation indicates negative dissipated energy at some pad eigenfrequencies predicted to be stable by the CEA. A transient nonlinear time domain analysis with no external excitation indicates that squeal could develop at these eigenfrequencies. Here, the acoustic radiation characteristics of those pad modes are determined by analysing the acoustic power levels and radiation efficiencies of simplified brake models in the form of a pad rubbing on a plate or on a disc using the acoustic boundary element method based on velocities extracted from the forced response analysis. Results show that unstable pad modes trigger unstable disc vibrations resulting in instantaneous mode squeal similar to those observed experimentally. Changes in the radiation efficiency with pressure variations are smaller than those with friction coefficient variations and are caused by the phase difference of the velocities out-of-plane vibration between the pad and the disc.
Quantitative analysis of numerical solvers for oscillatory biomolecular system models
Quo, Chang F; Wang, May D
2008-01-01
Background This article provides guidelines for selecting optimal numerical solvers for biomolecular system models. Because various parameters of the same system could have drastically different ranges from 10-15 to 1010, the ODEs can be stiff and ill-conditioned, resulting in non-unique, non-existing, or non-reproducible modeling solutions. Previous studies have not examined in depth how to best select numerical solvers for biomolecular system models, which makes it difficult to experimentally validate the modeling results. To address this problem, we have chosen one of the well-known stiff initial value problems with limit cycle behavior as a test-bed system model. Solving this model, we have illustrated that different answers may result from different numerical solvers. We use MATLAB numerical solvers because they are optimized and widely used by the modeling community. We have also conducted a systematic study of numerical solver performances by using qualitative and quantitative measures such as convergence, accuracy, and computational cost (i.e. in terms of function evaluation, partial derivative, LU decomposition, and "take-off" points). The results show that the modeling solutions can be drastically different using different numerical solvers. Thus, it is important to intelligently select numerical solvers when solving biomolecular system models. Results The classic Belousov-Zhabotinskii (BZ) reaction is described by the Oregonator model and is used as a case study. We report two guidelines in selecting optimal numerical solver(s) for stiff, complex oscillatory systems: (i) for problems with unknown parameters, ode45 is the optimal choice regardless of the relative error tolerance; (ii) for known stiff problems, both ode113 and ode15s are good choices under strict relative tolerance conditions. Conclusions For any given biomolecular model, by building a library of numerical solvers with quantitative performance assessment metric, we show that it is possible
Nonspinning numerical relativity waveform surrogates: Building the model
NASA Astrophysics Data System (ADS)
Galley, Chad
2015-04-01
Simulating binary black hole coalescences involves solving Einstein's equations with large-scale computing resources that can take months to complete for a single numerical solution. This engenders a computationally intractable problem for multiple-query applications related to parameter space exploration, data analysis for gravitational wave detectors like LIGO, and semi-analytical waveform fits. I discuss how reduced order modeling techniques are used to build accurate surrogates that can be evaluated quickly in place of numerically solving Einstein's equations for generating gravitational waveforms of nonspinning binary black hole coalescences. To within error, the surrogate can model all modes available from a numerical simulation including, for example, troublesome modes such as the (3,2) mode and memory modes. A companion talk will cover quantifying the best surrogate model's errors. The results of this work represent a significant advance by making it possible to use numerical relativity waveforms for multiple-query applications.
Numerical modelling of the solidification of ductile iron
NASA Astrophysics Data System (ADS)
Liu, J.; Elliott, R.
1998-01-01
Numerical calculations are presented describing the solidification of a ductile iron based on the Stefanescu macroscopic heat transfer-microscopic solidification kinetic model but using a different kinetic model than that used by Stefanescu. The results show that the kinetic model used influences the recalescence behaviour predicted by the modelling. Cooling curves calculated with the present model show reasonable agreement with experimentally measured cooling curves for four different cooling rates.
Numerical simulations of a reduced model for blood coagulation
NASA Astrophysics Data System (ADS)
Pavlova, Jevgenija; Fasano, Antonio; Sequeira, Adélia
2016-04-01
In this work, the three-dimensional numerical resolution of a complex mathematical model for the blood coagulation process is presented. The model was illustrated in Fasano et al. (Clin Hemorheol Microcirc 51:1-14, 2012), Pavlova et al. (Theor Biol 380:367-379, 2015). It incorporates the action of the biochemical and cellular components of blood as well as the effects of the flow. The model is characterized by a reduction in the biochemical network and considers the impact of the blood slip at the vessel wall. Numerical results showing the capacity of the model to predict different perturbations in the hemostatic system are discussed.
Quantitative comparisons of numerical models of brittle wedge dynamics
NASA Astrophysics Data System (ADS)
Buiter, Susanne
2010-05-01
Numerical and laboratory models are often used to investigate the evolution of deformation processes at various scales in crust and lithosphere. In both approaches, the freedom in choice of simulation method, materials and their properties, and deformation laws could affect model outcomes. To assess the role of modelling method and to quantify the variability among models, we have performed a comparison of laboratory and numerical experiments. Here, we present results of 11 numerical codes, which use finite element, finite difference and distinct element techniques. We present three experiments that describe shortening of a sand-like, brittle wedge. The material properties of the numerical ‘sand', the model set-up and the boundary conditions are strictly prescribed and follow the analogue setup as closely as possible. Our first experiment translates a non-accreting wedge with a stable surface slope of 20 degrees. In agreement with critical wedge theory, all models maintain the same surface slope and do not deform. This experiment serves as a reference that allows for testing against analytical solutions for taper angle, root-mean-square velocity and gravitational rate of work. The next two experiments investigate an unstable wedge in a sandbox-like setup, which deforms by inward translation of a mobile wall. The models accommodate shortening by formation of forward and backward shear zones. We compare surface slope, rate of dissipation of energy, root-mean-square velocity, and the location, dip angle and spacing of shear zones. We show that we successfully simulate sandbox-style brittle behaviour using different numerical modelling techniques and that we obtain the same styles of deformation behaviour in numerical and laboratory experiments at similar levels of variability. The GeoMod2008 Numerical Team: Markus Albertz, Michelle Cooke, Tony Crook, David Egholm, Susan Ellis, Taras Gerya, Luke Hodkinson, Boris Kaus, Walter Landry, Bertrand Maillot, Yury Mishin
Numerical modelling in biosciences using delay differential equations
NASA Astrophysics Data System (ADS)
Bocharov, Gennadii A.; Rihan, Fathalla A.
2000-12-01
Our principal purposes here are (i) to consider, from the perspective of applied mathematics, models of phenomena in the biosciences that are based on delay differential equations and for which numerical approaches are a major tool in understanding their dynamics, (ii) to review the application of numerical techniques to investigate these models. We show that there are prima facie reasons for using such models: (i) they have a richer mathematical framework (compared with ordinary differential equations) for the analysis of biosystem dynamics, (ii) they display better consistency with the nature of certain biological processes and predictive results. We analyze both the qualitative and quantitative role that delays play in basic time-lag models proposed in population dynamics, epidemiology, physiology, immunology, neural networks and cell kinetics. We then indicate suitable computational techniques for the numerical treatment of mathematical problems emerging in the biosciences, comparing them with those implemented by the bio-modellers.
Nonlinear oscillator metamaterial model: numerical and experimental verification.
Poutrina, E; Huang, D; Urzhumov, Y; Smith, D R
2011-04-25
We verify numerically and experimentally the accuracy of an analytical model used to derive the effective nonlinear susceptibilities of a varactor-loaded split ring resonator (VLSRR) magnetic medium. For the numerical validation, a nonlinear oscillator model for the effective magnetization of the metamaterial is applied in conjunction with Maxwell equations and the two sets of equations solved numerically in the time-domain. The computed second harmonic generation (SHG) from a slab of a nonlinear material is then compared with the analytical model. The computed SHG is in excellent agreement with that predicted by the analytical model, both in terms of magnitude and spectral characteristics. Moreover, experimental measurements of the power transmitted through a fabricated VLSRR metamaterial at several power levels are also in agreement with the model, illustrating that the effective medium techniques associated with metamaterials can accurately be transitioned to nonlinear systems. PMID:21643082
Sheet Hydroforming Process Numerical Model Improvement Through Experimental Results Analysis
NASA Astrophysics Data System (ADS)
Gabriele, Papadia; Antonio, Del Prete; Alfredo, Anglani
2010-06-01
The increasing application of numerical simulation in metal forming field has helped engineers to solve problems one after another to manufacture a qualified formed product reducing the required time [1]. Accurate simulation results are fundamental for the tooling and the product designs. The wide application of numerical simulation is encouraging the development of highly accurate simulation procedures to meet industrial requirements. Many factors can influence the final simulation results and many studies have been carried out about materials [2], yield criteria [3] and plastic deformation [4,5], process parameters [6] and their optimization. In order to develop a reliable hydromechanical deep drawing (HDD) numerical model the authors have been worked out specific activities based on the evaluation of the effective stiffness of the blankholder structure [7]. In this paper after an appropriate tuning phase of the blankholder force distribution, the experimental activity has been taken into account to improve the accuracy of the numerical model. In the first phase, the effective capability of the blankholder structure to transfer the applied load given by hydraulic actuators to the blank has been explored. This phase ended with the definition of an appropriate subdivision of the blankholder active surface in order to take into account the effective pressure map obtained for the given loads configuration. In the second phase the numerical results obtained with the developed subdivision have been compared with the experimental data of the studied model. The numerical model has been then improved, finding the best solution for the blankholder force distribution.
Numerical modelling of river morphodynamics: Latest developments and remaining challenges
NASA Astrophysics Data System (ADS)
Siviglia, Annunziato; Crosato, Alessandra
2016-07-01
Numerical morphodynamic models provide scientific frameworks for advancing our understanding of river systems. The research on involved topics is an important and socially relevant undertaking regarding our environment. Nowadays numerical models are used for different purposes, from answering questions about basic morphodynamic research to managing complex river engineering problems. Due to increasing computer power and the development of advanced numerical techniques, morphodynamic models are now more and more used to predict the bed patterns evolution to a broad spectrum of spatial and temporal scales. The development and the success of application of such models are based upon a wide range of disciplines from applied mathematics for the numerical solution of the equations to geomorphology for the physical interpretation of the results. In this light we organized this special issue (SI) soliciting multidisciplinary contributions which encompass any aspect needed for the development and applications of such models. Most of the papers in the SI stem from contributions to session HS9.5/GM7.11 on numerical modelling and experiments in river morphodynamics at the European Geosciences Union (EGU) General Assembly held in Vienna, April 27th to May 2nd 2014.
Nonspinning numerical relativity waveform surrogates: assessing the model
NASA Astrophysics Data System (ADS)
Field, Scott; Blackman, Jonathan; Galley, Chad; Scheel, Mark; Szilagyi, Bela; Tiglio, Manuel
2015-04-01
Recently, multi-modal gravitational waveform surrogate models have been built directly from data numerically generated by the Spectral Einstein Code (SpEC). I will describe ways in which the surrogate model error can be quantified. This task, in turn, requires (i) characterizing differences between waveforms computed by SpEC with those predicted by the surrogate model and (ii) estimating errors associated with the SpEC waveforms from which the surrogate is built. Both pieces can have numerous sources of numerical and systematic errors. We make an attempt to study the most dominant error sources and, ultimately, the surrogate model's fidelity. These investigations yield information about the surrogate model's uncertainty as a function of time (or frequency) and parameter, and could be useful in parameter estimation studies which seek to incorporate model error. Finally, I will conclude by comparing the numerical relativity surrogate model to other inspiral-merger-ringdown models. A companion talk will cover the building of multi-modal surrogate models.
NASA Astrophysics Data System (ADS)
Chatterjee, Krishnendu; Tuli, Suneet
2013-11-01
Lock-in thermography is increasingly becoming popular as a non-destructive testing technique for defect detection in composite materials for its low heating excitation. The experimental data is processed with Fourier transformation to produce phase and amplitude images. Phase images, though immune to surface emissivity variation, suffer from blind frequency effect, where a defect becomes invisible at a certain excitation frequency. There exists no analytical model to predict this 3-dimensional heat flow phenomenon. This paper presents a study of blind frequency using electro-thermal model based numerical simulation on a piece of thermally anisotropic carbon fibre composite. The performance of the simulator is optimized for spatial mesh size. Further the effect of paint layer, which is often applied to the sample surface for better thermal imaging, has been incorporated in the simulation. Finally, both experimental and simulation results are presented side-by-side for easy comparison.
Comparison of two numerical techniques for aerodynamic model identification
NASA Technical Reports Server (NTRS)
Verhaegen, M. H.
1987-01-01
An algorithm, called the Minimal Residual QR algorithm, is presented to solve subset regression problems. It is shown that this scheme can be used as a numerically reliable implementation of the stepwise regression technique, which is widely used to identify an aerodynamic model from flight test data. This capability as well as the numerical superiority of this scheme over the stepwise regression technique is demonstrated in an experimental simulation study.
Numerical Models of Broad Bandwidth Nanosecond Optical Parametric Oscillators
Bowers, M.S.; Gehr, R.J.; Smith, A.V.
1998-10-14
We describe results from three new methods of numerically modeling broad-bandwidth, nanosecond OPO's in the plane-wave approximate ion. They account for differences in group velocities among the three mixing waves, and also include a qutt~ttun noise model.
A survey of numerical models for wind prediction
NASA Technical Reports Server (NTRS)
Schonfeld, D.
1980-01-01
A literature review is presented of the work done in the numerical modeling of wind flows. Pertinent computational techniques are described, as well as the necessary assumptions used to simplify the governing equations. A steady state model is outlined, based on the data obtained at the Deep Space Communications complex at Goldstone, California.
Hydroforming Of Patchwork Blanks — Numerical Modeling And Experimental Validation
NASA Astrophysics Data System (ADS)
Lamprecht, Klaus; Merklein, Marion; Geiger, Manfred
2005-08-01
In comparison to the commonly applied technology of tailored blanks the concept of patchwork blanks offers a number of additional advantages. Potential application areas for patchwork blanks in automotive industry are e.g. local reinforcements of automotive closures, structural reinforcements of rails and pillars as well as shock towers. But even if there is a significant application potential for patchwork blanks in automobile production, industrial realization of this innovative technique is decelerated due to a lack of knowledge regarding the forming behavior and the numerical modeling of patchwork blanks. Especially for the numerical simulation of hydroforming processes, where one part of the forming tool is replaced by a fluid under pressure, advanced modeling techniques are required to ensure an accurate prediction of the blanks' forming behavior. The objective of this contribution is to provide an appropriate model for the numerical simulation of patchwork blanks' forming processes. Therefore, different finite element modeling techniques for patchwork blanks are presented. In addition to basic shell element models a combined finite element model consisting of shell and solid elements is defined. Special emphasis is placed on the modeling of the weld seam. For this purpose the local mechanical properties of the weld metal, which have been determined by means of Martens-hardness measurements and uniaxial tensile tests, are integrated in the finite element models. The results obtained from the numerical simulations are compared to experimental data from a hydraulic bulge test. In this context the focus is laid on laser- and spot-welded patchwork blanks.
Experimental, numerical and analytical models of mantle starting plumes
NASA Astrophysics Data System (ADS)
Coulliette, D. L.; Loper, D. E.
1995-12-01
The results of a combined experimental, numerical and analytical investigation of starting thermal plumes are described, to obtain a better perspective on plumes within the Earth's mantle. Thermal plumes were generated experimentally in a tank of corn syrup by means of an electrical heater. Viscosity ratios of 400, 30 000, and 10 8 were generated by varying the temperature of the tank. Plumes for the smaller ratios had the traditional 'balloon-on-astring' shape, but that at the highest ratio had a novel morphology. The plume heads in the first two cases were observed to rise at roughly a constant speed, in contrast to most previous studies which found the plume heads to accelerate. Loss of buoyancy from the plume head owing to heat loss is believed to be responsible for this difference. Starting plumes were simulated numerically using an axisymmetric, finite-element code to solve the Boussinesq equations at finite Prandtl numbers. The constant rise speed observed experimentally was confirmed by the numerical simulation for the viscosity ratios of 400 and 30 000, but numerical instability prevented simulation of the case with a viscosity ratio of 10 8. There was very good agreement between the experimental and numerical rise speeds. An analytical model was developed which reduces to previous models in limiting cases. This parameterization gives better agreement with the experimental and numerical results than does any previous model.
Development of a numerical simulation model of the cardiovascular system.
Geertsema, A A; Rakhorst, G; Mihaylov, D; Blanksma, P K; Verkerke, G J
1997-12-01
A numerical simulation model of the cardiovascular system has been developed. It consists of a model of the left atrium, the left ventricle, the coronary vascular system, the aorta, the arterial system, and the venous system. The input of the complete model is the elastance (pressure/volume ratio) developed by the left ventricle. The shape of this elastance is constant in different circumstances. Left ventricular (LV) myocardial oxygen consumption and the amount of oxygen offered to the left ventricle can be calculated with the model. The model has been validated using data from a patient suffering from coronary artery disease. The measured clinical hemodynamical waveforms could be fitted to those generated by the model. With the numerical simulation model, it is possible to predict the functioning of the left ventricle under different circumstances. This makes it possible to study in vitro various pathological clinical situations. PMID:9423983
Rapid installation of numerical models in multiple parent codes
Brannon, R.M.; Wong, M.K.
1996-10-01
A set of``model interface guidelines``, called MIG, is offered as a means to more rapidly install numerical models (such as stress-strain laws) into any parent code (hydrocode, finite element code, etc.) without having to modify the model subroutines. The model developer (who creates the model package in compliance with the guidelines) specifies the model`s input and storage requirements in a standardized way. For portability, database management (such as saving user inputs and field variables) is handled by the parent code. To date, NUG has proved viable in beta installations of several diverse models in vectorized and parallel codes written in different computer languages. A NUG-compliant model can be installed in different codes without modifying the model`s subroutines. By maintaining one model for many codes, MIG facilitates code-to-code comparisons and reduces duplication of effort potentially reducing the cost of installing and sharing models.
Ensemble-type numerical uncertainty information from single model integrations
Rauser, Florian Marotzke, Jochem; Korn, Peter
2015-07-01
We suggest an algorithm that quantifies the discretization error of time-dependent physical quantities of interest (goals) for numerical models of geophysical fluid dynamics. The goal discretization error is estimated using a sum of weighted local discretization errors. The key feature of our algorithm is that these local discretization errors are interpreted as realizations of a random process. The random process is determined by the model and the flow state. From a class of local error random processes we select a suitable specific random process by integrating the model over a short time interval at different resolutions. The weights of the influences of the local discretization errors on the goal are modeled as goal sensitivities, which are calculated via automatic differentiation. The integration of the weighted realizations of local error random processes yields a posterior ensemble of goal approximations from a single run of the numerical model. From the posterior ensemble we derive the uncertainty information of the goal discretization error. This algorithm bypasses the requirement of detailed knowledge about the models discretization to generate numerical error estimates. The algorithm is evaluated for the spherical shallow-water equations. For two standard test cases we successfully estimate the error of regional potential energy, track its evolution, and compare it to standard ensemble techniques. The posterior ensemble shares linear-error-growth properties with ensembles of multiple model integrations when comparably perturbed. The posterior ensemble numerical error estimates are of comparable size as those of a stochastic physics ensemble.
A general numerical model for wave rotor analysis
NASA Technical Reports Server (NTRS)
Paxson, Daniel W.
1992-01-01
Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.
Numerical models of laser fusion of intestinal tissues.
Pearce, John A
2009-01-01
Numerical models of continuous wave Tm:YAG thermal fusion in rat intestinal tissues were compared to experiment. Optical and thermal FDM models that included tissue damage based on Arrhenius kinetics were used to predict birefringence loss in collagen as the standard of comparison. The models also predicted collagen shrinkage, jellification and water loss. The inclusion of variable optical and thermal properties is essential to achieve favorable agreement between predicted and measured damage boundaries. PMID:19964349
Numerical models for the evaluation of geothermal systems
Bodvarsson, G.S.; Pruess, K.; Lippmann, M.J.
1986-08-01
We have carried out detailed simulations of various fields in the USA (Bada, New Mexico; Heber, California); Mexico (Cerro Prieto); Iceland (Krafla); and Kenya (Olkaria). These simulation studies have illustrated the usefulness of numerical models for the overall evaluation of geothermal systems. The methodology for modeling the behavior of geothermal systems, different approaches to geothermal reservoir modeling and how they can be applied in comprehensive evaluation work are discussed.
Numerical modeling of runback water on ice protected aircraft surfaces
NASA Technical Reports Server (NTRS)
Al-Khalil, Kamel M.; Keith, Theo G., Jr.; Dewitt, Kenneth J.
1992-01-01
A numerical simulation for 'running wet' aircraft anti-icing systems is developed. The model includes breakup of the water film, which exists in regions of direct impingement, into individual rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted in the numerical solution procedure. The solid wall is modeled as a multilayer structure and the anti-icing system used is of the thermal type utilizing hot air and/or electrical heating elements embedded with the layers. Details of the calculation procedure and the methods used are presented.
Feedbacks Between Numerical and Analytical Models in Hydrogeology
NASA Astrophysics Data System (ADS)
Zlotnik, V. A.; Cardenas, M. B.; Toundykov, D.; Cohn, S.
2012-12-01
Hydrogeology is a relatively young discipline which combines elements of Earth science and engineering. Mature fundamental disciplines (e.g., physics, chemistry, fluid mechanics) have centuries-long history of mathematical modeling even prior to discovery of Darcy's law. Thus, in hydrogeology, relatively few classic analytical models (such those by Theis, Polubarinova-Kochina, Philip, Toth, Henry, Dagan, Neuman) were developed by the early 1970's. The advent of computers and practical demands refocused mathematical models towards numerical techniques. With more diverse but less mathematically-oriented training, most hydrogeologists shifted from analytical methods to use of standardized computational software. Spatial variability in internal properties and external boundary conditions and geometry, and the added complexity of chemical and biological processes will remain major challenges for analytical modeling. Possibly, analytical techniques will play a subordinate role to numerical approaches in many applications. On the other hand, the rise of analytical element modeling of groundwater flow is a strong alternative to numerical models when data demand and computational efficiency is considered. The hallmark of analytical models - transparency and accuracy - will remain indispensable for scientific exploration of complex phenomena and for benchmarking numerical models. Therefore, there will always be feedbacks and complementarities between numerical and analytical techniques, as well as a certain ideological schism among various views to modeling. We illustrate the idea of feedbacks by reviewing evolution of Joszef Toth's analytical model of gravity driven flow systems. Toth's (1963) approach was to reduce the flow domain to a rectangle which allowed for closed-form solution of the governing equations. Succeeding numerical finite-element models by Freeze and Witherspoon (1966-1968) explored the effects of geometry and heterogeneity on regional groundwater flow
Development of Numerical Grids for UZ Flow and Transport Modeling
J. Hinds
2001-12-18
This Analysis/Model Report (AMR) describes the methods used to develop numerical grids of the unsaturated hydrogeologic system beneath Yucca Mountain. Numerical grid generation is an integral part of the development of a complex, three-dimensional (3-D) model, such as the Unsaturated-Zone Flow and Transport Model (UZ Model) of Yucca Mountain. The resulting numerical grids, developed using current geologic, hydrogeologic, and mineralogic data, provide the necessary framework to: (1) develop calibrated hydrogeologic property sets and flow fields, (2) test conceptual hypotheses of flow and transport, and (3) predict flow and transport behavior under a variety of climatic and thermal loading conditions. Revision 00 of the work described herein follows the planning and work direction outlined in the ''Development of Numerical Grids for UZ Flow and Transport Modeling'' (CRWMS M&O 1999c). The technical scope, content, and management of ICN 01 of this AMR is currently controlled by the planning document, ''Technical Work Plan for Unsaturated Zone (UZ) Flow and Transport Process Model Report'' (BSC 2001a). The scope for the TBV resolution actions in this ICN is described in the ''Technical Work Plan for: Integrated Management of Technical Product Input Department'' (BSC 2001 b, Addendum B, Section 4.1). The steps involved in numerical grid development include: (1) defining the location of important calibration features, (2) determining model grid layers and fault geometry based on the Geologic Framework Model (GFM), the Integrated Site Model (ISM), and definition of hydrogeologic units (HGUs), (3) analyzing and extracting GFM and ISM data pertaining to layer contacts and property distributions, (4) discretizing and refining the two-dimensional (2-D), plan-view numerical grid, (5) generating the 3-D grid with finer resolution at the repository horizon and within the Calico Hills nonwelded (CHn) hydrogeologic unit, and (6) formulating the dual-permeability mesh. The products
Numerical Modeling of Shear Bands and Dynamic Fracture in Metals
NASA Astrophysics Data System (ADS)
McAuliffe, Colin James
Understanding the failure of metals at high strain rate is of utmost importance in the design of a broad range of engineering systems. Numerical methods offer the ability to analyze such complex physics and aid the design of structural systems. The objective of this research will be to develop reliable finite element models for high strain rate failure modelling, incorporating shear bands and fracture. Shear band modelling is explored first, and the subsequent developments are extended to incorporate fracture. Mesh sensitivity, the spurious dependence of failure on the discretization, is a well known hurdle in achieving reliable numerical results for shear bands and fracture, or any other strain softening model. Mesh sensitivity is overcome by regularization, and while details of regularization techniques may differ, all are similar in that a length scale is introduced which serves as a localization limiter. This dissertation contains two main contributions, the first of which presents several developments in shear band modeling. The importance of using a monolithic nonlinear solver in combination with a PDE model accounting for thermal diffusion is demonstrated. In contrast, excluding one or both of these components leads to unreliable numerical results. The Pian-Sumihara stress interpolants are also employed in small and finite deformation and shown to significantly improve the computational cost of shear band modelling. This is partly due to the fact that fewer unknowns than an irreducible discretization result from the same mesh, and more significantly, the fact that convergence of numerical results upon mesh refinement is improved drastically. This means coarser meshes are adequate to resolve shear bands, alleviating some of the computational cost of numerical modelling, which are notoriously significant. Since extremely large deformations are present during shear banding, a mesh to mesh transfer algorithm is presented for the Pian Sumihara element and used as
Numerical modelling and verification of Polish ventricular assist device.
Milenin, Andrzej; Kopernik, Magdalena; Jurkojć, Dorota; Gawlikowski, Maciej; Rusin, Tomasz; Darłak, Maciej; Kustosz, Roman
2012-01-01
The developed multiscale model of blood chamber of POLVAD (Polish ventricular assist device) was introduced. The tension test for polymer and digital image correlation (DIC) were performed for verification of the strains and displacements obtained in the numerical model of POLVAD_EXT. The numerical simulations were carried out in conditions given in the experiment to compare the results obtained on external surfaces of blood chamber of the POLVAD_EXT. The examined polymer applied in the POLVADs is sensitive to changes of temperature and this observation is considered in all prepared numerical models. The comparison of experimental and numerical results shows acceptable coincidence. There are some heterogeneous distributions of strains in experiment with respect to analysis of computed parameters. The comparison of two versions of blood chambers (POLVAD and POLVAD_EXT) in numerical analysis shows that POLVAD_EXT construction is better with respect to analysis of strain and stress. The maximum values of computed parameters are located in the regions between connectors on the internal surfaces of blood chambers of POLVAD. PMID:23140381
Uncertainty evaluation in numerical modeling of complex devices
NASA Astrophysics Data System (ADS)
Cheng, X.; Monebhurrun, V.
2014-10-01
Numerical simulation is an efficient tool for exploring and understanding the physics of complex devices, e.g. mobile phones. For meaningful results, it is important to evaluate the uncertainty of the numerical simulation. Uncertainty quantification in specific absorption rate (SAR) calculation using a full computer-aided design (CAD) mobile phone model is a challenging task. Since a typical SAR numerical simulation is computationally expensive, the traditional Monte Carlo (MC) simulation method proves inadequate. The unscented transformation (UT) is an alternative and numerically efficient method herein investigated to evaluate the uncertainty in the SAR calculation using the realistic models of two commercially available mobile phones. The electromagnetic simulation process is modeled as a nonlinear mapping with the uncertainty in the inputs e.g. the relative permittivity values of the mobile phone material properties, inducing an uncertainty in the output, e.g. the peak spatial-average SAR value.The numerical simulation results demonstrate that UT may be a potential candidate for the uncertainty quantification in SAR calculations since only a few simulations are necessary to obtain results similar to those obtained after hundreds or thousands of MC simulations.
Fuller, Sam M; Butz, Daniel R; Vevang, Curt B; Makhlouf, Mansour V
2014-09-01
Three-dimensional printing is being rapidly incorporated in the medical field to produce external prosthetics for improved cosmesis and fabricated molds to aid in presurgical planning. Biomedically engineered products from 3-dimensional printers are also utilized as implantable devices for knee arthroplasty, airway orthoses, and other surgical procedures. Although at first expensive and conceptually difficult to construct, 3-dimensional printing is now becoming more affordable and widely accessible. In hand surgery, like many other specialties, new or customized instruments would be desirable; however, the overall production cost restricts their development. We are presenting our step-by-step experience in creating a bone reduction clamp for finger fractures using 3-dimensional printing technology. Using free, downloadable software, a 3-dimensional model of a bone reduction clamp for hand fractures was created based on the senior author's (M.V.M.) specific design, previous experience, and preferences for fracture fixation. Once deemed satisfactory, the computer files were sent to a 3-dimensional printing company for the production of the prototypes. Multiple plastic prototypes were made and adjusted, affording a fast, low-cost working model of the proposed clamp. Once a workable design was obtained, a printing company produced the surgical clamp prototype directly from the 3-dimensional model represented in the computer files. This prototype was used in the operating room, meeting the expectations of the surgeon. Three-dimensional printing is affordable and offers the benefits of reducing production time and nurturing innovations in hand surgery. This article presents a step-by-step description of our design process using online software programs and 3-dimensional printing services. As medical technology advances, it is important that hand surgeons remain aware of available resources, are knowledgeable about how the process works, and are able to take advantage of
Numerical Model Studies of the Martian Mesoscale Circulations
NASA Technical Reports Server (NTRS)
Segal, Moti; Arritt, Raymond W.
1997-01-01
The study objectives were to evaluate by numerical modeling various possible mesoscale circulation on Mars and related atmospheric boundary layer processes. The study was in collaboration with J. Tillman of the University of Washington (who supported the study observationally). Interaction has been made with J. Prusa of Iowa State University in numerical modeling investigation of dynamical effects of topographically-influenced flow. Modeling simulations included evaluations of surface physical characteristics on: (i) the Martian atmospheric boundary layer and (ii) their impact on thermally and dynamically forced mesoscale flows. Special model evaluations were made in support of selection of the Pathfinder landing sites. J. Tillman's finding of VL-2 inter-annual temperature difference was followed by model simulations attempting to point out the forcing for this feature. Publication of the results in the reviewed literature in pending upon completion of the manuscripts in preparation as indicated later.
Numerical characterization and modeling of adiabatic slot film cooling
NASA Astrophysics Data System (ADS)
Voegele, Andrew
Film cooling is a technique used to protect critical surfaces in combustors, thrust chambers, turbines and nozzles from hot, chemically reacting gases. Accurately predicting the film's performance is especially challenging in the vicinity of the wall and the film injection plane due to the complex interactions of two highly turbulent, shearing, boundary layer flows. Properly characterizing the streams at the inlet of a numerical simulation and the choice of turbulence model are crucial to accurately predicting the decay of the film. To address these issues, this study employs a RANS solver that is used to model a film cooled wall. Menter's baseline model, and shear-stress transport model and the Spalart-Allmaras model are employed to determine the effect on film cooling predictions. Several methods for prescribing the inlet planes are explored. These numerical studies are compared with experimental data obtained in a UMD film cooling wind tunnel.
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.
Development, validation and application of numerical space environment models
NASA Astrophysics Data System (ADS)
Honkonen, Ilja
2013-10-01
Currently the majority of space-based assets are located inside the Earth's magnetosphere where they must endure the effects of the near-Earth space environment, i.e. space weather, which is driven by the supersonic flow of plasma from the Sun. Space weather refers to the day-to-day changes in the temperature, magnetic field and other parameters of the near-Earth space, similarly to ordinary weather which refers to changes in the atmosphere above ground level. Space weather can also cause adverse effects on the ground, for example, by inducing large direct currents in power transmission systems. The performance of computers has been growing exponentially for many decades and as a result the importance of numerical modeling in science has also increased rapidly. Numerical modeling is especially important in space plasma physics because there are no in-situ observations of space plasmas outside of the heliosphere and it is not feasible to study all aspects of space plasmas in a terrestrial laboratory. With the increasing number of computational cores in supercomputers, the parallel performance of numerical models on distributed memory hardware is also becoming crucial. This thesis consists of an introduction, four peer reviewed articles and describes the process of developing numerical space environment/weather models and the use of such models to study the near-Earth space. A complete model development chain is presented starting from initial planning and design to distributed memory parallelization and optimization, and finally testing, verification and validation of numerical models. A grid library that provides good parallel scalability on distributed memory hardware and several novel features, the distributed cartesian cell-refinable grid (DCCRG), is designed and developed. DCCRG is presently used in two numerical space weather models being developed at the Finnish Meteorological Institute. The first global magnetospheric test particle simulation based on the