Avoiding Overfitting in Inverse Modeling of Chloride Migration in Concrete

Inverse analysis for reactive transport of chlorides through concrete in the presence of electric field is presented. The model is solved using MATLAB’s built-in solvers “pdepe.m” and “ode15s.m”. The results from the model are compared with experimental measurements from accelerated migration test and a function representing the lack of fit is formed. This function is optimised with respect to varying amount of key parameters defining the model. Levenberg-Marquardt trust-region optimisation approach is employed. The paper presents a method by which the degree of inter-dependency between parameters and sensitivity (significance) of each parameter towards model predictions can be studied on models with or without clearly defined governing equations. Eigen value analysis of the Hessian matrix was employed to investigate and avoid over-parametrisation in inverse analysis. We investigated simultaneous fitting of parameters for diffusivity, chloride binding as defined by Freundlich isotherm (thermodynamic) and binding rate (kinetic parameter). Fitting of more than 2 parameters, simultaneously, demonstrates a high degree of parameter inter-dependency. This finding is significant as mathematical models for representing chloride transport rely on several parameters for each mode of transport (i.e., diffusivity, binding, etc.), which combined may lead to unreliable simultaneous estimation of parameters.

Least Squares Approach to Inverse Problems in Musculoskeletal Biomechanics

Inverse simulations of musculoskeletal models computes the internal forces such as muscle and joint reaction forces, which are hard to measure, using the more easily measured motion and external forces as input data. Because of the difficulties of measuring muscle forces and joint reactions, simulations are hard to validate. One way of reducing errors for the simulations is to ensure that the mathematical problem is well-posed. This paper presents a study of regularity aspects for an inverse simulation method, often called forward dynamics or dynamical optimization, that takes into account both measurement errors and muscle dynamics. The simulation method is explained in detail. Regularity is examined for a test problem around the optimum using the approximated quadratic problem. The results shows improved rank by including a regularization term in the objective that handles the mechanical over-determinancy. Using the 3-element Hill muscle model the chosen regularization term is the norm of the activation. To make the problem full-rank only the excitation bounds should be included in the constraints. However, this results in small negative values of the activation which indicates that muscles are pushing and not pulling. Despite this unrealistic behavior the error maybe small enough to be accepted for specific applications. These results is a starting point start for achieving better results of inverse musculoskeletal simulations from a numerical point of view.

Laboratory performance evaluation of CIR-emulsion and its comparison against CIR-foam test results from Phase III, TR-578

Currently, no standard mix design procedure is available for CIR-emulsion in Iowa. The CIR-foam mix design process developed during the previous phase is applied for CIR-emulsion mixtures with varying emulsified asphalt contents. Dynamic modulus test, dynamic creep test, static creep test and raveling test were conducted to evaluate the short- and long-term performance of CIR-emulsion mixtures at various testing temperatures and loading conditions. A potential benefit of this research is a better understanding of CIR-emulsion material properties in comparison with those of CIR-foam material that would allow for the selection of the most appropriate CIR technology and the type and amount of the optimum stabilization material. Dynamic modulus, flow number and flow time of CIR-emulsion mixtures using CSS-1h were generally higher than those of HFMS-2p. Flow number and flow time of CIR-emulsion using RAP materials from Story County was higher than those from Clayton County. Flow number and flow time of CIR-emulsion with 0.5% emulsified asphalt was higher than CIR-emulsion with 1.0% or 1.5%. Raveling loss of CIR-emulsion with 1.5% emulsified was significantly less than those with 0.5% and 1.0%. Test results in terms of dynamic modulus, flow number, flow time and raveling loss of CIR-foam mixtures are generally better than those of CIR-emulsion mixtures. Given the limited RAP sources used for this study, it is recommended that the CIR-emulsion mix design procedure should be validated against several RAP sources and emulsion types.

Some Inverse Problems in Analysis and Geometry

Thesis (Ph.D.)--University of Washington, 2016-08

A domain decomposition based algorithm for non-linear 2D inverse heat conduction problems

Inverse heat conduction problems (IHCPs) appear in many important scientific and technological fields. Hence analysis, design, implementation and testing of inverse algorithms are also of great scientific and technological interest. The numerical simulation of 2-D and –D inverse (or even direct) problems involves a considerable amount of computation. Therefore, the investigation and exploitation of parallel properties of such algorithms are equally becoming very important. Domain decomposition (DD) methods are widely used to solve large scale engineering problems and to exploit their inherent ability for the solution of such problems.

The Inverse Jacobian Control Method Applied to a Four Degree of Freedom Confined Space Robot

Thesis (Master's)--University of Washington, 2016-06

Tensor Completion for Multidimensional Inverse Problems with Applications to Magnetic Resonance Relaxometry

This thesis deals with tensor completion for the solution of multidimensional inverse problems. We study the problem of reconstructing an approximately low rank tensor from a small number of noisy linear measurements. New recovery guarantees, numerical algorithms, non-uniform sampling strategies, and parameter selection algorithms are developed. We derive a fixed point continuation algorithm for tensor completion and prove its convergence. A restricted isometry property (RIP) based tensor recovery guarantee is proved. Probabilistic recovery guarantees are obtained for sub-Gaussian measurement operators and for measurements obtained by non-uniform sampling from a Parseval tight frame. We show how tensor completion can be used to solve multidimensional inverse problems arising in NMR relaxometry. Algorithms are developed for regularization parameter selection, including accelerated k-fold cross-validation and generalized cross-validation. These methods are validated on experimental and simulated data. We also derive condition number estimates for nonnegative least squares problems. Tensor recovery promises to significantly accelerate N-dimensional NMR relaxometry and related experiments, enabling previously impractical experiments. Our methods could also be applied to other inverse problems arising in machine learning, image processing, signal processing, computer vision, and other fields.

Hierarchical Reconstruction Method for Solving Ill-posed Linear Inverse Problems

We present a detailed analysis of the application of a multi-scale Hierarchical Reconstruction method for solving a family of ill-posed linear inverse problems. When the observations on the unknown quantity of interest and the observation operators are known, these inverse problems are concerned with the recovery of the unknown from its observations. Although the observation operators we consider are linear, they are inevitably ill-posed in various ways. We recall in this context the classical Tikhonov regularization method with a stabilizing function which targets the specific ill-posedness from the observation operators and preserves desired features of the unknown. Having studied the mechanism of the Tikhonov regularization, we propose a multi-scale generalization to the Tikhonov regularization method, so-called the Hierarchical Reconstruction (HR) method. First introduction of the HR method can be traced back to the Hierarchical Decomposition method in Image Processing. The HR method successively extracts information from the previous hierarchical residual to the current hierarchical term at a finer hierarchical scale. As the sum of all the hierarchical terms, the hierarchical sum from the HR method provides an reasonable approximate solution to the unknown, when the observation matrix satisfies certain conditions with specific stabilizing functions. When compared to the Tikhonov regularization method on solving the same inverse problems, the HR method is shown to be able to decrease the total number of iterations, reduce the approximation error, and offer self control of the approximation distance between the hierarchical sum and the unknown, thanks to using a ladder of finitely many hierarchical scales. We report numerical experiments supporting our claims on these advantages the HR method has over the Tikhonov regularization method.

Processing of oil enriched fish (Hypophthalmichthys molitrix) sausage using emulsion technology

In most countries along with various food products, fish sausage is supplied in different formulas. Unfortunately, in our country because of different reasons, production and supply of fish sausage in industrial level has not yet been successful and some efforts taken, has also been doomed to failure or not welcomed. Fat fish is a rich source of poly unsaturated fatty acids (PUFA) and co-3. In this research, efforts have been made to produce and enrich sausage with fish oil and maintenance of fatty acids has also been experimented using gas chromatography along with heating process. The stages of producing ground fish and fish sausage are as the following: Transferring and preparing fish, washing the cleared fish, filleting, separating fillet steak, washing and drying them, Refining meat, Producing and homogenizing mixture from basic ingredients in a cutter, filling, knotting and heat processing. The fish sausage produced by this method tried and welcomed by the subjects. In the product in which fish meat was used, the subjects was not recognized fish flavor and taste and when in addition to fish meat, fish oil was used during enrichment, the flavor and taste of fish was considered as highly acceptable. TVN measurement of the produced fish sausage was kept in the refrigerator in two month was at a maximum of 16.5, the amount of peroxide was at a maximum 1.5% after the period of two months. During this period the Colony count was at maximum of 19.5 x 104, the high maximum of the number of coliforms was 10/gr, and for mold and yeast 83/gr , but Escherichia coli, Staphylococcus aureus, Salmonella and Clostridium perfringens were not found. The protein of the resulting product was 15-18%, lipid at about 11-15% and moisture 60-65%. Comparing fatty acids, including unsaturated fatty acids in ground and oil fish used in producing fish sausage with those of fish sausage showed that the heat used in processing had the least effect on fatty acids of the meat and oil used here and the resulting fish sausage is considered as food for good health.

Utilisation d’un modèle de co-culture cellulaire pour l’évaluation in vitro du transport inverse du cholestérol

De nouveaux modèles cellulaires in vitro par transfert de milieu et par coculture ont été mis au point afin d’évaluer la capacité des HDL à éliminer l’excès de cholestérol des tissus périphériques et de le transporter vers le foie afin d’être excrété par le foie, un processus nommé le transport inverse du cholestérol (TIC). Le système cellulaire par transfert in vitro où des macrophages J774 sont gorgés de LDL acétylées et marqués au 3H-cholestérol a été préalablement établi afin de mesurer par scintillation l’efflux de cholestérol marqué vers le milieu de culture contenant des accepteurs de cholestérol. Ce milieu conditionné est transféré sur des cellules HepG2 afin d’étudier l’influx du cholestérol marqué. Ce dernier nous permet d’observer un transport de cholestérol de 25 % hors des J774 et un transport de 39 000 cpm dans les HepG2 en utilisant un milieu contenant 2 % de sérums humains mis en commun. Une stimulation des cellules J774 par l’AMPc augmente l’efflux et l’influx d’environ 45 %. Des tests de preuve de concept ont été effectués sur le système cellulaire par co-culture qui utilise des chambres de Boyden où les J774 sont localisées au fond d’un puits et les HepG2 dans un insert, et où le milieu est partagé entre les deux types cellulaires. On a déterminé qu’une confluence densité de 60 000 cellules/cm2 sur un insert constitué d’une membrane de polyester avec des pores de 3,0 μm, sans autre revêtement, permet d’observer un influx spécifique au sérum d’environ 6 000 cpm associés aux cellules HepG2, où 50 % des comptes radioactifs sont dans les cellules et l’autre moitié présente à la surface cellulaire.

Inverse identification of the bending stiffness of a braided polyethylene twine subject to large deformation: application to the identification of the mesh opening rigidity of fishing nets.

The evaluation of the mesh opening stiffness of fishing nets is an important issue in assessing the selectivity of trawls. It appeared that a larger bending rigidity of twines decreases the mesh opening and could reduce the escapement of fish. Nevertheless, netting structure is complex. A netting is made up of braided twines made of polyethylene or polyamide. These twines are tied with non-symmetrical knots. Thus, these assemblies develop contact-friction interactions. Moreover, the netting can be subject to large deformation. In this study, we investigate the responses of netting samples to different types of solicitations. Samples are loaded and unloaded with creep and relaxation stages, with different boundary conditions. Then, two models have been developed: an analytical model and a finite element model. The last one was used to assess, with an inverse identification algorithm, the bending stiffness of twines. In this paper, experimental results and a model for netting structures made up of braided twines are presented. During dry forming of a composite, for example, the matrix is not present or not active, and relative sliding can occur between constitutive fibres. So an accurate modelling of the mechanical behaviour of fibrous material is necessary. This study offers experimental data which could permit to improve current models of contact-friction interactions [4], to validate models for large deformation analysis of fibrous materials [1] on a new experimental case, then to improve the evaluation of the mesh opening stiffness of a fishing net

Inverse kinematics for upper limb compound movement estimation in exoskeleton-assisted rehabilitation

Robot-Assisted Rehabilitation (RAR) is relevant for treating patients affected by nervous system injuries (e.g., stroke and spinal cord injury) -- The accurate estimation of the joint angles of the patient limbs in RAR is critical to assess the patient improvement -- The economical prevalent method to estimate the patient posture in Exoskeleton-based RAR is to approximate the limb joint angles with the ones of the Exoskeleton -- This approximation is rough since their kinematic structures differ -- Motion capture systems (MOCAPs) can improve the estimations, at the expenses of a considerable overload of the therapy setup -- Alternatively, the Extended Inverse Kinematics Posture Estimation (EIKPE) computational method models the limb and Exoskeleton as differing parallel kinematic chains -- EIKPE has been tested with single DOFmovements of the wrist and elbow joints -- This paper presents the assessment of EIKPEwith elbow-shoulder compoundmovements (i.e., object prehension) -- Ground-truth for estimation assessment is obtained from an optical MOCAP (not intended for the treatment stage) -- The assessment shows EIKPE rendering a good numerical approximation of the actual posture during the compoundmovement execution, especially for the shoulder joint angles -- This work opens the horizon for clinical studies with patient groups, Exoskeleton models, and movements types --