Physics Based Approach to Multi-domain Multilevel System Simulation

Transcription

1 Physics Based Approach to Multi-domain Multilevel System Simulation 2011 REGIONAL CONFERENCES 1 ANSYS, Inc. Mark Christini

2 ANSYS & Process Engineering System Sub-System ANSYS - Simplorer Mixed-Signal Multi-Domain System Simulator Model Order Reduction Cosimulation Component ANSYS Workbench Electrical Magnetic Fluid Mechanical Thermal Acoustic 2

3 Application Examples Linear Time Invariant (LTI) model and Co-Simulation using CFD analysis for battery thermal analysis Parasitic Extraction (R,L,C) using Electromagnetic BEM for EMC/EMI analysis Co-Simulation with CFD for Gravity Separator analysis State Space Stiffness Matrix extraction using Mechanical FEA for drill string analysis Nonlinear Lookup Table Extraction and Co-Simulation using CFD for pipe analysis VHDL-AMS, a Multi-Physics and Multi-Domain modeling language 3

4 Linear Time Invariant (LTI) Model Extracted from CFD Analysis Li-Ion Battery Example 4

5 Linear Time Invariant (LTI) Model Extracted from CFD Analysis Battery 4 Battery 5 Battery 6 Battery 1 Battery 2 Battery 3 Inputs: heat source to each battery Outputs: battery volume average temperature 5

6 Impulse Response due to First Cell Self heating is much larger than cross heat (note the scale is different) The first cell has strong impact on the second and third cell Only one of the six sets of impulse responses are shown here. 6

7 Automatic Generation of Foster Network or State Space Model Module for Simplorer Some cross heating elements have negligible contribution (less than 0.1% compared with self heating) and thus no Foster network 7

8 Reduced Order Model for CFD Batteries LTI Model Fluid Flow Region 8 Voc = *exp(-35*(abs(IBatt.V/Vinit))) *(abs(IBatt.V/Vinit)) *(abs(IBatt.V/Vinit))^ *(abs(IBatt.V/Vinit))^ /30.0*(U1.Temp_block_1-273)

9 Simplorer Architecture C/C++ User Defined Model Matlab Real Time Workshop Matlab Simulink Co-Simulation ANSYS MBD ANSYS Maxwell Simulation Data Bus/Simulator Coupling Technology Blocks: P=I 2 *R Circuits: States: T R = f(t) 9 Model Extraction: P Equivalent >> T Circuit, Impulse Response Extracted LTI, Stiffness Matrix Electromagnetic (FEA) Mechanical (FEA) Thermal (FEA/CFD) Fluidic (CFD) VHDL-AMS IF (domain = quiescent_domain) V0 == init_v; ELSE Current == cap*voltage'dot; END USE;

10 Co-Simulation between System and CFD Power Temperature 10 Simplorer Battery Model Co-simulation CFD Battery Model

11 Parasitic Extraction (R,L,C) using Electromagnetic BEM EMC/EMI Analysis 11

12 IGBT Characterization Dynamic IGBT accurately captures the switching waveforms Static IGBT for fast system simulations 12

13 L,R,C Extraction EMI/EMC: Automatic L,R,C Extraction and Network Model Current Distribution The structure is meshed using automatic and adaptive meshing 13

14 Frequency Dependent Parameters Extract the resistance, inductance, capacitance and conductance (R,L,C,G) parameters of the entire package Frequency can have a significant impact on the design performance 14 Low Frequency High Frequency

15 Integration into Simplorer DGraphSel1 NIGBT71.IC m m 2DGraphCon FFT m GS_I m m k 3.00k 10.00k k 1.00Meg 15 Extract Power Loss

16 Radiated Emission from HFSS m m 2DGraphCon1 Ansoft HFSS m GS_I k 3.00k 10.00k k 1.00Meg Multiplied mage plots by Simplorer Freq. res. Emission Test MagE@10m by specified inputs Normalized S para. 16

17 HFSS 17 mag 100 MHz, Power = W Spectrum (MHz) Power (W) The E field is very localized close to the module even at 100 MHz However, the very high power can lead to large E field at 1m for 1000w Power values of E field even far Spectrum (MHz) (V/m) (W) from the module This design is fine at 110MHz. E field E field at at 1m 1m f (V/m) (V/m Power E field at m for 1000w Power E field E field at at 1m m f Spectrum (MHz) (W) Spectrum (MHz) (V/m) (W) (V/m) (V/m

18 Co-Simulation with CFD Gravity Separator Analysis 18

19 Demo Case: Gravity Separator Oil-Water Separator model in Fluent Oil+Water Mixture Oil Outlet 4 m Desired Interface level 1 m 0 m Coalescer Water Outlet 19

20 PID Controller in Simplorer Proportional-Integral-Derivative (PID) controller Operates on the error signal and sends a correction signal correction ( t) t 1 de( t) K p e( t) e( t) dt Td Ti dt 0 Simplorer 20

21 Fluent-Simplorer Co-simulation Fluent Outlet Pressure Interface level Simplorer Desired level 21

22 Results Initial Level = 0.5 m Final Level = 1.0 m 22

23 State Space Stiffness Matrix extraction using Mechanical FEA Drill String Analysis 23

24 Y1 [ rad] VM _RO TB1. O MEG A [ r ad_per _sec] Pipe Borehole Diagram Borehead Top: Motor Drive Curve Info VM _RO TB1. O M EG A TR Model: Torque Source Pipe Segment Curve Info SM _RO TB1. PHI - SM _RO TB2. PHI TR SM _RO TB1. PHI - SM _RO TB3. PHI TR Model: Structural behavior for pipe segments Pipe Connector Model: Wall impact and friction Borehead Bottom: Drill Head Model: Friction and Grinding 24

25 + S + S F S + F T From Lumped Elements to FEA Extracted Models PHI P_in Y_in X_in Z_in Y X Z Lumped Parameter Pipe Element P_in Y_in X_in Z_in P_out Y_out X_out Z_out F 0 Wall Element F_ROTB3 P_port 0 P_out Y_out X_out Z_out pipe_r_in uf_in X_port S_x FM_TRB1 F_TRB3 EQUBL1 S_TRB1 Wall EQUBL2 F Y_port F_TRB4 S_y EQUBL 0 0 EQUBL F FM_TRB2 0 25

26 Y1 [ r ad] VM _RO TB1. O MEG A [ r ad_per _sec] TR TR TR Curve Info VM _RO TB1. O M EG A Curve Info SM _RO TB1. PHI -SM _RO TB2. PHI SM _RO TB1. PHI -SM _RO TB3. PHI R1 Rh R2 R5 R3 R4 R SM _TRB3. S [ m m ] Curve Info SM _TRB3. S TR Curve Info TR R1 Curve Info TR R2 Curve Info TR R3 Curve Info R4 TR Curve Info TR R5 Curve Info TR R6 Curve Info TR Rh R2 R5 R6 System Results Radial Displacement Results Drill_Wall2_Pos 2_Simplorer_2blks2_subcir_longer_pipe1 ANSOFT Curve Info R2 TR Near Borehole Top Drill_Wall5_Pos ANSOFT Curve Info R5 TR Time [s] Drill_Wall6_Pos ANSOFT Curve Info TR R Time [s] 26 Near Borehole Bottom Time [s]

27 Nonlinear Lookup Table Extraction and Co-Simulation using CFD Pipe Analysis 27

28 PIPE1.P [kpascal] Level1 Level1 Piping System Pipe 1 Pressure A_LumpedElements ANSOFT PIPE PIPE5 Level1 Level1 PIPE4 0 Level1 PIPE7 PIPE PIPE3 PIPE1 Level1 PIPE2 Level1 Level1 P Q1 P Time [ms] MX1:

29 Simulate: Solving the Physics Correctly Co-Simulation Total Pressure Velocity Vectors 29

30 Verify: The Physics is Simulated Correctly in the System Workbench Integration 30

31 Level1 Level1 System & CFD Co-Simulation Less Than 1% Error Between Full CFD and System Behavioral Model System Simulation using Extracted Behavioral Model Lumped Circuit Element kpa kpa 13.9 k Pa PIPE6 QM1 CFD_PQ NL RHYD2 PIPE5 PIPE3 PIPE2 Level1 Level1 Level1 PIPE4 0 Level1 Level1 PIPE7 PIPE8 0 P2 0 Q1 P

32 VHDL-AMS, a Multi-Physics and Multi-Domain Modeling Language 32

33 VHDL-AMS Model Construct Entity Interface description of a subsystem or physical device Specifies input and output ports to the model Architecture Behavior description Can be dataflow, structural, procedural, etc Modeling can deal with both analog (continuous) and digital (discrete) domains input ports Entity Architecture 1 Architecture 2 Architecture 3 In-out ports output ports 33

34 VHDL-AMS Code for a Resistor Entity Interface description of a subsystem or physical device Specifies input and output ports to the model ENTITY resistor IS generic ( RNOM : real := 1.0e+03; ); port ( QUANTITY Temp : IN T:= 300.0; TERMINAL p : electrical; TERMINAL m : electrical ); END ENTITY resistor; Current Voltage p m resistor Resistance value Temperature The resistor model has one model constant, one input quantity, and two terminals 34

35 VHDL-AMS Code for a Resistor Architecture Description of the model and no solving information is required ARCHITECTURE Simple OF resistor IS QUANTITY voltage ACROSS current THROUGH p TO m; BEGIN voltage == current * RNOM; END ARCHITECTURE Simple; No solving information is needed!! ARCHITECTURE 2ndOrderT OF resistor IS QUANTITY voltage ACROSS current THROUGH p TO m; BEGIN voltage == current*rnom*( e- 2*(T-300.)+1.1e-3*(T-300.)**2); END ARCHITECTURE 2ndOrderT; A second architecture is possible!! 35

36 VHDL-AMS Code for a Capacitor Entity ENTITY cap IS generic ( capacitance : real := 1.0e-03; ); port ( TERMINAL p : electrical; TERMINAL m : electrical ); END ENTITY cap; Architecture ARCHITECTURE arch_cap OF cap IS quantity voltage across current through p to m; BEGIN current == capacitance * voltage'dot; END ARCHITECTURE cap; The capacitor model has one model constant and two terminals The model description essentially has one single line!! 36

37 What About PDEs Boundary Value Problems? Example: φ 1 φ 2 φ 3 φ 4 φ x=0 h L=1.0 Entity ENTITY steady_state_boundary_value IS generic ( phibc: real := 0.0 ); END ENTITY steady_state_boundary_value; Architecture ARCHITECTURE arch_steady_state_boundary_value OF steady_state_boundary_value IS quantity phi1,phi2,phi3,phi4,phi5 : real; constant h : real := 0.25; BEGIN phi1 == phibc; (phi3-phi1)/(2.0*h) == 1.0*h; (phi4-phi2)/(2.0*h) == 2.0*h; (phi5-phi3)/(2.0*h) == 3.0*h; (1.0*phi3-4.0*phi4+3.0*phi5)/(2.0*h) == 4.0*h; END ARCHITECTURE arch_steady_state_boundary_value; 37

38 CAD/FEA Depth Linkt to 0D-1D Knowledge Management Circuit/System Engineering Simulations for Realizing SDPD System Core Optimization Activity Order: 0D Order: 1D Component Order: 2D/3D Electrical Magnetic Fluid Mechanical Thermal 38 Process Compression HPC WB Adoption: Multi-Domain Breadth CAD, PLM, etc.

39 39