Calculation of the wetted area of a planing hull with a chine

Summary The extent to which a low-aspect-ratio flat ship with a chined hull is wetted when planing at infinite Froude number is investigated. A numerical method of solution for the wetted area, which is applicable to more general planing problems, is presented. The results obtained by this method are compared with those found by solving the inverse problem of determining the hull shape which produces a given waterplane shape and are shown to be in excellent agreement. Results are also presented which indicate that a vertical chine may be used to fix the shape of the wetted region.     

Toward modeling and simulation of critical national infrastructure interdependencies

Modern society's physical health depends vitally upon a number of real, interdependent, critical infrastructure networks that deliver power, petroleum, natural gas, water, and communications. Its economic health depends on a number of other infrastructure networks, some virtual and some real, that link residences, industries, commercial sectors, and transportation sectors. The continued prosperity and national security of the US depends on our ability to understand the vulnerabilities of and analyze the performance of both the individual infrastructures and the entire interconnected system of infrastructures. Only then can we respond to potential disruptions in a timely and effective manner. Collaborative efforts among Sandia, other government agencies, private industry, and academia have resulted in realistic models for many of the individual component infrastructures. In this paper, we propose an innovative modeling and analysis framework to study the entire system of physical and economic infrastructures. That framework uses the existing individual models together with system dynamics, functional models, and nonlinear optimization algorithms. We describe this framework and demonstrate its potential use to analyze, and propose a response for, a hypothetical disruption.     

The influence of gravity on the performance of planing vessels in calm water

Usually gravity can be neglected for planing vessels at very high planing speed. However, if the planing speed becomes lower, the influence of gravity must be considered. A 2D+t theory with gravity effects is applied to study the steady performance of planing vessels at moderate planing speeds. In the framework of potential theory, a computer program based on a boundary-element method (BEM) in two dimensions is first developed, in which a new numerical model for the jet flow is introduced. The spray evolving from the free surface is cut to avoid the plunging breaker to impact on the underlying water. Further, flow separation along a chine line can be simulated. The BEM program is verified by comparing with similarity solutions and validated by comparing with drop tests of V-shaped cylinders. Then the steady motion of prismatic planing vessels is studied by using the 2D+t theory. The numerical results are compared with the results by Savitsky’s empirical formula and the experiments by Troesch. Significant nonlinearities in the restoring force coefficients can be seen from the results. Three-dimensional effects are discussed to explain the difference between the numerical results and the experimental results. Finally, in the comparison of results at high planing speed and moderate planing speed, it is shown that the gravity not only affects the free-surface profile around the hull, but also influences the hydrodynamic force on the hull surface.     

Interstitial instrumentation for therapeutic ultrasonic heating:modeling the discrete blood vessels

Instrumentation for interstitial ultrasound (US) heating is an important emerging technology in thermotherapy of deep seated tumors or those hard to reach by external devices. The instrumentation has special significance in case of radio-and/or chemotherapy resistant lesions. Its efficacy strongly depends on local tissue properties, especially local blood vessels. We evaluate effects of the vessels on temperature distribution elevated from basal by deposition of ultrasound energy. In the proposed model, we take into account several micron diameter vessels in proximity to the US four-applicator array. At large distances from the array, the volume is assigned a modified effective thermal conductivity. Our Finite Element Analysis of the so-defined problem shows that modelling under the assumption of constant, basal temperature across the vessels' lumen leads to erroneous results. The simulations agree best with experiments if fixed nodal temperature is applied at 60% of the lumen. We specify requirements on the array to avoid local underheating that could lead to performance failure of the instrumentation     

Mathematical modeling of the lysis of clots in blood vessels

A mathematical model of the lysis of clots in blood vessels has been constructed on the basis of convectiondiffusion equations. The fibrin clot is considered as an immobile solid phase, and the plasminogen, plasmin, and plasminogen activators — as dissolved liquid phases. As a result of the numerical solution of the model, dependences predicting the development of the lysis process have been obtained. An important role of clot swelling in the process of lysis has been revealed.     

Interaction between an elastic structure and free-surface flows: experimental versus numerical comparisons using the PFEM

The paper aims to introduce new fluid–structure interaction (FSI) tests to compare experimental results with numerical ones. The examples have been chosen for a particular case for which experimental results are not much reported. This is the case of FSI including free surface flows. The possibilities of the Particle Finite Element Method (PFEM) [1] for the simulation of free surface flows is also tested. The simulations are run using the same scale as the experiment in order to minimize errors due to scale effects. Different scenarios are simulated by changing the boundary conditions for reproducing flows with the desired characteristics. Details of the input data for all the examples studied are given. The aim is to identifying benchmark problems for FSI including free surface flows for future comparisons between different numerical approaches.     

Rapid motions of free-surface avalanches down curved and twisted channels and their numerical simulation

This paper presents a new model and discussions about the motion of avalanches from initiation to run-out over moderately curved and twisted channels of complicated topography and its numerical simulations. The model is a generalization of a well established and widely used depth-averaged avalanche model of Savage & Hutter and is published with all its details in Pudasaini & Hutter (Pudasaini & Hutter 2003 J. Fluid Mech. 495, 193-208). The intention was to be able to describe the flow of a finite mass of snow, gravel, debris or mud, down a curved and twisted corrie of nearly arbitrary cross-sectional profile. The governing equations for the distribution of the avalanche thickness and the topography-parallel depth-averaged velocity components are a set of hyperbolic partial differential equations. They are solved for different topographic configurations, from simple to complex, by applying a high-resolution non-oscillatory central differencing scheme with total variation diminishing limiter. Here we apply the model to a channel with circular cross-section and helical talweg that merges into a horizontal channel which may or may not become flat in cross-section. We show that run-out position and geometry depend strongly on the curvature and twist of the talweg and cross-sectional geometry of the channel, and how the topography is shaped close to run-out zones.     

舰船辐射噪声模拟技术研究

本文对舰船辐射噪声的模拟技术进行了研究,给出了线谱,连续谱噪声,噪声调制的模拟方法,进行了计算机模拟,研制了可编程模拟器。通过对模拟器输出信号的处理和分析,结果表明该模拟器能很好地再现舰船辐射噪声的典型声学特征。     

Evaluation of numerical schemes using different simulation methods for the continuous phase modeling of cyclone separators

The steady and unsteady state simulations of Stairmand cyclone separator were carried out to investigate the performance of different interpolation schemes for discretization of pressure gradient and advection terms. The RSM turbulence model was revisited to explore its simulation capability of PVC phenomenon and fluctuating velocity profiles of cyclone separators. The combination of Presto, SO, standard and BFW schemes for discretization of pressure gradient and FOU, power law, SOU, QUICK and MUSCL schemes for discretization of advection terms were studied. The double precision solver of Fluent 6.3.26 and modified RSM turbulence model constants of Jiao et al. (Chem. Eng. Technol. 30 (2007) 15–20) were also verified for simulation of cyclone separators. The predicted mean and fluctuating velocity profiles and pressure drop inside the cyclone separator with steady and unsteady simulations have been compared to experimental results available in literature.The steady state simulation failed to predict velocity profiles and pressure drop inside cyclone separator accurately, whereas the unsteady state simulation predicted velocity profiles, pressure drop and PVC phenomenon close to experimental values. The prediction of fluctuating velocity profile was better than previously reported work in the core region compared to the off core region. The present study revealed that the SOU scheme for discretization of advection terms of momentum, kinetic energy and its dissipation rate equations and the FOU scheme for Reynolds stresses together with the Presto scheme for discretization of pressure gradient with unsteady simulation are the optimum choice for simulation of cyclone separators.     

Transillumination optical tomography of tissue-engineered blood vessels: a Monte Carlo simulation

A Monte Carlo technique has been developed to simulate the transillumination laser computed tomography of tissue-engineered blood vessels. The blood vessel was modeled as a single cylinder layer mounted on a tubular mandrel. Sequences of images were acquired while rotating the mandrel. The tomographic image was reconstructed by applying a standard Radon transform. Angular discrimination was applied to simulate a spatial filter, which was used to reject multiply scattered photons. The simulation results indicated that the scattering effect can be overcome with angular discrimination because of the thin tissue thickness. However, any refractive-index mismatch among the tissue, the surrounding media, and the mandrel could produce significant distortions in the reconstructed image.     

Toward a formalization of affordance modeling for engineering design

When developing an artifact, designers must first capture and represent user needs. These needs can then be transformed into system requirements or objectives. The contribution of this work is rooted in the formalization of the affordance-based approach for capturing user needs in the early stages of design. This formalization comes in three forms: the first affordance basis for engineering design (a defined set of affordances), a formal structure for affordance statements, and a new relational model structure. This formalization is intended to improve model quality and consistency, while managing model creation resources. Further, this affordance-based approach to capturing user needs imposes a level of abstraction that forces solution independence yet is capable of capturing the large range of user needs. As such, the approach provides a structured approach to problem abstraction—the process of specifying user needs without reference to specific solutions. This affordance-based problem representation relies on other design process tools to help develop the actual artifact, which is also discussed.     

Fuzzy modeling and simulation for lead removal using micellar-enhanced ultrafiltration (MEUF)

In the present paper, a three factor, three-level response surface design based on Box-Behnken design (BBD) was developed for maximizing lead removal from aqueous solution using micellar-enhanced ultrafiltration (MEUF). Due to extremely complexity and nonlinearity of membrane separation processes, fuzzy logic (FL) models have been driven to simulate MEUF process under a wide range of initial and hydrodynamic conditions. Instead of using mathematical model, fuzzy logic approach provides a simpler and easier approach to describe the relationships between the processing variables and the metal rejection and permeation flux. Statistical values, which quantify the degree of agreement between experimental observations and numerically calculated values, were found greater than 91% for all cases. The results show that predicted values obtained from the fuzzy model were in very good agreement with the reported experimental data.     

Wet process rotary cement kilns: modeling and simulation

The knowledge of wet process kilns design and operation in the cement industry has not been improved substantially during the last decades since the pre-calcination process has been developed to become the technological standard. In spite of the tendency to replace wet process rotary kilns for cement production by modern dry process kilns with pre-calcination, there are still a substantial number of wet process kilns in operation. Their efficient operation requires an individually designed chain system to enhance the kiln performance. Currently the design of chain systems for wet process rotary cement kilns depends strongly on experience. Systematic theoretical and experimental investigations for the dimensioning of chain systems are not known. Based on a one-dimensional cell model, a simulation program for wet process rotary cement kilns with chain systems has been developed. With the aid of this simulation tool, the influence of chain systems on solid composition and temperature profiles, clinker throughput and fuel requirement is predicted and the design of chain systems can be optimized.     

Cleavage fracture modeling of pressure vessels under transient thermo-mechanical loading

The next generation of fracture assessment procedures for nuclear reactor pressure vessels (RPVs) will combine non-linear analyses of crack front response with stochastic treatments of crack size, shape, orientation, location, material properties and thermal-pressure transients. The projected computational demands needed to support stochastic approaches with detailed 3-D, non-linear stress analyses of vessels containing defects appear well beyond current and near-term capabilities. In the interim, 2-D models become appealing to approximate certain classes of critical flaws in RPVs, and have computational demands within reach for stochastic frameworks. The present work focuses on the capability of 2-D models to provide values for the Weibull stress fracture parameter with accuracy comparable to those from very detailed 3-D models. Weibull stress approaches provide one route to connect non-linear vessel response with fracture toughness values measured using small laboratory specimens. The embedded axial flaw located in the RPV wall near the cladding-vessel interface emerges from current linear-elastic, stochastic investigations as a critical contributor to the conditional probability of initiation. Three different types of 2-D models reflecting this configuration are subjected to a thermal-pressure transient characteristic of a critical pressurized thermal shock event. The plane-strain, 2-D models include: the modified boundary layer (MBL) model, the middle tension (M(T)) model, and the 2-D RPV model. The 2-D MBL model provides a high quality estimate for the Weibull stress but only in crack front regions with a positive T-stress. For crack front locations with low constraint (T-stress < 0), the M(T) specimen provides very accurate Weibull stress values but only for pressure load acting alone on the RPV. For RPVs under a combined thermal-pressure transient, Weibull stresses computed from the 2-D RPV model demonstrate close agreement with those computed from the corresponding crack front locations in the 3-D RPV model having large negative T-stresses. Applications of this family of 2-D models provide Weibull stress values in excellent agreement with very detailed 3-D models while retaining practical levels of computational effort.     

Computational materials: Multi-scale modeling and simulation of nanostructured materials

The paper provides details on the current approach to multi-scale modeling and simulation of advanced materials for structural applications. Examples are given that illustrate the suggested approaches to predicting the behavior and influencing the design of nanostructured materials such as high-performance polymers, composites, and nanotube-reinforced polymers. Primary simulation and measurement methods applicable to multi-scale modeling are outlined. Key challenges including verification and validation are highlighted and discussed.     

基于Agent的水声对抗仿真系统建模与仿真分析

分析了水声对抗仿真系统的组成和特点,基于Agent方法进行了水声对抗仿真建模,建立了Agent实体的功能层、信息层和认知层的模型结构,给出了基于Agent的系统结构和个体Agent的建模过程,从一个新的视点上研究了水声对抗分布仿真系统,最后总结了这种仿真体系的特点和长处。     

薄膜生长模型与计算机模拟

阐述了计算机模拟薄膜生长研究在纳米材料设计中的重要意义。通过最近多篇文献报道，归娄介绍了常见的薄膜生长模型、原理及方法，并详细对比了分子动力学、蒙特卡罗和量子力学3类方法、原理及应用特点。报道了国内外近几年薄膜生长模型与模拟的一些研究成果。结合不同的方法与原理，以分子动力学结合蒙特卡罗模型为重点，简要评述了这些模型及模拟特点，并提出了作者的一些见解。分析了多元化合物薄膜生长模拟的技术难点，指出了量子力学计算理论与方法同分子动力学、蒙特卡罗模型相结合的新型模型是解决这一难题的发展方向，展望了今后国内外的发展趋势。     

A thermomechanical fracture modeling and simulation for functionally graded solids using a residual-strain formulation

This paper addresses mixed-mode crack growth in two-dimensional functionally graded solids under thermomechanical loads, and investigates the effect of mechanical and thermal loads as well as the T-stress on their crack growth behavior. A novel residual strain-based formulation in the interaction integral method is developed and used for the accurate evaluation of mixed-mode stress intensity factors and/or the T-stress. Simulation of mixed-mode crack propagation in functionally graded materials including solid oxide fuel cells under thermomechanical loads is performed by means of the finite element method and the generalized interaction integrals in conjunction with a remeshing algorithm. An iterative procedure is used for crack growth simulation including the calculation of mixed-mode stress intensity factors and/or the T-stress by means of the generalized interaction integral method, determination of crack growth direction and crack initiation condition based on selected fracture criteria, and local automatic remeshing along the crack path. The present approach employs a user-defined crack increment at the beginning of the simulation. Crack trajectories and fracture parameters obtained by the present simulation for thermomechanical loads are assessed for some numerical examples in comparison with those for mechanical loads.